Tangent normalization for somatic copy-number inference in cancer genome analysis.

MOTIVATION: Somatic copy-number alterations (SCNAs) play an important role in cancer development. Systematic noise in sequencing and array data present a significant challenge to the inference of SCNAs for cancer genome analyses. As part of The Cancer Genome Atlas, the Broad Institute Genome Characterization Center developed the Tangent normalization method to generate copy-number profiles using data from single-nucleotide polymorphism (SNP) arrays and whole-exome sequencing (WES) technologies for over 10 000 pairs of tumors and matched normal samples. Here, we describe the Tangent method, which uses a unique linear combination of normal samples as a reference for each tumor sample, to subtract systematic errors that vary across samples. We also describe a modification of Tangent, called Pseudo-Tangent, which enables denoising through comparisons between tumor profiles when few normal samples are available. RESULTS: Tangent normalization substantially increases signal-to-noise ratios (SNRs) compared to conventional normalization methods in both SNP array and WES analyses. Tangent and Pseudo-Tangent normalizations improve the SNR by reducing noise with minimal effect on signal and exceed the contribution of other steps in the analysis such as choice of segmentation algorithm. Tangent and Pseudo-Tangent are broadly applicable and enable more accurate inference of SCNAs from DNA sequencing and array data. AVAILABILITY AND IMPLEMENTATION: Tangent is available at https://github.com/broadinstitute/tangent and as a Docker image (https://hub.docker.com/r/broadinstitute/tangent). Tangent is also the normalization method for the copy-number pipeline in Genome Analysis Toolkit 4 (GATK4). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition.

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

Patterns of somatic structural variation in human cancer genomes.

A key mutational process in cancer is structural variation, in which rearrangements delete, amplify or reorder genomic segments that range in size from kilobases to whole chromosomes1-7. Here we develop methods to group, classify and describe somatic structural variants, using data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumour types8. Sixteen signatures of structural variation emerged. Deletions have a multimodal size distribution, assort unevenly across tumour types and patients, are enriched in late-replicating regions and correlate with inversions. Tandem duplications also have a multimodal size distribution, but are enriched in early-replicating regions-as are unbalanced translocations. Replication-based mechanisms of rearrangement generate varied chromosomal structures with low-level copy-number gains and frequent inverted rearrangements. One prominent structure consists of 2-7 templates copied from distinct regions of the genome strung together within one locus. Such cycles of templated insertions correlate with tandem duplications, and-in liver cancer-frequently activate the telomerase gene TERT. A wide variety of rearrangement processes are active in cancer, which generate complex configurations of the genome upon which selection can act.

Author Correction: Targeting wild-type KRAS-amplified gastroesophageal cancer through combined MEK and SHP2 inhibition.

In the Supplementary Information originally published with this article, a lane was missing in the ß-actin blot in Supplementary Fig. 2. The lane has been added. The error has been corrected in the Supplementary Information associated with this article.

SETD2 Is Recurrently Mutated in Whole-Exome Sequenced Canine Osteosarcoma.

Osteosarcoma is a debilitating bone cancer that affects humans, especially children and adolescents. A homologous form of osteosarcoma spontaneously occurs in dogs, and its differential incidence observed across breeds allows for the investigation of tumor mutations in the context of multiple genetic backgrounds. Using whole-exome sequencing and dogs from three susceptible breeds (22 golden retrievers, 21 Rottweilers, and 23 greyhounds), we found that osteosarcoma tumors show a high frequency of somatic copy-number alterations (SCNA), affecting key oncogenes and tumor-suppressor genes. The across-breed results are similar to what has been observed for human osteosarcoma, but the disease frequency and somatic mutation counts vary in the three breeds. For all breeds, three mutational signatures (one of which has not been previously reported) and 11 significantly mutated genes were identified. TP53 was the most frequently altered gene (83% of dogs have either mutations or SCNA in TP53), recapitulating observations in human osteosarcoma. The second most frequently mutated gene, histone methyltransferase SETD2, has known roles in multiple cancers, but has not previously been strongly implicated in osteosarcoma. This study points to the likely importance of histone modifications in osteosarcoma and highlights the strong genetic similarities between human and dog osteosarcoma, suggesting that canine osteosarcoma may serve as an excellent model for developing treatment strategies in both species.Significance: Canine osteosarcoma genomics identify SETD2 as a possible oncogenic driver of osteosarcoma, and findings establish the canine model as a useful comparative model for the corresponding human disease. Cancer Res; 78(13); 3421-31. ©2018 AACR.

Targeting wild-type KRAS-amplified gastroesophageal cancer through combined MEK and SHP2 inhibition.

The role of KRAS, when activated through canonical mutations, has been well established in cancer1. Here we explore a secondary means of KRAS activation in cancer: focal high-level amplification of the KRAS gene in the absence of coding mutations. These amplifications occur most commonly in esophageal, gastric and ovarian adenocarcinomas2-4. KRAS-amplified gastric cancer models show marked overexpression of the KRAS protein and are insensitive to MAPK blockade owing to their capacity to adaptively respond by rapidly increasing KRAS-GTP levels. Here we demonstrate that inhibition of the guanine-exchange factors SOS1 and SOS2 or the protein tyrosine phosphatase SHP2 can attenuate this adaptive process and that targeting these factors, both genetically and pharmacologically, can enhance the sensitivity of KRAS-amplified models to MEK inhibition in both in vitro and in vivo settings. These data demonstrate the relevance of copy-number amplification as a mechanism of KRAS activation, and uncover the therapeutic potential for targeting of these tumors through combined SHP2 and MEK inhibition.

Genomic and Functional Approaches to Understanding Cancer Aneuploidy.

Aneuploidy, whole chromosome or chromosome arm imbalance, is a near-universal characteristic of human cancers. In 10,522 cancer genomes from The Cancer Genome Atlas, aneuploidy was correlated with TP53 mutation, somatic mutation rate, and expression of proliferation genes. Aneuploidy was anti-correlated with expression of immune signaling genes, due to decreased leukocyte infiltrates in high-aneuploidy samples. Chromosome arm-level alterations show cancer-specific patterns, including loss of chromosome arm 3p in squamous cancers. We applied genome engineering to delete 3p in lung cells, causing decreased proliferation rescued in part by chromosome 3 duplication. This study defines genomic and phenotypic correlates of cancer aneuploidy and provides an experimental approach to study chromosome arm aneuploidy.

Erratum: Genomic landscape of high-grade meningiomas.

[This corrects the article DOI: 10.1038/s41525-017-0014-7.].

Genomic landscape of high-grade meningiomas.

High-grade meningiomas frequently recur and are associated with high rates of morbidity and mortality. To determine the factors that promote the development and evolution of these tumors, we analyzed the genomes of 134 high-grade meningiomas and compared this information with data from 587 previously published meningiomas. High-grade meningiomas had a higher mutation burden than low-grade meningiomas but did not harbor any statistically significant mutated genes aside from NF2. High-grade meningiomas also possessed significantly elevated rates of chromosomal gains and losses, especially among tumors with monosomy 22. Meningiomas previously treated with adjuvant radiation had significantly more copy number alterations than radiation-induced or radiation-naïve meningiomas. Across serial recurrences, genomic disruption preceded the emergence of nearly all mutations, remained largely uniform across time, and when present in low-grade meningiomas, correlated with subsequent progression to a higher grade. In contrast to the largely stable copy number alterations, mutations were strikingly heterogeneous across tumor recurrences, likely due to extensive geographic heterogeneity in the primary tumor. While high-grade meningiomas harbored significantly fewer overtly targetable alterations than low-grade meningiomas, they contained numerous mutations that are predicted to be neoantigens, suggesting that immunologic targeting may be of therapeutic value.

Corrigendum: Recurrent hormone-binding domain truncated ESR1 amplifications in primary endometrial cancers suggest their implication in hormone independent growth.

This corrects the article DOI: 10.1038/srep25521.

Somatic copy number alterations in gastric adenocarcinomas among Asian and Western patients.

Gastric cancer, a leading worldwide cause of cancer mortality, shows high geographic and ethnic variation in incidence rates, which are highest in East Asia. The anatomic locations and clinical behavior also differ by geography, leading to the controversial idea that Eastern and Western forms of the disease are distinct. In view of these differences, we investigated whether gastric cancers from Eastern and Western patients show distinct genomic profiles. We used high-density profiling of somatic copy-number aberrations to analyze the largest collection to date of gastric adenocarcinomas and utilized genotyping data to rigorously annotate ethnic status. The size of this collection allowed us to accurately identify regions of significant copy-number alteration and separately to evaluate tumors arising in Eastern and Western patients. Among molecular subtypes classified by The Cancer Genome Atlas, the frequency of gastric cancers showing chromosomal instability was modestly higher in Western patients. After accounting for this difference, however, gastric cancers arising in Easterners and Westerners have highly similar somatic copy-number patterns. Only one genomic event, focal deletion of the phosphatase gene PTPRD, was significantly enriched in Western cases, though also detected in Eastern cases. Thus, despite the different risk factors and clinical features, gastric cancer appears to be a fundamentally similar disease in both populations and the divergent clinical outcomes cannot be ascribed to different underlying structural somatic genetic aberrations.

Clinical targeted exome-based sequencing in combination with genome-wide copy number profiling: precision medicine analysis of 203 pediatric brain tumors.

BACKGROUND: Clinical genomics platforms are needed to identify targetable alterations, but implementation of these technologies and best practices in routine clinical pediatric oncology practice are not yet well established. METHODS: Profile is an institution-wide prospective clinical research initiative that uses targeted sequencing to identify targetable alterations in tumors. OncoPanel, a multiplexed targeted exome-sequencing platform that includes 300 cancer-causing genes, was used to assess single nucleotide variants and rearrangements/indels. Alterations were annotated (Tiers 1-4) based on clinical significance, with Tier 1 alterations having well-established clinical utility. OncoCopy, a clinical genome-wide array comparative genomic hybridization (aCGH) assay, was also performed to evaluate copy number alterations and better define rearrangement breakpoints. RESULTS: Cancer genomes of 203 pediatric brain tumors were profiled across histological subtypes, including 117 samples analyzed by OncoPanel, 146 by OncoCopy, and 60 tumors subjected to both methodologies. OncoPanel revealed clinically relevant alterations in 56% of patients (44 cancer mutations and 20 rearrangements), including BRAF alterations that directed the use of targeted inhibitors. Rearrangements in MYB-QKI, MYBL1, BRAF, and FGFR1 were also detected. Furthermore, while copy number profiles differed across histologies, the combined use of OncoPanel and OncoCopy identified subgroup-specific alterations in 89% (17/19) of medulloblastomas. CONCLUSION: The combination of OncoPanel and OncoCopy multiplex genomic assays can identify critical diagnostic, prognostic, and treatment-relevant alterations and represents an effective precision medicine approach for clinical evaluation of pediatric brain tumors.

Pan-Cancer Analysis Links PARK2 to BCL-XL-Dependent Control of Apoptosis.

Mutation of the PARK2 gene can promote both Parkinson's Disease and cancer, yet the underlying mechanisms of how PARK2 controls cellular physiology is incompletely understood. Here, we show that the PARK2 tumor suppressor controls the apoptotic regulator BCL-XL and modulates programmed cell death. Analysis of approximately 10,000 tumor genomes uncovers a striking pattern of mutual exclusivity between PARK2 genetic loss and amplification of BCL2L1, implicating these genes in a common pathway. PARK2 directly binds to and ubiquitinates BCL-XL. Inactivation of PARK2 leads to aberrant accumulation of BCL-XL both in vitro and in vivo, and cancer-specific mutations in PARK2 abrogate the ability of the ubiquitin E3 ligase to target BCL-XL for degradation. Furthermore, PARK2 modulates mitochondrial depolarization and apoptosis in a BCL-XL-dependent manner. Thus, like genes at the nodal points of growth arrest pathways such as p53, the PARK2 tumor suppressor is able to exert its antiproliferative effects by regulating both cell cycle progression and programmed cell death.

Landscape of Genomic Alterations in Pituitary Adenomas.

Purpose: Pituitary adenomas are the second most common primary brain tumor, yet their genetic profiles are incompletely understood.Experimental Design: We performed whole-exome sequencing of 42 pituitary macroadenomas and matched normal DNA. These adenomas included hormonally active and inactive tumors, ones with typical or atypical histology, and ones that were primary or recurrent.Results: We identified mutations, insertions/deletions, and copy-number alterations. Nearly one-third of samples (29%) had chromosome arm-level copy-number alterations across large fractions of the genome. Despite such widespread genomic disruption, these tumors had few focal events, which is unusual among highly disrupted cancers. The other 71% of tumors formed a distinct molecular class, with somatic copy number alterations involving less than 6% of the genome. Among the highly disrupted group, 75% were functional adenomas or atypical null-cell adenomas, whereas 87% of the less-disrupted group were nonfunctional adenomas. We confirmed this association between functional subtype and disruption in a validation dataset of 87 pituitary adenomas. Analysis of previously published expression data from an additional 50 adenomas showed that arm-level alterations significantly impacted transcript levels, and that the disrupted samples were characterized by expression changes associated with poor outcome in other cancers. Arm-level losses of chromosomes 1, 2, 11, and 18 were significantly recurrent. No significantly recurrent mutations were identified, suggesting no genes are altered by exonic mutations across large fractions of pituitary macroadenomas.Conclusions: These data indicate that sporadic pituitary adenomas have distinct copy-number profiles that associate with hormonal and histologic subtypes and influence gene expression. Clin Cancer Res; 23(7); 1841-51. ©2016 AACR.

Tumor-suppressor genes that escape from X-inactivation contribute to cancer sex bias.

There is a striking and unexplained male predominance across many cancer types. A subset of X-chromosome genes can escape X-inactivation, which would protect females from complete functional loss by a single mutation. To identify putative 'escape from X-inactivation tumor-suppressor' (EXITS) genes, we examined somatic alterations from >4,100 cancers across 21 tumor types for sex bias. Six of 783 non-pseudoautosomal region (PAR) X-chromosome genes (ATRX, CNKSR2, DDX3X, KDM5C, KDM6A, and MAGEC3) harbored loss-of-function mutations more frequently in males (based on a false discovery rate < 0.1), in comparison to zero of 18,055 autosomal and PAR genes (Fisher's exact P < 0.0001). Male-biased mutations in genes that escape X-inactivation were observed in combined analysis across many cancers and in several individual tumor types, suggesting a generalized phenomenon. We conclude that biallelic expression of EXITS genes in females explains a portion of the reduced cancer incidence in females as compared to males across a variety of tumor types.

Genomic evolution and chemoresistance in germ-cell tumours.

Germ-cell tumours (GCTs) are derived from germ cells and occur most frequently in the testes. GCTs are histologically heterogeneous and distinctly curable with chemotherapy. Gains of chromosome arm 12p and aneuploidy are nearly universal in GCTs, but specific somatic genomic features driving tumour initiation, chemosensitivity and progression are incompletely characterized. Here, using clinical whole-exome and transcriptome sequencing of precursor, primary (testicular and mediastinal) and chemoresistant metastatic human GCTs, we show that the primary somatic feature of GCTs is highly recurrent chromosome arm level amplifications and reciprocal deletions (reciprocal loss of heterozygosity), variations that are significantly enriched in GCTs compared to 19 other cancer types. These tumours also acquire KRAS mutations during the development from precursor to primary disease, and primary testicular GCTs (TGCTs) are uniformly wild type for TP53. In addition, by functional measurement of apoptotic signalling (BH3 profiling) of fresh tumour and adjacent tissue, we find that primary TGCTs have high mitochondrial priming that facilitates chemotherapy-induced apoptosis. Finally, by phylogenetic analysis of serial TGCTs that emerge with chemotherapy resistance, we show how TGCTs gain additional reciprocal loss of heterozygosity and that this is associated with loss of pluripotency markers (NANOG and POU5F1) in chemoresistant teratomas or transformed carcinomas. Our results demonstrate the distinct genomic features underlying the origins of this disease and associated with the chemosensitivity phenotype, as well as the rare progression to chemoresistance. These results identify the convergence of cancer genomics, mitochondrial priming and GCT evolution, and may provide insights into chemosensitivity and resistance in other cancers.