Paradoxical cancer cell proliferation after FGFR inhibition through decreased p21 signaling in FGFR1-amplified breast cancer cells.

Fibroblast growth factors (FGFs) control various cellular functions through fibroblast growth factor receptor (FGFR) activation, including proliferation, differentiation, migration, and survival. FGFR amplification in ER + breast cancer patients correlate with poor prognosis, and FGFR inhibitors are currently being tested in clinical trials. By comparing three-dimensional spheroid growth of ER + breast cancer cells with and without FGFR1 amplification, our research discovered that FGF2 treatment can paradoxically decrease proliferation in cells with FGFR1 amplification or overexpression. In contrast, FGF2 treatment in cells without FGFR1 amplification promotes classical FGFR proliferative signaling through the MAPK cascade. The growth inhibitory effect of FGF2 in FGFR1 amplified cells aligned with an increase in p21, a cell cycle inhibitor that hinders the G1 to S phase transition in the cell cycle. Additionally, FGF2 addition in FGFR1 amplified cells activated JAK-STAT signaling and promoted a stem cell-like state. FGF2-induced paradoxical effects were reversed by inhibiting p21 or the JAK-STAT pathway and with pan-FGFR inhibitors. Analysis of patient ER + breast tumor transcriptomes from the TCGA and METABRIC datasets demonstrated a strong positive association between expression of FGF2 and stemness signatures, which was further enhanced in tumors with high FGFR1 expression. Overall, our findings reveal a divergence in FGFR signaling, transitioning from a proliferative to stemness state driven by activation of JAK-STAT signaling and modulation of p21 levels. Activation of these divergent signaling pathways in FGFR amplified cancer cells and paradoxical growth effects highlight a challenge in the use of FGFR inhibitors in cancer treatment.

Knowledge-based classification of fine-grained immune cell types in single-cell RNA-Seq data.

Single-cell RNA sequencing (scRNA-Seq) is an emerging strategy for characterizing immune cell populations. Compared to flow or mass cytometry, scRNA-Seq could potentially identify cell types and activation states that lack precise cell surface markers. However, scRNA-Seq is currently limited due to the need to manually classify each immune cell from its transcriptional profile. While recently developed algorithms accurately annotate coarse cell types (e.g. T cells versus macrophages), making fine distinctions (e.g. CD8+ effector memory T cells) remains a difficult challenge. To address this, we developed a machine learning classifier called ImmClassifier that leverages a hierarchical ontology of cell type. We demonstrate that its predictions are highly concordant with flow-based markers from CITE-seq and outperforms other tools (+15% recall, +14% precision) in distinguishing fine-grained cell types with comparable performance on coarse ones. Thus, ImmClassifier can be used to explore more deeply the heterogeneity of the immune system in scRNA-Seq experiments.

ENDORSE: a prognostic model for endocrine therapy in estrogen-receptor-positive breast cancers.

Advanced and metastatic estrogen receptor-positive (ER+ ) breast cancers are often endocrine resistant. However, endocrine therapy remains the primary treatment for all advanced ER+ breast cancers. Treatment options that may benefit resistant cancers, such as add-on drugs that target resistance pathways or switching to chemotherapy, are only available after progression on endocrine therapy. Here we developed an endocrine therapy prognostic model for early and advanced ER+ breast cancers. The endocrine resistance (ENDORSE) model is composed of two components, each based on the empirical cumulative distribution function of ranked expression of gene signatures. These signatures include a feature set associated with long-term survival outcomes on endocrine therapy selected using lasso-regularized Cox regression and a pathway-based curated set of genes expressed in response to estrogen. We extensively validated ENDORSE in multiple ER+ clinical trial datasets and demonstrated superior and consistent performance of the model over clinical covariates, proliferation markers, and multiple published signatures. Finally, genomic and pathway analyses in patient data revealed possible mechanisms that may help develop rational stratification strategies for endocrine-resistant ER+ breast cancer patients.

Circulating immune cell phenotype dynamics reflect the strength of tumor-immune cell interactions in patients during immunotherapy.

The extent to which immune cell phenotypes in the peripheral blood reflect within-tumor immune activity prior to and early in cancer therapy is unclear. To address this question, we studied the population dynamics of tumor and immune cells, and immune phenotypic changes, using clinical tumor and immune cell measurements and single-cell genomic analyses. These samples were serially obtained from a cohort of advanced gastrointestinal cancer patients enrolled in a trial with chemotherapy and immunotherapy. Using an ecological population model, fitted to clinical tumor burden and immune cell abundance data from each patient, we find evidence of a strong tumor-circulating immune cell interaction in responder patients but not in those patients that progress on treatment. Upon initiation of therapy, immune cell abundance increased rapidly in responsive patients, and once the peak level is reached tumor burden decreases, similar to models of predator-prey interactions; these dynamic patterns were absent in nonresponder patients. To interrogate phenotype dynamics of circulating immune cells, we performed single-cell RNA sequencing at serial time points during treatment. These data show that peripheral immune cell phenotypes were linked to the increased strength of patients' tumor-immune cell interaction, including increased cytotoxic differentiation and strong activation of interferon signaling in peripheral T cells in responder patients. Joint modeling of clinical and genomic data highlights the interactions between tumor and immune cell populations and reveals how variation in patient responsiveness can be explained by differences in peripheral immune cell signaling and differentiation soon after the initiation of immunotherapy.

Topsentinol L Trisulfate, a Marine Natural Product That Targets Basal-like and Claudin-Low Breast Cancers.

Patients diagnosed with basal-like breast cancer suffer from poor prognosis and limited treatment options. There is an urgent need to identify new targets that can benefit patients with basal-like and claudin-low (BL-CL) breast cancers. We screened fractions from our Marine Invertebrate Compound Library (MICL) to identify compounds that specifically target BL-CL breast cancers. We identified a previously unreported trisulfated sterol, i.e., topsentinol L trisulfate (TLT), which exhibited increased efficacy against BL-CL breast cancers relative to luminal/HER2+ breast cancer. Biochemical investigation of the effects of TLT on BL-CL cell lines revealed its ability to inhibit activation of AMP-activated protein kinase (AMPK) and checkpoint kinase 1 (CHK1) and to promote activation of p38. The importance of targeting AMPK and CHK1 in BL-CL cell lines was validated by treating a panel of breast cancer cell lines with known small molecule inhibitors of AMPK (dorsomorphin) and CHK1 (Ly2603618) and recording the increased effectiveness against BL-CL breast cancers as compared with luminal/HER2+ breast cancer. Finally, we generated a drug response gene-expression signature and projected it against a human tumor panel of 12 different cancer types to identify other cancer types sensitive to the compound. The TLT sensitivity gene-expression signature identified breast and bladder cancer as the most sensitive to TLT, while glioblastoma multiforme was the least sensitive.

Exploiting collateral sensitivity controls growth of mixed culture of sensitive and resistant cells and decreases selection for resistant cells in a cell line model.

BACKGROUND: CDK4/6 inhibitors such as ribociclib are becoming widely used targeted therapies in hormone-receptor-positive (HR+) human epidermal growth factor receptor 2-negative (HER2-) breast cancer. However, cancers can advance due to drug resistance, a problem in which tumor heterogeneity and evolution are key features. METHODS: Ribociclib-resistant HR+/HER2- CAMA-1 breast cancer cells were generated through long-term ribociclib treatment. Characterization of sensitive and resistant cells were performed using RNA sequencing and whole exome sequencing. Lentiviral labeling with different fluorescent proteins enabled us to track the proliferation of sensitive and resistant cells under different treatments in a heterogeneous, 3D spheroid coculture system using imaging microscopy and flow cytometry. RESULTS: Transcriptional profiling of sensitive and resistant cells revealed the downregulation of the G2/M checkpoint in the resistant cells. Exploiting this acquired vulnerability; resistant cells exhibited collateral sensitivity for the Wee-1 inhibitor, adavosertib (AZD1775). The combination of ribociclib and adavosertib achieved additional antiproliferative effect exclusively in the cocultures compared to monocultures, while decreasing the selection for resistant cells. CONCLUSIONS: Our results suggest that optimal antiproliferative effects in heterogeneous cancers can be achieved via an integrative therapeutic approach targeting sensitive and resistant cancer cell populations within a tumor, respectively.

Pathway activity profiling of growth factor receptor network and stemness pathways differentiates metaplastic breast cancer histological subtypes.

BACKGROUND: Gene expression profiling of rare cancers has proven challenging due to limited access to patient materials and requirement of intact, non-degraded RNA for next-generation sequencing. We customized a gene expression panel compatible with degraded RNA from formalin-fixed, paraffin-embedded (FFPE) patient cancer samples and investigated its utility in pathway activity profiling in patients with metaplastic breast cancer (MpBC). METHODS: Activity of various biological pathways was profiled in samples from nineteen patients with MpBC and 8 patients with invasive ductal carcinoma with triple negative breast cancer (TNBC) phenotype using a custom gene expression-based assay of 345 genes. RESULTS: MpBC samples of mesenchymal (chondroid and/or osteoid) histology demonstrated increased SNAI1 and BCL2L11 pathway activity compared to samples with non-mesenchymal histology. Additionally, late cornified envelope and keratinization genes were downregulated in MpBC compared to TNBC, and epithelial-to-mesenchymal transition (EMT) and collagen genes were upregulated in MpBC. Patients with high activity of an invasiveness gene expression signature, as well as high expression of the mesenchymal marker and extracellular matrix glycoprotein gene SPARC, experienced worse outcomes than those with low invasiveness activity and low SPARC expression. CONCLUSIONS: This study demonstrates the utility of gene expression profiling of metaplastic breast cancer FFPE samples with a custom counts-based assay. Gene expression patterns identified by this assay suggest that, although often histologically triple negative, patients with MpBC have distinct pathway activation compared to patients with invasive ductal TNBC. Incorporation of targeted therapies may lead to improved outcome for MpBC patients, especially in those patients expressing increased activity of invasiveness pathways.

The value of genomics in dissecting the RAS-network and in guiding therapeutics for RAS-driven cancers.

The rise in genomic knowledge over the past decade has revealed the molecular etiology of many diseases, and has identified intricate signaling network activity in human cancers. Genomics provides the opportunity to determine genome structure and capture the activity of thousands of molecular events concurrently, which is important for deciphering highly complex genetic diseases such as cancer. In this review, we focus on genomic efforts directed towards one of cancer's most frequently mutated networks, the RAS pathway. Genomic tools such as gene expression signatures and assessment of mutations across the RAS network enable the capture of RAS signaling complexity. Due to this high level of interaction and cross-talk within the network, efforts to target RAS signaling in the clinic have generally failed, and we currently lack the ability to directly inhibit the RAS protein with high efficacy. We propose that the use of gene expression data can identify effective treatments that broadly inhibit the RAS network as this approach measures pathway activity independent of mutation status or any single mechanism of activation. Here, we review the genomic studies that map the complexity of the RAS network in cancer, and that show how genomic measurements of RAS pathway activation can identify effective RAS inhibition strategies. We also address the challenges and future directions for treating RAS-driven tumors. In summary, genomic assessment of RAS signaling provides a level of complexity necessary to accurately map the network that matches the intricacy of RAS pathway interactions in cancer.

Integrative analyses reveal signaling pathways underlying familial breast cancer susceptibility.

The signaling events that drive familial breast cancer (FBC) risk remain poorly understood. While the majority of genomic studies have focused on genetic risk variants, known risk variants account for at most 30% of FBC cases. Considering that multiple genes may influence FBC risk, we hypothesized that a pathway-based strategy examining different data types from multiple tissues could elucidate the biological basis for FBC. In this study, we performed integrated analyses of gene expression and exome-sequencing data from peripheral blood mononuclear cells and showed that cell adhesion pathways are significantly and consistently dysregulated in women who develop FBC. The dysregulation of cell adhesion pathways in high-risk women was also identified by pathway-based profiling applied to normal breast tissue data from two independent cohorts. The results of our genomic analyses were validated in normal primary mammary epithelial cells from high-risk and control women, using cell-based functional assays, drug-response assays, fluorescence microscopy, and Western blotting assays. Both genomic and cell-based experiments indicate that cell-cell and cell-extracellular matrix adhesion processes seem to be disrupted in non-malignant cells of women at high risk for FBC and suggest a potential role for these processes in FBC development.

A genomic strategy to elucidate modules of oncogenic pathway signaling networks.

Recent studies have emphasized the importance of pathway-specific interpretations for understanding the functional relevance of gene alterations in human cancers. Although signaling activities are often conceptualized as linear events, in reality, they reflect the activity of complex functional networks assembled from modules that each respond to input signals. To acquire a deeper understanding of this network structure, we developed an approach to deconstruct pathways into modules represented by gene expression signatures. Our studies confirm that they represent units of underlying biological activity linked to known biochemical pathway structures. Importantly, we show that these signaling modules provide tools to dissect the complexity of oncogenic states that define disease outcomes as well as response to pathway-specific therapeutics. We propose that this model of pathway structure constitutes a framework to study the processes by which information propogates through cellular networks and to elucidate the relationships of fundamental modules to cellular and clinical phenotypes.

Alternative preprocessing of RNA-Sequencing data in The Cancer Genome Atlas leads to improved analysis results.

MOTIVATION: The Cancer Genome Atlas (TCGA) RNA-Sequencing data are used widely for research. TCGA provides 'Level 3' data, which have been processed using a pipeline specific to that resource. However, we have found using experimentally derived data that this pipeline produces gene-expression values that vary considerably across biological replicates. In addition, some RNA-Sequencing analysis tools require integer-based read counts, which are not provided with the Level 3 data. As an alternative, we have reprocessed the data for 9264 tumor and 741 normal samples across 24 cancer types using the Rsubread package. We have also collated corresponding clinical data for these samples. We provide these data as a community resource. RESULTS: We compared TCGA samples processed using either pipeline and found that the Rsubread pipeline produced fewer zero-expression genes and more consistent expression levels across replicate samples than the TCGA pipeline. Additionally, we used a genomic-signature approach to estimate HER2 (ERBB2) activation status for 662 breast-tumor samples and found that the Rsubread data resulted in stronger predictions of HER2 pathway activity. Finally, we used data from both pipelines to classify 575 lung cancer samples based on histological type. This analysis identified various non-coding RNA that may influence lung-cancer histology. AVAILABILITY AND IMPLEMENTATION: The RNA-Sequencing and clinical data can be downloaded from Gene Expression Omnibus (accession number GSE62944). Scripts and code that were used to process and analyze the data are available from https://github.com/srp33/TCGA_RNASeq_Clinical. CONTACT: stephen_piccolo@byu.edu or andreab@genetics.utah.edu SUPPLEMENTARY INFORMATION: Supplementary material is available at Bioinformatics online.

ASSIGN: context-specific genomic profiling of multiple heterogeneous biological pathways.

MOTIVATION: Although gene-expression signature-based biomarkers are often developed for clinical diagnosis, many promising signatures fail to replicate during validation. One major challenge is that biological samples used to generate and validate the signature are often from heterogeneous biological contexts-controlled or in vitro samples may be used to generate the signature, but patient samples may be used for validation. In addition, systematic technical biases from multiple genome-profiling platforms often mask true biological variation. Addressing such challenges will enable us to better elucidate disease mechanisms and provide improved guidance for personalized therapeutics. RESULTS: Here, we present a pathway profiling toolkit, Adaptive Signature Selection and InteGratioN (ASSIGN), which enables robust and context-specific pathway analyses by efficiently capturing pathway activity in heterogeneous sets of samples and across profiling technologies. The ASSIGN framework is based on a flexible Bayesian factor analysis approach that allows for simultaneous profiling of multiple correlated pathways and for the adaptation of pathway signatures into specific disease. We demonstrate the robustness and versatility of ASSIGN in estimating pathway activity in simulated data, cell lines perturbed pathways and in primary tissues samples including The Cancer Genome Atlas breast carcinoma samples and liver samples exposed to genotoxic carcinogens. AVAILABILITY AND IMPLEMENTATION: Software for our approach is available for download at: http://www.bioconductor.org/packages/release/bioc/html/ASSIGN.html and https://github.com/wevanjohnson/ASSIGN.

Multiplatform single-sample estimates of transcriptional activation.

Over the past two decades, many biotechnology platforms have been developed for high-throughput gene expression profiling. However, because each platform is subject to technology-specific biases and produces distinct raw-data distributions, researchers have experienced difficulty in integrating data across platforms. Data integration is crucial to data-generating consortiums, researchers transitioning to newer profiling technologies, and individuals seeking to aggregate data across experiments. We address this need with our Universal exPression Code (UPC) approach, which corrects for platform-specific background noise using models that account for the genomic base composition and length of target regions; this approach also uses a mixture model to estimate whether a gene is active in a particular profiling sample. The latter produces standardized UPC values on a zero-to-one scale, so that they can be interpreted consistently, irrespective of profiling technology, thus enabling downstream analysis pipelines to be developed in a platform-agnostic manner. The UPC method can be applied to one- and two-channel expression microarrays and to next-generation sequencing data (RNA sequencing). Furthermore, UPCs are derived using information from within a given sample only--no ancillary samples are required at processing time. Thus, UPCs are suitable for personalized-medicine workflows where samples must be processed individually rather than in batches. In a variety of analyses and comparisons, UPCs perform comparably to other methods designed specifically for microarrays or RNA sequencing in most settings. Software for calculating UPCs is freely available at www.bioconductor.org/packages/release/bioc/html/SCAN.UPC.html.

Linking oncogenic pathways with therapeutic opportunities.

The accumulation of multiple mutations and alterations in the cancer genome underlies the complexity of cancer phenotypes. A consequence of these alterations is the deregulation of various cell-signalling pathways that control cell function. Molecular-profiling studies, particularly DNA microarray analyses, have the potential to describe this complexity. These studies also provide an opportunity to link pathway deregulation with potential therapeutic strategies. This approach, when coupled with other methods for identifying pathway activation, provides an opportunity to both match individual patients with the most appropriate therapeutic strategy and identify potential options for combination therapy.

Patient-derived xenografts of triple-negative breast cancer reproduce molecular features of patient tumors and respond to mTOR inhibition.

INTRODUCTION: Triple-negative breast cancer (TNBC) is aggressive and lacks targeted therapies. Phosphatidylinositide 3-kinase (PI3K)/mammalian target of rapamycin (mTOR) pathways are frequently activated in TNBC patient tumors at the genome, gene expression and protein levels, and mTOR inhibitors have been shown to inhibit growth in TNBC cell lines. We describe a panel of patient-derived xenografts representing multiple TNBC subtypes and use them to test preclinical drug efficacy of two mTOR inhibitors, sirolimus (rapamycin) and temsirolimus (CCI-779). METHODS: We generated a panel of seven patient-derived orthotopic xenografts from six primary TNBC tumors and one metastasis. Patient tumors and corresponding xenografts were compared by histology, immunohistochemistry, array comparative genomic hybridization (aCGH) and phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA) sequencing; TNBC subtypes were determined. Using a previously published logistic regression approach, we generated a rapamycin response signature from Connectivity Map gene expression data and used it to predict rapamycin sensitivity in 1,401 human breast cancers of different intrinsic subtypes, prompting in vivo testing of mTOR inhibitors and doxorubicin in our TNBC xenografts. RESULTS: Patient-derived xenografts recapitulated histology, biomarker expression and global genomic features of patient tumors. Two primary tumors had PIK3CA coding mutations, and five of six primary tumors showed flanking intron single nucleotide polymorphisms (SNPs) with conservation of sequence variations between primary tumors and xenografts, even on subsequent xenograft passages. Gene expression profiling showed that our models represent at least four of six TNBC subtypes. The rapamycin response signature predicted sensitivity for 94% of basal-like breast cancers in a large dataset. Drug testing of mTOR inhibitors in our xenografts showed 77 to 99% growth inhibition, significantly more than doxorubicin; protein phosphorylation studies indicated constitutive activation of the mTOR pathway that decreased with treatment. However, no tumor was completely eradicated. CONCLUSIONS: A panel of patient-derived xenograft models covering a spectrum of TNBC subtypes was generated that histologically and genomically matched original patient tumors. Consistent with in silico predictions, mTOR inhibitor testing in our TNBC xenografts showed significant tumor growth inhibition in all, suggesting that mTOR inhibitors can be effective in TNBC, but will require use with additional therapies, warranting investigation of optimal drug combinations.

A dynamic bronchial airway gene expression signature of chronic obstructive pulmonary disease and lung function impairment.

RATIONALE: Molecular phenotyping of chronic obstructive pulmonary disease (COPD) has been impeded in part by the difficulty in obtaining lung tissue samples from individuals with impaired lung function. OBJECTIVES: We sought to determine whether COPD-associated processes are reflected in gene expression profiles of bronchial airway epithelial cells obtained by bronchoscopy. METHODS: Gene expression profiling of bronchial brushings obtained from 238 current and former smokers with and without COPD was performed using Affymetrix Human Gene 1.0 ST Arrays. MEASUREMENTS AND MAIN RESULTS: We identified 98 genes whose expression levels were associated with COPD status, FEV1% predicted, and FEV1/FVC. In silico analysis identified activating transcription factor 4 (ATF4) as a potential transcriptional regulator of genes with COPD-associated airway expression, and ATF4 overexpression in airway epithelial cells in vitro recapitulates COPD-associated gene expression changes. Genes with COPD-associated expression in the bronchial airway epithelium had similarly altered expression profiles in prior studies performed on small-airway epithelium and lung parenchyma, suggesting that transcriptomic alterations in the bronchial airway epithelium reflect molecular events found at more distal sites of disease activity. Many of the airway COPD-associated gene expression changes revert toward baseline after therapy with the inhaled corticosteroid fluticasone in independent cohorts. CONCLUSIONS: Our findings demonstrate a molecular field of injury throughout the bronchial airway of active and former smokers with COPD that may be driven in part by ATF4 and is modifiable with therapy. Bronchial airway epithelium may ultimately serve as a relatively accessible tissue in which to measure biomarkers of disease activity for guiding clinical management of COPD.

A novel combination therapy for ER+ breast cancer suppresses drug resistance via an evolutionary double-bind.

Chemotherapy remains a commonly used and important treatment option for metastatic breast cancer. A majority of ER+ metastatic breast cancer patients ultimately develop resistance to chemotherapy, resulting in disease progression. We hypothesized that an "evolutionary double-bind", where treatment with one drug improves the response to a different agent, would improve the effectiveness and durability of responses to chemotherapy. This approach exploits vulnerabilities in acquired resistance mechanisms. Evolutionary models can be used in refractory cancer to identify alternative treatment strategies that capitalize on acquired vulnerabilities and resistance traits for improved outcomes. To develop and test these models, ER+ breast cancer cell lineages sensitive and resistant to chemotherapy are grown in spheroids with varied initial population frequencies to measure cross-sensitivity and efficacy of chemotherapy and add-on treatments such as disulfiram combination treatment. Different treatment schedules then assessed the best strategy for reducing the selection of resistant populations. We developed and parameterized a game-theoretic mathematical model from this in vitro experimental data, and used it to predict the existence of a double-bind where selection for resistance to chemotherapy induces sensitivity to disulfiram. The model predicts a dose-dependent re-sensitization (a double-bind) to chemotherapy for monotherapy disulfiram.

A single-sample microarray normalization method to facilitate personalized-medicine workflows.

Gene-expression microarrays allow researchers to characterize biological phenomena in a high-throughput fashion but are subject to technological biases and inevitable variabilities that arise during sample collection and processing. Normalization techniques aim to correct such biases. Most existing methods require multiple samples to be processed in aggregate; consequently, each sample's output is influenced by other samples processed jointly. However, in personalized-medicine workflows, samples may arrive serially, so renormalizing all samples upon each new arrival would be impractical. We have developed Single Channel Array Normalization (SCAN), a single-sample technique that models the effects of probe-nucleotide composition on fluorescence intensity and corrects for such effects, dramatically increasing the signal-to-noise ratio within individual samples while decreasing variation across samples. In various benchmark comparisons, we show that SCAN performs as well as or better than competing methods yet has no dependence on external reference samples and can be applied to any single-channel microarray platform.

Phylogenetic inference from single-cell RNA-seq data.

Tumors are comprised of subpopulations of cancer cells that harbor distinct genetic profiles and phenotypes that evolve over time and during treatment. By reconstructing the course of cancer evolution, we can understand the acquisition of the malignant properties that drive tumor progression. Unfortunately, recovering the evolutionary relationships of individual cancer cells linked to their phenotypes remains a difficult challenge. To address this need, we have developed PhylinSic, a method that reconstructs the phylogenetic relationships among cells linked to their gene expression profiles from single cell RNA-sequencing (scRNA-Seq) data. This method calls nucleotide bases using a probabilistic smoothing approach and then estimates a phylogenetic tree using a Bayesian modeling algorithm. We showed that PhylinSic identified evolutionary relationships underpinning drug selection and metastasis and was sensitive enough to identify subclones from genetic drift. We found that breast cancer tumors resistant to chemotherapies harbored multiple genetic lineages that independently acquired high K-Ras and ß-catenin, suggesting that therapeutic strategies may need to control multiple lineages to be durable. These results demonstrated that PhylinSic can reconstruct evolution and link the genotypes and phenotypes of cells across monophyletic tumors using scRNA-Seq.

ONC201/TIC10 enhances durability of mTOR inhibitor everolimus in metastatic ER+ breast cancer.

The mTOR inhibitor, everolimus, is an important clinical management component of metastatic ER+ breast cancer (BC). However, most patients develop resistance and progress on therapy, highlighting the need to discover strategies that increase mTOR inhibitor effectiveness. We developed ER+ BC cell lines, sensitive or resistant to everolimus, and discovered that combination treatment of ONC201/TIC10 with everolimus inhibited cell growth in 2D/3D in vitro studies. We confirmed increased therapeutic response in primary patient cells progressing on everolimus, supporting clinical relevance. We show that ONC201/TIC10 mechanism in metastatic ER+ BC cells involves oxidative phosphorylation inhibition and stress response activation. Transcriptomic analysis in everolimus resistant breast patient tumors and mitochondrial functional assays in resistant cell lines demonstrated increased mitochondrial respiration dependency, contributing to ONC201/TIC10 sensitivity. We propose that ONC201/TIC10 and modulation of mitochondrial function may provide an effective add-on therapy strategy for patients with metastatic ER+ BCs resistant to mTOR inhibitors.

Breast cancer is one of the most frequently diagnosed cancers globally, particularly among women. The most common type of breast cancer expresses a receptor for the hormone estrogen. Many treatments block the activity of estrogen and therefore slow or block the development and spread of this type of breast cancer. For patients with advanced breast cancer, hormone-blocking treatments work best in combination with other drugs, including one called everolimus. However, in many patients the cancer cells become resistant to these therapies, leading to disease progression and decreased survival. To explore treatment strategies that could enhance the effectiveness of existing therapies for breast cancer, Farmaki et al. studied how cancer cells which had become resistant to everolimus responded when treated with an experimental drug called ONC201/TIC10. A combination of everolimus and ONC201/TIC10 inhibited growth of resistant cancer cells that had been grown in a three-dimensional arrangement to mimic human tumors. Moreover, the drug combination effectively targeted breast cancer cells collected from patients whose cancer had progressed while being treated with everolimus, suggesting that ONC201/TIC10 could be relevant in a clinical setting. Finally, molecular and biochemical experiments revealed that the drug ONC201/TIC10 works by disrupting the pathways that everolimus-resistant cancer cells use to generate the energy required to grow and proliferate. Taken together these findings suggest that ONC201/TIC10 may provide an effective add-on therapy for patients with certain types of advanced breast cancer that are no longer responding to everolimus. Before this becomes a reality for patients, however, there will have to be more experimental testing of ONC201/TIC10 to determine optimal dosing and timing strategy for future clinical trials.