Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Genomic basis for RNA alterations in cancer

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Abstract SY10-02: Pan-cancer study of recurrent and heterogeneous RNA aberrations and association with whole-genome variants

Whole-exome sequencing studies have transformed our understanding of recurrent somatic mutations that contribute to cancer pathogenesis; however, these studies limit our ability to identify cancer-associated mutations to those that cause protein-coding changes. To more comprehensively catalogue cancer-associated gene alterations, we have extensively characterized tumor transcriptomes from 1,220 donors with matched whole-genome sequence data to identify recurrent RNA-level aberrations. Specifically, we created a unified RNA-Seq analysis pipeline including sequence alignment and quality control and subsequently identified gene alterations through outlier detection from estimated gene expression levels, alternative splicing, alternative transcription starts, and allele-specific expression and through identified RNA-edited sites and gene fusions. Our data represent an extensive catalog of RNA aberrations for each gene across 27 cancer types. We have also tested for genetic associations with these RNA phenotypes. Using an integrative analysis approach, we have mapped genome-wide cis and trans effects on individual RNA phenotypes, considering both common germline variants as well as somatic SNVs in gene promoters, enhancers, and intronic and other regions. Many of the regulatory associations we identify are not accessible by exome sequencing, underlining the importance of whole-genome sequence data. Utilizing this RNA-centric view, we have identified genes that are recurrently altered, yet have not been previously characterized as cancer genes or identified through DNA-level driver gene analysis. To identify further supporting evidence that these recurrent alterations are potential drivers, we identified genes with mutually exclusive RNA-level alterations. Our findings reveal new insights into selective advantages of somatic changes and molecular mechanisms of cancer. This work is by the Transcriptome Working Group of the Pan-Cancer Analysis of Whole Genomes (PCAWG) consortium and authors are listed in alphabetical order.

Citation Format: Samirkumar Amin, Philip Awadalla, Andrew Biankin, Paul Boutros, Alvis Brazma, Angela Norie Brooks, Claudia Calabrese, David Chang, Aurélien Chateigner, Ken Chen, Zechen Chong, Brian Craft, Chad Creighton, Deniz Demircioğlu, Nuno Fonseca, Milana Frenkel-Morgenstern, Gad Getz, Jonathan Göke, Mary Goldman, Liliana Greger, Syed Haider, Yao He, Katherine Hoadley, Yuan Ji, Andre Kahles, Ekta Khurana, Jan Korbel, Kjong Lehmann, Han Liang, Fenglin Liu, Maximillian Marin, Matthew Meyerson, Akinyemi Ojesina, Francis Ouellette, Chandra Pedamallu, Marc Perry, Gunnar Rätsch, Roland Schwarz, Yuichi Shiraishi, Cameron Soulette, Oliver Stegle, Patrick Tan, Alfonso Valencia, Linda Xiang, Christina Yung, Junjun Zhang, Fan Zhang, Zemin Zhang, Jingchun Zhu. Pan-cancer study of recurrent and heterogeneous RNA aberrations and association with whole-genome variants [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr SY10-02. doi:10.1158/1538-7445.AM2017-SY10-02

Subtype-specific micro-RNA expression signatures in breast cancer progression

Robust markers of invasiveness may help reduce the overtreatment of in situ carcinomas. Breast cancer is a heterogeneous disease and biological mechanisms for carcinogenesis vary between subtypes. Stratification by subtype is therefore necessary to identify relevant and robust signatures of invasive disease. We have identified microRNA (miRNA) alterations during breast cancer progression in two separate datasets and used stratification and external validation to strengthen the findings. We analyzed two separate datasets (METABRIC and AHUS) consisting of a total of 186 normal breast tissue samples, 18 ductal carcinoma in situ (DCIS) and 1,338 invasive breast carcinomas. Validation in a separate dataset and stratification by molecular subtypes based on immunohistochemistry, PAM50 and integrated cluster classifications were performed. We propose subtype-specific miRNA signatures of invasive carcinoma and a validated signature of DCIS. miRNAs included in the invasive signatures include downregulation of miR-139-5p in aggressive subtypes and upregulation of miR-29c-5p expression in the luminal subtypes. No miRNAs were differentially expressed in the transition from DCIS to invasive carcinomas on the whole, indicating the need for subtype stratification. A total of 27 miRNAs were included in our proposed DCIS signature. Significant alterations of expression included upregulation of miR-21-5p and the miR-200 family and downregulation of let-7 family members in DCIS samples. The signatures proposed here can form the basis for studies exploring DCIS samples with increased invasive potential and serum biomarkers for in situ and invasive breast cancer.

Variation in genomic landscape of clear cell renal cell carcinoma across Europe

The incidence of renal cell carcinoma (RCC) is increasing worldwide, and its prevalence is particularly high in some parts of Central Europe. Here we undertake whole-genome and transcriptome sequencing of clear cell RCC (ccRCC), the most common form of the disease, in patients from four different European countries with contrasting disease incidence to explore the underlying genomic architecture of RCC. Our findings support previous reports on frequent aberrations in the epigenetic machinery and PI3K/mTOR signalling, and uncover novel pathways and genes affected by recurrent mutations and abnormal transcriptome patterns including focal adhesion, components of extracellular matrix (ECM) and genes encoding FAT cadherins. Furthermore, a large majority of patients from Romania have an unexpected high frequency of A:T>T:A transversions, consistent with exposure to aristolochic acid (AA). These results show that the processes underlying ccRCC tumorigenesis may vary in different populations and suggest that AA may be an important ccRCC carcinogen in Romania, a finding with major public health implications.

Tandem RNA chimeras contribute to transcriptome diversity in human population and are associated with intronic genetic variants

Chimeric RNAs originating from two or more different genes are known to exist not only in cancer, but also in normal tissues, where they can play a role in human evolution. However, the exact mechanism of their formation is unknown. Here, we use RNA sequencing data from 462 healthy individuals representing 5 human populations to systematically identify and in depth characterize 81 RNA tandem chimeric transcripts, 13 of which are novel. We observe that 6 out of these 81 chimeras have been regarded as cancer-specific. Moreover, we show that a prevalence of long introns at the fusion breakpoint is associated with the chimeric transcripts formation. We also find that tandem RNA chimeras have lower abundances as compared to their partner genes. Finally, by combining our results with genomic data from the same individuals we uncover intronic genetic variants associated with the chimeric RNA formation. Taken together our findings provide an important insight into the chimeric transcripts formation and open new avenues of research into the role of intronic genetic variants in post-transcriptional processing events.

Transcriptome and genome sequencing uncovers functional variation in humans

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of mRNA and miRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project – the first uniformly processed RNA-seq data from multiple human populations with high-quality genome sequences. We discovered extremely widespread genetic variation affecting regulation of the majority of genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on cellular mechanisms of regulatory and loss-of-function variation, and allowed us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Identification of recurrent FGFR3 fusion genes in lung cancer through kinome-centred RNA sequencing

Oncogenic fusion genes that involve kinases have proven to be effective targets for therapy in a wide range of cancers. Unfortunately, the diagnostic approaches required to identify these events are struggling to keep pace with the diverse array of genetic alterations that occur in cancer. Diagnostic screening in solid tumours is particularly challenging, as many fusion genes occur with a low frequency. To overcome these limitations, we developed a capture enrichment strategy to enable high-throughput transcript sequencing of the human kinome. This approach provides a global overview of kinase fusion events, irrespective of the identity of the fusion partner. To demonstrate the utility of this system, we profiled 100 non-small cell lung cancers and identified numerous genetic alterations impacting fibroblast growth factor receptor 3 (FGFR3) in lung squamous cell carcinoma and a novel ALK fusion partner in lung adenocarcinoma.

Mucosal genome-wide methylation changes in inflammatory bowel disease

BACKGROUND

DNA methylation constitutes a key epigenetic mechanism by which cells regulate gene transcription. Among its roles are the dynamic regulation of gene expression, for example, as part of an evolving immune response, and cell differentiation in specialized tissues. Here our aim was to study the impact of differences in methylation patterns in the intestine with regard to inflammatory bowel disease (IBD) susceptibility and activity.

METHODS

Having extracted DNA from rectal biopsies, we conducted genome-wide methylation profiling using the HumanMethylation27 BeadChip microarray to identify genes showing evidence of differential methylation between cases of ulcerative colitis and Crohn's disease and healthy controls. Selected methylation signals were validated in an independent replication panel by pyrosequencing. Correlation with gene expression was sought by quantitative real-time polymerase chain reaction (RT-PCR).

RESULTS

Multiple genes showed significant evidence of differential methylation, several appearing in both ulcerative colitis and Crohn's disease comparisons including THRAP2, FANCC, GBGT1, DOK2, TNFSF4, TNFSF12, and FUT7. Many more than expected by chance overlapped with genes previously implicated as playing a role in IBD susceptibility in genome-wide association scans, including CARD9, ICAM3, and IL8RB (P < 0.001). Correlation between methylation and gene expression was identified for selected transcripts.

CONCLUSIONS

Consistent differences in DNA methylation between IBD cases and controls at regulatory sites within these genes suggest that their altered transcription contributes to IBD pathogenesis.

Sequential DNA methylation changes are associated with DNMT3B overexpression in colorectal neoplastic progression

BACKGROUND AND AIMS

Although aberrant methylation of key genes in the progression of colorectal neoplasia has been reported, no model-based analysis of the incremental changes through the intermediate adenoma stage has been described. In addition, the biological drivers for these methylation changes have yet to be defined. Linear mixed-effects modelling was used in this study to understand the onset and patterns of the methylation changes of SFRP2, IGF2 DMR0, H19, LINE-1 and a CpG island methylator phenotype (CIMP) marker panel, and they were correlated with DNA methyltransferase 3B (DNMT3B) levels of expression in a sample set representative of colorectal neoplastic progression.

METHODS

Methylation of the above CpG islands was measured using quantitative pyrosequencing assays in 261 tissue samples. This included a prospective collection of 44 colectomy specimens with concurrent normal mucosa, adenoma and invasive cancer tissues. Tissue microarrays from a subset of 64 cases were used for immunohistochemical analysis of DNMT3B expression.

RESULTS

It is shown that the onset and pattern of methylation changes during colorectal neoplastic progression are locus dependent. The CIMP marker RUNX3 was the earliest CpG island showing significant change, followed by the CIMP markers NEUROG1 and CACNA1G at the hyperplastic polyp stage. SFRP2 and IGF2 DMR0 showed significant methylation changes at the adenomatous polyp stage, followed by the CIMP markers CDKN2A and hMLH1 at the adenocarcinoma stage. DNMT3B levels of immunohistochemical expression increased significantly (p < 0.001) from normal to hyperplastic and from adenomatous polyps to carcinoma samples. DNMT3B expression correlated positively with SFRP2 methylation (r = 0.42, p < 0.001, 95% CI 0.25 to 0.56), but correlated negatively with IGF2 DMR0 methylation (r = 0.26, p = 0.01, 95% CI -0.45 to -0.05). A subset of the CIMP panel (NEUROG1, CACNA1G and CDKN2A) positively correlated with DNMT3B levels of expression (p < 0.05).

CONCLUSION

Hierarchical epigenetic alterations occur at transition points during colorectal neoplastic progression. These cumulative changes are closely correlated with a gain of DNMT3B expression, suggesting a causal relationship.