Search

INTRODUCTION: Most patients with HPV-positive oropharyngeal squamous cell carcinoma (OPSCC) have an excellent response to chemoradiation, and trials are now investigating de-escalated treatment. However, up to 25% of patients with HPV-positive OPSCC will experience recurrence, and up to 5% will even progress through primary treatment. Currently, there are no molecular markers to identify patients with poor prognosis who would be harmed by de-escalation. Herein we report the clinical and genomic characteristics of persistent HPV-positive OPSCC after definitive platinum-based chemoradiation therapy. METHODS: Patients with HPV-positive OPSCC treated with curative intent platinum-based chemoradiation between 2007 and 2017 at two institutions and with a persistent locoregional disease were included. We evaluated clinical characteristics, including smoking status, age, stage, treatment, and overall survival. A subset of five patients had tissue available for targeted exome DNA sequencing and RNA sequencing. Genomic analysis was compared to a previously published cohort of 47 treatment-responsive HPV+ OPSCC tumors after batch correction. Mutational landscape, pathway activation, and OncoGPS tumor states were employed to characterize these tumors. RESULTS: Ten patients met the inclusion criteria. The tumor and nodal stages ranged from T1 to T4 and N1 to N2 by AJCC 8th edition staging. All patients were p16-positive by immunohistochemistry, and eight with available in situ hybridization were confirmed to be HPV-positive. The 1-year overall survival from the time of diagnosis was 57%, and the 2-year overall survival was 17%. TP53 mutations were present in three of five (60%) persistent tumors compared to 2% (one of 47) of treatment-responsive HPV-positive tumors (p = 0.008). Other genes with recurrent mutations in persistent HPV-positive OPSCC tumors were NF1, KMT2D, PIK3C2B, and TFGBR2. Compared to treatment-responsive HPV-positive tumors, persistent tumors demonstrated activation of DNA Repair and p53, EMT, MYC, SRC, and TGF-beta signaling pathways, with post-treatment samples demonstrating significant activation of the PI3K-EMT-Stem pathways compared to pretreatment samples. CONCLUSION: Chemoradiation-resistant HPV-positive OPSCC occurs infrequently but portends a poor prognosis. These tumors demonstrate higher rates of p53 mutation and activation of MYC, SRC, and TGF-beta pathways. A comparison of tumors before and after treatment demonstrates PI3K-EMT-Stem pathways post-treatment in HPV-positive tumors with persistent disease after platinum-based chemoradiation.

β-arrestins play a key role in G protein-coupled receptor (GPCR) internalization, trafficking, and signaling. Whether β-arrestins act independently of G protein-mediated signaling has not been fully elucidated. Studies using genome-editing approaches revealed that whereas G proteins are essential for mitogen-activated protein kinase activation by GPCRs., β-arrestins play a more prominent role in signal compartmentalization. However, in the absence of G proteins, GPCRs may not activate β-arrestins, thereby limiting the ability to distinguish G protein from β-arrestin-mediated signaling events. We used β2-adrenergic receptor (β2AR) and its β2AR-C tail mutant expressed in human embryonic kidney 293 cells wildtype or CRISPR-Cas9 gene edited for Gαs, β-arrestin1/2, or GPCR kinases 2/3/5/6 in combination with arrestin conformational sensors to elucidate the interplay between Gαs and β-arrestins in controlling gene expression. We found that Gαs is not required for β2AR and β-arrestin conformational changes, β-arrestin recruitment, and receptor internalization, but that Gαs dictates the GPCR kinase isoforms involved in β-arrestin recruitment. By RNA-Seq analysis, we found that protein kinase A and mitogen-activated protein kinase gene signatures were activated by stimulation of β2AR in wildtype and β-arrestin1/2-KO cells but absent in Gαs-KO cells. These results were validated by re-expressing Gαs in the corresponding KO cells and silencing β-arrestins in wildtype cells. These findings were extended to cellular systems expressing endogenous levels of β2AR. Overall, our results support that Gs is essential for β2AR-promoted protein kinase A and mitogen-activated protein kinase gene expression signatures, whereas β-arrestins initiate signaling events modulating Gαs-driven nuclear transcriptional activity.

Pancreatic ductal adenocarcinoma (PDAC) is an insidious disease with a low 5-year survival rate. PDAC is characterized by infiltration of abundant tumor-associated macrophages (TAM), which promote immune tolerance and immunotherapeutic resistance. Here we report that macrophage spleen tyrosine kinase (Syk) promotes PDAC growth and metastasis. In orthotopic PDAC mouse models, genetic deletion of myeloid Syk reprogrammed macrophages into immunostimulatory phenotype, increased the infiltration, proliferation, and cytotoxicity of CD8+ T cells, and repressed PDAC growth and metastasis. Furthermore, gemcitabine (Gem) treatment induced an immunosuppressive microenvironment in PDAC by promoting protumorigenic polarization of macrophages. In contrast, treatment with the FDA-approved Syk inhibitor R788 (fostamatinib) remodeled the tumor immune microenvironment, "re-educated" protumorigenic macrophages towards an immunostimulatory phenotype and boosted CD8+ T-cell responses in Gem-treated PDAC in orthotopic mouse models and an ex vivo human pancreatic slice culture model. These findings illustrate the potential of Syk inhibition for enhancing the antitumor immune responses in PDAC and support the clinical evaluation of R788 either alone or together with Gem as a potential treatment strategy for PDAC.SignificanceSyk blockade induces macrophage polarization to an immunostimulatory phenotype, which enhances CD8+ T-cell responses and improves gemcitabine efficacy in pancreatic ductal adenocarcinoma, a clinically challenging malignancy.

Cell-cycle control is accomplished by cyclin-dependent kinases (CDKs), motivating extensive research into CDK targeting small-molecule drugs as cancer therapeutics. Here we use combinatorial CRISPR/Cas9 perturbations to uncover an extensive network of functional interdependencies among CDKs and related factors, identifying 43 synthetic-lethal and 12 synergistic interactions. We dissect CDK perturbations using single-cell RNAseq, for which we develop a novel computational framework to precisely quantify cell-cycle effects and diverse cell states orchestrated by specific CDKs. While pairwise disruption of CDK4/6 is synthetic-lethal, only CDK6 is required for normal cell-cycle progression and transcriptional activation. Multiple CDKs (CDK1/7/9/12) are synthetic-lethal in combination with PRMT5, independent of cell-cycle control. In-depth analysis of mRNA expression and splicing patterns provides multiple lines of evidence that the CDK-PRMT5 dependency is due to aberrant transcriptional regulation resulting in premature termination. These inter-dependencies translate to drug-drug synergies, with therapeutic implications in cancer and other diseases.

BackgroundNeuroblastoma (NB) is considered an immunologically cold tumor and is usually less responsive to immune checkpoint blockade (ICB). Tumor-associated macrophages (TAMs) are highly infiltrated in NB tumors and promote immune escape and resistance to ICB. Hence therapeutic strategies targeting immunosuppressive TAMs can improve responses to ICB in NB. We recently discovered that spleen tyrosine kinase (Syk) reprograms TAMs toward an immunostimulatory phenotype and enhances T-cell responses in the lung adenocarcinoma model. Here we investigated if Syk is an immune-oncology target in NB and tested whether a novel immunotherapeutic approach utilizing Syk inhibitor together with radiation and ICB could provide a durable anti-tumor immune response in an MYCN amplified murine model of NB.MethodsMyeloid Syk KO mice and syngeneic MYCN-amplified cell lines were used to elucidate the effect of myeloid Syk on the NB tumor microenvironment (TME). In addition, the effect of Syk inhibitor, R788, on anti-tumor immunity alone or in combination with anti-PDL1 mAb and radiation was also determined in murine NB models. The underlying mechanism of action of this novel therapeutic combination was also investigated.ResultsHerein, we report that Syk is a marker of NB-associated macrophages and plays a crucial role in promoting immunosuppression in the NB TME. We found that the blockade of Syk in NB-bearing mice markedly impairs tumor growth. This effect is facilitated by macrophages that become immunogenic in the absence of Syk, skewing the suppressive TME towards immunostimulation and activating anti-tumor immune responses. Moreover, combining FDA-approved Syk inhibitor, R788 (fostamatinib) along with anti-PDL1 mAb provides a synergistic effect leading to complete tumor regression and durable anti-tumor immunity in mice bearing small tumors (50 mm3) but not larger tumors (250 mm3). However, combining radiation to R788 and anti-PDL1 mAb prolongs the survival of mice bearing large NB9464 tumors.ConclusionCollectively, our findings demonstrate the central role of macrophage Syk in NB progression and demonstrate that Syk blockade can "reeducate" TAMs towards immunostimulatory phenotype, leading to enhanced T cell responses. These findings further support the clinical evaluation of fostamatinib alone or with radiation and ICB, as a novel therapeutic intervention in neuroblastoma.

BackgroundCDK12 inactivation leading to increased neoantigen burden has been hypothesized to sensitize tumors to immune checkpoint inhibition. Pan-cancer data regarding the frequency of CDK12 alterations are limited. We aimed to characterize CDK12 alterations across all cancer types through real-world clinical-grade sequencing.MethodsThis was a single-center retrospective analysis of 4994 cancer patients who underwent tissue or blood genomic profiling, including CDK12 assessment, conducted as part of routine care from December 2012 to January 2020. Prevalence, clinical characteristics, and treatment outcomes of patients with tumors with pathogenic CDK12 alterations were described.ResultsIn all, 39 (0.78%, n = 39/4994) patients had pathogenic CDK12 alterations. Among CDK12-altered tumors, the most common organ site was prostate (n = 9, 23.1%) followed by colorectal (n = 5, 12.8%). Adenocarcinoma was the most common histology (n = 26, 66.7%). Median follow-up from time of diagnosis was 4.02 years. Median overall survival from time of metastasis was 4.43 years (95% CI: 3.11-5.74). Ten patients with CDK12-altered tumors received at least one immune checkpoint inhibitor-containing regimen. The majority of patients (n = 6/10, 60%) experienced an objective response. Progression-free survival for patients who had metastatic disease and received a checkpoint inhibitor-containing regimen was 1.16 years (95% CI: 0.32-2.00).ConclusionCDK12 alterations are rare events across hematologic and solid tumor malignancies. They represent a clinically distinct molecular cancer subtype which may have increased responsiveness to checkpoint inhibition. Prospective studies are warranted to investigate checkpoint inhibition in CDK12-altered tumors.

BackgroundThe role of the inflammatory milieu in prostate cancer progression is not well understood. Differences in inflammatory signaling between localized and metastatic disease may point to opportunities for early intervention.MethodsWe modeled PCa disease progression by analyzing RNA-seq of localized vs. metastatic patient samples, followed by CIBERSORTx to assess their immune cell populations. The VHA CDW registry of PCa patients was analyzed for anti-TNF clinical outcomes.ResultsWe observed statistically significant opposing patterns of IL-6 and TNFα expression between localized and metastatic disease. IL-6 was robustly expressed in localized disease and downregulated in metastatic disease. The reverse was observed with TNFα expression. Metastatic disease was also characterized by downregulation of adhesion molecule E-selectin, matrix metalloproteinase ADAMTS-4 and a shift to M2 macrophages whereas localized disease demonstrated a preponderance of M1 macrophages. Treatment with anti-TNF agents was associated with earlier stage disease at diagnosis.ConclusionsOur data points to clearly different inflammatory contexts between localized and metastatic prostate cancer. Primary localized disease demonstrates local inflammation and adaptive immunity, whereas metastases are characterized by immune cold microenvironments and a shift towards resolution of inflammation and tissue repair. Therapies that interfere with these inflammatory networks may offer opportunities for early intervention in monotherapy or in combination with immunotherapies and anti-angiogenic approaches.

Natural killer (NK) cells are known to mediate killing of various cancer types, but tumor cells can develop resistance mechanisms to escape NK cell-mediated killing. Here, we use a "two cell type" whole genome CRISPR-Cas9 screening system to discover key regulators of tumor sensitivity and resistance to NK cell-mediated cytotoxicity in human glioblastoma stem cells (GSC). We identify CHMP2A as a regulator of GSC resistance to NK cell-mediated cytotoxicity and we confirm these findings in a head and neck squamous cells carcinoma (HNSCC) model. We show that deletion of CHMP2A activates NF-κB in tumor cells to mediate increased chemokine secretion that promotes NK cell migration towards tumor cells. In the HNSCC model we demonstrate that CHMP2A mediates tumor resistance to NK cells via secretion of extracellular vesicles (EVs) that express MICA/B and TRAIL. These secreted ligands induce apoptosis of NK cells to inhibit their antitumor activity. To confirm these in vitro studies, we demonstrate that deletion of CHMP2A in CAL27 HNSCC cells leads to increased NK cell-mediated killing in a xenograft immunodeficient mouse model. These findings illustrate a mechanism of tumor immune escape through EVs secretion and identify inhibition of CHMP2A and related targets as opportunities to improve NK cell-mediated immunotherapy.

With an increasing number of patients with degenerative hepatic diseases, such as liver fibrosis, and a limited supply of donor organs, there is an unmet need for therapies that can repair or regenerate damaged liver tissue. Treatment with macrophages that are capable of phagocytosis and anti-inflammatory activities such as secretion of matrix metalloproteinases (MMPs) provide an attractive cellular therapy approach. Human induced pluripotent stem cells (iPSCs) are capable of efficiently generating a large-scale, homogenous population of human macrophages using fully defined feeder- and serum-free differentiation protocol. Human iPSC-macrophages exhibit classical surface cell markers and phagocytic activity similar to peripheral blood-derived macrophages. Moreover, gene and cytokine expression analysis reveal that these macrophages can be efficiently polarized to pro-inflammatory M1 or anti-inflammatory M2 phenotypes in presence of LPS + IFN-γ and IL-4 + IL-13, respectively. M1 macrophages express high level of CD80, TNF-α, and IL-6 while M2 macrophages show elevated expression of CD206, CCL17, and CCL22. Here, we demonstrate that treatment of liver fibrosis with both human iPSC-derived macrophage populations and especially M2 subtype significantly reduces fibrogenic gene expression and disease associated histological markers including Sirius Red, αSMA and desmin in immunodeficient Rag2-/- γc-/- mice model, making this approach a promising cell-based avenue to ameliorate fibrosis.

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

MDSCs are immune cells of myeloid lineage that plays a key role in promoting tumor growth. The expansion of MDSCs in tumor-bearing hosts reduces the efficacy of checkpoint inhibitors and CAR-T therapies, and hence strategies that deplete or block the recruitment of MDSCs have shown benefit in improving responses to immunotherapy in various cancers, including NB. Ibrutinib, an irreversible molecular inhibitor of BTK, has been widely studied in B cell malignancies, and recently, this drug is repurposed for the treatment of solid tumors. Herein we report that BTK is highly expressed in both granulocytic and monocytic murine MDSCs isolated from mice bearing NB tumors, and its increased expression correlates with a poor relapse-free survival probability of NB patients. Moreover, in vitro treatment of murine MDSCs with ibrutinib altered NO production, decreased mRNA expression of Ido, Arg, Tgfβ, and displayed defects in T-cell suppression. Consistent with these findings, in vivo inhibition of BTK with ibrutinib resulted in reduced MDSC-mediated immune suppression, increased CD8+ T cell infiltration, decreased tumor growth, and improved response to anti-PDL1 checkpoint inhibitor therapy in a murine model of NB. These results demonstrate that ibrutinib modulates immunosuppressive functions of MDSC and can be used either alone or in combination with immunotherapy for augmenting antitumor immune responses in NB.

In adult tissue, stem and progenitor cells must tightly regulate the balance between proliferation and differentiation to sustain homeostasis. How this exquisite balance is achieved is an area of active investigation. Here, we show that epidermal genes, including ~30% of induced differentiation genes already contain stalled Pol II at the promoters in epidermal stem and progenitor cells which is then released into productive transcription elongation upon differentiation. Central to this process are SPT6 and PAF1 which are necessary for the elongation of these differentiation genes. Upon SPT6 or PAF1 depletion there is a loss of human skin differentiation and stratification. Unexpectedly, loss of SPT6 also causes the spontaneous transdifferentiation of epidermal cells into an intestinal-like phenotype due to the stalled transcription of the master regulator of epidermal fate P63. Our findings suggest that control of transcription elongation through SPT6 plays a prominent role in adult somatic tissue differentiation and the inhibition of alternative cell fate choices.

Medulloblastoma is among the most common malignant brain tumors in children. Recent studies have identified at least four subgroups of the disease that differ in terms of molecular characteristics and patient outcomes. Despite this heterogeneity, most patients with medulloblastoma receive similar therapies, including surgery, radiation, and intensive chemotherapy. Although these treatments prolong survival, many patients still die from the disease and survivors suffer severe long-term side effects from therapy. We hypothesize that each patient with medulloblastoma is sensitive to different therapies and that tailoring therapy based on the molecular and cellular characteristics of patients' tumors will improve outcomes. To test this, we assembled a panel of orthotopic patient-derived xenografts (PDX) and subjected them to DNA sequencing, gene expression profiling, and high-throughput drug screening. Analysis of DNA sequencing revealed that most medulloblastomas do not have actionable mutations that point to effective therapies. In contrast, gene expression and drug response data provided valuable information about potential therapies for every tumor. For example, drug screening demonstrated that actinomycin D, which is used for treatment of sarcoma but rarely for medulloblastoma, was active against PDXs representing Group 3 medulloblastoma, the most aggressive form of the disease. Functional analysis of tumor cells was successfully used in a clinical setting to identify more treatment options than sequencing alone. These studies suggest that it should be possible to move away from a one-size-fits-all approach and begin to treat each patient with therapies that are effective against their specific tumor. SIGNIFICANCE: These findings show that high-throughput drug screening identifies therapies for medulloblastoma that cannot be predicted by genomic or transcriptomic analysis.

The dominant paradigm for HPV carcinogenesis includes integration into the host genome followed by expression of E6 and E7 (E6/E7). We explored an alternative carcinogenic pathway characterized by episomal E2, E4, and E5 (E2/E4/E5) expression. Half of HPV positive cervical and pharyngeal cancers comprised a subtype with increase in expression of E2/E4/E5, as well as association with lack of integration into the host genome. Models of the E2/E4/E5 carcinogenesis show p53 dependent enhanced proliferation in vitro, as well as increased susceptibility to induction of cancer in vivo. Whole genomic expression analysis of the E2/E4/E5 pharyngeal cancer subtype is defined by activation of the fibroblast growth factor receptor (FGFR) pathway and this subtype is susceptible to combination FGFR and mTOR inhibition, with implications for targeted therapy.

In a substantial fraction of cancers TERT promoter (TERTp) mutations drive expression of the catalytic subunit of telomerase, contributing to their proliferative immortality. We conducted a pan-cancer analysis of cell lines and find a TERTp mutation expression signature dominated by epithelial-to-mesenchymal transition and MAPK signaling. These data indicate that TERTp mutants are likely to generate distinctive tumor microenvironments and intercellular interactions. Analysis of high-throughput screening tests of 546 small molecules on cell line growth indicated that TERTp mutants displayed heightened sensitivity to specific drugs, including RAS pathway inhibitors, and we found that inhibition of MEK1 and 2, key RAS/MAPK pathway effectors, inhibited TERT mRNA expression. Consistent with an enrichment of mesenchymal states in TERTp mutants, cell lines and some patient tumors displayed low expression of the central adherens junction protein E-cadherin, and we provide evidence that its expression in these cells is regulated by MEK1/2. Several mesenchymal transcription factors displayed elevated expression in TERTp mutants including ZEB1 and 2, TWIST1 and 2, and SNAI1. Of note, the developmental transcription factor SNAI2/SLUG was conspicuously elevated in a significant majority of TERTp-mutant cell lines, and knock-down experiments suggest that it promotes TERT expression. IMPLICATIONS: Cancers harboring TERT promoter mutations are often more lethal, but the basis for this higher mortality remains unknown. Our study identifies that TERTp mutants, as a class, associate with a distinct gene and protein expression signature likely to impact their biological and clinical behavior and provide new directions for investigating treatment approaches for these cancers.

PurposeHuman papillomavirus (HPV)-related head and neck squamous cell carcinoma (HNSCC) is associated with daily marijuana use and is also increasing in parallel with increased marijuana use in the United States. Our study is designed to define the interaction between cannabinoids and HPV-positive HNSCC.Experimental designThe expression of cannabinoid receptors CNR1 and CNR2 was analyzed using The Cancer Genome Atlas (TCGA) HNSCC data. We used agonists, antagonists, siRNAs, or shRNA-based models to explore the roles of CNR1 and CNR2 in HPV-positive HNSCC cell lines and animal models. Cannabinoid downstream pathways involved were determined by Western blotting and analyzed in a primary HPV HNSCC cohort with single-sample gene set enrichment analysis (ssGSEA) and the OncoGenome Positioning System (Onco-GPS).ResultsIn TCGA cohort, the expression of CNR1 and CNR2 was elevated in HPV-positive HNSCC compared with HPV-negative HNSCC, and knockdown of CNR1/CNR2 expression inhibited proliferation in HPV-positive HNSCC cell lines. Specific CNR1 and CNR2 activation as well as nonselective cannabinoid receptor activation in cell lines and animal models promoted cell growth, migration, and inhibited apoptosis through p38 MAPK pathway activation. CNR1/CNR2 antagonists suppressed cell proliferation and migration and induced apoptosis. Using whole-genome expression analysis in a primary HPV HNSCC cohort, we identified specific p38 MAPK pathway activation signature in tumors from HPV HNSCC patients with objective measurement of concurrent cannabinoid exposure.ConclusionsCannabinoids can promote progression of HPV-positive HNSCC through p38 MAPK pathway activation.

Macrophages (MΦ) play a critical role in tumor growth, immunosuppression, and inhibition of adaptive immune responses in cancer. Hence, targeting signaling pathways in MΦs that promote tumor immunosuppression will provide therapeutic benefit. PI3Kγ has been recently established by our group and others as a novel immuno-oncology target. Herein, we report that an MΦ Syk-PI3K axis drives polarization of immunosuppressive MΦs that establish an immunosuppressive tumor microenvironment in in vivo syngeneic tumor models. Genetic or pharmacologic inhibition of Syk and/or PI3Kγ in MΦs promotes a proinflammatory MΦ phenotype, restores CD8+ T-cell activity, destabilizes HIF under hypoxia, and stimulates an antitumor immune response. Assay for transposase-accessible Chromatin using Sequencing (ATAC-seq) analyses on the bone marrow-derived macrophages (BMDM) show that inhibition of Syk kinase promotes activation and binding of NF-κB motif in SykMC-KO BMDMs, thus stimulating immunostimulatory transcriptional programming in MΦs to suppress tumor growth. Finally, we have developed in silico the "first-in-class" dual Syk/PI3K inhibitor, SRX3207, for the combinatorial inhibition of Syk and PI3K in one small molecule. This chemotype demonstrates efficacy in multiple tumor models and represents a novel combinatorial approach to activate antitumor immunity.

Neuroblastoma is the most common extracranial solid tumor of childhood and accounts for 15% of all pediatric cancer-related deaths. New therapies are needed to improve outcomes for children with high-risk and relapsed tumors. Inhibitors of the RET kinase and the RAS-MAPK pathway have previously been shown to be effective against neuroblastoma, suggesting that combined inhibition may have increased efficacy. RXDX-105 is a small molecule inhibitor of multiple kinases, including the RET and BRAF kinases. We found that treatment of neuroblastoma cells with RXDX-105 resulted in a significant decrease in cell viability and proliferation in vitro and in tumor growth and tumor vascularity in vivo. Treatment with RXDX-105 inhibited RET phosphorylation and phosphorylation of the MEK and ERK kinases in neuroblastoma cells and xenograft tumors, and RXDX-105 treatment induced both apoptosis and cell cycle arrest. RXDX-105 also showed enhanced efficacy in combination with 13-cis-retinoic acid, which is currently a component of maintenance therapy for children with high-risk neuroblastoma. Our results demonstrate that RXDX-105 shows promise as a novel therapeutic agent for children with high-risk and relapsed neuroblastoma.

Recurrent mutations implicate several epigenetic regulators in the early molecular pathobiology of myelodysplastic syndromes (MDS). We hypothesized that MDS subtypes defined by DNA methylation (DNAm) patterns could enhance our understanding of MDS disease biology and identify patients with convergent epigenetic profiles. Bisulfite padlock probe sequencing was used to measure DNAm of ∼500 000 unique cytosine guanine dinucleotides covering 140 749 nonoverlapping regulatory regions across the genome in bone marrow DNA samples from 141 patients with MDS. Application of a nonnegative matrix factorization (NMF)-based decomposition of DNAm profiles identified 5 consensus clusters described by 5 NMF components as the most stable grouping solution. Each of the 5 NMF components identified by this approach correlated with specific genetic abnormalities and categorized patients into 5 distinct methylation clusters, each largely defined by a single NMF component. Methylation clusters displayed unique differentially methylated regulatory loci enriched for active and bivalent promoters and enhancers. Two clusters were enriched for samples with complex karyotypes, although only one had an increased number of TP53 mutations. Each of the 3 most frequently mutated splicing factors, SF3B1, U2AF1, and SRSF2, was enriched in different clusters. Mutations of ASXL1, EZH2, and RUNX1 were coenriched in the SRSF2-containing cluster. In multivariate analysis, methylation cluster membership remained independently associated with overall survival. Targeted DNAm profiles identify clinically relevant subtypes of MDS not otherwise distinguished by mutations or clinical features. Patients with diverse genetic lesions can converge on common DNAm states with shared pathogenic mechanisms and clinical outcomes.

Metformin may reduce the progression of head and neck squamous cell carcinoma (HNSCC); however, whether metformin acts by altering the host metabolism or targets cancer-initiating cells remains poorly understood. This gap in knowledge has prevented the stratification of patient populations who are most likely to benefit from metformin treatment. Here, we explored whether metformin acts directly on HNSCC cells to inhibit aberrant cell growth. To investigate the tumor cell autonomous effects of metformin, we engineered representative HPV- and HPV+ HNSCC cells harboring typical genetic alternations to express the yeast mitochondrial NADH dehydrogenase (NDI1) protein, which is insensitive to metformin. NDI1 expression rescued the inhibitory effects of metformin on mitochondrial complex I, abolished the ability of metformin to activate AMP-activated protein kinase, and inhibited mTOR signaling both in vitro and in vivo, and was sufficient to render metformin ineffective to prevent HNSCC tumor growth. This experimental system provided an opportunity to identify metformin-regulated transcriptional programs linked to cancer cell growth inhibition in the tumor microenvironment. Remarkably, computational analysis of the metformin-induced transcriptome revealed that metformin downregulated gene expression signatures associated with cancer stemness and epithelial-mesenchymal transition, concomitant with increased expression of squamous differentiation genes. These findings support that metformin may act directly on cancer-initiating cells to prevent their progression to HNSCC, which may inform the selection of patients at risk of developing HNSCC in future early-stage clinical trials. SIGNIFICANCE: Metformin's ability to directly target HNSCC-initiating cells instead of exerting cancer preventive activity based solely on its systemic effects may inform the selection of patients in future precision prevention trials.

Despite recent advances in targeted drugs and immunotherapy, cancer remains "the emperor of all maladies" due to almost inevitable emergence of resistance. Drug resistance is thought to be driven by genetic alterations and/or dynamic plasticity that deregulate pathway activities and regulatory programs of a highly heterogeneous tumour. In this study, we propose a modelling framework to simulate population dynamics of heterogeneous tumour cells with reversible drug resistance. Drug sensitivity of a tumour cell is determined by its internal states, which are demarcated by coordinated activities of multiple interconnected oncogenic pathways. Transitions between cellular states depend on the effects of targeted drugs and regulatory relations between the pathways. Under this framework, we build a simple model to capture drug resistance characteristics of BRAF-mutant melanoma, where two cell states are determined by two mutually inhibitory - main and alternative - pathways. We assume that cells with an activated main pathway are proliferative yet sensitive to the BRAF inhibitor, and cells with an activated alternative pathway are quiescent but resistant to the drug. We describe a dynamical process of tumour growth under various drug regimens using the explicit solutions of mean-field equations. Based on these solutions, we compare efficacy of three treatment strategies from simulated data: static treatments with continuous and constant dosages, periodic treatments with regular intermittent active phases and drug holidays, and treatments derived from optimal control theory (OCT). Periodic treatments outperform static treatments with a considerable margin, while treatments based on OCT outperform the best periodic treatment. Our results provide insights regarding optimal cancer treatment modalities for heterogeneous tumours, and may guide the development of optimal therapeutic strategies to circumvent plastic drug resistance. They can also be used to evaluate the efficacy of suboptimal treatments that may account for side effects of the treatment and the cost of its application.

G protein-coupled receptors (GPCRs) are the largest gene family of cell membrane-associated molecules mediating signal transmission, and their involvement in key physiological functions is well-established. The ability of GPCRs to regulate a vast array of fundamental biological processes, such as cardiovascular functions, immune responses, hormone and enzyme release from endocrine and exocrine glands, neurotransmission, and sensory perception (e.g. vision, odor, and taste), is largely due to the diversity of these receptors and the layers of their downstream signaling circuits. Dysregulated expression and aberrant functions of GPCRs have been linked to some of the most prevalent human diseases, which renders GPCRs one of the top targets for pharmaceutical drug development. However, the study of the role of GPCRs in tumor biology has only just begun to make headway. Recent studies have shown that GPCRs can contribute to the many facets of tumorigenesis, including proliferation, survival, angiogenesis, invasion, metastasis, therapy resistance, and immune evasion. Indeed, GPCRs are widely dysregulated in cancer and yet are underexploited in oncology. We present here a comprehensive analysis of GPCR gene expression, copy number variation, and mutational signatures in 33 cancer types. We also highlight the emerging role of GPCRs as part of oncocrine networks promoting tumor growth, dissemination, and immune evasion, and we stress the potential benefits of targeting GPCRs and their signaling circuits in the new era of precision medicine and cancer immunotherapies.

The original version of this Article contained an error in the author affiliations. Trey Ideker was incorrectly associated with 'Department of Medicine (Oncology), Stanford University School of Medicine, 875 Blake Wilbur Dr, Palo Alto, CA 94304, USA.' This has now been corrected in both the PDF and HTML versions of the Article.

Although promoter-associated CpG islands have been established as targets of DNA methylation changes in cancer, previous studies suggest that epigenetic dysregulation outside the promoter region may be more closely associated with transcriptional changes. Here we examine DNA methylation, chromatin marks, and transcriptional alterations to define the relationship between transcriptional modulation and spatial changes in chromatin structure. Using human papillomavirus-related oropharyngeal carcinoma as a model, we show aberrant enrichment of repressive H3K9me3 at the transcriptional start site (TSS) with methylation-associated, tumor-specific gene silencing. Further analysis identifies a hypermethylated subtype which shows a functional convergence on MYC targets and association with CREBBP/EP300 mutation. The tumor-specific shift to transcriptional repression associated with DNA methylation at TSSs was confirmed in multiple tumor types. Our data may show a common underlying epigenetic dysregulation in cancer associated with broad enrichment of repressive chromatin marks and aberrant DNA hypermethylation at TSSs in combination with MYC network activation.

Tumor-associated neo-antigens are mutated peptides that allow the immune system to recognize the affected cell as foreign. Cells carrying excessive mutation load often develop mechanisms of tolerance. PD-L1/PD-1 checkpoint immunotherapy is a highly promising approach to overcome these protective signals and induce tumor shrinkage. Yet, the nature of the neo-antigens driving those beneficial responses remains unclear. Here, we show that APOBEC-related mutagenesis - a mechanism at the crossroads between anti-viral immunity and endogenous nucleic acid editing - increases neo-peptide hydrophobicity (a feature of immunogenicity), as demonstrated by in silico computation and in the TCGA pan-cancer cohort, where APOBEC-related mutagenesis was also strongly associated with immune marker expression. Moreover, APOBEC-related mutagenesis correlated with immunotherapy response in a cohort of 99 patients with diverse cancers, and this correlation was independent of the tumor mutation burden (TMB). Combining APOBEC-related mutagenesis estimate and TMB resulted in greater predictive ability than either parameter alone. Based on these results, further investigation of APOBEC-related mutagenesis as a marker of response to anti-cancer checkpoint blockade is warranted.

Activating mutations in GNAQ/GNA11, encoding Gαq G proteins, are initiating oncogenic events in uveal melanoma (UM). However, there are no effective therapies for UM. Using an integrated bioinformatics pipeline, we found that PTK2, encoding focal adhesion kinase (FAK), represents a candidate synthetic lethal gene with GNAQ activation. We show that Gαq activates FAK through TRIO-RhoA non-canonical Gαq-signaling, and genetic ablation or pharmacological inhibition of FAK inhibits UM growth. Analysis of the FAK-regulated transcriptome demonstrated that GNAQ stimulates YAP through FAK. Dissection of the underlying mechanism revealed that FAK regulates YAP by tyrosine phosphorylation of MOB1, inhibiting core Hippo signaling. Our findings establish FAK as a potential therapeutic target for UM and other Gαq-driven pathophysiologies that involve unrestrained YAP function.

Background: The translation of genomic discoveries to the clinic is the cornerstone of precision medicine. However, incorporating next generation sequencing (NGS) of hematologic malignancies into clinical management remains limited. Methods: We describe 235 patients who underwent integrated NGS profiling (406 genes) and analyze the alterations and their potential actionability. Results: Overall, 227 patients (96.5%) had adequate tissue. Most common diagnoses included myelodysplastic syndrome (22.9%), chronic lymphocytic leukemia (17.2%), non-Hodgkin lymphoma (13.2%), acute myeloid leukemia (11%), myeloproliferative neoplasm (9.2%), acute lymphoblastic leukemia (8.8%), and multiple myeloma (7.5%). Most patients (N = 197/227 (87%)) harbored ≥1 genomic alteration(s); 170/227 (75%), ≥1 potentially actionable alteration(s) targetable by an FDA-approved (mostly off-label) or an investigational agent. Altogether, 546 distinct alterations were seen, most commonly involving TP53 (10.8%), TET2 (4.6%), and DNMT3A (4.2%). The median tumor mutational burden (TMB) was low (1.7 alterations/megabase); 12% of patients had intermediate or high TMB (higher TMB correlates with favorable response to anti-PD1/PDL1 inhibition in solid tumors). In conclusion, 96.5% of patients with hematologic malignancies have adequate tissue for comprehensive genomic profiling. Most patients had unique molecular signatures, and 75% had alterations that may be pharmacologically tractable with gene- or immune-targeted agents.

Background: The translation of genomic discoveries to the clinic is the cornerstone of precision medicine. However, incorporating next generation sequencing (NGS) of hematologic malignancies into clinical management remains limited. Methods: We describe 235 patients who underwent integrated NGS profiling (406 genes) and analyze the alterations and their potential actionability. Results: Overall, 227 patients (96.5%) had adequate tissue. Most common diagnoses included myelodysplastic syndrome (22.9%), chronic lymphocytic leukemia (17.2%), non-Hodgkin lymphoma (13.2%), acute myeloid leukemia (11%), myeloproliferative neoplasm (9.2%), acute lymphoblastic leukemia (8.8%), and multiple myeloma (7.5%). Most patients (N = 197/227 (87%)) harbored ≥1 genomic alteration(s); 170/227 (75%), ≥1 potentially actionable alteration(s) targetable by an FDA-approved (mostly off-label) or an investigational agent. Altogether, 546 distinct alterations were seen, most commonly involving TP53 (10.8%), TET2 (4.6%), and DNMT3A (4.2%). The median tumor mutational burden (TMB) was low (1.7 alterations/megabase); 12% of patients had intermediate or high TMB (higher TMB correlates with favorable response to anti-PD1/PDL1 inhibition in solid tumors). In conclusion, 96.5% of patients with hematologic malignancies have adequate tissue for comprehensive genomic profiling. Most patients had unique molecular signatures, and 75% had alterations that may be pharmacologically tractable with gene- or immune-targeted agents.

High-throughput single-cell gene expression profiling has enabled the definition of new cell types and developmental trajectories. Visualizing these datasets is crucial to biological interpretation, and a popular method is t-stochastic neighbor embedding (t-SNE), which visualizes local patterns well but distorts global structure, such as distances between clusters. We developed similarity weighted nonnegative embedding (SWNE), which enhances interpretation of datasets by embedding the genes and factors that separate cell states on the visualization alongside the cells and maintains fidelity when visualizing local and global structure for both developmental trajectories and discrete cell types. SWNE uses nonnegative matrix factorization to decompose the gene expression matrix into biologically relevant factors; embeds the cells, genes, and factors in a 2D visualization; and uses a similarity matrix to smooth the embeddings. We demonstrate SWNE on single-cell RNA-seq data from hematopoietic progenitors and human brain cells. SWNE is available as an R package at github.com/yanwu2014/swne.

Despite extensive efforts, oncogenic KRAS remains resistant to targeted therapy. Combined downstream RAL-TBK1 and MEK inhibition induces only transient lung tumor shrinkage in KRAS-driven genetically engineered mouse models (GEMMs). Using the sensitive KRAS;LKB1 (KL) mutant background, we identify YAP1 upregulation and a therapy-induced secretome as mediators of acquired resistance. This program is reversible, associated with H3K27 promoter acetylation, and suppressed by BET inhibition, resensitizing resistant KL cells to TBK1/MEK inhibition. Constitutive YAP1 signaling promotes intrinsic resistance in KRAS;TP53 (KP) mutant lung cancer. Intermittent treatment with the BET inhibitor JQ1 thus overcomes resistance to combined pathway inhibition in KL and KP GEMMs. Using potent and selective TBK1 and BET inhibitors we further develop an effective therapeutic strategy with potential translatability to the clinic.

There is a pressing need to identify therapeutic targets in tumors with low mutation rates such as the malignant pediatric brain tumor medulloblastoma. To address this challenge, we quantitatively profiled global proteomes and phospho-proteomes of 45 medulloblastoma samples. Integrated analyses revealed that tumors with similar RNA expression vary extensively at the post-transcriptional and post-translational levels. We identified distinct pathways associated with two subsets of SHH tumors, and found post-translational modifications of MYC that are associated with poor outcomes in group 3 tumors. We found kinases associated with subtypes and showed that inhibiting PRKDC sensitizes MYC-driven cells to radiation. Our study shows that proteomics enables a more comprehensive, functional readout, providing a foundation for future therapeutic strategies.

Mesenchymal tumor subpopulations secrete pro-tumorigenic cytokines and promote treatment resistance1-4. This phenomenon has been implicated in chemorefractory small cell lung cancer and resistance to targeted therapies5-8, but remains incompletely defined. Here, we identify a subclass of endogenous retroviruses (ERVs) that engages innate immune signaling in these cells. Stimulated 3 prime antisense retroviral coding sequences (SPARCS) are oriented inversely in 3' untranslated regions of specific genes enriched for regulation by STAT1 and EZH2. Derepression of these loci results in double-stranded RNA generation following IFN-γ exposure due to bi-directional transcription from the STAT1-activated gene promoter and the 5' long terminal repeat of the antisense ERV. Engagement of MAVS and STING activates downstream TBK1, IRF3, and STAT1 signaling, sustaining a positive feedback loop. SPARCS induction in human tumors is tightly associated with major histocompatibility complex class 1 expression, mesenchymal markers, and downregulation of chromatin modifying enzymes, including EZH2. Analysis of cell lines with high inducible SPARCS expression reveals strong association with an AXL/MET-positive mesenchymal cell state. While SPARCS-high tumors are immune infiltrated, they also exhibit multiple features of an immune-suppressed microenviroment. Together, these data unveil a subclass of ERVs whose derepression triggers pathologic innate immune signaling in cancer, with important implications for cancer immunotherapy.

Dysregulation of the Hippo signaling pathway and the consequent YAP1 activation is a frequent event in human malignancies, yet the underlying molecular mechanisms are still poorly understood. A pancancer analysis of core Hippo kinases and their candidate regulating molecules revealed few alterations in the canonical Hippo pathway, but very frequent genetic alterations in the FAT family of atypical cadherins. By focusing on head and neck squamous cell carcinoma (HNSCC), which displays frequent FAT1 alterations (29.8%), we provide evidence that FAT1 functional loss results in YAP1 activation. Mechanistically, we found that FAT1 assembles a multimeric Hippo signaling complex (signalome), resulting in activation of core Hippo kinases by TAOKs and consequent YAP1 inactivation. We also show that unrestrained YAP1 acts as an oncogenic driver in HNSCC, and that targeting YAP1 may represent an attractive precision therapeutic option for cancers harboring genomic alterations in the FAT1 tumor suppressor genes.

Functional redundancy shared by paralog genes may afford protection against genetic perturbations, but it can also result in genetic vulnerabilities due to mutual interdependency1-5. Here, we surveyed genome-scale short hairpin RNA and CRISPR screening data on hundreds of cancer cell lines and identified MAGOH and MAGOHB, core members of the splicing-dependent exon junction complex, as top-ranked paralog dependencies6-8. MAGOHB is the top gene dependency in cells with hemizygous MAGOH deletion, a pervasive genetic event that frequently occurs due to chromosome 1p loss. Inhibition of MAGOHB in a MAGOH-deleted context compromises viability by globally perturbing alternative splicing and RNA surveillance. Dependency on IPO13, an importin-β receptor that mediates nuclear import of the MAGOH/B-Y14 heterodimer9, is highly correlated with dependency on both MAGOH and MAGOHB. Both MAGOHB and IPO13 represent dependencies in murine xenografts with hemizygous MAGOH deletion. Our results identify MAGOH and MAGOHB as reciprocal paralog dependencies across cancer types and suggest a rationale for targeting the MAGOHB-IPO13 axis in cancers with chromosome 1p deletion.

Gene networks are rapidly growing in size and number, raising the question of which networks are most appropriate for particular applications. Here, we evaluate 21 human genome-wide interaction networks for their ability to recover 446 disease gene sets identified through literature curation, gene expression profiling, or genome-wide association studies. While all networks have some ability to recover disease genes, we observe a wide range of performance with STRING, ConsensusPathDB, and GIANT networks having the best performance overall. A general tendency is that performance scales with network size, suggesting that new interaction discovery currently outweighs the detrimental effects of false positives. Correcting for size, we find that the DIP network provides the highest efficiency (value per interaction). Based on these results, we create a parsimonious composite network with both high efficiency and performance. This work provides a benchmark for selection of molecular networks in human disease research.

We report a novel computational method, RegNetDriver, to identify tumorigenic drivers using the combined effects of coding and non-coding single nucleotide variants, structural variants, and DNA methylation changes in the DNase I hypersensitivity based regulatory network. Integration of multi-omics data from 521 prostate tumor samples indicated a stronger regulatory impact of structural variants, as they affect more transcription factor hubs in the tissue-specific network. Moreover, crosstalk between transcription factor hub expression modulated by structural variants and methylation levels likely leads to the differential expression of target genes. We report known prostate tumor regulatory drivers and nominate novel transcription factors (ERF, CREB3L1, and POU2F2), which are supported by functional validation.

African-American men have the highest incidence of and mortality from prostate cancer. Whether a biological basis exists for this disparity remains unclear. Exome sequencing (n = 102) and targeted validation (n = 90) of localized primary hormone-naïve prostate cancer in African-American men identified several gene mutations not previously observed in this context, including recurrent loss-of-function mutations in ERF, an ETS transcriptional repressor, in 5% of cases. Analysis of existing prostate cancer cohorts revealed ERF deletions in 3% of primary prostate cancers and mutations or deletions in ERF in 3% to 5% of lethal castration-resistant prostate cancers. Knockdown of ERF confers increased anchorage-independent growth and generates a gene expression signature associated with oncogenic ETS activation and androgen signaling. Together, these results suggest that ERF is a prostate cancer tumor-suppressor gene. More generally, our findings support the application of systematic cancer genomic characterization in settings of broader ancestral diversity to enhance discovery and, eventually, therapeutic applications.Significance: Systematic genomic sequencing of prostate cancer in African-American men revealed new insights into prostate cancer, including the identification of ERF as a prostate cancer gene; somatic copy-number alteration differences; and uncommon PIK3CA and PTEN alterations. This study highlights the importance of inclusion of underrepresented minorities in cancer sequencing studies. Cancer Discov; 7(9); 973-83. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 920.

The systematic sequencing of the cancer genome has led to the identification of numerous genetic alterations in cancer. However, a deeper understanding of the functional consequences of these alterations is necessary to guide appropriate therapeutic strategies. Here, we describe Onco-GPS (OncoGenic Positioning System), a data-driven analysis framework to organize individual tumor samples with shared oncogenic alterations onto a reference map defined by their underlying cellular states. We applied the methodology to the RAS pathway and identified nine distinct components that reflect transcriptional activities downstream of RAS and defined several functional states associated with patterns of transcriptional component activation that associates with genomic hallmarks and response to genetic and pharmacological perturbations. These results show that the Onco-GPS is an effective approach to explore the complex landscape of oncogenic cellular states across cancers, and an analytic framework to summarize knowledge, establish relationships, and generate more effective disease models for research or as part of individualized precision medicine paradigms.

Interactive analysis notebook environments promise to streamline genomics research through interleaving text, multimedia, and executable code into unified, sharable, reproducible "research narratives." However, current notebook systems require programming knowledge, limiting their wider adoption by the research community. We have developed the GenePattern Notebook environment (http://www.genepattern-notebook.org), to our knowledge the first system to integrate the dynamic capabilities of notebook systems with an investigator-focused, easy-to-use interface that provides access to hundreds of genomic tools without the need to write code.

Plasticity of the cell state has been proposed to drive resistance to multiple classes of cancer therapies, thereby limiting their effectiveness. A high-mesenchymal cell state observed in human tumours and cancer cell lines has been associated with resistance to multiple treatment modalities across diverse cancer lineages, but the mechanistic underpinning for this state has remained incompletely understood. Here we molecularly characterize this therapy-resistant high-mesenchymal cell state in human cancer cell lines and organoids and show that it depends on a druggable lipid-peroxidase pathway that protects against ferroptosis, a non-apoptotic form of cell death induced by the build-up of toxic lipid peroxides. We show that this cell state is characterized by activity of enzymes that promote the synthesis of polyunsaturated lipids. These lipids are the substrates for lipid peroxidation by lipoxygenase enzymes. This lipid metabolism creates a dependency on pathways converging on the phospholipid glutathione peroxidase (GPX4), a selenocysteine-containing enzyme that dissipates lipid peroxides and thereby prevents the iron-mediated reactions of peroxides that induce ferroptotic cell death. Dependency on GPX4 was found to exist across diverse therapy-resistant states characterized by high expression of ZEB1, including epithelial-mesenchymal transition in epithelial-derived carcinomas, TGFβ-mediated therapy-resistance in melanoma, treatment-induced neuroendocrine transdifferentiation in prostate cancer, and sarcomas, which are fixed in a mesenchymal state owing to their cells of origin. We identify vulnerability to ferroptic cell death induced by inhibition of a lipid peroxidase pathway as a feature of therapy-resistant cancer cells across diverse mesenchymal cell-state contexts.

Medulloblastoma is the most frequent malignant pediatric brain tumor and is divided into at least four subgroups known as WNT, SHH, Group 3, and Group 4. Here, we characterized gene regulation mechanisms in the most aggressive subtype, Group 3 tumors, through genome-wide chromatin and expression profiling. Our results show that most active distal sites in these tumors are occupied by the transcription factor OTX2. Highly active OTX2-bound enhancers are often arranged as clusters of adjacent peaks and are also bound by the transcription factor NEUROD1. These sites are responsive to OTX2 and NEUROD1 knockdown and could also be generated de novo upon ectopic OTX2 expression in primary cells, showing that OTX2 cooperates with NEUROD1 and plays a major role in maintaining and possibly establishing regulatory elements as a pioneer factor. Among OTX2 target genes, we identified the kinase NEK2, whose knockdown and pharmacologic inhibition decreased cell viability. Our studies thus show that OTX2 controls the regulatory landscape of Group 3 medulloblastoma through cooperative activity at enhancer elements and contributes to the expression of critical target genes.Significance: The gene regulation mechanisms that drive medulloblastoma are not well understood. Using chromatin profiling, we find that the transcription factor OTX2 acts as a pioneer factor and, in cooperation with NEUROD1, controls the Group 3 medulloblastoma active enhancer landscape. OTX2 itself or its target genes, including the mitotic kinase NEK2, represent attractive targets for future therapies. Cancer Discov; 7(3); 288-301. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 235.