We have profiled promoter DNA methylation alterations in 272 glioblastoma tumors in the context of The Cancer Genome Atlas (TCGA). We found that a distinct subset of samples displays concerted hypermethylation at a large number of loci, indicating the existence of a glioma-CpG island methylator phenotype (G-CIMP). We validated G-CIMP in a set of non-TCGA glioblastomas and low-grade gliomas. G-CIMP tumors belong to the proneural subgroup, are more prevalent among lower-grade gliomas, display distinct copy-number alterations, and are tightly associated with IDH1 somatic mutations. Patients with G-CIMP tumors are younger at the time of diagnosis and experience significantly improved outcome. These findings identify G-CIMP as a distinct subset of human gliomas on molecular and clinical grounds.

Contemporary high dimensional biological assays, such as mRNA expression microarrays, regularly involve multiple data processing steps, such as experimental processing, computational processing, sample selection, or feature selection (i.e. gene selection), prior to deriving any biological conclusions. These steps can dramatically change the interpretation of an experiment. Evaluation of processing steps has received limited attention in the literature. It is not straightforward to evaluate different processing methods and investigators are often unsure of the best method. We present a simple statistical tool, Standardized WithIn class Sum of Squares (SWISS), that allows investigators to compare alternate data processing methods, such as different experimental methods, normalizations, or technologies, on a dataset in terms of how well they cluster a priori biological classes. SWISS uses Euclidean distance to determine which method does a better job of clustering the data elements based on a priori classifications. We apply SWISS to three different gene expression applications. The first application uses four different datasets to compare different experimental methods, normalizations, and gene sets. The second application, using data from the MicroArray Quality Control (MAQC) project, compares different microarray platforms. The third application compares different technologies: a single Agilent two-color microarray versus one lane of RNA-Seq. These applications give an indication of the variety of problems that SWISS can be helpful in solving. The SWISS analysis of one-color versus two-color microarrays provides investigators who use two-color arrays the opportunity to review their results in light of a single-channel analysis, with all of the associated benefits offered by this design. Analysis of the MACQ data shows differential intersite reproducibility by array platform. SWISS also shows that one lane of RNA-Seq clusters data by biological phenotypes as well as a single Agilent two-color microarray.

The necessary infrastructure to carry out genomics-driven oncology is now widely available and has resulted in the exponential increase in characterized cancer genomes. While a subset of genomic alterations is clinically actionable, the majority of somatic events remain classified as variants of unknown significance and will require functional characterization. A careful cataloging of the genomic alterations and their response to therapeutic intervention should allow the compilation of an "actionability atlas" and the creation of a genomic taxonomy stratified by tumor type and oncogenic pathway activation. The next phase of genomic medicine will therefore require talented bioinformaticians, genomic navigators, and multidisciplinary approaches to decode complex cancer genomes and guide potential therapy. Equally important will be the ethical and interpretable return of results to practicing oncologists. Finally, the integration of genomics into clinical trials is likely to speed the development of predictive biomarkers of response to targeted therapy as well as define pathways to acquired resistance.

High-throughput sequencing technologies, including RNA-seq, have made it possible to move beyond gene expression analysis to study transcriptional events including alternative splicing and gene fusions. Furthermore, recent studies in cancer have suggested the importance of identifying transcriptionally altered loci as biomarkers for improved prognosis and therapy. While many statistical methods have been proposed for identifying novel transcriptional events with RNA-seq, nearly all rely on contrasting known classes of samples, such as tumor and normal. Few tools exist for the unsupervised discovery of such events without class labels. In this paper, we present SigFuge for identifying genomic loci exhibiting differential transcription patterns across many RNA-seq samples. SigFuge combines clustering with hypothesis testing to identify genes exhibiting alternative splicing, or differences in isoform expression. We apply SigFuge to RNA-seq cohorts of 177 lung and 279 head and neck squamous cell carcinoma samples from the Cancer Genome Atlas, and identify several cases of differential isoform usage including CDKN2A, a tumor suppressor gene known to be inactivated in a majority of lung squamous cell tumors. By not restricting attention to known sample stratifications, SigFuge offers a novel approach to unsupervised screening of genetic loci across RNA-seq cohorts. SigFuge is available as an R package through Bioconductor.

Germline mutations in LKB1 (STK11) are associated with the Peutz-Jeghers syndrome (PJS), which includes aberrant mucocutaneous pigmentation, and somatic LKB1 mutations occur in 10% of cutaneous melanoma. By somatically inactivating Lkb1 with K-Ras activation (±p53 loss) in murine melanocytes, we observed variably pigmented and highly metastatic melanoma with 100% penetrance. LKB1 deficiency resulted in increased phosphorylation of the SRC family kinase (SFK) YES, increased expression of WNT target genes, and expansion of a CD24(+) cell population, which showed increased metastatic behavior in vitro and in vivo relative to isogenic CD24(-) cells. These results suggest that LKB1 inactivation in the context of RAS activation facilitates metastasis by inducing an SFK-dependent expansion of a prometastatic, CD24(+) tumor subpopulation.

Identifying variants using high-throughput sequencing data is currently a challenge because true biological variants can be indistinguishable from technical artifacts. One source of technical artifact results from incorrectly aligning experimentally observed sequences to their true genomic origin ('mismapping') and inferring differences in mismapped sequences to be true variants. We developed BlackOPs, an open-source tool that simulates experimental RNA-seq and DNA whole exome sequences derived from the reference genome, aligns these sequences by custom parameters, detects variants and outputs a blacklist of positions and alleles caused by mismapping. Blacklists contain thousands of artifact variants that are indistinguishable from true variants and, for a given sample, are expected to be almost completely false positives. We show that these blacklist positions are specific to the alignment algorithm and read length used, and BlackOPs allows users to generate a blacklist specific to their experimental setup. We queried the dbSNP and COSMIC variant databases and found numerous variants indistinguishable from mapping errors. We demonstrate how filtering against blacklist positions reduces the number of potential false variants using an RNA-seq glioblastoma cell line data set. In summary, accounting for mapping-caused variants tuned to experimental setups reduces false positives and, therefore, improves genome characterization by high-throughput sequencing.

Medulloblastoma (MB) is the most common malignant brain tumor in children. Patients whose tumors exhibit overexpression or amplification of the MYC oncogene (c-MYC) usually have an extremely poor prognosis, but there are no animal models of this subtype of the disease. Here, we show that cerebellar stem cells expressing Myc and mutant Trp53 (p53) generate aggressive tumors following orthotopic transplantation. These tumors consist of large, pleiomorphic cells and resemble human MYC-driven MB at a molecular level. Notably, antagonists of PI3K/mTOR signaling, but not Hedgehog signaling, inhibit growth of tumor cells. These findings suggest that cerebellar stem cells can give rise to MYC-driven MB and identify a novel model that can be used to test therapies for this devastating disease.

Mutational inactivation of the SWI/SNF chromatin regulator ATRX occurs frequently in gliomas, the most common primary brain tumors. Whether and how ATRX deficiency promotes oncogenesis by epigenomic dysregulation remains unclear, despite its recent implication in both genomic instability and telomere dysfunction. Here we report that Atrx loss recapitulates characteristic disease phenotypes and molecular features in putative glioma cells of origin, inducing cellular motility although also shifting differentiation state and potential toward an astrocytic rather than neuronal histiogenic profile. Moreover, Atrx deficiency drives widespread shifts in chromatin accessibility, histone composition, and transcription in a distribution almost entirely restricted to genomic sites normally bound by the protein. Finally, direct gene targets of Atrx that mediate specific Atrx-deficient phenotypes in vitro exhibit similarly selective misexpression in ATRX-mutant human gliomas. These findings demonstrate that ATRX deficiency and its epigenomic sequelae are sufficient to induce disease-defining oncogenic phenotypes in appropriate cellular and molecular contexts.

Patients with laryngeal squamous cell carcinoma (LSCC) often fail radiation therapy (RT), when received as monotherapy or in combination with other treatment modalities. Mechanisms for RT failure are poorly understood. We hypothesized that tumors failing RT would have increased rates of somatic mutations in genes associated with radiation resistance, particularly in genes associated with the NFE2L2 oxidative stress pathway. Using targeted exome sequencing on pretreated LSCC tumors, we retrospectively compared somatic mutation profile with clinical data and response to treatment.

Bariatric surgery is a sustainable weight loss approach, including vertical sleeve gastrectomy (VSG). Obesity exacerbates tumor growth, while diet-induced weight loss impairs progression. It remains unknown how bariatric surgery-induced weight loss impacts cancer progression or alters response to therapy. Using a pre-clinical model of obesity followed by VSG or diet-induced weight loss, breast cancer progression and immune checkpoint blockade therapy were investigated. Weight loss by VSG or weight-matched dietary intervention before tumor engraftment protected against obesity-exacerbated tumor progression. However, VSG was not as effective as diet in reducing tumor burden despite achieving similar weight and adiposity loss. Leptin did not associate with changes in tumor burden; however, circulating IL-6 was elevated in VSG mice. Uniquely, VSG tumors displayed elevated inflammation and immune checkpoint ligand PD-L1+ myeloid and non-immune cells. VSG tumors also had reduced T lymphocytes and markers of cytolysis, suggesting an ineffective anti-tumor microenvironment which prompted investigation of immune checkpoint blockade. While obese mice were resistant to immune checkpoint blockade, anti-PD-L1 potently impaired tumor progression after VSG through improved anti-tumor immunity. Thus, in formerly obese mice, surgical weight loss followed by immunotherapy reduced breast cancer burden. Finally, we compared transcriptomic changes in adipose tissue after bariatric surgery from patients and mouse models. A conserved bariatric surgery-associated weight loss signature (BSAS) was identified which significantly associated with decreased tumor volume. Findings demonstrate conserved impacts of obesity and bariatric surgery-induced weight loss pathways associated with breast cancer progression.

Growth factors in tumor environments are regulators of cell survival and metastasis. Here, we reveal the dichotomy between TGF-β superfamily growth factors BMP and TGF-β/activin and their downstream SMAD effectors. Gene expression profiling uncovers SOX2 as a key contextual signaling node regulated in an opposing manner by BMP2, -4, and -9 and TGF-β and activin A to impact anchorage-independent cell survival. We find that SOX2 is repressed by BMPs, leading to a reduction in intraperitoneal tumor burden and improved survival of tumor-bearing mice. Repression of SOX2 is driven by SMAD1-dependent histone H3K27me3 recruitment and DNA methylation at SOX2's promoter. Conversely, TGF-β, which is elevated in patient ascites, and activin A can promote SOX2 expression and anchorage-independent survival by SMAD3-dependent histone H3K4me3 recruitment. Our findings identify SOX2 as a contextual and contrastingly regulated node downstream of TGF-β members controlling anchorage-independent survival and metastasis in ovarian cancers.

Salivary duct carcinoma with rhabdoid features (SDC-RF) is a rare form of salivary gland neoplasm that was recently described. We report a case of SDC-RF of the parotid gland with loss of E-cadherin and decreased β-catenin expression in a 73-year-old male who presented with right facial/neck swelling and intermittent pain. Morphologically, the tumor presented with a discohesive infiltrate of isolated and cords of pleomorphic round cells containing moderate amount of eosinophilic to fine-vacuolated cytoplasm and hyperchromatic nuclei infiltrating through fibroadipose tissue and salivary parenchyma. Immunophenotypically, the tumor was positive for Cytokeratins Oscar and 7, GATA3, GCDFP, HER2, and an androgen receptor but negative for CK20, S100, p40, Melan A, CDX2, TTF1, ER, SATB2, DOG1, synaptophysin, and chromogranin. Due to its diffuse infiltrating pattern, involvement of the parapharyngeal space, supraclavicular fat pad, dermis, and skin without a defined surgical target, the tumor was deemed unresectable. Anti-HER2 therapy (Herceptin and Pertuzumab) was utilized. At the last follow-up, the patient is alive, with complete locoregional control and brain metastases. An electronic search was performed in the following registries for papers published up to June 2023: PubMed, Embase, and Web of Science. For the database searches, the keywords searched were "salivary gland", "salivary duct carcinoma", and "salivary duct carcinoma with rhabdoid features". Our review of the literature identified 30 cases of SDC-RF that reveal there is a predilection for males (83%), parotid gland (72%), and patients older than the 6th decade of life (83%). Immunophenotypically, all SDC-RF cases except one were positive for AR and GCDFP (97%), 81% were positive for HER2, and loss or decreased expression of E-cadherin in 93% of cases. In conclusion, we described a rare case of SDF-RF of the parotid gland with no E-cadherin expression, decreased β-catenin expression, and its immunophenotypic profile.

Almost all Glioblastoma (GBM) are either intrinsically resistant to the chemotherapeutical drug temozolomide (TMZ) or acquire therapy-induced mutations that cause chemoresistance and recurrence. The genome maintenance mechanisms responsible for GBM chemoresistance and hypermutation are unknown. We show that the E3 ubiquitin ligase RAD18 (a proximal regulator of TLS) is activated in a Mismatch repair (MMR)-dependent manner in TMZ-treated GBM cells, promoting post-replicative gap-filling and survival. An unbiased CRISPR screen provides an aerial map of RAD18-interacting DNA damage response (DDR) pathways deployed by GBM to tolerate TMZ genotoxicity. Analysis of mutation signatures from TMZ-treated GBM reveals a role for RAD18 in error-free bypass of O6mG (the most toxic TMZ-induced lesion), and error-prone bypass of other TMZ-induced lesions. Our analyses of recurrent GBM patient samples establishes a correlation between low RAD18 expression and hypermutation. Taken together we define molecular underpinnings for the hallmark tumorigenic phenotypes of TMZ-treated GBM.

While aerobic glycolysis is linked to unconstrained proliferation in cancer, less is known about its physiological role. Why this metabolic program that promotes tumor growth is preserved in the genome has thus been unresolved. We tested the hypothesis that aerobic glycolysis derives from developmental processes that regulate rapid proliferation.

Background: Little is known about the prognostic significance of somatically mutated genes in metastatic melanoma (MM). We have employed a combined clinical and bioinformatics approach on tumor samples from cutaneous melanoma (SKCM) as part of The Cancer Genome Atlas project (TCGA) to identify mutated genes with potential clinical relevance. Methods: After limiting our DNA sequencing analysis to MM samples (n = 356) and to the CANCER CENSUS gene list, we filtered out mutations with low functional significance (snpEFF). We performed Cox analysis on 53 genes that were mutated in ≥3% of samples, and had ≥50% difference in incidence of mutations in deceased subjects versus alive subjects. Results: Four genes were potentially prognostic [RAC1, FGFR1, CARD11, CIITA; false discovery rate (FDR) < 0.2]. We identified 18 additional genes (e.g., SPEN, PDGFRB, GNAS, MAP2K1, EGFR, TSC2) that were less likely to have prognostic value (FDR < 0.4). Most somatic mutations in these 22 genes were infrequent (< 10%), associated with high somatic mutation burden, and were evenly distributed across all exons, except for RAC1 and MAP2K1. Mutations in only 9 of these 22 genes were also identified by RNA sequencing in >75% of the samples that exhibited corresponding DNA mutations. The low frequency, UV signature type and RNA expression of the 22 genes in MM samples were confirmed in a separate multi-institution validation cohort (n = 413). An underpowered analysis within a subset of this validation cohort with available patient follow-up (n = 224) showed that somatic mutations in SPEN and RAC1 reached borderline prognostic significance [log-rank favorable (p = 0.09) and adverse (p = 0.07), respectively]. Somatic mutations in SPEN, and to a lesser extent RAC1, were not associated with definite gene copy number or RNA expression alterations. High (>2+) nuclear plus cytoplasmic expression intensity for SPEN was associated with longer melanoma-specific overall survival (OS) compared to lower (≤ 2+) nuclear intensity (p = 0.048). We conclude that expressed somatic mutations in infrequently mutated genes beyond the well-characterized ones (e.g., BRAF, RAS, CDKN2A, PTEN, TP53), such as RAC1 and SPEN, may have prognostic significance in MM.

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smoking and/or human papillomavirus (HPV). SCCs harbor 3q, 5p, and other recurrent chromosomal copy-number alterations (CNAs), DNA mutations, and/or aberrant methylation of genes and microRNAs, which are correlated with the expression of multi-gene programs linked to squamous cell stemness, epithelial-to-mesenchymal differentiation, growth, genomic integrity, oxidative damage, death, and inflammation. Low-CNA SCCs tended to be HPV(+) and display hypermethylation with repression of TET1 demethylase and FANCF, previously linked to predisposition to SCC, or harbor mutations affecting CASP8, RAS-MAPK pathways, chromatin modifiers, and immunoregulatory molecules. We uncovered hypomethylation of the alternative promoter that drives expression of the ΔNp63 oncogene and embedded miR944. Co-expression of immune checkpoint, T-regulatory, and Myeloid suppressor cells signatures may explain reduced efficacy of immune therapy. These findings support possibilities for molecular classification and therapeutic approaches.

Breast cancer subtype can be classified using standard clinical markers (estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2)), supplemented with additional markers. However, automated biomarker scoring and classification schemes have not been standardized. The aim of this study was to optimize tumor classification using automated methods in order to describe subtype frequency in the African American Breast Cancer Epidemiology and Risk (AMBER) consortium.

Glioma is a unique neoplastic disease that develops exclusively in the central nervous system (CNS) and rarely metastasizes to other tissues. This feature strongly implicates the tumor-host CNS microenvironment in gliomagenesis and tumor progression. We investigated the differences and similarities in glioma biology as conveyed by transcriptomic patterns across four mammalian hosts: rats, mice, dogs, and humans. Given the inherent intra-tumoral molecular heterogeneity of human glioma, we focused this study on tumors with upregulation of the platelet-derived growth factor signaling axis, a common and early alteration in human gliomagenesis. The results reveal core neoplastic alterations in mammalian glioma, as well as unique contributions of the tumor host to neoplastic processes. Notable differences were observed in gene expression patterns as well as related biological pathways and cell populations known to mediate key elements of glioma biology, including angiogenesis, immune evasion, and brain invasion. These data provide new insights regarding mammalian models of human glioma, and how these insights and models relate to our current understanding of the human disease.

In-depth molecular pathogenesis of ground-glass nodular lung adenocarcinoma has not been well understood. The objectives of this study were to identify genomic alterations in ground-glass nodular lung adenocarcinomas and to investigate whether viral transcripts were detected in these tumors. Nine patients with pure (n = 4) and part-solid (n = 5) ground-glass nodular adenocarcinomas were included. Six were females with a median age of 58 years. We performed targeted exon sequencing and RNA sequencing. EGFR (n = 10), IDH2 (n = 2), TP53 (n = 1), PTEN (n = 1), EPHB4 (n = 1), and BRAF (n = 1) were identified as driver mutations by targeted exon sequencing. Vasculogenesis-associated genes including NOTCH4 and TGFBR3 expression were significantly downregulated in adenocarcinoma tissue versus normal tissue (adjusted P values < 0.001 for both NOTCH4 and TGFBR3). In addition, five novel fusion gene loci were identified in four lung adenocarcinomas. However, no significant virus-associated transcripts were detected in tumors. In conclusions, EGFR, IDH2, TP53, PTEN, EPHB4, and BRAF were identified as putative driver mutations of ground-glass nodular adenocarcinomas. Five novel fusion genes were also identified in four tumors. Viruses do not appear to be involved in the tumorigenesis of ground-glass nodular lung adenocarcinoma.

The C-terminal binding protein (CtBP) is a NADH-dependent transcriptional repressor that links carbohydrate metabolism to epigenetic regulation by recruiting diverse histone-modifying complexes to chromatin. Here global profiling of CtBP in breast cancer cells reveals that it drives epithelial-to-mesenchymal transition, stem cell pathways and genome instability. CtBP expression induces mesenchymal and stem cell-like features, whereas CtBP depletion or caloric restriction reverses gene repression and increases DNA repair. Multiple members of the CtBP-targeted gene network are selectively downregulated in aggressive breast cancer subtypes. Differential expression of CtBP-targeted genes predicts poor clinical outcome in breast cancer patients, and elevated levels of CtBP in patient tumours predict shorter median survival. Finally, both CtBP promoter targeting and gene repression can be reversed by small molecule inhibition. These findings define broad roles for CtBP in breast cancer biology and suggest novel chromatin-based strategies for pharmacologic and metabolic intervention in cancer.