Abstract 4340: Integrating whole genome and transcriptome analysis to inform treatment decisions in the metastatic cancer clinical setting

Genomic analysis is an established clinical tool for positioning patients to specific treatment rationales. However, given the vast range of genomic complexity that is observed between tumors the extent to which genomic characterization must be applied in order to generate a true picture of the actionable and informative genomic events in an individual's tumour is still a matter for intense debate. Historically, cost and the complexity of the bioinformatics analysis needed has limited the adoption of comprehensive genomic characterization. However with the emerging clinical significance of genomic signatures, mutation burden and immunogenomic profiling, information that can best be recovered from comprehensive genomic analyses, a better understanding of the utility of extensive genomic characterization of tumours in a clinical setting is needed. For the past 6 years, the Personalized Oncogenomics Project (POG) at the BC Cancer Agency has used paired tumour/normal Whole Genome and Transcriptome Analysis (WGTA) to comprehensively characterize tumour samples from patients with incurable metastatic cancer to inform treatment planning for each patient within a clinically relevant timeframe. To date tumours from over 700 patients have been genomically profiled using this approach. These data represent a cohort of tumours from the post-treatment metastatic setting that have been comprehensively characterized at a whole genome level augmented by comprehensive clinical annotation, including information about the primary diagnosis, progression of the disease and treatment history, as well as follow-up information including the response of the tumour to genome analysis-informed therapy. We will discuss findings from this powerful genomic cohort focusing on aspects that are particularly well represented by using a WGTA approach including; mutation signatures of biological and clinical significance, integration of expression data to complement genomic findings and highlighting the discovery of novel NRG1 gene fusions as a demonstration of how using a comprehensive multi-omic approach for genomic analysis to infer the oncogenic impact of novel genomic events can impact patient care.

Citation Format: Martin R. Jones, Yaoqing Shen, Erin Pleasance, Elisa Majounie, Laura Williamson, Eric Zhao, Eric Chuah, Karen L. Mungall, Andrew J. Mungall, Richard A. Moore, Yussanne Ma, Stephen Yip, Howard Lim, Daniel Renouf, Steven J. Jones, Janesssa Jaskin, Marco A. Marra. Integrating whole genome and transcriptome analysis to inform treatment decisions in the metastatic cancer clinical setting [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4340.

Successful targeting of the NRG1 pathway indicates novel treatment strategy for metastatic cancer

Background

NRG1 fusion-positive lung cancers have emerged as potentially actionable events in lung cancer, but clinical support is currently limited and no evidence of efficacy of this approach in cancers beyond lung has been shown.

Patients and methods

Here, we describe two patients with advanced cancers refractory to standard therapies. Patient 1 had lung adenocarcinoma and patient 2 cholangiocarcinoma. Whole-genome and transcriptome sequencing were carried out for these cases with select findings validated by fluorescence in situ hybridization.

Results

Both tumors were found to be positive for NRG1 gene fusions. In patient 1, an SDC4-NRG1 gene fusion was detected, similar gene fusions having been described in lung cancers previously. In patient 2, a novel ATP1B1-NRG1 gene fusion was detected. Cholangiocarcinoma is not a disease type in which NRG1 fusions had been described previously. Integrative genome analysis was used to assess the potential functional significance of the detected genomic events including the gene fusions, prioritizing therapeutic strategies targeting the HER-family of growth factor receptors. Both patients were treated with the pan HER-family kinase inhibitor afatinib and both displayed significant and durable response to treatment. Upon progression sites of disease were sequenced. The lack of obvious genomic events to describe the disease progression indicated that broad transcriptomic or epigenetic mechanisms could be attributed to the lack of prolonged response to afatinib.

Conclusion

These observations lend further support to the use of pan HER-tyrosine kinase inhibitors for the treatment of NRG1 fusion-positive in both cancers of lung and hepatocellular origin and indicate more broadly that cancers found to be NRG1 fusion-positive may benefit from such a clinical approach regardless of their site of origin.

Clinical trial information

Personalized Oncogenomics (POG) Program of British Columbia: Utilization of Genomic Analysis to Better Understand Tumour Heterogeneity and Evolution (NCT02155621).

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

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.

Personalized oncogenomic analysis of metastatic adenoid cystic carcinoma: using whole-genome sequencing to inform clinical decision-making

Metastatic adenoid cystic carcinomas (ACCs) can cause significant morbidity and mortality. Because of their slow growth and relative rarity, there is limited evidence for systemic therapy regimens. Recently, molecular profiling studies have begun to reveal the genetic landscape of these poorly understood cancers, and new treatment possibilities are beginning to emerge. The objective is to use whole-genome and transcriptome sequencing and analysis to better understand the genetic alterations underlying the pathology of metastatic and rare ACCs and determine potentially actionable therapeutic targets. We report five cases of metastatic ACC, not originating in the salivary glands, in patients enrolled in the Personalized Oncogenomics (POG) Program at the BC Cancer Agency. Genomic workup included whole-genome and transcriptome sequencing, detailed analysis of tumor alterations, and integration with existing knowledge of drug–target combinations to identify potential therapeutic targets. Analysis reveals low mutational burden in these five ACC cases, and mutation signatures that are commonly observed in multiple cancer types. Notably, the only recurrent structural aberration identified was the well-described MYB-NFIB fusion that was present in four of five cases, and one case exhibited a closely related MYBL1-NFIB fusion. Recurrent mutations were also identified in BAP1 and BCOR, with additional mutations in individual samples affecting NOTCH1 and the epigenetic regulators ARID2, SMARCA2, and SMARCB1. Copy changes were rare, and they included amplification of MYC and homozygous loss of CDKN2A in individual samples. Genomic analysis revealed therapeutic targets in all five cases and served to inform a therapeutic choice in three of the cases to date.

Whole genome and whole transcriptome genomic profiling of a metastatic eccrine porocarcinoma

Eccrine porocarcinomas (EPs) are rare malignant tumours of the intraepidermic sweat gland duct and most often arise from benign eccrine poromas. Some recurrent somatic genomic events have been identified in these malignancies, but very little is known about the complexity of their molecular pathophysiology. We describe the whole genome and whole transcriptome genomic profiling of a metastatic EP in a 66-year-old male patient with a previous history of localized porocarcinoma of the scalp. Whole genome and whole transcriptome genomic profiling was performed on the metastatic EP. Whole genome sequencing was performed on blood-derived DNA in order to allow a comparison between germline and somatic events. We found somatic copy losses of several tumour suppressor genes including APC, PTEN and CDKN2A, CDKN2B and CDKN1A. We identified a somatic hemizygous CDKN2A pathogenic splice site variant. De novo transcriptome assembly revealed abnormal splicing of CDKN2A p14ARF and p16INK4a. Elevated expression of oncogenes EGFR and NOTCH1 was noted and no somatic mutations were found in these genes. Wnt pathway somatic alterations were also observed. In conclusion, our results suggest that the molecular pathophysiology of malignant EP features high complexity and subtle interactions of multiple key genes. Cell cycle dysregulation and CDKN2A loss of function was found to be a new potential driver in EP tumourigenesis. Moreover, the combination of somatic copy number variants and abnormal gene expression perhaps partly related to epigenetic mechanisms, all likely contribute to the development of this rare malignancy in our patient.

The Genome of the Beluga Whale (Delphinapterus leucas)

The beluga whale is a cetacean that inhabits arctic and subarctic regions, and is the only living member of the genus Delphinapterus. The genome of the beluga whale was determined using DNA sequencing approaches that employed both microfluidic partitioning library and non-partitioned library construction. The former allowed for the construction of a highly contiguous assembly with a scaffold N50 length of over 19 Mbp and total reconstruction of 2.32 Gbp. To aid our understanding of the functional elements, transcriptome data was also derived from brain, duodenum, heart, lung, spleen, and liver tissue. Assembled sequence and all of the underlying sequence data are available at the National Center for Biotechnology Information (NCBI) under the Bioproject accession number PRJNA360851A.

The Genome of the Northern Sea Otter (Enhydra lutris kenyoni)

The northern sea otter inhabits coastal waters of the northern Pacific Ocean and is the largest member of the Mustelidae family. DNA sequencing methods that utilize microfluidic partitioned and non-partitioned library construction were used to establish the sea otter genome. The final assembly provided 2.426 Gbp of highly contiguous assembled genomic sequences with a scaffold N50 length of over 38 Mbp. We generated transcriptome data derived from a lymphoma to aid in the determination of functional elements. The assembled genome sequence and underlying sequence data are available at the National Center for Biotechnology Information (NCBI) under the BioProject accession number PRJNA388419.

The whole-genome landscape of medulloblastoma subtypes

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and ‘enhancer hijacking’ events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Genomic analysis of 491 medulloblastoma samples, including methylation profiling of 1,256 cases, effectively assigns candidate drivers to most tumours across all molecular subgroups.

Medulloblastomas are highly malignant brain tumours that develop during childhood. Paul Northcott and colleagues analysed the whole-genome sequences of 491 medulloblastomas in order to characterize the genomic landscape across tumours and identify new drivers and mutational signatures. Their integrative genomic analyses, including methylation profiling of 1,256 medulloblastomas, identifies subgroup-specific driver mutations and suggests additional tumour subtypes. The authors assign driver mutations to a high proportion of the less well characterized Group 3 and Group 4, which together contribute to more than 60% of all medulloblastomas.

Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

Genomic profiling of pelvic genital type leiomyosarcoma in a woman with a germline CHEK2:c.1100delC mutation and a concomitant diagnosis of metastatic invasive ductal breast carcinoma

We describe a woman with the known pathogenic germline variant CHEK2:c.1100delC and synchronous diagnoses of both pelvic genital type leiomyosarcoma (LMS) and metastatic invasive ductal breast carcinoma. CHEK2 (checkpoint kinase 2) is a tumor-suppressor gene encoding a serine/threonine-protein kinase (CHEK2) involved in double-strand DNA break repair and cell cycle arrest. The CHEK2:c.1100delC variant is a moderate penetrance allele resulting in an approximately twofold increase in breast cancer risk. Whole-genome and whole-transcriptome sequencing were performed on the leiomyosarcoma and matched blood-derived DNA. Despite the presence of several genomic hits within the double-strand DNA damage pathway (CHEK2 germline variant and multiple RAD51B somatic structural variants), tumor profiling did not show an obvious DNA repair deficiency signature. However, even though the LMS displayed clear malignant features, its genomic profiling revealed several characteristics classically associated with leiomyomas including a translocation, t(12;14), with one breakpoint disrupting RAD51B and the other breakpoint upstream of HMGA2 with very high expression of HMGA2 and PLAG1. This is the first report of LMS genomic profiling in a patient with the germline CHEK2:c.1100delC variant and an additional diagnosis of metastatic invasive ductal breast carcinoma. We also describe a possible mechanistic relationship between leiomyoma and LMS based on genomic and transcriptome data. Our findings suggest that RAD51B translocation and HMGA2 overexpression may play an important role in LMS oncogenesis.

Genomic and Cytogenetic Characterization of a Balanced Translocation Disrupting NUP98

A 41-year-old Asian woman with bilateral renal angiomyolipomas (AML) was incidentally identified to have a balanced translocation, 46,XX,t(11;12)(p15.4;q15). She had no other features or family history to suggest a diagnosis of tuberous sclerosis. Her healthy daughter had the same translocation and no renal AML at the age of 3 years. Whole-genome sequencing was performed on genomic maternal DNA isolated from blood. A targeted de novo assembly was then conducted with ABySS for chromosomes 11 and 12. Sanger sequencing was used to validate the translocation breakpoints. As a result, genomic characterization of chromosomes 11 and 12 revealed that the 11p breakpoint disrupted the NUP98 gene in intron 1, causing a separation of the promoter and transcription start site from the rest of the gene. The translocation breakpoint on chromosome 12q was located in a gene desert. NUP98 has not yet been associated with renal AML pathogenesis, but somatic NUP98 alterations are recurrently implicated in hematological malignancies, most often following a gene fusion event. We also found evidence for complex structural events involving chromosome 12, which appear to disrupt the TDG gene. We identified a TDGP1 partially processed pseudogene at 12p12.1, which adds complexity to the de novo assembly. In conclusion, this is the first report of a germline constitutional structural chromosome rearrangement disrupting NUP98 that occurred in a generally healthy woman with bilateral renal AML.

Integrative Analysis Identifies Four Molecular and Clinical Subsets in Uveal Melanoma

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.

Comprehensive whole genome sequence analyses yields novel genetic and structural insights for Intellectual Disability

BACKGROUND

Intellectual Disability (ID) is among the most common global disorders, yet etiology is unknown in ~30% of patients despite clinical assessment. Whole genome sequencing (WGS) is able to interrogate the entire genome, providing potential to diagnose idiopathic patients.

METHODS

We conducted WGS on eight children with idiopathic ID and brain structural defects, and their normal parents; carrying out an extensive data analyses, using standard and discovery approaches.

RESULTS

We verified de novo pathogenic single nucleotide variants (SNV) in ARID1B c.1595delG and PHF6 c.820C > T, potentially causative de novo two base indels in SQSTM1 c.115_116delinsTA and UPF1 c.1576_1577delinsA, and de novo SNVs in CACNB3 c.1289G > A, and SPRY4 c.508 T > A, of uncertain significance. We report results from a large secondary control study of 2081 exomes probing the pathogenicity of the above genes. We analyzed structural variation by four different algorithms including de novo genome assembly. We confirmed a likely contributory 165 kb de novo heterozygous 1q43 microdeletion missed by clinical microarray. The de novo assembly resulted in unmasking hidden genome instability that was missed by standard re-alignment based algorithms. We also interrogated regulatory sequence variation for known and hypothesized ID genes and present useful strategies for WGS data analyses for non-coding variation.

CONCLUSION

This study provides an extensive analysis of WGS in the context of ID, providing genetic and structural insights into ID and yielding diagnoses.

Loh M, P. Hunger S, Mullighan CG. Genomic analyses identify recurrent MEF2D fusions in acute lymphoblastic leukaemia

Chromosomal rearrangements are initiating events in acute lymphoblastic leukaemia (ALL). Here using RNA sequencing of 560 ALL cases, we identify rearrangements between MEF2D (myocyte enhancer factor 2D) and five genes (BCL9, CSF1R, DAZAP1, HNRNPUL1 and SS18) in 22 B progenitor ALL (B-ALL) cases with a distinct gene expression profile, the most common of which is MEF2D-BCL9. Examination of an extended cohort of 1,164 B-ALL cases identified 30 cases with MEF2D rearrangements, which include an additional fusion partner, FOXJ2; thus, MEF2D-rearranged cases comprise 5.3% of cases lacking recurring alterations. MEF2D-rearranged ALL is characterized by a distinct immunophenotype, DNA copy number alterations at the rearrangement sites, older diagnosis age and poor outcome. The rearrangements result in enhanced MEF2D transcriptional activity, lymphoid transformation, activation of HDAC9 expression and sensitive to histone deacetylase inhibitor treatment. Thus, MEF2D-rearranged ALL represents a distinct form of high-risk leukaemia, for which new therapeutic approaches should be considered.

Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways

Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.

MEN1 mutations in Hürthle cell (oncocytic) thyroid carcinoma

CONTEXT AND OBJECTIVE

Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation, and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management.

DESIGN AND PATIENTS

Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line, and five Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in the multiple endocrine neoplasia I (MEN1) gene.

RESULTS

Here we report the identification of MEN1 loss-of-function mutations in 4% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumors.

CONCLUSION

Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of MEN1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis, and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

Abstract 3087: Whole genome sequencing of rhabdoid tumors of the kidney

Rhabdoid tumors (RT) of the kidney (RTK) are aggressive pediatric solid tumors that predominantly affect infants. There is no effective chemotherapy and the overall 4-year survival rate is 23%. RT has a characteristic loss of SMARCB1 function, found in >90% of the patients. SMARCB1 is a conserved core subunit of the SWI/SNF chromatin-remodeling complex, which in turn is responsible for proper chromatin assembly and dynamic regulation of gene expression. Previous studies showed a remarkable paucity of mutations in coding regions of genomes and a highly penetrant cancer susceptibility in a conditional knockout mouse model. These findings support the interpretation that SMARCB1 is a tumor suppressor whose inactivation is the primary driver in RT and that RT follows a tumorigenesis model in which cancer is driven by aberrant epigenetic regulation and gene expression instead of accumulation of somatic mutations. Characterizing interplays of mutations, gene expression and epigenetic regulation will be important in understanding RT development and biology.

To achieve this goal, we will comprehensively characterize genetic and epigenetic aberrations in RT using HiSeq sequencing technology. Our research efforts include profiling whole genome, whole transcriptome, promoter methylation and histone modification in 40 primary RTK samples. Here, we report preliminary results from whole genome analyses.

Using an amplification-free library construction method, we sequenced whole genomes of 40 RTK and matched normal cases to an average haploid coverage of 39.4X. The RTK genomes were mostly diploid, but we found 35 loci that are either recurrently focally amplified or deleted using GISTIC 2.0 at FDR ≤0.05. Using the Trans-ABySS de novo short-read assembler, we assembled the RT cases’ whole genomes and identified a total of 19 genes that were recurrently rearranged in 8 out of 40 cases. Eleven of the genes were either known tumor suppressors (e.g. CABIN1, BCR) or associated with developmental or neurodegenerative diseases (e.g. UPB1, SPECC1L). The genome-wide single nucleotide variant and indel analysis showed an average somatic mutation rate at 0.37 per Mb in RTK, comparable to the previous finding of 0.19 per Mb in AT/RT. Approximately 99% of the somatic mutations occurred in non-genic regions. SMARCB1 had homozygous loss of function in 83% of cases by somatic homozygous deletion, or heterozygous deletion or truncating point mutations followed by loss of heterozygosity. The remaining cases appeared to have at least 1 copy of the gene unaffected and analyses are ongoing to investigate the inactivation mechanism in these cases.

Citation Format: Hye-Jung E. Chun, Kelsey Zhu, Jenny Q. Qian, Karen L. Mungall, Yussanne Ma, Yong-Jun Zhao, Andrew J. Mungall, Richard A. Moore, Jacquie E. Schein, Daniela S. Gerhard, Elizabeth J. Perlman, Marco A. Marra. Whole genome sequencing of rhabdoid tumors of the kidney. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3087. doi:10.1158/1538-7445.AM2014-3087

Recurrent targets of aberrant somatic hypermutation in lymphoma

Somatic hypermutation (SHM) in the variable region of immunoglobulin genes (IGV) naturally occurs in a narrow window of B cell development to provide high-affinity antibodies. However, SHM can also aberrantly target proto-oncogenes and cause genome instability. The role of aberrant SHM (aSHM) has been widely studied in various non-Hodgkin's lymphoma particularly in diffuse large B-cell lymphoma (DLBCL). Although, it has been speculated that aSHM targets a wide range of genome loci so far only twelve genes have been identified as targets of aSHM through the targeted sequencing of selected genes. A genome-wide study aiming at identifying a comprehensive set of aSHM targets recurrently occurring in DLBCL has not been previously undertaken. Here, we present a comprehensive assessment of the somatic hypermutated genes in DLBCL identified through an analysis of genomic and transcriptome data derived from 40 DLBCL patients. Our analysis verifies that there are indeed many genes that are recurrently affected by aSHM. In particular, we have identified 32 novel targets that show same or higher level of aSHM activity than genes previously reported. Amongst these novel targets, 22 genes showed a significant correlation between mRNA abundance and aSHM.

Complete genomic landscape of a recurring sporadic parathyroid carcinoma

Parathyroid carcinoma is a rare endocrine malignancy with an estimated incidence of less than 1 per million population. Excessive secretion of parathyroid hormone, extremely high serum calcium level, and the deleterious effects of hypercalcaemia are the clinical manifestations of the disease. Up to 60% of patients develop multiple disease recurrences and although long-term survival is possible with palliative surgery, permanent remission is rarely achieved. Molecular drivers of sporadic parathyroid carcinoma have remained largely unknown. Previous studies, mostly based on familial cases of the disease, suggested potential roles for the tumour suppressor MEN1 and proto-oncogene RET in benign parathyroid tumourigenesis, while the tumour suppressor HRPT2 and proto-oncogene CCND1 may also act as drivers in parathyroid cancer. Here, we report the complete genomic analysis of a sporadic and recurring parathyroid carcinoma. Mutational landscapes of the primary and recurrent tumour specimens were analysed using high-throughput sequencing technologies. Such molecular profiling allowed for identification of somatic mutations never previously identified in this malignancy. These included single nucleotide point mutations in well-characterized cancer genes such as mTOR, MLL2, CDKN2C, and PIK3CA. Comparison of acquired mutations in patient-matched primary and recurrent tumours revealed loss of PIK3CA activating mutation during the evolution of the tumour from the primary to the recurrence. Structural variations leading to gene fusions and regions of copy loss and gain were identified at a single-base resolution. Loss of the short arm of chromosome 1, along with somatic missense and truncating mutations in CDKN2C and THRAP3, respectively, provides new evidence for the potential role of these genes as tumour suppressors in parathyroid cancer. The key somatic mutations identified in this study can serve as novel diagnostic markers as well as therapeutic targets.

Barnacle: detecting and characterizing tandem duplications and fusions in transcriptome assemblies

Background

Chimeric transcripts, including partial and internal tandem duplications (PTDs, ITDs) and gene fusions, are important in the detection, prognosis, and treatment of human cancers.

Results

We describe Barnacle, a production-grade analysis tool that detects such chimeras in de novo assemblies of RNA-seq data, and supports prioritizing them for review and validation by reporting the relative coverage of co-occurring chimeric and wild-type transcripts. We demonstrate applications in large-scale disease studies, by identifying PTDs in MLL, ITDs in FLT3, and reciprocal fusions between PML and RARA, in two deeply sequenced acute myeloid leukemia (AML) RNA-seq datasets.

Conclusions

Our analyses of real and simulated data sets show that, with appropriate filter settings, Barnacle makes highly specific predictions for three types of chimeric transcripts that are important in a range of cancers: PTDs, ITDs, and fusions. High specificity makes manual review and validation efficient, which is necessary in large-scale disease studies. Characterizing an extended range of chimera types will help generate insights into progression, treatment, and outcomes for complex diseases.

The genetic landscape of high-risk neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 affected individuals (cases) using a combination of whole-exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per Mb (0.48 nonsilent) and notably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, and an additional 7.1% had focal deletions), MYCN (1.7%, causing a recurrent p.Pro44Leu alteration) and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1 and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies that rely on frequently altered oncogenic drivers.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-β signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

TBL1XR1/TP63: a novel recurrent gene fusion in B-cell non-Hodgkin lymphoma

Recently, the landscape of single base mutations in diffuse large B-cell lymphoma (DLBCL) was described. Here we report the discovery of a gene fusion between TBL1XR1 and TP63, the only recurrent somatic novel gene fusion identified in our analysis of transcriptome data from 96 DLBCL cases. Based on this cohort and a further 157 DLBCL cases analyzed by FISH, the incidence in de novo germinal center B cell-like (GCB) DLBCL is 5% (6 of 115). The fusion appears exclusive to GCB and was not seen in 138 non-GCB cases examined (P = .008, Fisher exact test) but was present at low incidence in follicular lymphoma (1 of 81). In all 7 cases identified, the 3' end of the fusion consists of exons 4 and onwards of TP63. The recurrence, subtype enrichment, and the remarkably conserved nature of the TP63 portion of the fusion suggest an important functional role in the lymphomas that harbor this event.