Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells.

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity.

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Systematic evolution of initiation factor 3 and the ribosomal protein uS12 optimizes Escherichia coli growth with an unconventional initiator tRNA.

The anticodon stem of initiator tRNA (i-tRNA) possesses the characteristic three consecutive GC base pairs (G29:C41, G30:C40, and G31:C39 abbreviated as GC/GC/GC or 3GC pairs) crucial to commencing translation. To understand the importance of this highly conserved element, we isolated two fast-growing suppressors of Escherichia coli sustained solely on an unconventional i-tRNA (i-tRNAcg/GC/cg ) having cg/GC/cg sequence instead of the conventional GC/GC/GC. Both suppressors have the common mutation of V93A in initiation factor 3 (IF3), and additional mutations of either V32L (Sup-1) or H76L (Sup-2) in small subunit ribosomal protein 12 (uS12). The V93A mutation in IF3 was necessary for relaxed fidelity of i-tRNA selection to sustain on i-tRNAcg/GC/cg though with a retarded growth. Subsequent mutations in uS12 salvaged the retarded growth by enhancing the fidelity of translation. The H76L mutation in uS12 showed better fidelity of i-tRNA selection. However, the V32L mutation compensated for the deficient fidelity of i-tRNA selection by ensuring an efficient fidelity check by ribosome recycling factor (RRF). We reveal unique genetic networks between uS12, IF3 and i-tRNA in initiation and between uS12, elongation factor-G (EF-G), RRF, and Pth (peptidyl-tRNA hydrolase) which, taken together, govern the fidelity of translation in bacteria.

A mutation in the ribosomal protein uS12 reveals novel functions of its universally conserved PNSA loop.

The ribosomal protein uS12 is conserved across all domains of life. Recently, a heterozygous spontaneous mutation in human uS12 (corresponding to R49K mutation immediately downstream of the universally conserved 44 PNSA47 loop in Escherichia coli uS12) was identified for causing ribosomopathy, highlighting the importance of the PNSA loop. To investigate the effects of a similar mutation in the absence of any wild-type alleles, we mutated the rpsL gene (encoding uS12) in E. coli. Consistent with its pathology (in humans), we were unable to generate the R49K mutation in E. coli in the absence of a support plasmid. However, we were able to generate the L48K mutation in its immediate vicinity. The L48K mutation resulted in a cold sensitive phenotype and ribosome biogenesis defect in the strain. We show that the L48K mutation impacts the steps of initiation and elongation. Furthermore, the genetic interactions of the L48K mutation with RRF and Pth suggest a novel role of the PNSA loop in ribosome recycling. Our studies reveal new functions of the PNSA loop in uS12, which has so far been studied in the context of translation elongation.

Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools.

Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.

Emergent dynamics of underlying regulatory network links EMT and androgen receptor-dependent resistance in prostate cancer.

Advanced prostate cancer patients initially respond to hormone therapy, be it in the form of androgen deprivation therapy or second-generation hormone therapies, such as abiraterone acetate or enzalutamide. However, most men with prostate cancer eventually develop hormone therapy resistance. This resistance can arise through multiple mechanisms, such as through genetic mutations, epigenetic mechanisms, or through non-genetic pathways, such as lineage plasticity along epithelial-mesenchymal or neuroendocrine-like axes. These mechanisms of hormone therapy resistance often co-exist within a single patient's tumor and can overlap within a single cell. There exists a growing need to better understand how phenotypic heterogeneity and plasticity results from emergent dynamics of the regulatory networks governing androgen independence. Here, we investigated the dynamics of a regulatory network connecting the drivers of androgen receptor (AR) splice variant-mediated androgen independence and those of epithelial-mesenchymal transition. Model simulations for this network revealed four possible phenotypes: epithelial-sensitive (ES), epithelial-resistant (ER), mesenchymal-resistant (MR) and mesenchymal-sensitive (MS), with the latter phenotype occurring rarely. We observed that well-coordinated "teams" of regulators working antagonistically within the network enable these phenotypes. These model predictions are supported by multiple transcriptomic datasets both at single-cell and bulk levels, including in vitro EMT induction models and clinical samples. Further, our simulations reveal spontaneous stochastic switching between the ES and MR states. Addition of the immune checkpoint molecule, PD-L1, to the network was able to capture the interactions between AR, PD-L1, and the mesenchymal marker SNAIL, which was also confirmed through quantitative experiments. This systems-level understanding of the driver of androgen independence and EMT could aid in understanding non-genetic transitions and progression of such cancers and help in identifying novel therapeutic strategies or targets.

An androgen receptor regulated gene score is associated with epithelial to mesenchymal transition features in triple negative breast cancers.

BACKGROUND: Androgen receptor (AR) is considered a marker of better prognosis in hormone receptor positive breast cancers (BC), however, its role in triple negative breast cancer (TNBC) is controversial. This may be attributed to intrinsic molecular differences or scoring methods for AR positivity. We derived AR regulated gene score and examined its utility in BC subtypes. METHODS: AR regulated genes were derived by applying a bioinformatic pipeline on publicly available microarray data sets of AR+ BC cell lines and gene score was calculated as average expression of six AR regulated genes. Tumors were divided into AR high and low based on gene score and associations with clinical parameters, circulating androgens, survival and epithelial to mesenchymal transition (EMT) markers were examined, further evaluated in invitro models and public datasets. RESULTS: 53% (133/249) tumors were classified as AR gene score high and were associated with significantly better clinical parameters, disease-free survival (86.13 vs 72.69 months, log rank p = 0.032) when compared to AR low tumors. 36% of TNBC (N = 66) were AR gene score high with higher expression of EMT markers (p = 0.024) and had high intratumoral levels of 5α-reductase, enzyme involved in intracrine androgen metabolism. In MDA-MB-453 treated with dihydrotestosterone, SLUG expression increased, E-cadherin decreased with increase in migration and these changes were reversed with bicalutamide. Similar results were obtained in public datasets. CONCLUSION: Deciphering the role of AR in BC is difficult based on AR protein levels alone. Our results support the context dependent function of AR in driving better prognosis in ER positive tumors and EMT features in TNBC tumors.

Multi-modal transcriptomic analysis unravels enrichment of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer.

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.

Quantitative landscapes reveal trajectories of cell-state transitions associated with drug resistance in melanoma.

Drug resistance and tumor relapse in patients with melanoma is attributed to a combination of genetic and non-genetic mechanisms. Dedifferentiation, a common mechanism of non-genetic resistance in melanoma is characterized by the loss of melanocytic markers. While various molecular attributes of de-differentiation have been identified, the transition dynamics remain poorly understood. Here, we construct cell-state transition landscapes, to quantify the stochastic dynamics driving phenotypic switching in melanoma based on its underlying regulatory network. These landscapes reveal the existence of multiple alternative paths to resistance-de-differentiation and transition to a hyper-pigmented phenotype. Finally, by visualizing the changes in the landscape during in silico molecular perturbations, we identify combinatorial strategies that can lead to the most optimal outcome-a landscape with the minimum occupancy of the two drug-resistant states. Therefore, we present these landscapes as platforms to screen possible therapeutic interventions in terms of their ability to lead to the most favorable patient outcomes.

Mapping phenotypic heterogeneity in melanoma onto the epithelial-hybrid-mesenchymal axis.

Epithelial to mesenchymal transition (EMT) is a well-studied hallmark of epithelial-like cancers that is characterized by loss of epithelial markers and gain of mesenchymal markers. Melanoma, which is derived from melanocytes of the skin, also undergo phenotypic plasticity toward mesenchymal-like phenotypes under the influence of various micro-environmental cues. Our study connects EMT to the phenomenon of de-differentiation (i.e., transition from proliferative to more invasive phenotypes) observed in melanoma cells during drug treatment. By analyzing 78 publicly available transcriptomic melanoma datasets, we found that de-differentiation in melanoma is accompanied by upregulation of mesenchymal genes, but not necessarily a concomitant loss of an epithelial program, suggesting a more "one-dimensional" EMT that leads to a hybrid epithelial/mesenchymal phenotype. Samples lying in the hybrid epithelial/mesenchymal phenotype also correspond to the intermediate phenotypes in melanoma along the proliferative-invasive axis - neural crest and transitory ones. As melanoma cells progress along the invasive axis, the mesenchymal signature does not increase monotonically. Instead, we observe a peak in mesenchymal scores followed by a decline, as cells further de-differentiate. This biphasic response recapitulates the dynamics of melanocyte development, suggesting close interactions among genes controlling differentiation and mesenchymal programs in melanocytes. Similar trends were noted for metabolic changes often associated with EMT in carcinomas in which progression along mesenchymal axis correlates with the downregulation of oxidative phosphorylation, while largely maintaining glycolytic capacity. Overall, these results provide an explanation for how EMT and de-differentiation axes overlap with respect to their transcriptional and metabolic programs in melanoma.

Systems-level network modeling deciphers the master regulators of phenotypic plasticity and heterogeneity in melanoma.

Phenotypic (i.e. non-genetic) heterogeneity in melanoma drives dedifferentiation, recalcitrance to targeted therapy and immunotherapy, and consequent tumor relapse and metastasis. Various markers or regulators associated with distinct phenotypes in melanoma have been identified, but, how does a network of interactions among these regulators give rise to multiple "attractor" states and phenotypic switching remains elusive. Here, we inferred a network of transcription factors (TFs) that act as master regulators for gene signatures of diverse cell-states in melanoma. Dynamical simulations of this network predicted how this network can settle into different "attractors" (TF expression patterns), suggesting that TF network dynamics drives the emergence of phenotypic heterogeneity. These simulations can recapitulate major phenotypes observed in melanoma and explain de-differentiation trajectory observed upon BRAF inhibition. Our systems-level modeling framework offers a platform to understand trajectories of phenotypic transitions in the landscape of a regulatory TF network and identify novel therapeutic strategies targeting melanoma plasticity.