Massively parallel identification of functionally consequential noncoding genetic variants in undiagnosed rare disease patients.

Clinical whole genome sequencing has enabled the discovery of potentially pathogenic noncoding variants in the genomes of rare disease patients with a prior history of negative genetic testing. However, interpreting the functional consequences of noncoding variants and distinguishing those that contribute to disease etiology remains a challenge. Here we address this challenge by experimentally profiling the functional consequences of rare noncoding variants detected in a cohort of undiagnosed rare disease patients at scale using a massively parallel reporter assay. We demonstrate that this approach successfully identifies rare noncoding variants that alter the regulatory capacity of genomic sequences. In addition, we describe an integrative analysis that utilizes genomic features alongside patient clinical data to further prioritize candidate variants with an increased likelihood of pathogenicity. This work represents an important step towards establishing a framework for the functional interpretation of clinically detected noncoding variants.

Novel temporal and spatial patterns of metastatic colonization from breast cancer rapid-autopsy tumor biopsies.

BACKGROUND: Metastatic breast cancer is a deadly disease with a low 5-year survival rate. Tracking metastatic spread in living patients is difficult and thus poorly understood. METHODS: Via rapid autopsy, we have collected 30 tumor samples over 3 timepoints and across 8 organs from a triple-negative metastatic breast cancer patient. The large number of sites sampled, together with deep whole-genome sequencing and advanced computational analysis, allowed us to comprehensively reconstruct the tumor's evolution at subclonal resolution. RESULTS: The most unique, previously unreported aspect of the tumor's evolution that we observed in this patient was the presence of "subclone incubators," defined as metastatic sites where substantial tumor evolution occurs before colonization of additional sites and organs by subclones that initially evolved at the incubator site. Overall, we identified four discrete waves of metastatic expansions, each of which resulted in a number of new, genetically similar metastasis sites that also enriched for particular organs (e.g., abdominal vs bone and brain). The lung played a critical role in facilitating metastatic spread in this patient: the lung was the first site of metastatic escape from the primary breast lesion, subclones at this site were likely the source of all four subsequent metastatic waves, and multiple sites in the lung acted as subclone incubators. Finally, functional annotation revealed that many known drivers or metastasis-promoting tumor mutations in this patient were shared by some, but not all metastatic sites, highlighting the need for more comprehensive surveys of a patient's metastases for effective clinical intervention. CONCLUSIONS: Our analysis revealed the presence of substantial tumor evolution at metastatic incubator sites in a patient, with potentially important clinical implications. Our study demonstrated that sampling of a large number of metastatic sites affords unprecedented detail for studying metastatic evolution.

Tepoxalin increases chemotherapy efficacy in drug-resistant breast cancer cells overexpressing the multidrug transporter gene ABCB1.

Effective cancer chemotherapy treatment requires both therapy delivery and retention by malignant cells. Cancer cell overexpression of the multidrug transmembrane transporter gene ABCB1 (MDR1, multi-drug resistance protein 1) thwarts therapy retention, leading to a drug-resistant phenotype. We explored the phenotypic impact of ABCB1 overexpression in normal human mammary epithelial cells (HMECs) via acute adenoviral delivery and in breast cancer cell lines with stable integration of inducible ABCB1 expression. One hundred sixty-two genes were differentially expressed between ABCB1-expressing and GFP-expressing HMECs, including the gene encoding the cyclooxygenase-2 protein, PTGS2. Several breast cancer cell lines with inducible ABCB1 expression demonstrated sensitivity to the 5-lipoxygenase, cyclooxygenase-1/2 inhibitor tepoxalin in two-dimensional drug response assays, and combination treatment of tepoxalin either with chemotherapies or with histone deacetylase (HDAC) inhibitors improved therapeutic efficacy in these lines. Moreover, selection for the ABCB1-expressing cell population was reduced in three-dimensional co-cultures of ABCB1-expressing and GFP-expressing cells when chemotherapy was given in combination with tepoxalin. Further study is warranted to ascertain the clinical potential of tepoxalin, an FDA-approved therapeutic for use in domesticated mammals, to restore chemosensitivity and improve drug response in patients with ABCB1-overexpressing drug-resistant breast cancers.

A `one-two punch' therapy strategy to target chemoresistance in estrogen receptor positive breast cancer.

Cancer cell phenotypes evolve during a tumor's treatment. In some cases, tumor cells acquire cancer stem cell-like (CSL) traits such as resistance to chemotherapy and diminished differentiation; therefore, targeting these cells may be therapeutically beneficial. In this study we show that in progressive estrogen receptor positive (ER+) metastatic breast cancer tumors, resistant subclones that emerge following chemotherapy have increased CSL abundance. Further, in vitro organoid growth of ER+ patient cancer cells also shows that chemotherapy treatment leads to increased abundance of ALDH+/CD44+ CSL cells. Chemotherapy induced CSL abundance is blocked by treatment with a pan-HDAC inhibitor, belinostat. Belinostat treatment diminished both mammosphere formation and size following chemotherapy, indicating a decrease in progenitor CSL traits. HDAC inhibitors specific to class IIa (HDAC4, HDAC5) and IIb (HDAC6) were shown to primarily reverse the chemo-resistant CSL state. Single-cell RNA sequencing analysis with patient samples showed that HDAC targets and MYC signaling were promoted by chemotherapy and inhibited upon HDAC inhibitor treatment. In summary, HDAC inhibition can block chemotherapy-induced drug resistant phenotypes with 'one-two punch' strategy in refractory breast cancer cells.

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.

Co-stimulatory and co-inhibitory immune markers in solid tumors with MET alterations.

The implication of MET alterations in solid tumors and the immune microenvironment remains elusive. Formalin-fixed, paraffin-embedded samples of 21 patients with solid tumors harboring MET alterations were used for immunohistochemical staining. Extracted RNA was analyzed with the NanoString nCounter human PanCancer immune profiling panel (NanoString Technologies, Inc., WA, USA). Patients were diagnosed with lung (n = 10), breast (n = 5), genitourinary (n = 3) or colorectal cancer (n = 3). Eleven had a MET missense mutation, four had an exon 14 splice site mutation and six had MET amplification. CD6, CCL19, CD40LG, XCR1, MAGEA1, ATM and CCL19 genes were significantly differentially expressed in MET-altered cancers. MET alterations may have a role in various solid tumors as potential therapeutic targets and combination therapy candidates with immune checkpoint inhibitors.

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.

Publisher Correction: Combating subclonal evolution of resistant cancer phenotypes.

The originally published version of this Article contained an error in Figure 4. In panel a, grey boxes surrounding the subclones associated with patients #2 and #4 obscured adjacent portions of the heatmap. This error has now been corrected in both the PDF and HTML versions of the Article.

Combating subclonal evolution of resistant cancer phenotypes.

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves.

Mechanisms and clinical implications of tumor heterogeneity and convergence on recurrent phenotypes.

Tumor heterogeneity has been identified at various -omic levels. The tumor genome, transcriptome, proteome, and phenome can vary widely across cells in patient tumors and are influenced by tumor cell interactions with heterogeneous physical conditions and cellular components of the tumor microenvironment. Here, we explore the concept that while variation exists at multiple -omic levels, changes at each of these levels converge on the same pathways and lead to convergent phenotypes in tumors that can provide common drug targets. These phenotypes include cellular growth and proliferation, sustained oncogenic signaling, and immune avoidance, among others. Tumor heterogeneity complicates treatment of patient cancers as it leads to varied response to therapies. Identification of convergent cellular phenotypes arising in patient cancers and targeted therapies that reverse them has the potential to transform the way clinicians treat these cancers and to improve patient outcome.

Activity of distinct growth factor receptor network components in breast tumors uncovers two biologically relevant subtypes.

BACKGROUND: The growth factor receptor network (GFRN) plays a significant role in driving key oncogenic processes. However, assessment of global GFRN activity is challenging due to complex crosstalk among GFRN components, or pathways, and the inability to study complex signaling networks in patient tumors. Here, pathway-specific genomic signatures were used to interrogate GFRN activity in breast tumors and the consequent phenotypic impact of GRFN activity patterns. METHODS: Novel pathway signatures were generated in human primary mammary epithelial cells by overexpressing key genes from GFRN pathways (HER2, IGF1R, AKT1, EGFR, KRAS (G12V), RAF1, BAD). The pathway analysis toolkit Adaptive Signature Selection and InteGratioN (ASSIGN) was used to estimate pathway activity for GFRN components in 1119 breast tumors from The Cancer Genome Atlas (TCGA) and across 55 breast cancer cell lines from the Integrative Cancer Biology Program (ICBP43). These signatures were investigated for their relationship to pro- and anti-apoptotic protein expression and drug response in breast cancer cell lines. RESULTS: Application of these signatures to breast tumor gene expression data identified two novel discrete phenotypes characterized by concordant, aberrant activation of either the HER2, IGF1R, and AKT pathways ("the survival phenotype") or the EGFR, KRAS (G12V), RAF1, and BAD pathways ("the growth phenotype"). These phenotypes described a significant amount of the variability in the total expression data across breast cancer tumors and characterized distinctive patterns in apoptosis evasion and drug response. The growth phenotype expressed lower levels of BIM and higher levels of MCL-1 proteins. Further, the growth phenotype was more sensitive to common chemotherapies and targeted therapies directed at EGFR and MEK. Alternatively, the survival phenotype was more sensitive to drugs inhibiting HER2, PI3K, AKT, and mTOR, but more resistant to chemotherapies. CONCLUSIONS: Gene expression profiling revealed a bifurcation pattern in GFRN activity represented by two discrete phenotypes. These phenotypes correlate to unique mechanisms of apoptosis and drug response and have the potential of pinpointing targetable aberration(s) for more effective breast cancer treatments.

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.