Interplay of Darwinian Selection, Lamarckian Induction and Microvesicle Transfer on Drug Resistance in Cancer

Metabolic reprogramming of poly(morpho)nuclear giant cells determines glioblastoma recovery from doxorubicin-induced stress

BACKGROUND

Multi-drug resistance of poly(morpho)nuclear giant cells (PGCs) determines their cytoprotective and generative potential in cancer ecosystems. However, mechanisms underlying the involvement of PGCs in glioblastoma multiforme (GBM) adaptation to chemotherapeutic regimes remain largely obscure. In particular, metabolic reprogramming of PGCs has not yet been considered in terms of GBM recovery from doxorubicin (DOX)-induced stress.

METHODS

Long-term proteomic and metabolic cell profiling was applied to trace the phenotypic dynamics of GBM populations subjected to pulse DOX treatment in vitro, with a particular focus on PGC formation and its metabolic background. The links between metabolic reprogramming, drug resistance and drug retention capacity of PGCs were assessed, along with their significance for GBM recovery from DOX-induced stress.

RESULTS

Pulse DOX treatment triggered the transient formation of PGCs, followed by the appearance of small expanding cell (SEC) clusters. Development of PGCs was accompanied by the mobilization of their metabolic proteome, transient induction of oxidative phosphorylation (OXPHOS), and differential intracellular accumulation of NADH, NADPH, and ATP. The metabolic background of PGC formation was confirmed by the attenuation of GBM recovery from DOX-induced stress following the chemical inhibition of GSK-3β, OXPHOS, and the pentose phosphate pathway. Concurrently, the mobilization of reactive oxygen species (ROS) scavenging systems and fine-tuning of NADPH-dependent ROS production systems in PGCs was observed. These processes were accompanied by perinuclear mobilization of ABCB1 and ABCG2 transporters and DOX retention in the perinuclear PGC compartments.

CONCLUSIONS

These data demonstrate the cooperative pattern of GBM recovery from DOX-induced stress and the crucial role of metabolic reprogramming of PGCs in this process. Metabolic reprogramming enhances the efficiency of self-defense systems and increases the DOX retention capacity of PGCs, potentially reducing DOX bioavailability in the proximity of SECs. Consequently, the modulation of PGC metabolism is highlighted as a potential target for intervention in glioblastoma treatment.

Impact of genetic variants in the solute carrier (SLC) genes encoding drug uptake transporters on the response to anticancer chemotherapy

Cancer drug resistance constitutes a severe limitation for the satisfactory outcome of these patients. This is a complex problem due to the co-existence in cancer cells of multiple and synergistic mechanisms of chemoresistance (MOC). These mechanisms are accounted for by the expression of a set of genes included in the so-called resistome, whose effectiveness often leads to a lack of response to pharmacological treatment. Additionally, genetic variants affecting these genes further increase the complexity of the question. This review focuses on a set of genes encoding members of the transportome involved in drug uptake, which have been classified into the MOC-1A subgroup of the resistome. These proteins belong to the solute carrier (SLC) superfamily. More precisely, we have considered here several members of families SLC2, SLC7, SLC19, SLC22, SLCO, SLC28, SLC29, SLC31, SLC46, and SLC47 due to the impact of their expression and genetic variants in anticancer drug uptake by tumor cells or, in some cases, general bioavailability. Changes in their expression levels and the appearance of genetic variants can contribute to the Darwinian selection of more resistant clones and, hence, to the development of a more malignant phenotype. Accordingly, to address this issue in future personalized medicine, it is necessary to characterize both changes in resistome genes that can affect their function. It is also essential to consider the time-dependent dimension of these features, as the genetic expression and the appearance of genetic variants can change during tumor progression and in response to treatment.

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

Acquired resistance to KRAS G12C small-molecule inhibitors via genetic/nongenetic mechanisms in lung cancer

Inherent or acquired resistance to sotorasib poses a substantialt challenge for NSCLC treatment. Here, we demonstrate that acquired resistance to sotorasib in isogenic cells correlated with increased expression of integrin β4 (ITGB4), a component of the focal adhesion complex. Silencing ITGB4 in tolerant cells improved sotorasib sensitivity, while overexpressing ITGB4 enhanced tolerance to sotorasib by supporting AKT-mTOR bypass signaling. Chronic treatment with sotorasib induced WNT expression and activated the WNT/β-catenin signaling pathway. Thus, silencing both ITGB4 and β-catenin significantly improved sotorasib sensitivity in tolerant, acquired, and inherently resistant cells. In addition, the proteasome inhibitor carfilzomib (CFZ) exhibited synergism with sotorasib by down-regulating ITGB4 and β-catenin expression. Furthermore, adagrasib phenocopies the combination effect of sotorasib and CFZ by suppressing KRAS activity and inhibiting cell cycle progression in inherently resistant cells. Overall, our findings unveil previously unrecognized nongenetic mechanisms underlying resistance to sotorasib and propose a promising treatment strategy to overcome resistance.

Stochastic Fluctuations Drive Non-genetic Evolution of Proliferation in Clonal Cancer Cell Populations

Evolutionary dynamics allows us to understand many changes happening in a broad variety of biological systems, ranging from individuals to complete ecosystems. It is also behind a number of remarkable organizational changes that happen during the natural history of cancers. These reflect tumour heterogeneity, which is present at all cellular levels, including the genome, proteome and phenome, shaping its development and interrelation with its environment. An intriguing observation in different cohorts of oncological patients is that tumours exhibit an increased proliferation as the disease progresses, while the timescales involved are apparently too short for the fixation of sufficient driver mutations to promote explosive growth. Here, we discuss how phenotypic plasticity, emerging from a single genotype, may play a key role and provide a ground for a continuous acceleration of the proliferation rate of clonal populations with time. We address this question by combining the analysis of real-time growth of non-small-cell lung carcinoma cells (N-H460) together with stochastic and deterministic mathematical models that capture proliferation trait heterogeneity in clonal populations to elucidate the contribution of phenotypic transitions on tumour growth dynamics.

Single-Cell DNA Methylation Analysis in Cancer

Morphological, transcriptomic, and genomic defects are well-explored parameters of cancer biology. In more recent years, the impact of epigenetic influences, such as DNA methylation, is becoming more appreciated. Aberrant DNA methylation has been implicated in many types of cancers, influencing cell type, state, transcriptional regulation, and genomic stability to name a few. Traditionally, large populations of cells from the tissue of interest are coalesced for analysis, producing averaged methylome data. Considering the inherent heterogeneity of cancer, analysing populations of cells as a whole denies the ability to discover novel aberrant methylation patterns, identify subpopulations, and trace cell lineages. Due to recent advancements in technology, it is now possible to obtain methylome data from single cells. This has both research and clinical implications, ranging from the identification of biomarkers to improved diagnostic tools. As with all emerging technologies, distinct experimental, bioinformatic, and practical challenges present themselves. This review begins with exploring the potential impact of single-cell sequencing on understanding cancer biology and how it could eventually benefit a clinical setting. Following this, the techniques and experimental approaches which made this technology possible are explored. Finally, the present challenges currently associated with single-cell DNA methylation sequencing are described.

An Epigenetic LINE-1-Based Mechanism in Cancer

In the last fifty years, large efforts have been deployed in basic research, clinical oncology, and clinical trials, yielding an enormous amount of information regarding the molecular mechanisms of cancer and the design of effective therapies. The knowledge that has accumulated underpins the complexity, multifactoriality, and heterogeneity of cancer, disclosing novel landscapes in cancer biology with a key role of genome plasticity. Here, we propose that cancer onset and progression are determined by a stress-responsive epigenetic mechanism, resulting from the convergence of upregulation of LINE-1 (long interspersed nuclear element 1), the largest family of human retrotransposons, genome damage, nuclear lamina fragmentation, chromatin remodeling, genome reprogramming, and autophagy activation. The upregulated expression of LINE-1 retrotransposons and their protein products plays a key role in these processes, yielding an increased plasticity of the nuclear architecture with the ensuing reprogramming of global gene expression, including the reactivation of embryonic transcription profiles. Cancer phenotypes would thus emerge as a consequence of the unscheduled reactivation of embryonic gene expression patterns in an inappropriate context, triggering de-differentiation and aberrant proliferation in differentiated cells. Depending on the intensity of the stressing stimuli and the level of LINE-1 response, diverse degrees of malignity would be generated.

Mathematical characterization of population dynamics in breast cancer cells treated with doxorubicin

The development of chemoresistance remains a significant cause of treatment failure in breast cancer. We posit that a mathematical understanding of chemoresistance could assist in developing successful treatment strategies. Towards that end, we have developed a model that describes the cytotoxic effects of the standard chemotherapeutic drug doxorubicin on the MCF-7 breast cancer cell line. We assume that treatment with doxorubicin induces a compartmentalization of the breast cancer cell population into surviving cells, which continue proliferating after treatment, and irreversibly damaged cells, which gradually transition from proliferating to treatment-induced death. The model is fit to experimental data including variations in drug concentration, inter-treatment interval, and number of doses. Our model recapitulates tumor cell dynamics in all these scenarios (as quantified by the concordance correlation coefficient, CCC > 0.95). In particular, superior tumor control is observed with higher doxorubicin concentrations, shorter inter-treatment intervals, and a higher number of doses (p < 0.05). Longer inter-treatment intervals require adapting the model parameterization after each doxorubicin dose, suggesting the promotion of chemoresistance. Additionally, we propose promising empirical formulas to describe the variation of model parameters as functions of doxorubicin concentration (CCC > 0.78). Thus, we conclude that our mathematical model could deepen our understanding of the cytotoxic effects of doxorubicin and could be used to explore practical drug regimens achieving optimal tumor control.

In Vitro Characterization of Renal Drug Transporter Activity in Kidney Cancer

The activity of drug transporters is central to the secretory function of the kidneys and a defining feature of renal proximal tubule epithelial cells (RPTECs). The expression, regulation, and function of these membrane-bound proteins is well understood under normal renal physiological conditions. However, the impact of drug transporters on the pathophysiology of kidney cancer is still elusive. In the present study, we employed different renal cell carcinoma (RCC) cell lines and a prototypical non-malignant RPTEC cell line to characterize the activity, expression, and potential regulatory mechanisms of relevant renal drug transporters in RCC in vitro. An analysis of the uptake and efflux activity, the expression of drug transporters, and the evaluation of cisplatin cytotoxicity under the effects of methylation or epidermal growth factor receptor (EGFR) inhibition showed that the RCC cells retained substantial drug transport activity. In RCC cells, P-glycoprotein was localized in the nucleus and its pharmacological inhibition enhanced cisplatin toxicity in non-malignant RPTECs. On the other hand, methylation inhibition enhanced cisplatin toxicity by upregulating the organic cation uptake activity in RCC cells. Differential effects of methylation and EGFR were observed in transporter expression, showing regulatory heterogeneity in these cells. Interestingly, the non-malignant RPTEC cell line that was used lacked the machinery responsible for organic cation transport, which reiterates the functional losses that renal cells undergo in vitro.

Dynamic Phenotypic Switching and Group Behavior Help Non-Small Cell Lung Cancer Cells Evade Chemotherapy

Drug resistance, a major challenge in cancer therapy, is typically attributed to mutations and genetic heterogeneity. Emerging evidence suggests that dynamic cellular interactions and group behavior also contribute to drug resistance. However, the underlying mechanisms remain poorly understood. Here, we present a new mathematical approach with game theoretical underpinnings that we developed to model real-time growth data of non-small cell lung cancer (NSCLC) cells and discern patterns in response to treatment with cisplatin. We show that the cisplatin-sensitive and cisplatin-tolerant NSCLC cells, when co-cultured in the absence or presence of the drug, display dynamic group behavior strategies. Tolerant cells exhibit a 'persister-like' behavior and are attenuated by sensitive cells; they also appear to 'educate' sensitive cells to evade chemotherapy. Further, tolerant cells can switch phenotypes to become sensitive, especially at low cisplatin concentrations. Finally, switching treatment from continuous to an intermittent regimen can attenuate the emergence of tolerant cells, suggesting that intermittent chemotherapy may improve outcomes in lung cancer.

Global analysis of a cancer model with drug resistance due to Lamarckian induction and microvesicle transfer

Development of resistance to chemotherapy in cancer patients strongly effects the outcome of the treatment. Due to chemotherapeutic agents, resistance can emerge by Darwinian evolution. Besides this, acquired drug resistance may arise via changes in gene expression. A recent discovery in cancer research uncovered a third possibility, indicating that this phenotype conversion can occur through the transfer of microvesicles from resistant to sensitive cells, a mechanism resembling the spread of an infectious agent. We present a model describing the evolution of sensitive and resistant tumour cells considering Darwinian selection, Lamarckian induction and microvesicle transfer. We identify three threshold parameters which determine the existence and stability of the three possible equilibria. Using a simple Dulac function, we give a complete description of the dynamics of the model depending on the three threshold parameters. We also establish an agent based model as a spatial version of the ODE model and compare the outputs of the two models. We find that although the ODE model does not provide spatial information about the structure of the tumour, it is capable to determine the outcome in terms of tumour size and distribution of cell types. We demonstrate the possible effects of increasing drug concentration, and characterize the possible bifurcation sequences. Our results show that the presence of microvesicle transfer cannot ruin a therapy that otherwise leads to extinction, however it may doom a partially successful therapy to failure.

Multi-omics data integration considerations and study design for biological systems and disease

With the advancement of next-generation sequencing and mass spectrometry, there is a growing need for the ability to merge biological features in order to study a system as a whole. Features such as the transcriptome, methylome, proteome, histone post-translational modifications and the microbiome all influence the host response to various diseases and cancers. Each of these platforms have technological limitations due to sample preparation steps, amount of material needed for sequencing, and sequencing depth requirements. These features provide a snapshot of one level of regulation in a system. The obvious next step is to integrate this information and learn how genes, proteins, and/or epigenetic factors influence the phenotype of a disease in context of the system. In recent years, there has been a push for the development of data integration methods. Each method specifically integrates a subset of omics data using approaches such as conceptual integration, statistical integration, model-based integration, networks, and pathway data integration. In this review, we discuss considerations of the study design for each data feature, the limitations in gene and protein abundance and their rate of expression, the current data integration methods, and microbiome influences on gene and protein expression. The considerations discussed in this review should be regarded when developing new algorithms for integrating multi-omics data.

Single-cell RNA sequencing in cancer research

Single-cell RNA sequencing (scRNA-seq), a technology that analyzes transcriptomes of complex tissues at single-cell levels, can identify differential gene expression and epigenetic factors caused by mutations in unicellular genomes, as well as new cell-specific markers and cell types. scRNA-seq plays an important role in various aspects of tumor research. It reveals the heterogeneity of tumor cells and monitors the progress of tumor development, thereby preventing further cellular deterioration. Furthermore, the transcriptome analysis of immune cells in tumor tissue can be used to classify immune cells, their immune escape mechanisms and drug resistance mechanisms, and to develop effective clinical targeted therapies combined with immunotherapy. Moreover, this method enables the study of intercellular communication and the interaction of tumor cells and non-malignant cells to reveal their role in carcinogenesis. scRNA-seq provides new technical means for further development of tumor research and is expected to make significant breakthroughs in this field. This review focuses on the principles of scRNA-seq, with an emphasis on the application of scRNA-seq in tumor heterogeneity, pathogenesis, and treatment.

A mesoscopic simulator to uncover heterogeneity and evolutionary dynamics in tumors

Increasingly complex in silico modeling approaches offer a way to simultaneously access cancerous processes at different spatio-temporal scales. High-level models, such as those based on partial differential equations, are computationally affordable and allow large tumor sizes and long temporal windows to be studied, but miss the discrete nature of many key underlying cellular processes. Individual-based approaches provide a much more detailed description of tumors, but have difficulties when trying to handle full-sized real cancers. Thus, there exists a trade-off between the integration of macroscopic and microscopic information, now widely available, and the ability to attain clinical tumor sizes. In this paper we put forward a stochastic mesoscopic simulation framework that incorporates key cellular processes during tumor progression while keeping computational costs to a minimum. Our framework captures a physical scale that allows both the incorporation of microscopic information, tracking the spatio-temporal emergence of tumor heterogeneity and the underlying evolutionary dynamics, and the reconstruction of clinically sized tumors from high-resolution medical imaging data, with the additional benefit of low computational cost. We illustrate the functionality of our modeling approach for the case of glioblastoma, a paradigm of tumor heterogeneity that remains extremely challenging in the clinical setting.

Evolutionary dynamics at the tumor edge reveal metabolic imaging biomarkers

Human cancers are biologically and morphologically heterogeneous. A variety of clonal populations emerge within these neoplasms and their interaction leads to complex spatiotemporal dynamics during tumor growth. We studied the reshaping of metabolic activity in human cancers by means of continuous and discrete mathematical models and matched the results to positron emission tomography (PET) imaging data. Our models revealed that the location of increasingly active proliferative cellular spots progressively drifted from the center of the tumor to the periphery, as a result of the competition between gradually more aggressive phenotypes. This computational finding led to the development of a metric, normalized distance from 18F-fluorodeoxyglucose (18F-FDG) hotspot to centroid (NHOC), based on the separation from the location of the activity (proliferation) hotspot to the tumor centroid. The NHOC metric can be computed for patients using 18F-FDG PET-computed tomography (PET/CT) images where the voxel of maximum uptake (standardized uptake value [SUV]max) is taken as the activity hotspot. Two datasets of 18F-FDG PET/CT images were collected, one from 61 breast cancer patients and another from 161 non-small-cell lung cancer patients. In both cohorts, survival analyses were carried out for the NHOC and for other classical PET/CT-based biomarkers, finding that the former had a high prognostic value, outperforming the latter. In summary, our work offers additional insights into the evolutionary mechanisms behind tumor progression, provides a different PET/CT-based biomarker, and reveals that an activity hotspot closer to the tumor periphery is associated to a worst patient outcome.

Engineering in Medicine To Address the Challenge of Cancer Drug Resistance: From Micro- and Nanotechnologies to Computational and Mathematical Modeling

Drug resistance has profoundly limited the success of cancer treatment, driving relapse, metastasis, and mortality. Nearly all anticancer drugs and even novel immunotherapies, which recalibrate the immune system for tumor recognition and destruction, have succumbed to resistance development. Engineers have emerged across mechanical, physical, chemical, mathematical, and biological disciplines to address the challenge of drug resistance using a combination of interdisciplinary tools and skill sets. This review explores the developing, complex, and under-recognized role of engineering in medicine to address the multitude of challenges in cancer drug resistance. Looking through the "lens" of intrinsic, extrinsic, and drug-induced resistance (also referred to as "tolerance"), we will discuss three specific areas where active innovation is driving novel treatment paradigms: (1) nanotechnology, which has revolutionized drug delivery in desmoplastic tissues, harnessing physiochemical characteristics to destroy tumors through photothermal therapy and rationally designed nanostructures to circumvent cancer immunotherapy failures, (2) bioengineered tumor models, which have benefitted from microfluidics and mechanical engineering, creating a paradigm shift in physiologically relevant environments to predict clinical refractoriness and enabling platforms for screening drug combinations to thwart resistance at the individual patient level, and (3) computational and mathematical modeling, which blends in silico simulations with molecular and evolutionary principles to map mutational patterns and model interactions between cells that promote resistance. On the basis that engineering in medicine has resulted in discoveries in resistance biology and successfully translated to clinical strategies that improve outcomes, we suggest the proliferation of multidisciplinary science that embraces engineering.

A Non-genetic Mechanism Involving the Integrin β4/Paxillin Axis Contributes to Chemoresistance in Lung Cancer

Tumor heterogeneity and cisplatin resistance are major causes of tumor relapse and poor survival. Here, we show that in lung cancer, interaction between paxillin (PXN) and integrin β4 (ITGB4), components of the focal adhesion (FA) complex, contributes to cisplatin resistance. Knocking down PXN and ITGB4 attenuated cell growth and improved cisplatin sensitivity, both in 2D and 3D cultures. PXN and ITGB4 independently regulated expression of several genes. In addition, they also regulated expression of common genes including USP1 and VDAC1, which are required for maintaining genomic stability and mitochondrial function, respectively. Mathematical modeling suggested that bistability could lead to stochastic phenotypic switching between cisplatin-sensitive and resistant states in these cells. Consistently, purified subpopulations of sensitive and resistant cells re-created the mixed parental population when cultured separately. Altogether, these data point to an unexpected role of the FA complex in cisplatin resistance and highlight a novel non-genetic mechanism.

Identification of a transient state during the acquisition of temozolomide resistance in glioblastoma

Drug resistance limits the therapeutic efficacy in cancers and leads to tumor recurrence through ill-defined mechanisms. Glioblastoma (GBM) are the deadliest brain tumors in adults. GBM, at diagnosis or after treatment, are resistant to temozolomide (TMZ), the standard chemotherapy. To better understand the acquisition of this resistance, we performed a longitudinal study, using a combination of mathematical models, RNA sequencing, single cell analyses, functional and drug assays in a human glioma cell line (U251). After an initial response characterized by cell death induction, cells entered a transient state defined by slow growth, a distinct morphology and a shift of metabolism. Specific genes expression associated to this population revealed chromatin remodeling. Indeed, the histone deacetylase inhibitor trichostatin (TSA), specifically eliminated this population and thus prevented the appearance of fast growing TMZ-resistant cells. In conclusion, we have identified in glioblastoma a population with tolerant-like features, which could constitute a therapeutic target.

Radiation protraction schedules for low-grade gliomas: a comparison between different mathematical models

We optimize radiotherapy (RT) administration strategies for treating low-grade gliomas. Specifically, we consider different tumour growth laws, both with and without spatial effects. In each scenario, we find the optimal treatment in the sense of maximizing the overall survival time of a virtual low-grade glioma patient, whose tumour progresses according to the examined growth laws. We discover that an extreme protraction therapeutic strategy, which amounts to substantially extending the time interval between RT sessions, may lead to better tumour control. The clinical implications of our results are also presented.

Targeting Oncogenic BRAF: Past, Present, and Future

Identifying recurrent somatic genetic alterations of, and dependency on, the kinase BRAF has enabled a "precision medicine" paradigm to diagnose and treat BRAF-driven tumors. Although targeted kinase inhibitors against BRAF are effective in a subset of mutant BRAF tumors, resistance to the therapy inevitably emerges. In this review, we discuss BRAF biology, both in wild-type and mutant settings. We discuss the predominant BRAF mutations and we outline therapeutic strategies to block mutant BRAF and cancer growth. We highlight common mechanistic themes that underpin different classes of resistance mechanisms against BRAF-targeted therapies and discuss tumor heterogeneity and co-occurring molecular alterations as a potential source of therapy resistance. We outline promising therapy approaches to overcome these barriers to the long-term control of BRAF-driven tumors and emphasize how an extensive understanding of these themes can offer more pre-emptive, improved therapeutic strategies.