Adenocarcinoma of the lung is the leading cause of cancer death worldwide. Here we report molecular profiling of 230 resected lung adenocarcinomas using messenger RNA, microRNA and DNA sequencing integrated with copy number, methylation and proteomic analyses. High rates of somatic mutation were seen (mean 8.9 mutations per megabase). Eighteen genes were statistically significantly mutated, including RIT1 activating mutations and newly described loss-of-function MGA mutations which are mutually exclusive with focal MYC amplification. EGFR mutations were more frequent in female patients, whereas mutations in RBM10 were more common in males. Aberrations in NF1, MET, ERBB2 and RIT1 occurred in 13% of cases and were enriched in samples otherwise lacking an activated oncogene, suggesting a driver role for these events in certain tumours. DNA and mRNA sequence from the same tumour highlighted splicing alterations driven by somatic genomic changes, including exon 14 skipping in MET mRNA in 4% of cases. MAPK and PI(3)K pathway activity, when measured at the protein level, was explained by known mutations in only a fraction of cases, suggesting additional, unexplained mechanisms of pathway activation. These data establish a foundation for classification and further investigations of lung adenocarcinoma molecular pathogenesis.

Summary We report a comprehensive molecular characterization of pheochromocytomas and paragangliomas (PCC/PGLs), a rare tumor type. Multi-platform integration revealed that PCC/PGLs are driven by diverse alterations affecting multiple genes and pathways. Pathogenic germline mutations occurred in eight PCC/PGL susceptibility genes. We identified CSDE1 as a somatically-mutated driver gene, complementing four known drivers (HRAS, RET, EPAS1, NF1). We also discovered fusion genes in PCC/PGL, involving MAML3, BRAF, NGFR and NF1. Integrated analysis classified PCC/PGLs into four molecularly-defined groups: a kinase signaling subtype, a pseudohypoxia subtype, a Wnt-altered subtype, driven by MAML3 and CSDE1, and a cortical admixture subtype. Correlates of metastatic PCC/PGL included the MAML3 fusion gene. This integrated molecular characterization provides a comprehensive foundation for developing PCC/PGL precision medicine.

SUMMARY Comprehensive multiplatform analysis of 80 uveal melanomas (UM) identifies four molecularly distinct, clinically relevant subtypes: two associated with poor-prognosis monosomy 3 (M3) and two with better-prognosis disomy 3 (D3). We show that BAP1 loss follows M3 occurrence and correlates with a global DNA methylation state that is distinct from D3-UM. Poor-prognosis M3-UM divide into subsets with divergent genomic aberrations, transcriptional features, and clinical outcomes. We report change-of-function SRSF2 mutations. Within D3-UM, EIF1AX- and SRSF2/SF3B1 -mutant tumors have distinct somatic copy number alterations and DNA methylation profiles, providing insight into the biology of these low-versus intermediate-risk clinical mutation subtypes.

Previous studies have established that a subset of head and neck tumors contains human papillomavirus (HPV) sequences and that HPV-driven head and neck cancers display distinct biological and clinical features. HPV is known to drive cancer by the actions of the E6 and E7 oncoproteins, but the molecular architecture of HPV infection and its interaction with the host genome in head and neck cancers have not been comprehensively described. We profiled a cohort of 279 head and neck cancers with next generation RNA and DNA sequencing and show that 35 (12.5%) tumors displayed evidence of high-risk HPV types 16, 33, or 35. Twentyfive cases had integration of the viral genome into one or more locations in the human genome with statistical enrichment for genic regions. Integrations had a marked impact on the human genome and were associated with alterations in DNA copy number, mRNA transcript abundance and splicing, and both inter-and intrachromosomal rearrangements. Many of these events involved genes with documented roles in cancer. Cancers with integrated vs. nonintegrated HPV displayed different patterns of DNA methylation and both human and viral gene expressions. Together, these data provide insight into the mechanisms by which HPV interacts with the human genome beyond expression of viral oncoproteins and suggest that specific integration events are an integral component of viral oncogenesis.cancer | head and neck | papilloma virus | genome rearrangement | integration sites H ead and neck cancer (HNC) is a heterogeneous group of tumors characterized by a common anatomic origin, and most such tumors develop from within the mucosa and are classified as head and neck squamous cell carcinomas (HNSCCs) (1). HNSCC, the sixth most common cancer diagnosed worldwide and the eighth most common cause of cancer death (2), is frequently associated with human papillomavirus (HPV) infection (3, 4). Depending on the anatomic site of the tumor, HPV prevalence is estimated at 23-36% (5). HPV-positive HNSCCs form a distinct subset of HNCs that differs from HPV-negative HNSCCs in tumor biology and clinical characteristics, including superior clinical outcomes (6-9).The molecular pathogenesis of HPV-driven HNSCC also seems distinct from HPV-negative tumors, with previous studies showing a divergent spectrum of alterations in gene expression, mutations, amplifications, and deletions as well as distinct epigenome alterations (10-15). HPV is known to drive tumorigenesis through the actions of its major oncoproteins E6 and E7, which target numerous cellular pathways, including inactivation of p53 and the retinoblastoma (Rb) protein (16-18). Together with E5, they also play an important role in immune evasion, being involved in both innate and adaptive immunity (19,20).Initially after infection, HPV is identified in circular extrachromosomal particles or episomes. A critical step in progression to cancer is the integration of viral DNA into the host cell Significance A significant proportion of head and neck cancer is driven by human papil...

In the originally published version of this paper, there were two instances in which a paper by Raab et al. was cited, but the full reference was accidentally omitted from the reference list. First, in the subsection titled ''Arid1a Haploinsufficiency Alters Global Chromatin Occupancy and Metastasis Genes,'' paragraph 4, the authors wrote, ''To determine if these genes have direct interactions with ARID1A and SWI/SNF components in human HCC, we examined ChIP-seq experiments from HepG2 cells performed by Raab et al.'' Second, in the legend for Figure 7G, the authors wrote, ''ChIP-seq data showing binding of EMILIN1, MAT1A, LCN2, and IL1R1 loci by ARID1A and SNF5 over input in human HepG2 hepatoma cells (data from Raab et al.).''

Thymic epithelial tumors (TETs) are one of the rarest adult malignancies. Among TETs, thymoma is the most predominant, characterized by a unique association with autoimmune diseases, followed by thymic carcinoma, which is less common but more clinically aggressive. Using multi-platform omics analyses on 117 TETs, we define four subtypes of these tumors defined by genomic hallmarks and an association with survival and World Health Organization histological subtype. We further demonstrate a marked prevalence of a thymoma-specific mutated oncogene, GTF2I, and explore its biological effects on multi-platform analysis. We further observe enrichment of mutations in HRAS, NRAS, and TP53. Last, we identify a molecular link between thymoma and the autoimmune disease myasthenia gravis, characterized by tumoral overexpression of muscle autoantigens, and increased aneuploidy.

Summary Cluster analysis has proved to be an invaluable tool for the exploratory and unsupervised analysis of high dimensional datasets. Among methods for clustering, hierarchical approaches have enjoyed substantial popularity in genomics and other fields for their ability to simultaneously uncover multiple layers of clustering structure. A critical and challenging question in cluster analysis is whether the identified clusters represent important underlying structure or are artifacts of natural sampling variation. Few approaches have been proposed for addressing this problem in the context of hierarchical clustering, for which the problem is further complicated by the natural tree structure of the partition, and the multiplicity of tests required to parse the layers of nested clusters. In this paper, we propose a Monte Carlo based approach for testing statistical significance in hierarchical clustering which addresses these issues. The approach is implemented as a sequential testing procedure guaranteeing control of the family-wise error rate. Theoretical justification is provided for our approach, and its power to detect true clustering structure is illustrated through several simulation studies and applications to two cancer gene expression datasets.

Background In high-throughput studies, hundreds to millions of hypotheses are typically tested. Statistical methods that control the false discovery rate (FDR) have emerged as popular and powerful tools for error rate control. While classic FDR methods use only p values as input, more modern FDR methods have been shown to increase power by incorporating complementary information as informative covariates to prioritize, weight, and group hypotheses. However, there is currently no consensus on how the modern methods compare to one another. We investigate the accuracy, applicability, and ease of use of two classic and six modern FDR-controlling methods by performing a systematic benchmark comparison using simulation studies as well as six case studies in computational biology. Results Methods that incorporate informative covariates are modestly more powerful than classic approaches, and do not underperform classic approaches, even when the covariate is completely uninformative. The majority of methods are successful at controlling the FDR, with the exception of two modern methods under certain settings. Furthermore, we find that the improvement of the modern FDR methods over the classic methods increases with the informativeness of the covariate, total number of hypothesis tests, and proportion of truly non-null hypotheses. Conclusions Modern FDR methods that use an informative covariate provide advantages over classic FDR-controlling procedures, with the relative gain dependent on the application and informativeness of available covariates. We present our findings as a practical guide and provide recommendations to aid researchers in their choice of methods to correct for false discoveries. Electronic supplementary material The online version of this article (10.1186/s13059-019-1716-1) contains supplementary material, which is available to authorized users.