1
|
Xiao H, Chen H, Zhang L, Duolikun M, Zhen B, Kuerban S, Li X, Wang Y, Chen L, Lin J. Cytoskeletal gene alterations linked to sorafenib resistance in hepatocellular carcinoma. World J Surg Oncol 2024; 22:152. [PMID: 38849867 PMCID: PMC11157844 DOI: 10.1186/s12957-024-03417-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2024] [Accepted: 05/20/2024] [Indexed: 06/09/2024] Open
Abstract
BACKGROUND Although sorafenib has been consistently used as a first-line treatment for advanced hepatocellular carcinoma (HCC), most patients will develop resistance, and the mechanism of resistance to sorafenib needs further study. METHODS Using KAS-seq technology, we obtained the ssDNA profiles within the whole genome range of SMMC-7721 cells treated with sorafenib for differential analysis. We then intersected the differential genes obtained from the analysis of hepatocellular carcinoma patients in GSE109211 who were ineffective and effective with sorafenib treatment, constructed a PPI network, and obtained hub genes. We then analyzed the relationship between the expression of these genes and the prognosis of hepatocellular carcinoma patients. RESULTS In this study, we identified 7 hub ERGs (ACTB, CFL1, ACTG1, ACTN1, WDR1, TAGLN2, HSPA8) related to drug resistance, and these genes are associated with the cytoskeleton. CONCLUSIONS The cytoskeleton is associated with sorafenib resistance in hepatocellular carcinoma. Using KAS-seq to analyze the early changes in tumor cells treated with drugs is feasible for studying the drug resistance of tumors, which provides reference significance for future research.
Collapse
Affiliation(s)
- Hong Xiao
- Key Laboratory of Tropical Biological Resources of Ministry of Education, School of Pharmaceutical Sciences, Hainan University, Hainan, China
| | - Hangyu Chen
- Department of Pharmacy, Peking University Third Hospital, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China
| | - Lei Zhang
- Department of Pharmacy, Peking University Third Hospital, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China
| | - Maimaitiyasen Duolikun
- Key Laboratory of Tropical Biological Resources of Ministry of Education, School of Pharmaceutical Sciences, Hainan University, Hainan, China
| | - Baixin Zhen
- Department of Pharmacology, Xinjiang Medical University, Urumqi, China
| | - Subinuer Kuerban
- Department of Pharmacology, Xinjiang Medical University, Urumqi, China
| | - Xuehui Li
- Department of Pharmacology, Xinjiang Medical University, Urumqi, China
| | - Yuxi Wang
- Department of Pharmacy, Peking University Third Hospital, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China
| | - Long Chen
- Department of Pharmacy, Peking University Third Hospital, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China.
- Peking University, Third Hospital Cancer Center, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China.
| | - Jian Lin
- Key Laboratory of Tropical Biological Resources of Ministry of Education, School of Pharmaceutical Sciences, Hainan University, Hainan, China.
- Department of Pharmacy, Peking University Third Hospital, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China.
- Synthetic and Functional Biomolecules Center, Beijing National Laboratory for Molecular Sciences, Peking University, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China.
- Peking University, Third Hospital Cancer Center, 49 Huayuan North Rd, Haidian District, Beijing, 100191, China.
| |
Collapse
|
2
|
Hu A, Ojwang' AME, Olumoyin KD, Rejniak KA. LinG3D: visualizing the spatio-temporal dynamics of clonal evolution. BMC Bioinformatics 2024; 25:201. [PMID: 38802748 PMCID: PMC11131251 DOI: 10.1186/s12859-024-05813-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2023] [Accepted: 05/16/2024] [Indexed: 05/29/2024] Open
Abstract
BACKGROUND Cancers are spatially heterogenous, thus their clonal evolution, especially following anti-cancer treatments, depends on where the mutated cells are located within the tumor tissue. For example, cells exposed to different concentrations of drugs, such as cells located near the vessels in contrast to those residing far from the vasculature, can undergo a different evolutionary path. However, classical representations of cell lineage trees do not account for this spatial component of emerging cancer clones. Here, we propose routines to trace spatial and temporal clonal evolution in computer simulations of the tumor evolution models. RESULTS The LinG3D (Lineage Graphs in 3D) is an open-source collection of routines (in MATLAB, Python, and R) that enables spatio-temporal visualization of clonal evolution in a two-dimensional tumor slice from computer simulations of the tumor evolution models. These routines draw traces of tumor clones in both time and space, and may include a projection of a selected microenvironmental factor, such as the drug or oxygen distribution within the tumor, if such a microenvironmental factor is used in the tumor evolution model. The utility of LinG3D has been demonstrated through examples of simulated tumors with different number of clones and, additionally, in experimental colony growth assay. CONCLUSIONS This routine package extends the classical lineage trees, that show cellular clone relationships in time, by adding the space component to show the locations of cellular clones within the 2D tumor tissue patch from computer simulations of tumor evolution models.
Collapse
Affiliation(s)
- Anjun Hu
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- Department of Chemistry, College of Arts and Sciences, University of South Florida, Tampa, FL, 33612, USA
| | - Awino Maureiq E Ojwang'
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Kayode D Olumoyin
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Katarzyna A Rejniak
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA.
- Department of Oncologic Sciences, Morsani College of Medicine, University of South Florida, Tampa, FL, 33612, USA.
| |
Collapse
|
3
|
Bishop RT, Miller AK, Froid M, Nerlakanti N, Li T, Frieling JS, Nasr MM, Nyman KJ, Sudalagunta PR, Canevarolo RR, Silva AS, Shain KH, Lynch CC, Basanta D. The bone ecosystem facilitates multiple myeloma relapse and the evolution of heterogeneous drug resistant disease. Nat Commun 2024; 15:2458. [PMID: 38503736 PMCID: PMC10951361 DOI: 10.1038/s41467-024-46594-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2022] [Accepted: 03/04/2024] [Indexed: 03/21/2024] Open
Abstract
Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse.
Collapse
Affiliation(s)
- Ryan T Bishop
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Anna K Miller
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Matthew Froid
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Niveditha Nerlakanti
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Tao Li
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Jeremy S Frieling
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Mostafa M Nasr
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Karl J Nyman
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Praneeth R Sudalagunta
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Rafael R Canevarolo
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Ariosto Siqueira Silva
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Kenneth H Shain
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- Department of Malignant Hematology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Conor C Lynch
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA.
| | - David Basanta
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA.
| |
Collapse
|
4
|
Ojwang' AME, Bazargan S, Johnson JO, Pilon-Thomas S, Rejniak KA. Histology-guided mathematical model of tumor oxygenation: sensitivity analysis of physical and computational parameters. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.05.583363. [PMID: 38496532 PMCID: PMC10942376 DOI: 10.1101/2024.03.05.583363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
A hybrid off-lattice agent-based model has been developed to reconstruct the tumor tissue oxygenation landscape based on histology images and simulated interactions between vasculature and cells with microenvironment metabolites. Here, we performed a robustness sensitivity analysis of that model's physical and computational parameters. We found that changes in the domain boundary conditions, the initial conditions, and the Michaelis constant are negligible and, thus, do not affect the model outputs. The model is also not sensitive to small perturbations of the vascular influx or the maximum consumption rate of oxygen. However, the model is sensitive to large perturbations of these parameters and changes in the tissue boundary condition, emphasizing an imperative aim to measure these parameters experimentally.
Collapse
Affiliation(s)
- Awino Maureiq E Ojwang'
- Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Sarah Bazargan
- Department of Immunology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Joseph O Johnson
- Analytic Microscopy Core, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Shari Pilon-Thomas
- Department of Immunology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Katarzyna A Rejniak
- Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
- Department of Oncologic Sciences, Morsani College of Medicine, University of South Florida, Tampa, FL, USA
| |
Collapse
|
5
|
Hu A, Ojwang' AME, Olumoyin KD, Rejniak KA. Visualizing the Spatio-Temporal Dynamics of Clonal Evolution with LinG3D software. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.05.583631. [PMID: 38496472 PMCID: PMC10942425 DOI: 10.1101/2024.03.05.583631] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
Cancer clonal evolution, especially following anti-cancer treatments, depends on the locations of the mutated cells within the tumor tissue. Cells near the vessels, exposed to higher concentrations of drugs, will undergo a different evolutionary path than cells residing far from the vasculature in the areas of lower drug levels. However, classical representations of cell lineage trees do not account for this spatial component of emerging cancer clones. Here, we propose the LinG3D (Lineage Graphs in 3D) algorithms to trace clonal evolution in space and time. These are an open-source collection of routines (in MATLAB, Python, and R) that enables spatio-temporal visualization of clonal evolution in a two-dimensional tumor slice from computer simulations of the tumor evolution models. These routines draw traces of tumor clones in both time and space, with an option to include a projection of a selected microenvironmental factor, such as the drug or oxygen distribution within the tumor. The utility of LinG3D has been demonstrated through examples of simulated tumors with different number of clones and, additionally, in experimental colony growth assay. This routine package extends the classical lineage trees, that show cellular clone relationships in time, by adding the space component to show the locations of cellular clones within the 2D tumor tissue patch from computer simulations of tumor evolution models.
Collapse
Affiliation(s)
- Anjun Hu
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL 33612, USA
| | - Awino Maureiq E Ojwang'
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL 33612, USA
| | - Kayode D Olumoyin
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL 33612, USA
| | - Katarzyna A Rejniak
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL 33612, USA
| |
Collapse
|
6
|
Gallaher J, Strobl M, West J, Gatenby R, Zhang J, Robertson-Tessi M, Anderson AR. Intermetastatic and Intrametastatic Heterogeneity Shapes Adaptive Therapy Cycling Dynamics. Cancer Res 2023; 83:2775-2789. [PMID: 37205789 PMCID: PMC10425736 DOI: 10.1158/0008-5472.can-22-2558] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 03/11/2023] [Accepted: 05/17/2023] [Indexed: 05/21/2023]
Abstract
Adaptive therapies that alternate between drug applications and drug-free vacations can exploit competition between sensitive and resistant cells to maximize the time to progression. However, optimal dosing schedules depend on the properties of metastases, which are often not directly measurable in clinical practice. Here, we proposed a framework for estimating features of metastases through tumor response dynamics during the first adaptive therapy treatment cycle. Longitudinal prostate-specific antigen (PSA) levels in 16 patients with metastatic castration-resistant prostate cancer undergoing adaptive androgen deprivation treatment were analyzed to investigate relationships between cycle dynamics and clinical variables such as Gleason score, the change in the number of metastases over a cycle, and the total number of cycles over the course of treatment. The first cycle of adaptive therapy, which consists of a response period (applying therapy until 50% PSA reduction), and a regrowth period (removing treatment until reaching initial PSA levels), delineated several features of the computational metastatic system: larger metastases had longer cycles; a higher proportion of drug-resistant cells slowed the cycles; and a faster cell turnover rate sped up drug response time and slowed regrowth time. The number of metastases did not affect cycle times, as response dynamics were dominated by the largest tumors rather than the aggregate. In addition, systems with higher intermetastasis heterogeneity responded better to continuous therapy and correlated with dynamics from patients with high or low Gleason scores. Conversely, systems with higher intrametastasis heterogeneity responded better to adaptive therapy and correlated with dynamics from patients with intermediate Gleason scores. SIGNIFICANCE Multiscale mathematical modeling combined with biomarker dynamics during adaptive therapy helps identify underlying features of metastatic cancer to inform treatment decisions.
Collapse
Affiliation(s)
- Jill Gallaher
- Department of Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida
| | - Maximilian Strobl
- Department of Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida
| | - Jeffrey West
- Department of Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida
| | - Robert Gatenby
- Department of Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida
- Department of Radiology, Moffitt Cancer Center, Tampa, Florida
| | - Jingsong Zhang
- Department of Genitourinary Oncology, Moffitt Cancer Center, Tampa, Florida
| | - Mark Robertson-Tessi
- Department of Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida
| | | |
Collapse
|
7
|
Jørgensen ACS, Hill CS, Sturrock M, Tang W, Karamched SR, Gorup D, Lythgoe MF, Parrinello S, Marguerat S, Shahrezaei V. Data-driven spatio-temporal modelling of glioblastoma. ROYAL SOCIETY OPEN SCIENCE 2023; 10:221444. [PMID: 36968241 PMCID: PMC10031411 DOI: 10.1098/rsos.221444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.
Collapse
Affiliation(s)
| | - Ciaran Scott Hill
- Department of Neurosurgery, The National Hospital for Neurology and Neurosurgery, London WC1N 3BG, UK
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Marc Sturrock
- Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, Dublin D02 YN77, Ireland
| | - Wenhao Tang
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| | - Saketh R. Karamched
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Dunja Gorup
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Mark F. Lythgoe
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Simona Parrinello
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Samuel Marguerat
- Genomics Translational Technology Platform, UCL Cancer Institute, University College London, London WC1E 6DD, UK
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
8
|
Jørgensen ACS, Ghosh A, Sturrock M, Shahrezaei V. Efficient Bayesian inference for stochastic agent-based models. PLoS Comput Biol 2022; 18:e1009508. [PMID: 36197919 PMCID: PMC9576090 DOI: 10.1371/journal.pcbi.1009508] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 10/17/2022] [Accepted: 09/21/2022] [Indexed: 11/14/2022] Open
Abstract
The modelling of many real-world problems relies on computationally heavy simulations of randomly interacting individuals or agents. However, the values of the parameters that underlie the interactions between agents are typically poorly known, and hence they need to be inferred from macroscopic observations of the system. Since statistical inference rests on repeated simulations to sample the parameter space, the high computational expense of these simulations can become a stumbling block. In this paper, we compare two ways to mitigate this issue in a Bayesian setting through the use of machine learning methods: One approach is to construct lightweight surrogate models to substitute the simulations used in inference. Alternatively, one might altogether circumvent the need for Bayesian sampling schemes and directly estimate the posterior distribution. We focus on stochastic simulations that track autonomous agents and present two case studies: tumour growths and the spread of infectious diseases. We demonstrate that good accuracy in inference can be achieved with a relatively small number of simulations, making our machine learning approaches orders of magnitude faster than classical simulation-based methods that rely on sampling the parameter space. However, we find that while some methods generally produce more robust results than others, no algorithm offers a one-size-fits-all solution when attempting to infer model parameters from observations. Instead, one must choose the inference technique with the specific real-world application in mind. The stochastic nature of the considered real-world phenomena poses an additional challenge that can become insurmountable for some approaches. Overall, we find machine learning approaches that create direct inference machines to be promising for real-world applications. We present our findings as general guidelines for modelling practitioners.
Collapse
Affiliation(s)
| | | | - Marc Sturrock
- Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, Dublin, Ireland
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London, United Kingdom
| |
Collapse
|
9
|
Cess CG, Finley SD. Multiscale modeling of tumor adaption and invasion following anti‐angiogenic therapy. COMPUTATIONAL AND SYSTEMS ONCOLOGY 2022. [DOI: 10.1002/cso2.1032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Affiliation(s)
- Colin G. Cess
- Department of Biomedical Engineering University of Southern California Los Angeles California USA
| | - Stacey D. Finley
- Department of Biomedical Engineering University of Southern California Los Angeles California USA
- Department of Quantitative and Computational Biology University of Southern California Los Angeles California USA
- Mork Family Department of Chemical Engineering and Materials Science University of Southern California Los Angeles California USA
| |
Collapse
|
10
|
Kuosmanen T, Cairns J, Noble R, Beerenwinkel N, Mononen T, Mustonen V. Drug-induced resistance evolution necessitates less aggressive treatment. PLoS Comput Biol 2021; 17:e1009418. [PMID: 34555024 PMCID: PMC8491903 DOI: 10.1371/journal.pcbi.1009418] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Revised: 10/05/2021] [Accepted: 09/03/2021] [Indexed: 12/24/2022] Open
Abstract
Increasing body of experimental evidence suggests that anticancer and antimicrobial therapies may themselves promote the acquisition of drug resistance by increasing mutability. The successful control of evolving populations requires that such biological costs of control are identified, quantified and included to the evolutionarily informed treatment protocol. Here we identify, characterise and exploit a trade-off between decreasing the target population size and generating a surplus of treatment-induced rescue mutations. We show that the probability of cure is maximized at an intermediate dosage, below the drug concentration yielding maximal population decay, suggesting that treatment outcomes may in some cases be substantially improved by less aggressive treatment strategies. We also provide a general analytical relationship that implicitly links growth rate, pharmacodynamics and dose-dependent mutation rate to an optimal control law. Our results highlight the important, but often neglected, role of fundamental eco-evolutionary costs of control. These costs can often lead to situations, where decreasing the cumulative drug dosage may be preferable even when the objective of the treatment is elimination, and not containment. Taken together, our results thus add to the ongoing criticism of the standard practice of administering aggressive, high-dose therapies and motivate further experimental and clinical investigation of the mutagenicity and other hidden collateral costs of therapies.
Collapse
Affiliation(s)
- Teemu Kuosmanen
- Organismal and Evolutionary Biology Research Programme, Department of Computer Science, University of Helsinki, Helsinki, Finland
| | - Johannes Cairns
- Organismal and Evolutionary Biology Research Programme, Department of Computer Science, University of Helsinki, Helsinki, Finland
| | - Robert Noble
- Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland
- Present address: Department of Mathematics, City, University of London, London, United Kingdom
| | - Niko Beerenwinkel
- Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland
| | - Tommi Mononen
- Organismal and Evolutionary Biology Research Programme, Department of Computer Science, University of Helsinki, Helsinki, Finland
| | - Ville Mustonen
- Organismal and Evolutionary Biology Research Programme, Department of Computer Science, University of Helsinki, Helsinki, Finland
- Institute of Biotechnology, Helsinki Institute for Information Technology, University of Helsinki, Helsinki, Finland
| |
Collapse
|
11
|
Bridging cell-scale simulations and radiologic images to explain short-time intratumoral oxygen fluctuations. PLoS Comput Biol 2021; 17:e1009206. [PMID: 34310608 PMCID: PMC8341701 DOI: 10.1371/journal.pcbi.1009206] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Revised: 08/05/2021] [Accepted: 06/22/2021] [Indexed: 11/19/2022] Open
Abstract
Radiologic images provide a way to monitor tumor development and its response to therapies in a longitudinal and minimally invasive fashion. However, they operate on a macroscopic scale (average value per voxel) and are not able to capture microscopic scale (cell-level) phenomena. Nevertheless, to examine the causes of frequent fast fluctuations in tissue oxygenation, models simulating individual cells’ behavior are needed. Here, we provide a link between the average data values recorded for radiologic images and the cellular and vascular architecture of the corresponding tissues. Using hybrid agent-based modeling, we generate a set of tissue morphologies capable of reproducing oxygenation levels observed in radiologic images. We then use these in silico tissues to investigate whether oxygen fluctuations can be explained by changes in vascular oxygen supply or by modulations in cellular oxygen absorption. Our studies show that intravascular changes in oxygen supply reproduce the observed fluctuations in tissue oxygenation in all considered regions of interest. However, larger-magnitude fluctuations cannot be recreated by modifications in cellular absorption of oxygen in a biologically feasible manner. Additionally, we develop a procedure to identify plausible tissue morphologies for a given temporal series of average data from radiology images. In future applications, this approach can be used to generate a set of tissues comparable with radiology images and to simulate tumor responses to various anti-cancer treatments at the tissue-scale level. Low levels of oxygen, called hypoxia, are observable in many solid tumors. They are associated with more aggressive malignant cells that are resistant to chemo-, radio-, and immunotherapies. Recently developed imaging techniques provide a way to measure the magnitude of frequent short-term oxygen fluctuations, but they operate on a macro-scale voxel level. To examine the possible causes of rapid oxygen fluctuations at the cell level, we developed a hybrid agent-based mathematical model. We tested two different mechanisms that may be responsible for these cyclic effects on tissue oxygenation: temporal variations in vascular influx of oxygen and modulations in cellular oxygen absorption. Additionally, we developed a procedure to identify plausible tissue morphologies from data collected from radiological images. This can provide a bridge between the micro-scale simulations with individual cells and the longitudinal medical images containing average values. In future applications, this approach can be used to generate a set of tissues compatible with radiology images and to simulate tumor responses to various anticancer treatments at the cell-scale level.
Collapse
|
12
|
Cunningham JJ, Bukkuri A, Brown JS, Gillies RJ, Gatenby RA. Coupled Source-Sink Habitats Produce Spatial and Temporal Variation of Cancer Cell Molecular Properties as an Alternative to Branched Clonal Evolution and Stem Cell Paradigms. Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.676071] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Intratumoral molecular cancer cell heterogeneity is conventionally ascribed to the accumulation of random mutations that occasionally generate fitter phenotypes. This model is built upon the “mutation-selection” paradigm in which mutations drive ever-fitter cancer cells independent of environmental circumstances. An alternative model posits spatio-temporal variation (e.g., blood flow heterogeneity) drives speciation by selecting for cancer cells adapted to each different environment. Here, spatial genetic variation is the consequence rather than the cause of intratumoral evolution. In nature, spatially heterogenous environments are frequently coupled through migration. Drawing from ecological models, we investigate adjacent well-perfused and poorly-perfused tumor regions as “source” and “sink” habitats, respectively. The source habitat has a high carrying capacity resulting in more emigration than immigration. Sink habitats may support a small (“soft-sink”) or no (“hard-sink”) local population. Ecologically, sink habitats can reduce the population size of the source habitat so that, for example, the density of cancer cells directly around blood vessels may be lower than expected. Evolutionarily, sink habitats can exert a selective pressure favoring traits different from those in the source habitat so that, for example, cancer cells adjacent to blood vessels may be suboptimally adapted for that habitat. Soft sinks favor a generalist cancer cell type that moves between the environment but can, under some circumstances, produce speciation events forming source and sink habitat specialists resulting in significant molecular variation in cancer cells separated by small distances. Finally, sink habitats, with limited blood supply, may receive reduced concentrations of systemic drug treatments; and local hypoxia and acidosis may further decrease drug efficacy allowing cells to survive treatment and evolve resistance. In such cases, the sink transforms into the source habitat for resistant cancer cells, leading to treatment failure and tumor progression. We note these dynamics will result in spatial variations in molecular properties as an alternative to the conventional branched evolution model and will result in cellular migration as well as variation in cancer cell phenotype and proliferation currently described by the stem cell paradigm.
Collapse
|
13
|
Suveges S, Chamseddine I, Rejniak KA, Eftimie R, Trucu D. Collective Cell Migration in a Fibrous Environment: A Hybrid Multiscale Modelling Approach. FRONTIERS IN APPLIED MATHEMATICS AND STATISTICS 2021; 7:680029. [PMID: 34322539 PMCID: PMC8315487 DOI: 10.3389/fams.2021.680029] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The specific structure of the extracellular matrix (ECM), and in particular the density and orientation of collagen fibres, plays an important role in the evolution of solid cancers. While many experimental studies discussed the role of ECM in individual and collective cell migration, there are still unanswered questions about the impact of nonlocal cell sensing of other cells on the overall shape of tumour aggregation and its migration type. There are also unanswered questions about the migration and spread of tumour that arises at the boundary between different tissues with different collagen fibre orientations. To address these questions, in this study we develop a hybrid multi-scale model that considers the cells as individual entities and ECM as a continuous field. The numerical simulations obtained through this model match experimental observations, confirming that tumour aggregations are not moving if the ECM fibres are distributed randomly, and they only move when the ECM fibres are highly aligned. Moreover, the stationary tumour aggregations can have circular shapes or irregular shapes (with finger-like protrusions), while the moving tumour aggregations have elongate shapes (resembling to clusters, strands or files). We also show that the cell sensing radius impacts tumour shape only when there is a low ratio of fibre to non-fibre ECM components. Finally, we investigate the impact of different ECM fibre orientations corresponding to different tissues, on the overall tumour invasion of these neighbouring tissues.
Collapse
Affiliation(s)
| | - Ibrahim Chamseddine
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa Florida, USA
| | - Katarzyna A. Rejniak
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa Florida, USA
- Department of Oncologic Sciences, Morsani College of Medicine, University of South Florida, Tampa Florida, USA
| | - Raluca Eftimie
- Laboratoire Mathématiques de Besançon, UMR-CNRS 6623, Université de Bourgogne Franche-Comté, 16 Route de Gray, Besançon, France
| | - Dumitru Trucu
- Department of Mathematics, University of Dundee, Dundee, UK
| |
Collapse
|
14
|
Wu Z, Lawrence PJ, Ma A, Zhu J, Xu D, Ma Q. Single-Cell Techniques and Deep Learning in Predicting Drug Response. Trends Pharmacol Sci 2020; 41:1050-1065. [PMID: 33153777 PMCID: PMC7669610 DOI: 10.1016/j.tips.2020.10.004] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 10/04/2020] [Accepted: 10/09/2020] [Indexed: 12/19/2022]
Abstract
Rapidly developing single-cell sequencing analyses produce more comprehensive profiles of the genomic, transcriptomic, and epigenomic heterogeneity of tumor subpopulations than do traditional bulk sequencing analyses. Moreover, single-cell techniques allow the response of a tumor to drug exposure to be more thoroughlyinvestigated. Deep learning (DL) models have successfully extracted features from complex bulk sequence data to predict drug responses. We review recent innovations in single-cell technologies and DL-based approaches related to drug sensitivity predictions. We believe that, by using insights from bulk sequencedata, deep transfer learning (DTL) can facilitate the use of single-cell data for training superior DL-based drug prediction models.
Collapse
Affiliation(s)
- Zhenyu Wu
- Department of Biomedical Informatics, The Ohio State University, Columbus, OH 43210, USA
| | - Patrick J Lawrence
- Department of Biomedical Informatics, The Ohio State University, Columbus, OH 43210, USA
| | - Anjun Ma
- Department of Biomedical Informatics, The Ohio State University, Columbus, OH 43210, USA
| | - Jian Zhu
- Department of Pathology, The Ohio State University, Columbus, OH 43210, USA
| | - Dong Xu
- Department of Electrical Engineering and Computer Science, and Christopher S. Bond Life Sciences Center, University of Missouri, Columbia, MO 65211, USA
| | - Qin Ma
- Department of Biomedical Informatics, The Ohio State University, Columbus, OH 43210, USA.
| |
Collapse
|
15
|
The CCL5/CCR5 Axis in Cancer Progression. Cancers (Basel) 2020; 12:cancers12071765. [PMID: 32630699 PMCID: PMC7407580 DOI: 10.3390/cancers12071765] [Citation(s) in RCA: 188] [Impact Index Per Article: 47.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 06/25/2020] [Accepted: 06/30/2020] [Indexed: 02/07/2023] Open
Abstract
Tumor cells can “hijack” chemokine networks to support tumor progression. In this context, the C-C chemokine ligand 5/C-C chemokine receptor type 5 (CCL5/CCR5) axis is gaining increasing attention, since abnormal expression and activity of CCL5 and its receptor CCR5 have been found in hematological malignancies and solid tumors. Numerous preclinical in vitro and in vivo studies have shown a key role of the CCL5/CCR5 axis in cancer, and thus provided the rationale for clinical trials using the repurposed drug maraviroc, a CCR5 antagonist used to treat HIV/AIDS. This review summarizes current knowledge on the role of the CCL5/CCR5 axis in cancer. First, it describes the involvement of the CCL5/CCR5 axis in cancer progression, including autocrine and paracrine tumor growth, ECM (extracellular matrix) remodeling and migration, cancer stem cell expansion, DNA damage repair, metabolic reprogramming, and angiogenesis. Then, it focuses on individual hematological and solid tumors in which CCL5 and CCR5 have been studied preclinically. Finally, it discusses clinical trials of strategies to counteract the CCL5/CCR5 axis in different cancers using maraviroc or therapeutic monoclonal antibodies.
Collapse
|
16
|
Berrouet C, Dorilas N, Rejniak KA, Tuncer N. Comparison of Drug Inhibitory Effects ([Formula: see text]) in Monolayer and Spheroid Cultures. Bull Math Biol 2020; 82:68. [PMID: 32495209 PMCID: PMC9773863 DOI: 10.1007/s11538-020-00746-7] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Accepted: 05/06/2020] [Indexed: 12/25/2022]
Abstract
Traditionally, the monolayer (two-dimensional) cell cultures are used for initial evaluation of the effectiveness of anticancer drugs. In particular, these experiments provide the [Formula: see text] curves that determine drug concentration that can inhibit growth of a tumor colony by half when compared to the cells grown with no exposure to the drug. Low [Formula: see text] value means that the drug is effective at low concentrations, and thus will show lower systemic toxicity when administered to the patient. However, in these experiments cells are grown in a monolayer, all well exposed to the drug, while in vivo tumors expand as three-dimensional multicellular masses, where inner cells have a limited access to the drug. Therefore, we performed computational studies to compare the [Formula: see text] curves for cells grown as a two-dimensional monolayer and a cross section through a three-dimensional spheroid. Our results identified conditions (drug diffusivity, drug action mechanisms and cell proliferation capabilities) under which these [Formula: see text] curves differ significantly. This will help experimentalists to better determine drug dosage for future in vivo experiments and clinical trials.
Collapse
Affiliation(s)
- Catherine Berrouet
- Department of Mathematical Sciences, Charles E. Schmidt College of Science, Florida Atlantic University, Boca Raton, FL, USA
| | - Naika Dorilas
- Department of Mathematical Sciences, Charles E. Schmidt College of Science, Florida Atlantic University, Boca Raton, FL, USA
| | - Katarzyna A. Rejniak
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
- Department of Oncologic Sciences, Morsani College of Medicine, University of South Florida, Tampa, FL, USA
| | - Necibe Tuncer
- Department of Mathematical Sciences, Charles E. Schmidt College of Science, Florida Atlantic University, Boca Raton, FL, USA
| |
Collapse
|