ExomeAI: detection of recurrent allelic imbalance in tumors using whole-exome sequencing data

Germline Missense Variants in CDC20 Result in Aberrant Mitotic Progression and Familial Cancer

CDC20 is a coactivator of the anaphase promoting complex/cyclosome (APC/C) and is essential for mitotic progression. APC/CCDC20 is inhibited by the spindle assembly checkpoint (SAC), which prevents premature separation of sister chromatids and aneuploidy in daughter cells. Although overexpression of CDC20 is common in many cancers, oncogenic mutations have never been identified in humans. Using whole-exome sequencing, we identified heterozygous missense CDC20 variants (L151R and N331K) that segregate with ovarian germ cell tumors in two families. Functional characterization showed these mutants retain APC/C activation activity but have impaired binding to BUBR1, a component of the SAC. Expression of L151R and N331K variants promoted mitotic slippage in HeLa cells and primary skin fibroblasts derived from carriers. Generation of mice carrying the N331K variant using CRISPR-Cas9 showed that, although homozygous N331K mice were nonviable, heterozygotes displayed accelerated oncogenicity of Myc-driven cancers. These findings highlight an unappreciated role for CDC20 variants as tumor-promoting genes.

SIGNIFICANCE

Two germline CDC20 missense variants that segregate with cancer in two families compromise the spindle assembly checkpoint and lead to aberrant mitotic progression, which could predispose cells to transformation. See related commentary by Villarroya-Beltri and Malumbres, p. 3432.

Common clonal origin of chronic myelomonocytic leukemia and B-cell acute lymphoblastic leukemia in a patient with a germline CHEK2 variant

Hematological malignancies are broadly divided into myeloid and lymphoid neoplasms, reflecting the two major cellular lineages of the hematopoietic system. It is generally rare for hematological malignancies to spontaneously progress with a switch from myeloid to lymphoid lineage. We describe the exceptional case of a patient who sequentially developed myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML), and B-cell acute lymphoblastic leukemia (B-ALL), as well as our investigation into the underlying pathogenesis. Using whole-exome sequencing (WES) performed on sorted CMML and B-ALL cell fractions, we identified both common and unique potential driver mutations, suggesting a branching clonal evolution giving rise to both diseases. Interestingly, we also identified a germline variant in the cancer susceptibility gene CHEK2 We validated that this variant (c.475T > C; p.Y159H), located in the forkhead-associated (FHA) domain, impairs its capacity to bind BRCA1 in cellulo. This unique case provides novel insight into the genetics of complex hematological diseases and highlights the possibility that such patients may carry inherited predispositions.

Outcome and molecular analysis of young children with choroid plexus carcinoma treated with non-myeloablative therapy: results from the SJYC07 trial

Background

Choroid plexus carcinoma (CPC) is a rare and aggressive tumor of infancy without a clear treatment strategy. This study describes the outcomes of children with CPC treated on the multi-institutional phase 2 SJYC07 trial and reports on the significance of clinical and molecular characteristics.

Methods

Eligible children <3 years-old with CPC were postoperatively stratified to intermediate-risk (IR) stratum if disease was localized or high-risk (HR) stratum, if metastatic. All received high-dose methotrexate-containing induction chemotherapy. IR-stratum patients received focal irradiation as consolidation whereas HR-stratum patients received additional chemotherapy. Consolidation was followed by oral antiangiogenic maintenance regimen. Survival rates and potential prognostic factors were analyzed.

Results

Thirteen patients (median age: 1.41 years, range: 0.21-2.93) were enrolled; 5 IR, 8 HR. Gross-total resection or near-total resection was achieved in ten patients and subtotal resection in 3. Seven patients had TP53-mutant tumors, including 4 who were germline carriers. Five patients experienced progression and died of disease; 8 (including 5 HR) are alive without progression. The 5-year progression-free survival (PFS) and overall survival rates were 61.5 ± 13.5% and 68.4 ± 13.1%. Patients with TP53-wild-type tumors had a 5-year PFS of 100% as compared to 28.6 ± 17.1% for TP53-mutant tumors (P = .012). Extent of resection, metastatic status, and use of radiation therapy were not significantly associated with survival.

Conclusions

Non-myeloablative high-dose methotrexate-containing therapy with maximal surgical resection resulted in long-term PFS in more than half of patients with CPC. TP53-mutational status was the only significant prognostic variable and should form the basis of risk-stratification in future trials.

DGCR8 microprocessor defect characterizes familial multinodular goiter with schwannomatosis

BACKGROUNDDICER1 is the only miRNA biogenesis component associated with an inherited tumor syndrome, featuring multinodular goiter (MNG) and rare pediatric-onset lesions. Other susceptibility genes for familial forms of MNG likely exist.METHODSWhole-exome sequencing of a kindred with early-onset MNG and schwannomatosis was followed by investigation of germline pathogenic variants that fully segregated with the disease. Genome-wide analyses were performed on 13 tissue samples from familial and nonfamilial DGCR8-E518K-positive tumors, including MNG, schwannomas, papillary thyroid cancers (PTCs), and Wilms tumors. miRNA profiles of 4 tissue types were compared, and sequencing of miRNA, pre-miRNA, and mRNA was performed in a subset of 9 schwannomas, 4 of which harbor DGCR8-E518K.RESULTSWe identified c.1552G>A;p.E518K in DGCR8, a microprocessor component located in 22q, in the kindred. The variant identified is a somatic hotspot in Wilms tumors and has been identified in 2 PTCs. Copy number loss of chromosome 22q, leading to loss of heterozygosity at the DGCR8 locus, was found in all 13 samples harboring c.1552G>A;p.E518K. miRNA profiling of PTCs, MNG, schwannomas, and Wilms tumors revealed a common profile among E518K hemizygous tumors. In vitro cleavage demonstrated improper processing of pre-miRNA by DGCR8-E518K. MicroRNA and RNA profiling show that this variant disrupts precursor microRNA production, impacting populations of canonical microRNAs and mirtrons.CONCLUSIONWe identified DGCR8 as the cause of an unreported autosomal dominant mendelian tumor susceptibility syndrome: familial multinodular goiter with schwannomatosis.FUNDINGCanadian Institutes of Health Research, Compute Canada, Alex's Lemonade Stand Foundation, the Mia Neri Foundation for Childhood Cancer, Cassa di Sovvenzioni e Risparmio fra il Personale della Banca d'Italia, and the KinderKrebsInitiative Buchholz/Holm-Seppensen.

Poorly differentiated thyroid carcinoma of childhood and adolescence: a distinct entity characterized by DICER1 mutations

Poorly differentiated thyroid carcinomas (PDTC) in young individuals are rare and their clinical and histopathologic features, genetic mechanisms, and outcomes remain largely unknown. Here, we report a detailed characterization of a series of six PDTC in patients ≤21 years old defined by Turin diagnostic criteria studied for mutations and gene fusions characteristic of thyroid cancer using targeted next-generation sequencing (NGS) and whole-exome sequencing (WES). All tumors had solid, insular, or trabecular growth pattern and high mitotic rate, and five out of six tumors showed tumor necrosis. Targeted NGS assay identified somatic mutations in the DICER1 gene in five of six (83%) tumors, all of which were "hotspot" mutations encoding the metal-ion binding sites of the RNase IIIb domain of DICER1. WES was performed in five cases which confirmed all hotspot mutations and detected two tumors with additional inactivating DICER1 alterations. Of these two, one was a germline pathogenic DICER1 variant and the other had loss of heterozygosity for DICER1. No other mutations or gene fusions characteristic of adult well-differentiated thyroid cancer and PDTC (BRAF, RAS, TERT, RET/PTC, and other) were detected. On follow-up, available for five patients, three patients died of disease 8-24 months after diagnosis, whereas two were alive with no disease. The results of our study demonstrate that childhood- and adolescent-onset PDTC are genetically distinct from adult-onset PDTC in that they are strongly associated with DICER1 mutations and may herald DICER1 syndrome in a minority. As such, all young persons with PDTC may benefit from genetic counseling. Furthermore, their clinically aggressive behavior contrasts sharply with the indolent nature of the great majority of thyroid tumors with DICER1 mutations reported to date.

Atypical tuberous sclerosis complex presenting as familial renal cell carcinoma with leiomyomatous stroma

We report an atypical tuberous sclerosis complex (TSC) phenotype presenting as familial multiple renal cell carcinomas (RCCs) with (angio)leiomyomatous stroma (RCCLS) (5/7 familial RCCs) on a background of multiple angiomyolipomas, hypopigmented skin macules, and absence of neurological anomalies. In the index case and three relatives, germline genetic testing identified a heterozygous TSC2 missense pathogenic variant [c.2714 G > A, (p.Arg905Gln)], a rare TSC‐associated alteration which has previously been associated with a milder TSC phenotype. Whole‐exome sequencing of five RCCs from the index case and one RCC from his mother demonstrated either unique tumour‐specific deleterious second hits in TSC2 or significant allelic imbalance at the TSC2 gene locus (5/6 RCCs). This study confirms the key tumourigenic role of tumour‐specific TSC2 second hits in TSC‐associated RCCs and supports the notion that RCCLS may be strongly related to abnormalities of the mTOR pathway.

A landscape of germ line mutations in a cohort of inherited bone marrow failure patients

Bone marrow (BM) failure (BMF) in children and young adults is often suspected to be inherited, but in many cases diagnosis remains uncertain. We studied a cohort of 179 patients (from 173 families) with BMF of suspected inherited origin but unresolved diagnosis after medical evaluation and Fanconi anemia exclusion. All patients had cytopenias, and 12.0% presented ≥5% BM blast cells. Median age at genetic evaluation was 11 years; 20.7% of patients were aged ≤2 years and 36.9% were ≥18 years. We analyzed genomic DNA from skin fibroblasts using whole-exome sequencing, and were able to assign a causal or likely causal germ line mutation in 86 patients (48.0%), involving a total of 28 genes. These included genes in familial hematopoietic disorders (GATA2, RUNX1), telomeropathies (TERC, TERT, RTEL1), ribosome disorders (SBDS, DNAJC21, RPL5), and DNA repair deficiency (LIG4). Many patients had an atypical presentation, and the mutated gene was often not clinically suspected. We also found mutations in genes seldom reported in inherited BMF (IBMF), such as SAMD9 and SAMD9L (N = 16 of the 86 patients, 18.6%), MECOM/EVI1 (N = 6, 7.0%), and ERCC6L2 (N = 7, 8.1%), each of which was associated with a distinct natural history; SAMD9 and SAMD9L patients often experienced transient aplasia and monosomy 7, whereas MECOM patients presented early-onset severe aplastic anemia, and ERCC6L2 patients, mild pancytopenia with myelodysplasia. This study broadens the molecular and clinical portrait of IBMF syndromes and sheds light on newly recognized disease entities. Using a high-throughput sequencing screen to implement precision medicine at diagnosis can improve patient management and family counseling.

Characterizing temporal genomic heterogeneity in pediatric high-grade gliomas

Pediatric high-grade gliomas (pHGGs) are aggressive neoplasms representing approximately 20% of brain tumors in children. Current therapies offer limited disease control, and patients have a poor prognosis. Empiric use of targeted therapy, especially at progression, is increasingly practiced despite a paucity of data regarding temporal and therapy-driven genomic evolution in pHGGs. To study the genetic landscape of pHGGs at recurrence, we performed whole exome and methylation analyses on matched primary and recurrent pHGGs from 16 patients. Tumor mutational profiles identified three distinct subgroups. Group 1 (n = 7) harbored known hotspot mutations in Histone 3 (H3) (K27M or G34V) or IDH1 (H3/IDH1 mutants) and co-occurring TP53 or ACVR1 mutations in tumor pairs across the disease course. Group 2 (n = 7), H3/IDH1 wildtype tumor pairs, harbored novel mutations in chromatin modifiers (ZMYND11, EP300 n = 2), all associated with TP53 alterations, or had BRAF V600E mutations (n = 2) conserved across tumor pairs. Group 3 included 2 tumors with NF1 germline mutations. Pairs from primary and relapsed pHGG samples clustered within the same DNA methylation subgroup. ATRX mutations were clonal and retained in H3G34V and H3/IDH1 wildtype tumors, while different genetic alterations in this gene were observed at diagnosis and recurrence in IDH1 mutant tumors. Mutations in putative drug targets (EGFR, ERBB2, PDGFRA, PI3K) were not always shared between primary and recurrence samples, indicating evolution during progression. Our findings indicate that specific key driver mutations in pHGGs are conserved at recurrence and are prime targets for therapeutic development and clinical trials (e.g. H3 post-translational modifications, IDH1, BRAF V600E). Other actionable mutations are acquired or lost, indicating that re-biopsy at recurrence will provide better guidance for effective targeted therapy of pHGGs.

Functionally Null RAD51D Missense Mutation Associates Strongly with Ovarian Carcinoma

RAD51D is a key player in DNA repair by homologous recombination (HR), and RAD51D truncating variant carriers have an increased risk for ovarian cancer. However, the contribution of nontruncating RAD51D variants to cancer predisposition remains uncertain. Using deep sequencing and case-control genotyping studies, we show that in French Canadians, the missense RAD51D variant c.620C>T;p.S207L is highly prevalent and is associated with a significantly increased risk for ovarian high-grade serous carcinoma (HGSC; 3.8% cases vs. 0.2% controls). The frequency of the p.S207L variant did not significantly differ from that of controls in breast, endometrial, pancreas, or colorectal adenocarcinomas. Functionally, we show that this mutation impairs HR by disrupting the RAD51D-XRCC2 interaction and confers PARP inhibitor sensitivity. These results highlight the importance of a functional RAD51D-XRCC2 interaction to promote HR and prevent the development of HGSC. This study identifies c.620C>T;p.S207L as the first bona fide pathogenic RAD51D missense cancer susceptibility allele and supports the use of targeted PARP-inhibitor therapies in ovarian cancer patients carrying deleterious missense RAD51D variants. Cancer Res; 77(16); 4517-29. ©2017 AACR.

A Survey of Computational Tools to Analyze and Interpret Whole Exome Sequencing Data

Whole Exome Sequencing (WES) is the application of the next-generation technology to determine the variations in the exome and is becoming a standard approach in studying genetic variants in diseases. Understanding the exomes of individuals at single base resolution allows the identification of actionable mutations for disease treatment and management. WES technologies have shifted the bottleneck in experimental data production to computationally intensive informatics-based data analysis. Novel computational tools and methods have been developed to analyze and interpret WES data. Here, we review some of the current tools that are being used to analyze WES data. These tools range from the alignment of raw sequencing reads all the way to linking variants to actionable therapeutics. Strengths and weaknesses of each tool are discussed for the purpose of helping researchers make more informative decisions on selecting the best tools to analyze their WES data.

Molecular analyses reveal close similarities between small cell carcinoma of the ovary, hypercalcemic type and atypical teratoid/rhabdoid tumor

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is the most common undifferentiated ovarian malignancy diagnosed in women under age 40. We and others recently determined that germline and/or somatic deleterious mutations in SMARCA4 characterize SCCOHT. Alterations in this gene, or the related SWI/SNF chromatin remodeling gene SMARCB1, have been previously reported in atypical teratoid/rhabdoid tumors (ATRTs) and malignant rhabdoid tumors (MRTs). To further describe the somatic landscape of SCCOHT, we performed whole exome sequencing on 14 tumors and their matched normal tissues and compared their genomic alterations with those in ATRT and ovarian high grade serous carcinoma (HGSC). We confirmed that SMARCA4 is the only recurrently mutated gene in SCCOHT, and show that recurrent allelic imbalance is observed exclusively on chromosome 19p, where SMARCA4 resides. By comparing genomic alterations between SCCOHT, ATRT and HGSC, we demonstrate that SCCOHTs, like ATRTs, have a remarkably simple genome and harbor significantly fewer somatic protein-coding mutations and chromosomal alterations than HGSC. Furthermore, a comparison of global DNA methylation profiles of 45 SCCOHTs, 65 ATRTs, and 92 HGSCs demonstrates a strong epigenetic correlation between SCCOHT and ATRT. Our results further confirm that the genomic and epigenomic signatures of SCCOHT are more similar to those of ATRT than HGSC, supporting our previous hypothesis that SCCOHT is a rhabdoid tumor and should be renamed MRT of the ovary. Furthermore, we conclude that SMARCA4 inactivation is the main cause of SCCOHT, and that new distinct therapeutic approaches should be developed to specifically target this devastating tumor.

Germline and somatic FGFR1 abnormalities in dysembryoplastic neuroepithelial tumors

Dysembryoplastic neuroepithelial tumor (DNET) is a benign brain tumor associated with intractable drug-resistant epilepsy. In order to identify underlying genetic alterations and molecular mechanisms, we examined three family members affected by multinodular DNETs as well as 100 sporadic tumors from 96 patients, which had been referred to us as DNETs. We performed whole-exome sequencing on 46 tumors and targeted sequencing for hotspot FGFR1 mutations and BRAF p.V600E was used on the remaining samples. FISH, copy number variation assays and Sanger sequencing were used to validate the findings. By whole-exome sequencing of the familial cases, we identified a novel germline FGFR1 mutation, p.R661P. Somatic activating FGFR1 mutations (p.N546K or p.K656E) were observed in the tumor samples and further evidence for functional relevance was obtained by in silico modeling. The FGFR1 p.K656E mutation was confirmed to be in cis with the germline p.R661P variant. In 43 sporadic cases, in which the diagnosis of DNET could be confirmed on central blinded neuropathology review, FGFR1 alterations were also frequent and mainly comprised intragenic tyrosine kinase FGFR1 duplication and multiple mutants in cis (25/43; 58.1 %) while BRAF p.V600E alterations were absent (0/43). In contrast, in 53 cases, in which the diagnosis of DNET was not confirmed, FGFR1 alterations were less common (10/53; 19 %; p < 0.0001) and hotspot BRAF p.V600E (12/53; 22.6 %) (p < 0.001) prevailed. We observed overexpression of phospho-ERK in FGFR1 p.R661P and p.N546K mutant expressing HEK293 cells as well as FGFR1 mutated tumor samples, supporting enhanced MAP kinase pathway activation under these conditions. In conclusion, constitutional and somatic FGFR1 alterations and MAP kinase pathway activation are key events in the pathogenesis of DNET. These findings point the way towards existing targeted therapies.

Spatial and temporal homogeneity of driver mutations in diffuse intrinsic pontine glioma

Diffuse Intrinsic Pontine Gliomas (DIPGs) are deadly paediatric brain tumours where needle biopsies help guide diagnosis and targeted therapies. To address spatial heterogeneity, here we analyse 134 specimens from various neuroanatomical structures of whole autopsy brains from nine DIPG patients. Evolutionary reconstruction indicates histone 3 (H3) K27M--including H3.2K27M--mutations potentially arise first and are invariably associated with specific, high-fidelity obligate partners throughout the tumour and its spread, from diagnosis to end-stage disease, suggesting mutual need for tumorigenesis. These H3K27M ubiquitously-associated mutations involve alterations in TP53 cell-cycle (TP53/PPM1D) or specific growth factor pathways (ACVR1/PIK3R1). Later oncogenic alterations arise in sub-clones and often affect the PI3K pathway. Our findings are consistent with early tumour spread outside the brainstem including the cerebrum. The spatial and temporal homogeneity of main driver mutations in DIPG implies they will be captured by limited biopsies and emphasizes the need to develop therapies specifically targeting obligate oncohistone partnerships.

Whole-exome sequencing as a diagnostic tool: current challenges and future opportunities

Whole-exome sequencing (WES) represents a significant breakthrough in the field of human genetics. This technology has largely contributed to the identification of new disease-causing genes and is now entering clinical laboratories. WES represents a powerful tool for diagnosis and could reduce the 'diagnostic odyssey' for many patients. In this review, we present a technical overview of WES analysis, variants annotation and interpretation in a clinical setting. We evaluate the usefulness of clinical WES in different clinical indications, such as rare diseases, cancer and complex diseases. Finally, we discuss the efficacy of WES as a diagnostic tool and the impact on patient management.