CONSERTING: integrating copy-number analysis with structural-variation detection

Biologic and Clinical Analysis of Childhood Gamma Delta T-ALL Identifies LMO2/STAG2 Rearrangements as Extremely High Risk

Acute lymphoblastic leukemia expressing the gamma delta T-cell receptor (γδ T-ALL) is a poorly understood disease. We studied 200 children with γδ T-ALL from 13 clinical study groups to understand the clinical and genetic features of this disease. We found age and genetic drivers were significantly associated with outcome. γδ T-ALL diagnosed in children under 3 years of age was extremely high-risk and enriched for genetic alterations that result in both LMO2 activation and STAG2 inactivation. Mechanistically, using patient samples and isogenic cell lines, we show that inactivation of STAG2 profoundly perturbs chromatin organization by altering enhancer-promoter looping, resulting in deregulation of gene expression associated with T-cell differentiation. High-throughput drug screening identified a vulnerability in DNA repair pathways arising from STAG2 inactivation, which can be targeted by poly(ADP-ribose) polymerase inhibition. These data provide a diagnostic framework for classification and risk stratification of pediatric γδ T-ALL. Significance: Patients with acute lymphoblastic leukemia expressing the gamma delta T-cell receptor under 3 years old or measurable residual disease ≥1% at end of induction showed dismal outcomes and should be classified as having high-risk disease. The STAG2/LMO2 subtype was enriched in this very young age group. STAG2 inactivation may perturb chromatin conformation and cell differentiation and confer vulnerability to poly(ADP-ribose) polymerase inhibition.

A new genomic framework to categorize pediatric acute myeloid leukemia

Recent studies on pediatric acute myeloid leukemia (pAML) have revealed pediatric-specific driver alterations, many of which are underrepresented in the current classification schemas. To comprehensively define the genomic landscape of pAML, we systematically categorized 887 pAML into 23 mutually distinct molecular categories, including new major entities such as UBTF or BCL11B, covering 91.4% of the cohort. These molecular categories were associated with unique expression profiles and mutational patterns. For instance, molecular categories characterized by specific HOXA or HOXB expression signatures showed distinct mutation patterns of RAS pathway genes, FLT3 or WT1, suggesting shared biological mechanisms. We show that molecular categories were strongly associated with clinical outcomes using two independent cohorts, leading to the establishment of a new prognostic framework for pAML based on these updated molecular categories and minimal residual disease. Together, this comprehensive diagnostic and prognostic framework forms the basis for future classification of pAML and treatment strategies.

Genetic and epigenetic features of bilateral Wilms tumor predisposition in patients from the Children's Oncology Group AREN18B5-Q

Developing synchronous bilateral Wilms tumor suggests an underlying (epi)genetic predisposition. Here, we evaluate this predisposition in 68 patients using whole exome or genome sequencing (n = 85 tumors from 61 patients with matched germline blood DNA), RNA-seq (n = 99 tumors), and DNA methylation analysis (n = 61 peripheral blood, n = 29 non-diseased kidney, n = 99 tumors). We determine the predominant events for bilateral Wilms tumor predisposition: 1)pre-zygotic germline genetic variants readily detectable in blood DNA [WT1 (14.8%), NYNRIN (6.6%), TRIM28 (5%), and BRCA-related genes (5%)] or 2)post-zygotic epigenetic hypermethylation at 11p15.5 H19/ICR1 that may require analysis of multiple tissue types for diagnosis. Of 99 total tumor specimens, 16 (16.1%) have 11p15.5 normal retention of imprinting, 25 (25.2%) have 11p15.5 copy neutral loss of heterozygosity, and 58 (58.6%) have 11p15.5 H19/ICR1 epigenetic hypermethylation (loss of imprinting). Here, we ascertain the epigenetic and genetic modes of bilateral Wilms tumor predisposition.

T cell immune deficiency rather than chromosome instability predisposes patients with short telomere syndromes to squamous cancers

Patients with short telomere syndromes (STS) are predisposed to developing cancer, believed to stem from chromosome instability in neoplastic cells. We tested this hypothesis in a large cohort assembled over the last 20 years. We found that the only solid cancers to which patients with STS are predisposed are squamous cell carcinomas of the head and neck, anus, or skin, a spectrum reminiscent of cancers seen in patients with immunodeficiency. Whole-genome sequencing showed no increase in chromosome instability, such as translocations or chromothripsis. Moreover, STS-associated cancers acquired telomere maintenance mechanisms, including telomerase reverse transcriptase (TERT) promoter mutations. A detailed study of the immune status of patients with STS revealed a striking T cell immunodeficiency at the time of cancer diagnosis. A similar immunodeficiency that impaired tumor surveillance was documented in mice with short telomeres. We conclude that STS patients’ predisposition to solid cancers is due to T cell exhaustion rather than autonomous defects in the neoplastic cells themselves.

Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities

Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to NOTCH inhibitors, which corresponded to high expression of NOTCH receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

The Genetic and Epigenetic Features of Bilateral Wilms Tumor Predisposition: A Report from the Children's Oncology Group AREN18B5-Q Study

This study comprehensively evaluated the landscape of genetic and epigenetic events that predispose to synchronous bilateral Wilms tumor (BWT). We performed whole exome or whole genome sequencing, total-strand RNA-seq, and DNA methylation analysis using germline and/or tumor samples from 68 patients with BWT from St. Jude Children's Research Hospital and the Children's Oncology Group. We found that 25/61 (41%) of patients evaluated harbored pathogenic or likely pathogenic germline variants, with WT1 (14.8%), NYNRIN (6.6%), TRIM28 (5%) and the BRCA-related genes (5%) BRCA1, BRCA2, and PALB2 being most common. Germline WT1 variants were strongly associated with somatic paternal uniparental disomy encompassing the 11p15.5 and 11p13/WT1 loci and subsequent acquired pathogenic CTNNB1 variants. Somatic coding variants or genome-wide copy number alterations were almost never shared between paired synchronous BWT, suggesting that the acquisition of independent somatic variants leads to tumor formation in the context of germline or early embryonic, post-zygotic initiating events. In contrast, 11p15.5 status (loss of heterozygosity, loss or retention of imprinting) was shared among paired synchronous BWT in all but one case. The predominant molecular events for BWT predisposition include pathogenic germline variants or post-zygotic epigenetic hypermethylation at the 11p15.5 H19/ICR1 locus (loss of imprinting). This study demonstrates that post-zygotic somatic mosaicism for 11p15.5 hypermethylation/loss of imprinting is the single most common initiating molecular event predisposing to BWT. Evidence of somatic mosaicism for 11p15.5 loss of imprinting was detected in leukocytes of a cohort of BWT patients and long-term survivors, but not in unilateral Wilms tumor patients and long-term survivors or controls, further supporting the hypothesis that post-zygotic 11p15.5 alterations occurred in the mesoderm of patients who go on to develop BWT. Due to the preponderance of BWT patients with demonstrable germline or early embryonic tumor predisposition, BWT exhibits a unique biology when compared to unilateral Wilms tumor and therefore warrants continued refinement of its own treatment-relevant biomarkers which in turn may inform directed treatment strategies in the future.

RNAseqCNV: analysis of large-scale copy number variations from RNA-seq data

Transcriptome sequencing (RNA-seq) is widely used to detect gene rearrangements and quantitate gene expression in acute lymphoblastic leukemia (ALL), but its utility and accuracy in identifying copy number variations (CNVs) has not been well described. CNV information inferred from RNA-seq can be highly informative to guide disease classification and risk stratification in ALL due to the high incidence of aneuploid subtypes within this disease. Here we describe RNAseqCNV, a method to detect large scale CNVs from RNA-seq data. We used models based on normalized gene expression and minor allele frequency to classify arm level CNVs with high accuracy in ALL (99.1% overall and 98.3% for non-diploid chromosome arms, respectively), and the models were further validated with excellent performance in acute myeloid leukemia (accuracy 99.8% overall and 99.4% for non-diploid chromosome arms). RNAseqCNV outperforms alternative RNA-seq based algorithms in calling CNVs in the ALL dataset, especially in samples with a high proportion of CNVs. The CNV calls were highly concordant with DNA-based CNV results and more reliable than conventional cytogenetic-based karyotypes. RNAseqCNV provides a method to robustly identify copy number alterations in the absence of DNA-based analyses, further enhancing the utility of RNA-seq to classify ALL subtype.

Integrated Genomic Analysis Identifies UBTF Tandem Duplications as a Recurrent Lesion in Pediatric Acute Myeloid Leukemia

The genetics of relapsed pediatric acute myeloid leukemia (AML) has yet to be comprehensively defined. Here, we present the spectrum of genomic alterations in 136 relapsed pediatric AMLs. We identified recurrent exon 13 tandem duplications (TD) in upstream binding transcription factor (UBTF) in 9% of relapsed AML cases. UBTF-TD AMLs commonly have normal karyotype or trisomy 8 with cooccurring WT1 mutations or FLT3-ITD but not other known oncogenic fusions. These UBTF-TD events are stable during disease progression and are present in the founding clone. In addition, we observed that UBTF-TD AMLs account for approximately 4% of all de novo pediatric AMLs, are less common in adults, and are associated with poor outcomes and MRD positivity. Expression of UBTF-TD in primary hematopoietic cells is sufficient to enhance serial clonogenic activity and to drive a similar transcriptional program to UBTF-TD AMLs. Collectively, these clinical, genomic, and functional data establish UBTF-TD as a new recurrent mutation in AML.

SIGNIFICANCE

We defined the spectrum of mutations in relapsed pediatric AML and identified UBTF-TDs as a new recurrent genetic alteration. These duplications are more common in children and define a group of AMLs with intermediate-risk cytogenetic abnormalities, FLT3-ITD and WT1 alterations, and are associated with poor outcomes. See related commentary by Hasserjian and Nardi, p. 173. This article is highlighted in the In This Issue feature, p. 171.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

Comprehensive molecular characterization of pediatric radiation-induced high-grade glioma

Radiation-induced high-grade gliomas (RIGs) are an incurable late complication of cranial radiation therapy. We performed DNA methylation profiling, RNA-seq, and DNA sequencing on 32 RIG tumors and an in vitro drug screen in two RIG cell lines. We report that based on DNA methylation, RIGs cluster primarily with the pediatric receptor tyrosine kinase I high-grade glioma subtype. Common copy-number alterations include Chromosome (Ch.) 1p loss/1q gain, and Ch. 13q and Ch. 14q loss; focal alterations include PDGFRA and CDK4 gain and CDKN2A and BCOR loss. Transcriptomically, RIGs comprise a stem-like subgroup with lesser mutation burden and Ch. 1p loss and a pro-inflammatory subgroup with greater mutation burden and depleted DNA repair gene expression. Chromothripsis in several RIG samples is associated with extrachromosomal circular DNA-mediated amplification of PDGFRA and CDK4. Drug screening suggests microtubule inhibitors/stabilizers, DNA-damaging agents, MEK inhibition, and, in the inflammatory subgroup, proteasome inhibitors, as potentially effective therapies.

Radiation-induced high-grade gliomas (RIGs) are an incurable late complication of cranial radiation therapy. In the largest study to date, we report the results of DNA methylation profiling, RNA-Seq and genomic sequencing of 32 RIG tumors, and an in vitro drug screen in two RIG cell lines.

Chemotherapy and mismatch repair deficiency cooperate to fuel TP53 mutagenesis and ALL relapse

Chemotherapy is a standard treatment for pediatric acute lymphoblastic leukemia (ALL), which sometimes relapses with chemoresistant features. However, whether acquired drug-resistance mutations in relapsed ALL pre-exist or are induced by treatment remains unknown. Here we provide direct evidence of a specific mechanism by which chemotherapy induces drug-resistance-associated mutations leading to relapse. Using genomic and functional analysis of relapsed ALL we show that thiopurine treatment in mismatch repair (MMR)-deficient leukemias induces hotspot TP53 R248Q mutations through a specific mutational signature (thio-dMMR). Clonal evolution analysis reveals sequential MMR inactivation followed by TP53 mutation in some patients with ALL. Acquired TP53 R248Q mutations are associated with on-treatment relapse, poor treatment response and resistance to multiple chemotherapeutic agents, which could be reversed by pharmacological p53 reactivation. Our findings indicate that TP53 R248Q in relapsed ALL originates through synergistic mutagenesis from thiopurine treatment and MMR deficiency and suggest strategies to prevent or treat TP53-mutant relapse.

Retinoblastoma from human stem cell-derived retinal organoids

Retinoblastoma is a childhood cancer of the developing retina that initiates with biallelic inactivation of the RB1 gene. Children with germline mutations in RB1 have a high likelihood of developing retinoblastoma and other malignancies later in life. Genetically engineered mouse models of retinoblastoma share some similarities with human retinoblastoma but there are differences in their cellular differentiation. To develop a laboratory model of human retinoblastoma formation, we make induced pluripotent stem cells (iPSCs) from 15 participants with germline RB1 mutations. Each of the stem cell lines is validated, characterized and then differentiated into retina using a 3-dimensional organoid culture system. After 45 days in culture, the retinal organoids are dissociated and injected into the vitreous of eyes of immunocompromised mice to support retinoblastoma tumor growth. Retinoblastomas formed from retinal organoids made from patient-derived iPSCs have molecular, cellular and genomic features indistinguishable from human retinoblastomas. This model of human cancer based on patient-derived iPSCs with germline cancer predisposing mutations provides valuable insights into the cellular origins of this debilitating childhood disease as well as the mechanism of tumorigenesis following RB1 gene inactivation.

Retinoblastoma is a heritable pediatric cancer driven by mutations in RB1. Here, the authors demonstrate the first patient derived model of retinoblastoma using iPSCs from patients with germline mutations in RB1.

Genomes for Kids: The scope of pathogenic mutations in pediatric cancer revealed by comprehensive DNA and RNA sequencing

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole genome (WGS), exome, and RNA sequencing, to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically-relevant (25%), and/or cancer predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors) and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers.

Patient-derived models recapitulate heterogeneity of molecular signatures and drug response in pediatric high-grade glioma

Pediatric high-grade glioma (pHGG) is a major contributor to cancer-related death in children. In vitro and in vivo disease models reflecting the intimate connection between developmental context and pathogenesis of pHGG are essential to advance understanding and identify therapeutic vulnerabilities. Here we report establishment of 21 patient-derived pHGG orthotopic xenograft (PDOX) models and eight matched cell lines from diverse groups of pHGG. These models recapitulate histopathology, DNA methylation signatures, mutations and gene expression patterns of the patient tumors from which they were derived, and include rare subgroups not well-represented by existing models. We deploy 16 new and existing cell lines for high-throughput screening (HTS). In vitro HTS results predict variable in vivo response to PI3K/mTOR and MEK pathway inhibitors. These unique new models and an online interactive data portal for exploration of associated detailed molecular characterization and HTS chemical sensitivity data provide a rich resource for pediatric brain tumor research.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Integrative reconstruction of cancer genome karyotypes using InfoGenomeR

Annotation of structural variations (SVs) and base-level karyotyping in cancer cells remains challenging. Here, we present Integrative Framework for Genome Reconstruction (InfoGenomeR)-a graph-based framework that can reconstruct individual SVs into karyotypes based on whole-genome sequencing data, by integrating SVs, total copy number alterations, allele-specific copy numbers, and haplotype information. Using whole-genome sequencing data sets of patients with breast cancer, glioblastoma multiforme, and ovarian cancer, we demonstrate the analytical potential of InfoGenomeR. We identify recurrent derivative chromosomes derived from chromosomes 11 and 17 in breast cancer samples, with homogeneously staining regions for CCND1 and ERBB2, and double minutes and breakage-fusion-bridge cycles in glioblastoma multiforme and ovarian cancer samples, respectively. Moreover, we show that InfoGenomeR can discriminate private and shared SVs between primary and metastatic cancer sites that could contribute to tumour evolution. These findings indicate that InfoGenomeR can guide targeted therapies by unravelling cancer-specific SVs on a genome-wide scale.

Network-based systems pharmacology reveals heterogeneity in LCK and BCL2 signaling and therapeutic sensitivity of T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia' drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

The acquisition of molecular drivers in pediatric therapy-related myeloid neoplasms

Pediatric therapy-related myeloid neoplasms (tMN) occur in children after exposure to cytotoxic therapy and have a dismal prognosis. The somatic and germline genomic alterations that drive these myeloid neoplasms in children and how they arise have yet to be comprehensively described. We use whole exome, whole genome, and/or RNA sequencing to characterize the genomic profile of 84 pediatric tMN cases (tMDS: n = 28, tAML: n = 56). Our data show that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements are frequent somatic drivers, and we identify cases with aberrant MECOM expression secondary to enhancer hijacking. Unlike adults with tMN, we find no evidence of pre-existing minor tMN clones (including those with TP53 mutations), but rather the majority of cases are unrelated clones arising as a consequence of cytotoxic therapy. These studies also uncover rare cases of lineage switch disease rather than true secondary neoplasms.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.

St. Jude Cloud: A Pediatric Cancer Genomic Data-Sharing Ecosystem

Effective data sharing is key to accelerating research to improve diagnostic precision, treatment efficacy, and long-term survival in pediatric cancer and other childhood catastrophic diseases. We present St. Jude Cloud (https://www.stjude.cloud), a cloud-based data-sharing ecosystem for accessing, analyzing, and visualizing genomic data from >10,000 pediatric patients with cancer and long-term survivors, and >800 pediatric sickle cell patients. Harmonized genomic data totaling 1.25 petabytes are freely available, including 12,104 whole genomes, 7,697 whole exomes, and 2,202 transcriptomes. The resource is expanding rapidly, with regular data uploads from St. Jude's prospective clinical genomics programs. Three interconnected apps within the ecosystem-Genomics Platform, Pediatric Cancer Knowledgebase, and Visualization Community-enable simultaneously performing advanced data analysis in the cloud and enhancing the Pediatric Cancer knowledgebase. We demonstrate the value of the ecosystem through use cases that classify 135 pediatric cancer subtypes by gene expression profiling and map mutational signatures across 35 pediatric cancer subtypes. SIGNIFICANCE: To advance research and treatment of pediatric cancer, we developed St. Jude Cloud, a data-sharing ecosystem for accessing >1.2 petabytes of raw genomic data from >10,000 pediatric patients and survivors, innovative analysis workflows, integrative multiomics visualizations, and a knowledgebase of published data contributed by the global pediatric cancer community.This article is highlighted in the In This Issue feature, p. 995.

Novel somatic alterations underlie Chinese papillary thyroid carcinoma

To characterize the somatic alterations of papillary thyroid carcinomas (PTC) in Chinese patients, we performed the next-generation-sequencing (NGS) study of the tumor-normal pairs of DNA and RNA samples extracted from 16 Chinese PTC patients. The whole genome sequencing (WGS) and transcriptome sequencing (RNA-seq) were conducted for 6 patients who were either current or former smokers and the whole exome sequencing (WES) and RNA-seq were conducted for another 10 patients who were never smokers. The NGS data were analyzed to identify somatic alteration events that may underlie PTC in Chinese patients. We identified a number of PTC driver genes harboring somatic driver mutations with significant functional impact such as COL11A1, TP53, PLXNA4, UBA1, AHNAK, CSMD2 and TTLL5 etc. Significant driver pathways underlying PTC were found, namely, the metabolic pathway, the pathway in cancer, the olfactory transduction pathway and the calcium signaling pathway. In addition, this study revealed genes with significant somatic copy number aberrations and corresponding somatic gene expression changes in PTC tumors, the most promising ones being BRD9, TRIP13, FZD3, and TFDP1 etc. We also identified several structural variants of PTCs, especially the novel in-frame fusion proteins such as TRNAU1AP-RCC1, RAB3GAP1-R3HDM1, and ENAH-ZSWIM5. Our study provided a list of novel PTC candidate genes with somatic alterations that may function as biomarkers for PTC in Chinese patients. The follow-up mechanism studies may be conducted based on the findings from this study.

Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome data from 702 neuroblastoma samples. Forty percent of samples harbor at least one recurrent driver gene alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years. COSMIC mutational signature 18, previously associated with reactive oxygen species, is the most common cause of driver point mutations in neuroblastoma, including most ALK and Ras-activating variants. Signature 18 appears early and is continuous throughout disease evolution. Signature 18 is enriched in neuroblastomas with MYCN amplification, 17q gain, and increased expression of mitochondrial ribosome and electron transport-associated genes. Recurrent FGFR1 variants in six patients, and ALK N-terminal structural alterations in five samples, identify additional patients potentially amenable to precision therapy.

FPGS relapse-specific mutations in relapsed childhood acute lymphoblastic leukemia

Although the cure rate for childhood acute lymphoblastic leukemia (ALL) has exceeded 80% with contemporary therapy, relapsed ALL remains a leading cause of cancer-related death in children. Relapse-specific mutations can be identified by comprehensive genome sequencing and might have clinical significance. Applying whole-exome sequencing to eight triplicate samples, we identified in one patient relapse-specific mutations in the folylpolyglutamate synthetase (FPGS) gene, whose product catalyzes the addition of multiple glutamate residues (polyglutamation) to methotrexate upon their entry into the cells. To determine the prevalence of mutations of the FPGS mutations, and those of two important genes in the thiopurine pathway, NT5C2 and PRPS1, we studied 299 diagnostic and 73 relapsed samples in 372 patients. Three more FPGS mutants were identified in two patients, NT5C2 mutations in six patients, and PRPS1 mutants in two patients. One patient had both NT5C2 and PRPS1 mutants. None of these alterations were detected at diagnosis with a sequencing depth of 1000X, suggesting that treatment pressure led to increased prevalence of mutations during therapy. Functional characterization of the FPGS mutants showed that they directly resulted in decreased enzymatic activity, leading to significant reduction in methotrexate polyglutamation, and therefore likely contributed to drug resistance and relapse in these cases. Thus, besides genomic alterations in thiopurine metabolizing enzymes, the relapse-specific mutations of FPGS represent another critical mechanism of acquired antimetabolite drug resistance in relapsed childhood ALL.

Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X

We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.

MYCN amplification and ATRX mutations are incompatible in neuroblastoma

Aggressive cancers often have activating mutations in growth-controlling oncogenes and inactivating mutations in tumor-suppressor genes. In neuroblastoma, amplification of the MYCN oncogene and inactivation of the ATRX tumor-suppressor gene correlate with high-risk disease and poor prognosis. Here we show that ATRX mutations and MYCN amplification are mutually exclusive across all ages and stages in neuroblastoma. Using human cell lines and mouse models, we found that elevated MYCN expression and ATRX mutations are incompatible. Elevated MYCN levels promote metabolic reprogramming, mitochondrial dysfunction, reactive-oxygen species generation, and DNA-replicative stress. The combination of replicative stress caused by defects in the ATRX-histone chaperone complex, and that induced by MYCN-mediated metabolic reprogramming, leads to synthetic lethality. Therefore, ATRX and MYCN represent an unusual example, where inactivation of a tumor-suppressor gene and activation of an oncogene are incompatible. This synthetic lethality may eventually be exploited to improve outcomes for patients with high-risk neuroblastoma.

Therapy-induced mutations drive the genomic landscape of relapsed acute lymphoblastic leukemia

To study the mechanisms of relapse in acute lymphoblastic leukemia (ALL), we performed whole-genome sequencing of 103 diagnosis-relapse-germline trios and ultra-deep sequencing of 208 serial samples in 16 patients. Relapse-specific somatic alterations were enriched in 12 genes (NR3C1, NR3C2, TP53, NT5C2, FPGS, CREBBP, MSH2, MSH6, PMS2, WHSC1, PRPS1, and PRPS2) involved in drug response. Their prevalence was 17% in very early relapse (<9 months from diagnosis), 65% in early relapse (9-36 months), and 32% in late relapse (>36 months) groups. Convergent evolution, in which multiple subclones harbor mutations in the same drug resistance gene, was observed in 6 relapses and confirmed by single-cell sequencing in 1 case. Mathematical modeling and mutational signature analysis indicated that early relapse resistance acquisition was frequently a 2-step process in which a persistent clone survived initial therapy and later acquired bona fide resistance mutations during therapy. In contrast, very early relapses arose from preexisting resistant clone(s). Two novel relapse-specific mutational signatures, one of which was caused by thiopurine treatment based on in vitro drug exposure experiments, were identified in early and late relapses but were absent from 2540 pan-cancer diagnosis samples and 129 non-ALL relapses. The novel signatures were detected in 27% of relapsed ALLs and were responsible for 46% of acquired resistance mutations in NT5C2, PRPS1, NR3C1, and TP53. These results suggest that chemotherapy-induced drug resistance mutations facilitate a subset of pediatric ALL relapses.

Patient-derived orthotopic xenografts of pediatric brain tumors: a St. Jude resource

Pediatric brain tumors are the leading cause of cancer-related death in children. Patient-derived orthotopic xenografts (PDOX) of childhood brain tumors have recently emerged as a biologically faithful vehicle for testing novel and more effective therapies. Herein, we provide the histopathological and molecular analysis of 37 novel PDOX models generated from pediatric brain tumor patients treated at St. Jude Children's Research Hospital. Using a combination of histopathology, whole-genome and whole-exome sequencing, RNA-sequencing, and DNA methylation arrays, we demonstrate the overall fidelity and inter-tumoral molecular heterogeneity of pediatric brain tumor PDOX models. These models represent frequent as well as rare childhood brain tumor entities, including medulloblastoma, ependymoma, atypical teratoid rhabdoid tumor, and embryonal tumor with multi-layer rosettes. PDOX models will be valuable platforms for evaluating novel therapies and conducting pre-clinical trials to accelerate progress in the treatment of brain tumors in children. All described PDOX models and associated datasets can be explored using an interactive web-based portal and will be made freely available to the research community upon request.

Population-wide copy number variation calling using variant call format files from 6,898 individuals

Copy number variants (CNVs) play an important role in a number of human diseases, but the accurate calling of CNVs remains challenging. Most current approaches to CNV detection use raw read alignments, which are computationally intensive to process. We use a regression tree-based approach to call germline CNVs from whole-genome sequencing (WGS, >18x) variant call sets in 6,898 samples across four European cohorts, and describe a rich large variation landscape comprising 1,320 CNVs. Eighty-one percent of detected events have been previously reported in the Database of Genomic Variants. Twenty-three percent of high-quality deletions affect entire genes, and we recapitulate known events such as the GSTM1 and RHD gene deletions. We test for association between the detected deletions and 275 protein levels in 1,457 individuals to assess the potential clinical impact of the detected CNVs. We describe complex CNV patterns underlying an association with levels of the CCL3 protein (MAF = 0.15, p = 3.6x10^-12 ) at the CCL3L3 locus, and a novel cis-association between a low-frequency NOMO1 deletion and NOMO1 protein levels (MAF = 0.02, p = 2.2x10^-7 ). This study demonstrates that existing population-wide WGS call sets can be mined for germline CNVs with minimal computational overhead, delivering insight into a less well-studied, yet potentially impactful class of genetic variant.

VCF2CNA: A tool for efficiently detecting copy-number alterations in VCF genotype data and tumor purity

VCF2CNA is a tool (Linux commandline or web-interface) for copy-number alteration (CNA) analysis and tumor purity estimation of paired tumor-normal VCF variant file formats. It operates on whole genome and whole exome datasets. To benchmark its performance, we applied it to 46 adult glioblastoma and 146 pediatric neuroblastoma samples sequenced by Illumina and Complete Genomics (CGI) platforms respectively. VCF2CNA was highly consistent with a state-of-the-art algorithm using raw sequencing data (mean F1-score = 0.994) in high-quality whole genome glioblastoma samples and was robust to uneven coverage introduced by library artifacts. In the whole genome neuroblastoma set, VCF2CNA identified MYCN high-level amplifications in 31 of 32 clinically validated samples compared to 15 found by CGI’s HMM-based CNA model. Moreover, VCF2CNA achieved highly consistent CNA profiles between WGS and WXS platforms (mean F1 score 0.97 on a set of 15 rhabdomyosarcoma samples). In addition, VCF2CNA provides accurate tumor purity estimates for samples with sufficient CNAs. These results suggest that VCF2CNA is an accurate, efficient and platform-independent tool for CNA and tumor purity analyses without accessing raw sequence data.

Germline deletion of ETV6 in familial acute lymphoblastic leukemia

Recent studies have identified germline mutations in TP53, PAX5, ETV6, and IKZF1 in kindreds with familial acute lymphoblastic leukemia (ALL), but the genetic basis of ALL in many kindreds is unknown despite mutational analysis of the exome. Here, we report a germline deletion of ETV6 identified by linkage and structural variant analysis of whole-genome sequencing data segregating in a kindred with thrombocytopenia, B-progenitor acute lymphoblastic leukemia, and diffuse large B-cell lymphoma. The 75-nt deletion removed the ETV6 exon 7 splice acceptor, resulting in exon skipping and protein truncation. The ETV6 deletion was also identified by optimal structural variant analysis of exome sequencing data. These findings identify a new mechanism of germline predisposition in ALL and implicate ETV6 germline variation in predisposition to lymphoma. Importantly, these data highlight the importance of germline structural variant analysis in the search for germline variants predisposing to familial leukemia.

Analysis of error profiles in deep next-generation sequencing data

BACKGROUND

Sequencing errors are key confounding factors for detecting low-frequency genetic variants that are important for cancer molecular diagnosis, treatment, and surveillance using deep next-generation sequencing (NGS). However, there is a lack of comprehensive understanding of errors introduced at various steps of a conventional NGS workflow, such as sample handling, library preparation, PCR enrichment, and sequencing. In this study, we use current NGS technology to systematically investigate these questions.

RESULTS

By evaluating read-specific error distributions, we discover that the substitution error rate can be computationally suppressed to 10-5 to 10-4, which is 10- to 100-fold lower than generally considered achievable (10-3) in the current literature. We then quantify substitution errors attributable to sample handling, library preparation, enrichment PCR, and sequencing by using multiple deep sequencing datasets. We find that error rates differ by nucleotide substitution types, ranging from 10-5 for A>C/T>G, C>A/G>T, and C>G/G>C changes to 10-4 for A>G/T>C changes. Furthermore, C>T/G>A errors exhibit strong sequence context dependency, sample-specific effects dominate elevated C>A/G>T errors, and target-enrichment PCR led to ~ 6-fold increase of overall error rate. We also find that more than 70% of hotspot variants can be detected at 0.1 ~ 0.01% frequency with the current NGS technology by applying in silico error suppression.

CONCLUSIONS

We present the first comprehensive analysis of sequencing error sources in conventional NGS workflows. The error profiles revealed by our study highlight new directions for further improving NGS analysis accuracy both experimentally and computationally, ultimately enhancing the precision of deep sequencing.

Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas

Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents, and ranges on the spectrum of malignancy from low grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF^1,2. However, in approximately 50% of cases no genetic driver has been established². Clinical whole-genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine-threonine kinase that activates MEK^3,4. The patient, who had exhausted all other therapeutic options, was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations of MAP3K8 in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements-uncovered by comprehensive clinical sequencing of a single case-are the most common genetic event in spitzoid melanoma, are present in adult melanomas and could be amenable to MEK inhibition.

Insights of Acute Lymphoblastic Leukemia with Development of Genomic Investigation

Treatment outcomes for acute lymphoblastic leukemia (ALL), especially pediatric ALL, have greatly improved due to the risk-adapted therapy. Combination of drug development, clinical practice, as well as basic genetic researches has brought the survival rate of ALL from less than 10% to more than 90% today, not only increasing the treatment efficacy but also limiting adverse drug reactions (ADRs). In this review, we summarized the landscape identification of ALL genetic alterations, which provided the opportunity to increase the survival rate and especially minimize the relapse risk of ALL, and highlighted the importance of the development of new technologies of genomic investigation for translational medicine.

The Clonal Evolution of Metastatic Osteosarcoma as Shaped by Cisplatin Treatment

To investigate the genomic evolution of metastatic pediatric osteosarcoma, we performed whole-genome and targeted deep sequencing on 14 osteosarcoma metastases and two primary tumors from four patients (two to eight samples per patient). All four patients harbored ancestral (truncal) somatic variants resulting in TP53 inactivation and cell-cycle aberrations, followed by divergence into relapse-specific lineages exhibiting a cisplatin-induced mutation signature. In three of the four patients, the cisplatin signature accounted for >40% of mutations detected in the metastatic samples. Mutations potentially acquired during cisplatin treatment included NF1 missense mutations of uncertain significance in two patients and a KIT G565R activating mutation in one patient. Three of four patients demonstrated widespread ploidy differences between samples from the sample patient. Single-cell seeding of metastasis was detected in most metastatic samples. Cross-seeding between metastatic sites was observed in one patient, whereas in another patient a minor clone from the primary tumor seeded both metastases analyzed. These results reveal extensive clonal heterogeneity in metastatic osteosarcoma, much of which is likely cisplatin-induced. IMPLICATIONS: The extent and consequences of chemotherapy-induced damage in pediatric cancers is unknown. We found that cisplatin treatment can potentially double the mutational burden in osteosarcoma, which has implications for optimizing therapy for recurrent, chemotherapy-resistant disease.

Structure and evolution of double minutes in diagnosis and relapse brain tumors

Double minute chromosomes are extrachromosomal circular DNA fragments frequently found in brain tumors. To understand their evolution, we characterized the double minutes in paired diagnosis and relapse tumors from a pediatric high-grade glioma and four adult glioblastoma patients. We determined the full structures of the major double minutes using a novel approach combining multiple types of supporting genomic evidence. Among the double minutes identified in the pediatric patient, only one carrying EGFR was maintained at high abundance in both samples, whereas two others were present in only trace amounts at diagnosis but abundant at relapse, and the rest were found either in the relapse sample only or in the diagnosis sample only. For the EGFR-carrying double minutes, we found a secondary somatic deletion in all copies at relapse, after erlotinib treatment. However, the somatic mutation was present at very low frequency at diagnosis, suggesting potential resistance to the EGFR inhibitor. This mutation caused an in-frame RNA transcript to skip exon 16, a novel transcript isoform absent in EST database, as well as about 700 RNA-seq of normal brains that we reviewed. We observed similar patterns involving longitudinal copy number shift of double minutes in another four pairs (diagnosis/relapse) of adult glioblastoma. Overall, in three of five paired tumor samples, we found that although the same oncogenes were amplified at diagnosis and relapse, they were amplified on different double minutes. Our results suggest that double minutes readily evolve, increasing tumor heterogeneity rapidly. Understanding patterns of double minute evolution can shed light on future therapeutic solutions to brain tumors carrying such variants.

Clinical cancer genomic profiling by three-platform sequencing of whole genome, whole exome and transcriptome

To evaluate the potential of an integrated clinical test to detect diverse classes of somatic and germline mutations relevant to pediatric oncology, we performed three-platform whole-genome (WGS), whole exome (WES) and transcriptome (RNA-Seq) sequencing of tumors and normal tissue from 78 pediatric cancer patients in a CLIA-certified, CAP-accredited laboratory. Our analysis pipeline achieves high accuracy by cross-validating variants between sequencing types, thereby removing the need for confirmatory testing, and facilitates comprehensive reporting in a clinically-relevant timeframe. Three-platform sequencing has a positive predictive value of 97–99, 99, and 91% for somatic SNVs, indels and structural variations, respectively, based on independent experimental verification of 15,225 variants. We report 240 pathogenic variants across all cases, including 84 of 86 known from previous diagnostic testing (98% sensitivity). Combined WES and RNA-Seq, the current standard for precision oncology, achieved only 78% sensitivity. These results emphasize the critical need for incorporating WGS in pediatric oncology testing.

Clinical oncology is rapidly adopting next-generation sequencing technology for nucleotide variant and indel detection. Here the authors present a three-platform approach (whole-genome, whole-exome, and whole-transcriptome) in pediatric patients for the detection of diverse types of germline and somatic variants.

SEG - A Software Program for Finding Somatic Copy Number Alterations in Whole Genome Sequencing Data of Cancer

As next-generation sequencing technology advances and the cost decreases, whole genome sequencing (WGS) has become the preferred platform for the identification of somatic copy number alteration (CNA) events in cancer genomes. To more effectively decipher these massive sequencing data, we developed a software program named SEG, shortened from the word “segment”. SEG utilizes mapped read or fragment density for CNA discovery. To reduce CNA artifacts arisen from sequencing and mapping biases, SEG first normalizes the data by taking the log₂-ratio of each tumor density against its matching normal density. SEG then uses dynamic programming to find change-points among a contiguous log₂-ratio data series along a chromosome, dividing the chromosome into different segments. SEG finally identifies those segments having CNA. Our analyses with both simulated and real sequencing data indicate that SEG finds more small CNAs than other published software tools.

Prognostic significance of frequent CLDN18-ARHGAP26/6 fusion in gastric signet-ring cell cancer

Signet-ring cell carcinoma (SRCC) has specific epidemiology and oncogenesis in gastric cancer, however, with no systematical investigation for prognostic genomic features. Here we report a systematic investigation conducted in 1868 Chinese gastric cancer patients indicating that signet-ring cells content was related to multiple clinical characteristics and treatment outcomes. We thus perform whole-genome sequencing on 32 pairs of SRC samples, and identify frequent CLDN18-ARHGAP26/6 fusion (25%). With 797 additional patients for validation, prevalence of CLDN18-ARHGAP26/6 fusion is noticed to be associated with signet-ring cell content, age at diagnosis, female/male ratio, and TNM stage. Importantly, patients with CLDN18-ARHGAP26/6 fusion have worse survival outcomes, and get no benefit from oxaliplatin/fluoropyrimidines-based chemotherapy, which is consistent with the fact of chemo-drug resistance acquired in CLDN18-ARHGAP26 introduced cell lines. Overall, this study provides insights into the clinical and genomic features of SRCC, and highlights the importance of frequent CLDN18-ARHGAP26/6 fusions in chemotherapy response for SRCC.

Germline Genetic IKZF1 Variation and Predisposition to Childhood Acute Lymphoblastic Leukemia

Somatic genetic alterations of IKZF1, which encodes the lymphoid transcription factor IKAROS, are common in high-risk B-progenitor acute lymphoblastic leukemia (ALL) and are associated with poor prognosis. Such alterations result in the acquisition of stem cell-like features, overexpression of adhesion molecules causing aberrant cell-cell and cell-stroma interaction, and decreased sensitivity to tyrosine kinase inhibitors. Here we report coding germline IKZF1 variation in familial childhood ALL and 0.9% of presumed sporadic B-ALL, identifying 28 unique variants in 45 children. The majority of variants adversely affected IKZF1 function and drug responsiveness of leukemic cells. These results identify IKZF1 as a leukemia predisposition gene, and emphasize the importance of germline genetic variation in the development of both familial and sporadic ALL.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

The genomic landscape of pediatric myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are uncommon in children and have a poor prognosis. In contrast to adult MDS, little is known about the genomic landscape of pediatric MDS. Here, we describe the somatic and germline changes of pediatric MDS using whole exome sequencing, targeted amplicon sequencing, and/or RNA-sequencing of 46 pediatric primary MDS patients. Our data show that, in contrast to adult MDS, Ras/MAPK pathway mutations are common in pediatric MDS (45% of primary cohort), while mutations in RNA splicing genes are rare (2% of primary cohort). Surprisingly, germline variants in SAMD9 or SAMD9L were present in 17% of primary MDS patients, and these variants were routinely lost in the tumor cells by chromosomal deletions (e.g., monosomy 7) or copy number neutral loss of heterozygosity (CN-LOH). Our data confirm that adult and pediatric MDS are separate diseases with disparate mechanisms, and that SAMD9/SAMD9L mutations represent a new class of MDS predisposition.

Myelodysplastic syndromes (MDS) are uncommon in children and have poor prognosis. Here, the authors interrogate the genomic landscape of MDS, confirming adult and paediatric MDS are separate diseases with disparate mechanisms, and highlighting that SAMD9/SAMD9L mutations represent a new class of MDS predisposition.

GRIDSS: sensitive and specific genomic rearrangement detection using positional de Bruijn graph assembly

The identification of genomic rearrangements with high sensitivity and specificity using massively parallel sequencing remains a major challenge, particularly in precision medicine and cancer research. Here, we describe a new method for detecting rearrangements, GRIDSS (Genome Rearrangement IDentification Software Suite). GRIDSS is a multithreaded structural variant (SV) caller that performs efficient genome-wide break-end assembly prior to variant calling using a novel positional de Bruijn graph-based assembler. By combining assembly, split read, and read pair evidence using a probabilistic scoring, GRIDSS achieves high sensitivity and specificity on simulated, cell line, and patient tumor data, recently winning SV subchallenge #5 of the ICGC-TCGA DREAM8.5 Somatic Mutation Calling Challenge. On human cell line data, GRIDSS halves the false discovery rate compared to other recent methods while matching or exceeding their sensitivity. GRIDSS identifies nontemplate sequence insertions, microhomologies, and large imperfect homologies, estimates a quality score for each breakpoint, stratifies calls into high or low confidence, and supports multisample analysis.

Clear cell sarcoma of kidney involving a horseshoe kidney and harboring EGFR internal tandem duplication

Clear cell sarcoma of kidney (CCSK) is a rare renal malignancy, previously unreported in horseshoe kidney (HSK). B-cell lymphoma 6 corepressor (BCOR) gene internal tandem duplication (ITD) was identified as a recurrent somatic alteration in approximately 85% of CCSKs. This and the YWHAE-NUTM2B/E fusion, the second most common recurrent molecular alteration in CCSK (10%), are considered to be mutually exclusive. However, there is a subset of CCSKs that do not harbor either the BCOR-ITD or YWHAE-NUTM2 translocation and lack known molecular alterations. Herein, we report the first case of CCSK arising in HSK and harboring epidermal growth factor receptor ITD.

The genomic landscape of core-binding factor acute myeloid leukemias

Acute myeloid leukemia (AML) comprises a heterogeneous group of leukemias frequently defined by recurrent cytogenetic abnormalities, including rearrangements involving the core-binding factor (CBF) transcriptional complex. To better understand the genomic landscape of CBF-AMLs, we analyzed both pediatric (n = 87) and adult (n = 78) samples, including cases with RUNX1-RUNX1T1 (n = 85) or CBFB-MYH11 (n = 80) rearrangements, by whole-genome or whole-exome sequencing. In addition to known mutations in the Ras pathway, we identified recurrent stabilizing mutations in CCND2, suggesting a previously unappreciated cooperating pathway in CBF-AML. Outside of signaling alterations, RUNX1-RUNX1T1 and CBFB-MYH11 AMLs demonstrated remarkably different spectra of cooperating mutations, as RUNX1-RUNX1T1 cases harbored recurrent mutations in DHX15 and ZBTB7A, as well as an enrichment of mutations in epigenetic regulators, including ASXL2 and the cohesin complex. This detailed analysis provides insights into the pathogenesis and development of CBF-AML, while highlighting dramatic differences in the landscapes of cooperating mutations for these related AML subtypes.

Comparison of genome sequencing and clinical genotyping for pharmacogenes

We compared whole exome sequencing (WES, n = 176 patients) and whole genome sequencing (WGS, n = 68) and clinical genotyping (DMET array-based approach) for interrogating 13 genes with Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines. We focused on 127 CPIC important variants: 103 single nucleotide variations (SNV), 21 insertion/deletions (Indel), HLA-B alleles, and two CYP2D6 structural variations. WES and WGS provided interrogation of nonoverlapping sets of 115 SNV/Indels with call rate >98%. Among 68 loci interrogated by both WES and DMET, 64 loci (94.1%, confidence interval [CI]: 85.6-98.4%) showed no discrepant genotyping calls. Among 66 loci interrogated by both WGS and DMET, 63 loci (95.5%, CI: 87.2-99.0%) showed no discrepant genotyping calls. In conclusion, even without optimization to interrogate pharmacogenetic variants, WES and WGS displayed potential to provide reliable interrogation of most pharmacogenes and further validation of genome sequencing in a clinical lab setting is warranted.

Allele-Specific Quantification of Structural Variations in Cancer Genomes

Aneuploidy and structural variations (SVs) generate cancer genomes containing a mixture of rearranged genomic segments with extensive somatic copy number alterations. However, existing methods can identify either SVs or allele-specific copy number alterations, but not both simultaneously, which provides a limited view of cancer genome structure. Here we introduce Weaver, an algorithm for the quantification and analysis of allele-specific copy numbers of SVs. Weaver uses a Markov Random Field to estimate joint probabilities of allele-specific copy number of SVs and their inter-connectivity based on paired-end whole-genome sequencing data. Weaver also predicts the timing of SVs relative to chromosome amplifications. We demonstrate the accuracy of Weaver using simulations and findings from whole-genome Optical Mapping. We apply Weaver to generate allele-specific copy numbers of SVs for MCF-7 and HeLa cell lines, and identify recurrent SV patterns in 44 TCGA ovarian cancer whole-genome sequencing datasets. Our approach provides a more complete assessment of the complex genomic architectures inherent to many cancer genomes.

The clinical utility of molecular genetic cancer profiling

INTRODUCTION

Next-Generation-Sequencing (NGS) has enabled gene mutation profiling - cataloguing sequence variants and modifications in clinical assays encompassing tens to thousands of genes in tumors and in germlines. The clinical benefit of applying multi-gene NGS to diverse applications in various malignancies remains to be demonstrated.

AREAS COVERED

Applications of gene mutation profiling in oncology include screening cancer-prone families, classification of malignancies, treatment selection, and monitoring the response to treatment of solid tumors (the 'liquid biopsy'). Google Scholar was used to search PubMed for the period 2011-2016 using combinations of the following search terms: 'clinical utility', NGS, 'molecular diagnostics'. Expert commentary: Clinical studies are in progress pairing mutation profiling with streamlined new trial designs to speed identification of promising drug-target combinations and to see if genotype-informed treatment selection will improve outcome across a spectrum of histologies. The analytical advantages and falling cost of NGS make focused gene panels likely to become the dominant modality in molecular diagnostic testing even if trials eventually discourage use of large panels to test all malignancies.