An in vitro model of tumor heterogeneity resolves genetic, epigenetic, and stochastic sources of cell state variability.

Tumor heterogeneity is a primary cause of treatment failure and acquired resistance in cancer patients. Even in cancers driven by a single mutated oncogene, variability in response to targeted therapies is well known. The existence of additional genomic alterations among tumor cells can only partially explain this variability. As such, nongenetic factors are increasingly seen as critical contributors to tumor relapse and acquired resistance in cancer. Here, we show that both genetic and nongenetic factors contribute to targeted drug response variability in an experimental model of tumor heterogeneity. We observe significant variability to epidermal growth factor receptor (EGFR) inhibition among and within multiple versions and clonal sublines of PC9, a commonly used EGFR mutant nonsmall cell lung cancer (NSCLC) cell line. We resolve genetic, epigenetic, and stochastic components of this variability using a theoretical framework in which distinct genetic states give rise to multiple epigenetic "basins of attraction," across which cells can transition driven by stochastic noise. Using mutational impact analysis, single-cell differential gene expression, and correlations among Gene Ontology (GO) terms to connect genomics to transcriptomics, we establish a baseline for genetic differences driving drug response variability among PC9 cell line versions. Applying the same approach to clonal sublines, we conclude that drug response variability in all but one of the sublines is due to epigenetic differences; in the other, it is due to genetic alterations. Finally, using a clonal drug response assay together with stochastic simulations, we attribute subclonal drug response variability within sublines to stochastic cell fate decisions and confirm that one subline likely contains genetic resistance mutations that emerged in the absence of drug treatment.

Thunor: visualization and analysis of high-throughput dose-response datasets.

High-throughput cell proliferation assays to quantify drug-response are becoming increasingly common and powerful with the emergence of improved automation and multi-time point analysis methods. However, pipelines for analysis of these datasets that provide reproducible, efficient, and interactive visualization and interpretation are sorely lacking. To address this need, we introduce Thunor, an open-source software platform to manage, analyze, and visualize large, dose-dependent cell proliferation datasets. Thunor supports both end-point and time-based proliferation assays as input. It provides a simple, user-friendly interface with interactive plots and publication-quality images of cell proliferation time courses, dose-response curves, and derived dose-response metrics, e.g. IC50, including across datasets or grouped by tags. Tags are categorical labels for cell lines and drugs, used for aggregation, visualization and statistical analysis, e.g. cell line mutation or drug class/target pathway. A graphical plate map tool is included to facilitate plate annotation with cell lines, drugs and concentrations upon data upload. Datasets can be shared with other users via point-and-click access control. We demonstrate the utility of Thunor to examine and gain insight from two large drug response datasets: a large, publicly available cell viability database and an in-house, high-throughput proliferation rate dataset. Thunor is available from www.thunor.net.

Systems-level network modeling of Small Cell Lung Cancer subtypes identifies master regulators and destabilizers.

Adopting a systems approach, we devise a general workflow to define actionable subtypes in human cancers. Applied to small cell lung cancer (SCLC), the workflow identifies four subtypes based on global gene expression patterns and ontologies. Three correspond to known subtypes (SCLC-A, SCLC-N, and SCLC-Y), while the fourth is a previously undescribed ASCL1+ neuroendocrine variant (NEv2, or SCLC-A2). Tumor deconvolution with subtype gene signatures shows that all of the subtypes are detectable in varying proportions in human and mouse tumors. To understand how multiple stable subtypes can arise within a tumor, we infer a network of transcription factors and develop BooleaBayes, a minimally-constrained Boolean rule-fitting approach. In silico perturbations of the network identify master regulators and destabilizers of its attractors. Specific to NEv2, BooleaBayes predicts ELF3 and NR0B1 as master regulators of the subtype, and TCF3 as a master destabilizer. Since the four subtypes exhibit differential drug sensitivity, with NEv2 consistently least sensitive, these findings may lead to actionable therapeutic strategies that consider SCLC intratumoral heterogeneity. Our systems-level approach should generalize to other cancer types.

An unbiased metric of antiproliferative drug effect in vitro.

In vitro cell proliferation assays are widely used in pharmacology, molecular biology, and drug discovery. Using theoretical modeling and experimentation, we show that current metrics of antiproliferative small molecule effect suffer from time-dependent bias, leading to inaccurate assessments of parameters such as drug potency and efficacy. We propose the drug-induced proliferation (DIP) rate, the slope of the line on a plot of cell population doublings versus time, as an alternative, time-independent metric.

A Nonquiescent "Idling" Population State in Drug-Treated, BRAF-Mutated Melanoma.

Targeted therapy is an effective standard of care in BRAF-mutated malignant melanoma. However, the duration of tumor remission varies unpredictably among patients, and relapse is almost inevitable. Here, we examine the responses of several BRAF-mutated melanoma cell lines (including isogenic subclones) to BRAF inhibitors. We observe complex response dynamics across cell lines, with short-term responses (<100 h) varying from cell line to cell line. In the long term, however, we observe equilibration of all drug-treated populations into a nonquiescent state characterized by a balanced rate of death and division, which we term the "idling" state, and to our knowledge, this state has not been previously reported. Using mathematical modeling, we propose that the observed population-level dynamics are the result of cells transitioning between basins of attraction within a drug-modified phenotypic landscape. Each basin is associated with a drug-induced proliferation rate, a recently introduced metric of an antiproliferative drug effect. The idling population state represents a new dynamic equilibrium in which cells are distributed across the landscape such that the population achieves zero net growth. By fitting our model to experimental drug-response data, we infer the phenotypic landscapes of all considered melanoma cell lines and provide a unifying view of how BRAF-mutated melanomas respond to BRAF inhibition. We hypothesize that the residual disease observed in patients after targeted therapy is composed of a significant number of idling cells. Thus, defining molecular determinants of the phenotypic landscape that idling populations occupy may lead to "targeted landscaping" therapies based on rational modification of the landscape to favor basins with greater drug susceptibility.

Altered TGF-α/ß signaling drives cooperation between breast cancer cell populations.

The role of tumor heterogeneity in regulating disease progression is poorly understood. We hypothesized that interactions between subpopulations of cancer cells can affect the progression of tumors selecting for a more aggressive phenotype. We developed an in vivo assay based on the immortalized nontumorigenic breast cell line MCF10A and its Ras-transformed derivatives AT1 (mildly tumorigenic) and CA1d (highly tumorigenic). CA1d cells outcompeted MCF10A, forming invasive tumors. AT1 grafts were approximately 1% the size of CA1d tumors when initiated using identical cell numbers. In contrast, CA1d/AT1 mixed tumors were larger than tumors composed of AT1 alone (100-fold) or CA1d (3-fold), suggesting cooperation in tumor growth. One of the mechanisms whereby CA1d and AT1 were found to cooperate was by modulation of TGF-α and TGF-ß signaling. Both of these molecules were sufficient to induce changes in AT1 proliferative potential in vitro. Reisolation of AT1 tumor-derived (AT1-TD) cells from these mixed tumors revealed that AT1-TD cells grew in vivo, forming tumors as large as tumorigenic CA1d cells. Cooperation between subpopulations of cancer epithelium is an understudied mechanism of tumor growth and invasion that may have implications on tumor resistance to current therapies.-Franco, O. E., Tyson, D. R., Konvinse, K. C., Udyavar, A. R., Estrada, L., Quaranta, V., Crawford, S. E., Hayward, S. W. Altered TGF-α/ß signaling drives cooperation between breast cancer cell populations.

Quantifying heterogeneity and dynamics of clonal fitness in response to perturbation.

The dynamics of heterogeneous clonal lineages within a cell population, in aggregate, shape both normal and pathological biological processes. Studies of clonality typically relate the fitness of clones to their relative abundance, thus requiring long-term experiments and limiting conclusions about the heterogeneity of clonal fitness in response to perturbation. We present, for the first time, a method that enables a dynamic, global picture of clonal fitness within a mammalian cell population. This novel assay allows facile comparison of the structure of clonal fitness in a cell population across many perturbations. By utilizing high-throughput imaging, our methodology provides ample statistical power to define clonal fitness dynamically and to visualize the structure of perturbation-induced clonal fitness within a cell population. We envision that this technique will be a powerful tool to investigate heterogeneity in biological processes involving cell proliferation, including development and drug response.

Fractional proliferation: a method to deconvolve cell population dynamics from single-cell data.

We present an integrated method that uses extended time-lapse automated imaging to quantify the dynamics of cell proliferation. Cell counts are fit with a quiescence-growth model that estimates rates of cell division, entry into quiescence and death. The model is constrained with rates extracted experimentally from the behavior of tracked single cells over time. We visualize the output of the analysis in fractional proliferation graphs, which deconvolve dynamic proliferative responses to perturbations into the relative contributions of dividing, quiescent (nondividing) and dead cells. The method reveals that the response of 'oncogene-addicted' human cancer cells to tyrosine kinase inhibitors is a composite of altered rates of division, death and entry into quiescence, a finding that challenges the notion that such cells simply die in response to oncogene-targeted therapy.

Purinergic Ca2+ signaling as a novel mechanism of drug tolerance in BRAF mutant melanoma.

Drug tolerance is a major cause of relapse after cancer treatment. In spite of intensive efforts1-9, its molecular basis remains poorly understood, hampering actionable intervention. We report a previously unrecognized signaling mechanism supporting drug tolerance in BRAF-mutant melanoma treated with BRAF inhibitors that could be of general relevance to other cancers. Its key features are cell-intrinsic intracellular Ca2+ signaling initiated by P2X7 receptors (purinergic ligand-gated cation channels), and an enhanced ability for these Ca2+ signals to reactivate ERK1/2 in the drug-tolerant state. Extracellular ATP, virtually ubiquitous in living systems, is the ligand that can initiate Ca2+ spikes via P2X7 channels. ATP is abundant in the tumor microenvironment and is released by dying cells, ironically implicating treatment-initiated cancer cell death as a source of trophic stimuli that leads to ERK reactivation and drug tolerance. Such a mechanism immediately offers an explanation of the inevitable relapse after BRAFi treatment in BRAF-mutant melanoma, and points to actionable strategies to overcome it.

Purinergic Ca2+ Signaling as a Novel Mechanism of Drug Tolerance in BRAF-Mutant Melanoma.

Drug tolerance is a major cause of relapse after cancer treatment. Despite intensive efforts, its molecular basis remains poorly understood, hampering actionable intervention. We report a previously unrecognized signaling mechanism supporting drug tolerance in BRAF-mutant melanoma treated with BRAF inhibitors that could be of general relevance to other cancers. Its key features are cell-intrinsic intracellular Ca2+ signaling initiated by P2X7 receptors (purinergic ligand-gated cation channels) and an enhanced ability for these Ca2+ signals to reactivate ERK1/2 in the drug-tolerant state. Extracellular ATP, virtually ubiquitous in living systems, is the ligand that can initiate Ca2+ spikes via P2X7 channels. ATP is abundant in the tumor microenvironment and is released by dying cells, ironically implicating treatment-initiated cancer cell death as a source of trophic stimuli that leads to ERK reactivation and drug tolerance. Such a mechanism immediately offers an explanation of the inevitable relapse after BRAFi treatment in BRAF-mutant melanoma and points to actionable strategies to overcome it.

DNA damage response signatures are associated with frontline chemotherapy response and routes of tumor evolution in extensive stage small cell lung cancer.

Introduction: A hallmark of small cell lung cancer (SCLC) is its recalcitrance to therapy. While most SCLCs respond to frontline therapy, resistance inevitably develops. Identifying phenotypes potentiating chemoresistance and immune evasion is a crucial unmet need. Previous reports have linked upregulation of the DNA damage response (DDR) machinery to chemoresistance and immune evasion across cancers. However, it is unknown if SCLCs exhibit distinct DDR phenotypes. Methods: To study SCLC DDR phenotypes, we developed a new DDR gene analysis method and applied it to SCLC clinical samples, in vitro , and in vivo model systems. We then investigated how DDR regulation is associated with SCLC biology, chemotherapy response, and tumor evolution following therapy. Results: Using multi-omic profiling, we demonstrate that SCLC tumors cluster into three DDR phenotypes with unique molecular features. Hallmarks of these DDR clusters include differential expression of DNA repair genes, increased replication stress, and heightened G2/M cell cycle arrest. SCLCs with elevated DDR phenotypes exhibit increased neuroendocrine features and decreased "inflamed" biomarkers, both within and across SCLC subtypes. Treatment naive DDR status identified SCLC patients with different responses to frontline chemotherapy. Tumors with initial DDR Intermediate and DDR High phenotypes demonstrated greater tendency for subtype switching and emergence of heterogeneous phenotypes following treatment. Conclusions: We establish that SCLC can be classified into one of three distinct, clinically relevant DDR clusters. Our data demonstrates that DDR status plays a key role in shaping SCLC phenotypes, chemotherapy response, and patterns of tumor evolution. Future work targeting DDR specific phenotypes will be instrumental in improving patient outcomes.

Involvement of Epithelial-Mesenchymal Transition Genes in Small Cell Lung Cancer Phenotypic Plasticity.

Small cell lung cancer (SCLC) is an aggressive cancer recalcitrant to treatment, arising predominantly from epithelial pulmonary neuroendocrine (NE) cells. Intratumor heterogeneity plays critical roles in SCLC disease progression, metastasis, and treatment resistance. At least five transcriptional SCLC NE and non-NE cell subtypes were recently defined by gene expression signatures. Transition from NE to non-NE cell states and cooperation between subtypes within a tumor likely contribute to SCLC progression by mechanisms of adaptation to perturbations. Therefore, gene regulatory programs distinguishing SCLC subtypes or promoting transitions are of great interest. Here, we systematically analyze the relationship between SCLC NE/non-NE transition and epithelial to mesenchymal transition (EMT)-a well-studied cellular process contributing to cancer invasiveness and resistance-using multiple transcriptome datasets from SCLC mouse tumor models, human cancer cell lines, and tumor samples. The NE SCLC-A2 subtype maps to the epithelial state. In contrast, SCLC-A and SCLC-N (NE) map to a partial mesenchymal state (M1) that is distinct from the non-NE, partial mesenchymal state (M2). The correspondence between SCLC subtypes and the EMT program paves the way for further work to understand gene regulatory mechanisms of SCLC tumor plasticity with applicability to other cancer types.

Addressing persistent challenges in digital image analysis of cancerous tissues.

The National Cancer Institute (NCI) supports many research programs and consortia, many of which use imaging as a major modality for characterizing cancerous tissue. A trans-consortia Image Analysis Working Group (IAWG) was established in 2019 with a mission to disseminate imaging-related work and foster collaborations. In 2022, the IAWG held a virtual hackathon focused on addressing challenges of analyzing high dimensional datasets from fixed cancerous tissues. Standard image processing techniques have automated feature extraction, but the next generation of imaging data requires more advanced methods to fully utilize the available information. In this perspective, we discuss current limitations of the automated analysis of multiplexed tissue images, the first steps toward deeper understanding of these limitations, what possible solutions have been developed, any new or refined approaches that were developed during the Image Analysis Hackathon 2022, and where further effort is required. The outstanding problems addressed in the hackathon fell into three main themes: 1) challenges to cell type classification and assessment, 2) translation and visual representation of spatial aspects of high dimensional data, and 3) scaling digital image analyses to large (multi-TB) datasets. We describe the rationale for each specific challenge and the progress made toward addressing it during the hackathon. We also suggest areas that would benefit from more focus and offer insight into broader challenges that the community will need to address as new technologies are developed and integrated into the broad range of image-based modalities and analytical resources already in use within the cancer research community.

Inhibition of phosphodiesterase-4 (PDE4) activity triggers luminal apoptosis and AKT dephosphorylation in a 3-D colonic-crypt model.

BACKGROUND: We previously established a three-dimensional (3-D) colonic crypt model using HKe3 cells which are human colorectal cancer (CRC) HCT116 cells with a disruption in oncogenic KRAS, and revealed the crucial roles of oncogenic KRAS both in inhibition of apoptosis and in disruption of cell polarity; however, the molecular mechanism of KRAS-induced these 3-D specific biological changes remains to be elucidated. RESULTS: Among the genes that were upregulated by oncogenic KRAS in this model, we focused on the phosphodiesterase 4B (PDE4B) of which expression levels were found to be higher in clinical tumor samples from CRC patients in comparison to those from healthy control in the public datasets of gene expression analysis. PDE4B2 was specifically overexpressed among other PDE4 isoforms, and re-expression of oncogenic KRAS in HKe3 cells resulted in PDE4B overexpression. Furthermore, the inhibition of PDE4 catalytic activity using rolipram reverted the disorganization of HCT116 cells into the normal physiologic state of the epithelial cell polarity by inducing the apical assembly of ZO-1 (a tight junction marker) and E-cadherin (an adherens junction marker) and by increasing the activity of caspase-3 (an apoptosis marker) in luminal cavities. Notably, rolipram reduced the AKT phosphorylation, which is known to be associated with the disruption of luminal cavity formation and CRC development. Similar results were also obtained using PDE4B2-shRNAs. In addition, increased expression of PDE4B mRNA was found to be correlated with relapsed CRC in a public datasets of gene expression analysis. CONCLUSIONS: These results collectively suggested that PDE4B is upregulated by oncogenic KRAS, and also that the inhibition of PDE4 catalytic activity can induce both epithelial cell polarity and luminal apoptosis in CRC, thus highlighting the utility of our 3-D culture (3 DC) model for the KRAS-induced development of CRC in 3-D microenvironment. Indeed, using this model, we found that PDE4B is a promising candidate for a therapeutic target as well as prognostic molecular marker in CRC. Further elucidation of the signaling network of PDE4B2 in 3 DC would provide a better understanding of CRC in vivo.

The contribution of age structure to cell population responses to targeted therapeutics.

Cells grown in culture act as a model system for analyzing the effects of anticancer compounds, which may affect cell behavior in a cell cycle position-dependent manner. Cell synchronization techniques have been generally employed to minimize the variation in cell cycle position. However, synchronization techniques are cumbersome and imprecise and the agents used to synchronize the cells potentially have other unknown effects on the cells. An alternative approach is to determine the age structure in the population and account for the cell cycle positional effects post hoc. Here we provide a formalism to use quantifiable lifespans from live cell microscopy experiments to parameterize an age-structured model of cell population response.

Real-time luminescence enables continuous drug-response analysis in adherent and suspension cell lines.

The drug-induced proliferation (DIP) rate is a metric of in vitro drug response that avoids inherent biases in commonly used metrics such as 72 h viability. However, DIP rate measurements rely on direct cell counting over time, a laborious task that is subject to numerous challenges, including the need to fluorescently label cells and automatically segment nuclei. Moreover, it is incredibly difficult to directly count cells and accurately measure DIP rates for cell populations in suspension. As an alternative, we use real-time luminescence measurements derived from the cellular activity of NAD(P)H oxidoreductase to efficiently estimate drug response in both adherent and suspension cell populations to a panel of known anticancer agents. For the adherent cell lines, we collect both luminescence reads and direct cell counts over time simultaneously to assess their congruency. Our results demonstrate that the proposed approach significantly speeds up data collection, avoids the need for cellular labels and image segmentation, and opens the door to significant advances in high-throughput screening of anticancer drugs.

Utilizing Three-Dimensional Culture Methods to Improve High-Throughput Drug Screening in Anaplastic Thyroid Carcinoma.

Anaplastic thyroid carcinoma (ATC) is the most aggressive endocrine neoplasm, with a median survival of just four to six months post-diagnosis. Even with surgical and chemotherapeutic interventions, the five-year survival rate is less than 5%. Although combination dabrafenib/trametinib therapy was recently approved for treatment of the ~25% of ATCs harboring BRAFV600E mutations, there are no approved, effective treatments for BRAF-wildtype disease. Herein, we perform a screen of 1525 drugs and evaluate therapeutic candidates using monolayer cell lines and four corresponding spheroid models of anaplastic thyroid carcinoma. We utilize three-dimensional culture methods, as they have been shown to more accurately recapitulate tumor responses in vivo. These three-dimensional cultures include four distinct ATC spheroid lines representing unique morphology and mutational drivers to provide drug prioritization that will be more readily translatable to the clinic. Using this screen, we identify three exceptionally potent compounds (bortezomib, cabazitaxel, and YM155) that have established safety profiles and could potentially be moved into clinical trial for the treatment of anaplastic thyroid carcinoma, a disease with few treatment options.

A community-based approach to image analysis of cells, tissues and tumors.

Emerging multiplexed imaging platforms provide an unprecedented view of an increasing number of molecular markers at subcellular resolution and the dynamic evolution of tumor cellular composition. As such, they are capable of elucidating cell-to-cell interactions within the tumor microenvironment that impact clinical outcome and therapeutic response. However, the rapid development of these platforms has far outpaced the computational methods for processing and analyzing the data they generate. While being technologically disparate, all imaging assays share many computational requirements for post-collection data processing. As such, our Image Analysis Working Group (IAWG), composed of researchers in the Cancer Systems Biology Consortium (CSBC) and the Physical Sciences - Oncology Network (PS-ON), convened a workshop on "Computational Challenges Shared by Diverse Imaging Platforms" to characterize these common issues and a follow-up hackathon to implement solutions for a selected subset of them. Here, we delineate these areas that reflect major axes of research within the field, including image registration, segmentation of cells and subcellular structures, and identification of cell types from their morphology. We further describe the logistical organization of these events, believing our lessons learned can aid others in uniting the imaging community around self-identified topics of mutual interest, in designing and implementing operational procedures to address those topics and in mitigating issues inherent in image analysis (e.g., sharing exemplar images of large datasets and disseminating baseline solutions to hackathon challenges through open-source code repositories).

Archetype tasks link intratumoral heterogeneity to plasticity and cancer hallmarks in small cell lung cancer.

Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC. A record of this paper's Transparent Peer Review process is included in the supplemental information.

Composite thin film and electrospun biomaterials for urologic tissue reconstruction.

A replacement material for autologous grafts for urinary tract reconstruction would dramatically reduce the complications of surgery for these procedures. However, acellular materials have not proven to work sufficiently well, and cell-seeded materials are technically challenging and time consuming to generate. An important function of the urinary tract is to prevent urine leakage into the surrounding tissue--a function usually performed by the urothelium. We hypothesize that by providing an impermeable barrier in the acellular graft material, urine leakage would be minimized, as the urothelium forms in vivo. However, since urothelial cells require access to nutrients from the supporting vasculature, the impermeable barrier must degrade over time. Here we present the development of a novel biomaterial composed of the common degradable polymers, poly(ε-caprolactone) and poly(L-lactic acid) and generated by electrospinning directly onto spin-coated thin films. The composite scaffolds with thin films on the luminal surface were compared to their electrospun counterparts and commercially available small intestinal submucosa by surface analysis using scanning electron microscopy and by analysis of permeability to small molecules. In addition, the materials were examined for their ability to support urothelial cell adhesion, proliferation, and multilayered urothelium formation. We provide evidence that these unique composite scaffolds provide significant benefit over commonly used acellular materials in vitro and suggest that they be further examined in vivo.