ABC of the future Pesonen, Henri; Simola, Umberto; Köhn‐Luque, Alvaro ...
International statistical review,
August 2023, 2023-08-00, 20230801, Letnik:
91, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Summary
Approximate Bayesian computation (ABC) has advanced in two decades from a seminal idea to a practically applicable inference tool for simulator‐based statistical models, which are becoming ...increasingly popular in many research domains. The computational feasibility of ABC for practical applications has been recently boosted by adopting techniques from machine learning to build surrogate models for the approximate likelihood or posterior and by the introduction of a general‐purpose software platform with several advanced features, including automated parallelisation. Here we demonstrate the strengths of the advances in ABC by going beyond the typical benchmark examples and considering real applications in astronomy, infectious disease epidemiology, personalised cancer therapy and financial prediction. We anticipate that the emerging success of ABC in producing actual added value and quantitative insights in the real world will continue to inspire a plethora of further applications across different fields of science, social science and technology.
Our understanding of how speed and persistence of cell migration affects the growth rate and size of tumors remains incomplete. To address this, we developed a mathematical model wherein cells ...migrate in two-dimensional space, divide, die or intravasate into the vasculature. Exploring a wide range of speed and persistence combinations, we find that tumor growth positively correlates with increasing speed and higher persistence. As a biologically relevant example, we focused on Golgi fragmentation, a phenomenon often linked to alterations of cell migration. Golgi fragmentation was induced by depletion of Giantin, a Golgi matrix protein, the downregulation of which correlates with poor patient survival. Applying the experimentally obtained migration and invasion traits of Giantin depleted breast cancer cells to our mathematical model, we predict that loss of Giantin increases the number of intravasating cells. This prediction was validated, by showing that circulating tumor cells express significantly less Giantin than primary tumor cells. Altogether, our computational model identifies cell migration traits that regulate tumor progression and uncovers a role of Giantin in breast cancer progression.
Estrogen receptor positive breast cancer is frequently treated with anti-hormonal treatment such as aromatase inhibitors (AI). Interestingly, a high body mass index has been shown to have a negative ...impact on AI efficacy, most likely due to disturbances in steroid metabolism and adipokine production. Here, we propose a mathematical model based on a system of ordinary differential equations to investigate the effect of high-fat diet on tumor growth. We inform the model with data from mouse experiments, where the animals are fed with high-fat or control (normal) diet. By incorporating AI treatment with drug resistance into the model and by solving optimal control problems we found differential responses for control and high-fat diet. To the best of our knowledge, this is the first attempt to model optimal anti-hormonal treatment for breast cancer in the presence of drug resistance. Our results underline the importance of considering high-fat diet and obesity as factors influencing clinical outcomes during anti-hormonal therapies in breast cancer patients.
Mathematical modeling and simulation is a promising approach to personalized cancer medicine. Yet, the complexity, heterogeneity and multi‐scale nature of cancer pose significant computational ...challenges. Coupling discrete cell‐based models with continuous models using hybrid cellular automata (CA) is a powerful approach for mimicking biological complexity and describing the dynamical exchange of information across different scales. However, when clinically relevant cancer portions are taken into account, such models become computationally very expensive. While efficient parallelization techniques for continuous models exist, their coupling with discrete models, particularly CA, necessitates more elaborate solutions. Building upon FEniCS, a popular and powerful scientific computing platform for solving partial differential equations, we developed parallel algorithms to link stochastic CA with differential equations (https://bitbucket.org/HTasken/cansim). The algorithms minimize the communication between processes that share CA neighborhood values while also allowing for reproducibility during stochastic updates. We demonstrated the potential of our solution on a complex hybrid cellular automaton model of breast cancer treated with combination chemotherapy. On a single‐core processor, we obtained nearly linear scaling with an increasing problem size, whereas weak parallel scaling showed moderate growth in solving time relative to increase in problem size. Finally, we applied the algorithm to a problem that is 500 times larger than previous work, allowing us to run personalized therapy simulations based on heterogeneous cell density and tumor perfusion conditions estimated from magnetic resonance imaging data on an unprecedented scale.
We developed an efficient and parallelized solver for a stochastic hybrid cellular automaton model of breast cancer treated with combination chemotherapy. We extended the parallelization capability of FEniCS by minimizing communication between processes that shared cellular automata neighborhood values while maintaining reproducibility during stochastic updates. We demonstrated the solver's scaling ability by applying to a problem that is 500 times larger than previous work.
Tumor heterogeneity is a complex and widely recognized trait that poses significant challenges in developing effective cancer therapies. In particular, many tumors harbor a variety of subpopulations ...with distinct therapeutic response characteristics. Characterizing this heterogeneity by determining the subpopulation structure within a tumor enables more precise and successful treatment strategies. In our prior work, we developed PhenoPop, a computational framework for unravelling the drug-response subpopulation structure within a tumor from bulk high-throughput drug screening data. However, the deterministic nature of the underlying models driving PhenoPop restricts the model fit and the information it can extract from the data. As an advancement, we propose a stochastic model based on the linear birth-death process to address this limitation. Our model can formulate a dynamic variance along the horizon of the experiment so that the model uses more information from the data to provide a more robust estimation. In addition, the newly proposed model can be readily adapted to situations where the experimental data exhibits a positive time correlation. We test our model on simulated data (in silico) and experimental data (in vitro), which supports our argument about its advantages.
Mechanistic learning refers to the synergistic combination of mechanistic mathematical modeling and data-driven machine or deep learning. This emerging field finds increasing applications in ...(mathematical) oncology. This review aims to capture the current state of the field and provides a perspective on how mechanistic learning may progress in the oncology domain. We highlight the synergistic potential of mechanistic learning and point out similarities and differences between purely data-driven and mechanistic approaches concerning model complexity, data requirements, outputs generated, and interpretability of the algorithms and their results. Four categories of mechanistic learning (sequential, parallel, extrinsic, intrinsic) of mechanistic learning are presented with specific examples. We discuss a range of techniques including physics-informed neural networks, surrogate model learning, and digital twins. Example applications address complex problems predominantly from the domain of oncology research such as longitudinal tumor response predictions or time-to-event modeling. As the field of mechanistic learning advances, we aim for this review and proposed categorization framework to foster additional collaboration between the data- and knowledge-driven modeling fields. Further collaboration will help address difficult issues in oncology such as limited data availability, requirements of model transparency, and complex input data which are embraced in a mechanistic learning framework.
The purpose of the present study is to investigate if consumption and supply hypoxia (CSH) MR-imaging can depict breast cancer hypoxia, using the CSH-method initially developed for prostate cancer. ...Furthermore, to develop a generalized pan-cancer application of the CSH-method that doesn't require a hypoxia reference standard for training the CSH-parameters. In a cohort of 69 breast cancer patients, we generated, based on the principles of intravoxel incoherent motion modelling, images reflecting cellular density (apparent diffusion coefficient; ADC) and vascular density (perfusion fraction; fp). Combinations of the information in these images were compared to a molecular hypoxia score made from gene expression data, aiming to identify a way to apply the CSH-methodology in breast cancer. Attempts to adapt previously proposed models for prostate cancer included direct transfers and model parameter rescaling. A novel approach, based on rescaling ADC and fp data to give more nuanced response in the relevant physiologic range, was also introduced. The new CSH-method was validated in a prostate cancer cohort with known hypoxia status. The proposed CSH-method gave estimates of hypoxia that was strongly correlated to the molecular hypoxia score in breast cancer, and hypoxia as measured in pathology slices stained with pimonidazole in prostate cancer. The generalized approach to CSH-imaging depicted hypoxia in both breast and prostate cancers and requires no model training. It is easy to implement using readily available technology and encourages further investigation of CSH-imaging in other cancer entities and in other settings, with the goal being to overcome hypoxia-induced resistance to treatment.
The role of dysfunction of the single ventricle in Fontan failure is incompletely understood. We aimed to evaluate hemodynamic responses to preload increase in Fontan circulation, to determine ...whether circulatory limitations in different locations identified by experimental preload increase are associated with cardiorespiratory fitness (CRF), and to assess the impact of left versus right ventricular morphology. In 38 consecutive patients (median age = 16.6 years, 16 females), heart catheterization was supplemented with a rapid 5-mL/kg body weight volume expansion. Central venous pressure (CVP), ventricular end-diastolic pressure (VEDP), and peak systolic pressure were averaged for 15‒30 s, 45‒120 s, and 4‒6 min (steady state), respectively. CRF was assessed by peak oxygen consumption (VO
2peak
) and ventilatory threshold (VT). Median CVP increased from 13 mmHg at baseline to 14.5 mmHg (
p
< 0.001) at steady state. CVP increased by more than 20% in eight patients. Median VEDP increased from 10 mmHg at baseline to 11.5 mmHg (
p
< 0.001). Ten patients had elevated VEDP at steady state, and in 21, VEDP increased more than 20%. The transpulmonary pressure difference (CVP‒VEDP) and CVP were consistently higher in patients with right ventricular morphology across repeated measurements. CVP at any stage was associated with VO
2peak
and VT. VEDP after volume expansion was associated with VT. Preload challenge demonstrates the limitations beyond baseline measurements. Elevation of both CVP and VEDP are associated with impaired CRF. Transpulmonary flow limitation was more pronounced in right ventricular morphology. Ventricular dysfunction may contribute to functional impairment after Fontan operation in young adulthood.
ClinicalTrials.gov
identifier
NCT02378857
Prediction models should be externally validated to assess their performance before implementation. Several prediction models for coronavirus disease-19 (COVID-19) have been published. This ...observational cohort study aimed to validate published models of severity for hospitalized patients with COVID-19 using clinical and laboratory predictors.
Prediction models fitting relevant inclusion criteria were chosen for validation. The outcome was either mortality or a composite outcome of mortality and ICU admission (severe disease). 1295 patients admitted with symptoms of COVID-19 at Kings Cross Hospital (KCH) in London, United Kingdom, and 307 patients at Oslo University Hospital (OUH) in Oslo, Norway were included. The performance of the models was assessed in terms of discrimination and calibration.
We identified two models for prediction of mortality (referred to as Xie and Zhang1) and two models for prediction of severe disease (Allenbach and Zhang2). The performance of the models was variable. For prediction of mortality Xie had good discrimination at OUH with an area under the receiver-operating characteristic (AUROC) 0.87 95% confidence interval (CI) 0.79-0.95 and acceptable discrimination at KCH, AUROC 0.79 0.76-0.82. In prediction of severe disease, Allenbach had acceptable discrimination (OUH AUROC 0.81 0.74-0.88 and KCH AUROC 0.72 0.68-0.75). The Zhang models had moderate to poor discrimination. Initial calibration was poor for all models but improved with recalibration.
The performance of the four prediction models was variable. The Xie model had the best discrimination for mortality, while the Allenbach model had acceptable results for prediction of severe disease.