Genetic determinants of the molecular portraits of epithelial cancers

The ability to characterize and predict tumor phenotypes is crucial to precision medicine. In this study, we present an integrative computational approach using a genome-wide association analysis and an Elastic Net prediction method to analyze the relationship between DNA copy number alterations and an archive of gene expression signatures. Across breast cancers, we are able to quantitatively predict many gene signatures levels within individual tumors with high accuracy based upon DNA copy number features alone, including proliferation status and Estrogen-signaling pathway activity. We can also predict many other key phenotypes, including intrinsic molecular subtypes, estrogen receptor status, and TP53 mutation. This approach is also applied to TCGA Pan-Cancer, which identify repeatedly predictable signatures across tumor types including immune features in lung squamous and basal-like breast cancers. These Elastic Net DNA predictors could also be called from DNA-based gene panels, thus facilitating their use as biomarkers to guide therapeutic decision making

Amplification of SOX4 promotes PI3K/Akt signaling in human breast cancer

Purpose: The PI3K/Akt signaling axis contributes to the dysregulation of many dominant features in breast cancer including cell proliferation, survival, metabolism, motility, and genomic instability. While multiple studies have demonstrated that basal-like or triple-negative breast tumors have uniformly high PI3K/Akt activity, genomic alterations that mediate dysregulation of this pathway in this subset of highly aggressive breast tumors remain to be determined. Methods: In this study, we present an integrated genomic analysis based on the use of a PI3K gene expression signature as a framework to analyze orthogonal genomic data from human breast tumors, including RNA expression, DNA copy number alterations, and protein expression. In combination with data from a genome-wide RNA-mediated interference screen in human breast cancer cell lines, we identified essential genetic drivers of PI3K/Akt signaling. Results: Our in silico analyses identified SOX4 amplification as a novel modulator of PI3K/Akt signaling in breast cancers and in vitro studies confirmed its role in regulating Akt phosphorylation. Conclusions: Taken together, these data establish a role for SOX4-mediated PI3K/Akt signaling in breast cancer and suggest that SOX4 may represent a novel therapeutic target and/or biomarker for current PI3K family therapies

Virus expression detection reveals RNA-sequencing contamination in TCGA

Background: Contamination of reagents and cross contamination across samples is a long-recognized issue in molecular biology laboratories. While often innocuous, contamination can lead to inaccurate results. Cantalupo et al., for example, found HeLa-derived human papillomavirus 18 (H-HPV18) in several of The Cancer Genome Atlas (TCGA) RNA-sequencing samples. This work motivated us to assess a greater number of samples and determine the origin of possible contaminations using viral sequences. To detect viruses with high specificity, we developed the publicly available workflow, VirDetect, that detects virus and laboratory vector sequences in RNA-seq samples. We applied VirDetect to 9143 RNA-seq samples sequenced at one TCGA sequencing center (28/33 cancer types) over 5 years. Results: We confirmed that H-HPV18 was present in many samples and determined that viral transcripts from H-HPV18 significantly co-occurred with those from xenotropic mouse leukemia virus-related virus (XMRV). Using laboratory metadata and viral transcription, we determined that the likely contaminant was a pool of cell lines known as the "common reference", which was sequenced alongside TCGA RNA-seq samples as a control to monitor quality across technology transitions (i.e. microarray to GAII to HiSeq), and to link RNA-seq to previous generation microarrays that standardly used the "common reference". One of the cell lines in the pool was a laboratory isolate of MCF-7, which we discovered was infected with XMRV; another constituent of the pool was likely HeLa cells. Conclusions: Altogether, this indicates a multi-step contamination process. First, MCF-7 was infected with an XMRV. Second, this infected cell line was added to a pool of cell lines, which contained HeLa. Finally, RNA from this pool of cell lines contaminated several TCGA tumor samples most-likely during library construction. Thus, these human tumors with H-HPV or XMRV reads were likely not infected with H-HPV 18 or XMRV

Prognostic value of B cells in cutaneous melanoma

Background: Measures of the adaptive immune response have prognostic and predictive associations in melanoma and other cancer types. Specifically, intratumoral T cell density and function have considerable prognostic and predictive value in skin cutaneous melanoma (SKCM). Less is known about the significance of tumor-infiltrating B cells in SKCM. Our goal was to understand the prognostic and predictive value of B cell phenotypic subsets in SKCM using RNA sequencing. Methods: We used our previously published algorithm, V'DJer, to assemble B cell receptor (BCR) repertoires and estimate diversity from short-read RNA sequencing (RNA-seq). We applied machine learning-based cellular phenotype classifiers to measure relative similarity of bulk tumor sample gene expression profiles and different B cell phenotypes. We assessed these aspects of B cell biology in 473 SKCM from the Cancer Genome Atlas Project (TCGA) as well as in RNA-seq data corresponding to tumor samples procured from patients who received CTLA-4 and PD-1 inhibitors for metastatic SKCM. Results: We found that the BCR repertoire was associated with different clinical factors, such as tumor tissue site and sex. However, increased clonality of the BCR repertoire was favorably prognostic in SKCM and was prognostic even after first conditioning on various clinical factors. Mutation burden was not correlated with any BCR measurement, and no specific mutation had an altered BCR repertoire. Lack of an assembled BCR in pre-treatment tumor tissues was associated with a lack of anti-tumor response to a CTLA-4 inhibitor in metastatic SKCM. Conclusions: These findings suggest an important prognostic and predictive role for B cell characteristics in SKCM. This has implications for melanoma immunobiology and potential development of immunogenomics features to predict survival and response to immunotherapy

FOXM1 Deubiquitination by USP21 Regulates Cell Cycle Progression and Paclitaxel Sensitivity in Basal-like Breast Cancer

The cell cycle transcription factor FOXM1 is activated in basal-like breast cancer (BLBC) and associated with therapeutic resistance and poor patient outcomes. Arceci et al. show USP21 antagonizes FOXM1 degradation, thereby promoting proliferation and paclitaxel resistance. USP21 is catalytically active and recurrently overexpressed in BLBC, representing a potential therapeutic target. © 2019 The Author(s)The transcription factor FOXM1 contributes to cell cycle progression and is significantly upregulated in basal-like breast cancer (BLBC). Despite its importance in normal and cancer cell cycles, we lack a complete understanding of mechanisms that regulate FOXM1. We identified USP21 in an RNAi-based screen for deubiquitinases that control FOXM1 abundance. USP21 increases the stability of FOXM1, and USP21 binds and deubiquitinates FOXM1 in vivo and in vitro, indicating a direct enzyme-substrate relationship. Depleting USP21 downregulates the FOXM1 transcriptional network and causes a significant delay in cell cycle progression. Significantly, USP21 depletion sensitized BLBC cell lines and mouse xenograft tumors to paclitaxel, an anti-mitotic, frontline therapy in BLBC treatment. USP21 is the most frequently amplified deubiquitinase in BLBC patient tumors, and its amplification co-occurs with the upregulation of FOXM1 protein. Altogether, these data suggest a role for USP21 in the proliferation and potentially treatment of FOXM1-high, USP21-high BLBC

Subtyping sub-Saharan esophageal squamous cell carcinoma by comprehensive molecular analysis

Esophageal squamous cell carcinoma (ESCC) is endemic in regions of sub-Saharan Africa (SSA), where it is the third most common cancer. Here, we describe whole-exome tumor/normal sequencing and RNA transcriptomic analysis of 59 patients with ESCC in Malawi. We observed similar genetic aberrations as reported in Asian and North American cohorts, including mutations of TP53, CDKN2A, NFE2L2, CHEK2, NOTCH1, FAT1, and FBXW7. Analyses for nonhuman sequences did not reveal evidence for infection with HPV or other occult pathogens. Mutational signature analysis revealed common signatures associated with aging, cytidine deaminase activity (APOBEC), and a third signature of unknown origin, but signatures of inhaled tobacco use, aflatoxin and mismatch repair were notably absent. Based on RNA expression analysis, ESCC could be divided into 3 distinct subtypes, which were distinguished by their expression of cell cycle and neural transcripts. This study demonstrates discrete subtypes of ESCC in SSA, and suggests that the endemic nature of this disease reflects exposure to a carcinogen other than tobacco and oncogenic viruses

PAM50 molecular intrinsic subtypes in the nurses' health Study cohorts

Background: Modified median and subgroup-specific gene subtypes by PAM50 and IHC surrogates improved to fair centering are two essential preprocessing methods to assign when Luminal subtypes were grouped together. Using the breast cancer molecular subtypes by PAM50. We evaluated the modified median method, our study consisted of 46% PAM50 subtypes derived from both methods in a subset of Luminal A, 18% Luminal B, 14% HER2-enriched, 15% Nurses' Health Study (NHS) and NHSII participants; correlat-Basal-like, and 8% Normal-like subtypes; 53% of tumor-ed tumor subtypes by PAM50 with IHC surrogates; and adjacent tissues were Normal-like. Women with the Basal-characterized the PAM50 subtype distribution, proliferation like subtype had a higher rate of relapse within 5 years. scores, and risk of relapse with proliferation and tumor size HER2-enriched subtypes had poorer outcomes prior to weighted (ROR-PT) scores in the NHS/NHSII. 1999. Methods: PAM50 subtypes, proliferation scores, and Conclusions: Either preprocessing method may be uti-ROR-PT scores were calculated for 882 invasive breast tumors lized to derive PAM50 subtypes for future studies. The and 695 histologically normal tumor-adjacent tissues. Cox majority of NHS/NHSII tumor and tumor-adjacent tissues proportional hazards models evaluated the relationship were classified as Luminal A and Normal-like, respectively. between PAM50 subtypes or ROR-PT scores/groups with Impact: Preprocessing methods are important for the recurrence-free survival (RFS) or distant RFS. accurate assignment of PAM50 subtypes. These data provide Results: PAM50 subtypes were highly comparable evidence that either preprocessing method can be used in between the two methods. The agreement between tumor epidemiologic studies

DNA Damage Repair Classifier Defines Distinct Groups in Hepatocellular Carcinoma

DNA repair pathways have been associated with variability in hepatocellular carcinoma (HCC) clinical outcomes, but the mechanism through which DNA repair varies as a function of liver regeneration and other HCC characteristics is poorly understood. We curated a panel of 199 genes representing 15 DNA repair pathways to identify DNA repair expression classes and evaluate their associations with liver features and clinicopathologic variables in The Cancer Genome Atlas (TCGA) HCC study. We identified two groups in HCC, defined by low or high expression across all DNA repair pathways. The low-repair group had lower grade and retained the expression of classical liver markers, whereas the high-repair group had more clinically aggressive features, increased p53 mutant-like gene expression, and high liver regenerative gene expression. These pronounced features overshadowed the variation in the low-repair subset, but when considered separately, the low-repair samples included three subgroups: L1, L2, and L3. L3 had high DNA repair expression with worse progression-free (HR 1.24, 95% CI 0.81–1.91) and overall (HR 1.63, 95% CI 0.98–2.71) survival. High-repair outcomes were also significantly worse compared with the L1 and L2 groups. HCCs vary in DNA repair expression, and a subset of tumors with high regeneration profoundly disrupts liver biology and poor prognosis

SCISSOR: a framework for identifying structural changes in RNA transcripts

High-throughput sequencing protocols such as RNA-seq have made it possible to interrogate the sequence, structure and abundance of RNA transcripts at higher resolution than previous microarray and other molecular techniques. While many computational tools have been proposed for identifying mRNA variation through differential splicing/alternative exon usage, challenges in its analysis remain. Here, we propose a framework for unbiased and robust discovery of aberrant RNA transcript structures using short read sequencing data based on shape changes in an RNA-seq coverage profile. Shape changes in selecting sample outliers in RNA-seq, SCISSOR, is a series of procedures for transforming and normalizing base-level RNA sequencing coverage data in a transcript independent manner, followed by a statistical framework for its analysis (https://github.com/hyochoi/SCISSOR). The resulting high dimensional object is amenable to unsupervised screening of structural alterations across RNA-seq cohorts with nearly no assumption on the mutational mechanisms underlying abnormalities. This enables SCISSOR to independently recapture known variants such as splice site mutations in tumor suppressor genes as well as novel variants that are previously unrecognized or difficult to identify by any existing methods including recurrent alternate transcription start sites and recurrent complex deletions in 3′ UTRs

The prognostic significance of low-frequency somatic mutations in metastatic cutaneous melanoma

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) 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