SJPedPanel: A pan-cancer gene panel for childhood malignancies.

Background: Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers. Results: SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission. Conclusions: SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Malignant Progression of an Ancestral Bone Marrow Clone Harboring a CIC-NUTM2A Fusion in Isolated Myeloid Sarcoma.

Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow, which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase (RTK) signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis. IMPLICATIONS: This study illustrates molecular involvement of phenotypically normal bone marrow in myeloid sarcoma, which has significant implications in clinical care. Further, it extends the definition of CIC-rearrangements to include hematologic malignancies and shows evidence of RTK activation that may be exploited therapeutically in cancer(s) driven by these fusions.

Design and Implementation of Improved SARS-CoV-2 Diagnostic Assays To Mitigate the Impact of Genomic Mutations on Target Failure: the Xpert Xpress SARS-CoV-2 Experience.

In 2020, the U.S. Food and Drug Administration (FDA) enabled manufacturers to request emergency use authorization (EUA) to facilitate the rapid authorization of in vitro diagnostic (IVD) platforms for the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Uncommon SARS-CoV-2 point mutations could cause nucleocapsid (N) gene target failure (NGTF) when using first-generation Xpert Xpress assays, so improvements were designed and implemented. In response to NGTF reports and with consideration of viral genomic information in public databases, the Xpress assays were redesigned to mitigate the impact of SARS-CoV-2 mutations on qualitative assay performance. The second-generation assays include a third gene target (RNA-dependent RNA polymerase [RdRp]) and redundant oligonucleotide probes for the N2 target. First- and second-generation assay performances were evaluated using a challenge set of samples. A second-generation assay with updated oligonucleotide chemistry received FDA EUA in September 2021. A prototype assay with oligonucleotide chemistry similar to that of the second-generation assay with FDA EUA successfully detected all three gene targets (N2, envelope [E], and RdRp) in all challenge samples (100%; 50/50), including variants with N gene mutations (g.29197C>T or g.29200C>T), which caused NGTF in the first-generation assays. Investigation and reporting of IVD target failures, public sharing of viral genomic sequence data, and the FDA EUA pathway were essential components in facilitating a short cycle time from the identification of mutations that impact the performance of an IVD assay to the design and implementation of an improved IVD assay. IMPORTANCE The SARS-CoV-2 genome has mutated during the coronavirus disease 2019 (COVID-19) pandemic. Some of these mutations have impacted the performance of nucleic acid amplification tests like PCR, which are commonly used as diagnostic tools to detect an infection. The U.S. Food and Drug Administration (FDA) emergency use authorization (EUA) process enables the rapid reformulation and regulatory authorization of improved PCRs. In our experience, the identification of SARS-CoV-2 mutations that impact PCR performance, the subsequent development of improved PCR chemistry, and the use of the FDA EUA regulatory pathway led to improved diagnostic performance during the SARS-CoV-2 pandemic that is able to keep pace with the rapidly evolving genome of SARS-CoV-2.

Antigen cross-presentation in young tumor-bearing hosts promotes CD8+ T cell terminal differentiation.

The immune system undergoes a progressive functional remodeling with age. Understanding how age bias shapes antitumor immunity is essential in designing effective immunotherapies, especially for pediatric patients. Here, we explore antitumor CD8+ T cell responses generated in young (prepubescent) and adult (presenescent) mice. Using an MHCI-deficient tumor model, we observed that tumor-reactive CD8+ T cells expanded in young tumor-bearing (TB) mice acquired a terminally differentiated phenotype characterized by overexpression of inhibitory receptors and the transcription factor Tox1. Furthermore, tumor-infiltrating CD8+ T cells from young tumors yielded a poor cytokine response compared with CD8+ T cells infiltrating adult tumors. Young migratory dendritic cells (migDCs) from the draining lymph nodes (dLNs), and mononuclear phagocytic cells (MPCs) infiltrating young tumors, were more competent in capturing and cross-presenting tumor antigen, leading to enhanced priming of CD8+ T cells in dLNs and their subsequent terminal differentiation in the tumors. Single-cell transcriptional profiling of tumor-infiltrating MPCs demonstrated that young MPCs are polarized toward an inflammatory, effector phenotype. Consistent with our observations in young versus adult TB mice, analysis of immune infiltrates from pediatric solid tumors showed a correlation between tumor-infiltrating CD8+ T cells with an exhaustion phenotype and the frequency of PD-L1-expressing monocytes/macrophages. Collectively, these data indicate that a young tissue microenvironment contributes to the generation of an immune response skewed toward a less pliable terminal effector state, thus narrowing the window for immunotherapeutic interventions.

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 sequencing (WGS), whole-exome sequencing (WES), and RNA sequencing (RNA-seq), 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. SIGNIFICANCE: Pediatric cancers are driven by diverse genomic lesions, and sequencing has proven useful in evaluating high-risk and relapsed/refractory cases. We show that combined WGS, WES, and RNA-seq of tumor and paired normal tissues enables identification and characterization of genetic drivers across the full spectrum of pediatric cancers. This article is highlighted in the In This Issue feature, p. 2945.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data.

BACKGROUND: There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations. RESULTS: We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket. CONCLUSIONS: Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

Lipoblastomas presenting in older children and adults: analysis of 22 cases with identification of novel PLAG1 fusion partners.

Lipoblastomas are benign neoplasms of embryonal white fat that typically present in the first 3 years of life and show a lobular arrangement of maturing adipocytes with variable degrees of myxoid change. We systematically studied the clinicopathologic and genetic features of lipoblastomas arising in older children and adults. Cases with a diagnosis of lipoblastoma or maturing lipoblastoma in patients >3 years of age were retrieved from our archives. Immunostaining for CD34 and desmin and molecular studies (FISH, RNA sequencing) were performed. Twenty-two cases (8F; 14M) were identified in patients ranging from 4 to 44 years of age (median 10 years). Sites included extremity (n = 15), head and neck (n = 4), and trunk (n = 3) with tumor sizes varying from 1.6 to 17.5 cm (median 5). Only three tumors had histologic features of "conventional" lipoblastoma. The majority of tumors (n = 14) were composed of variably sized lobules of mature adipose tissue partitioned by thin fibrous septa ("maturing"). The remaining five cases consisted predominantly of bland spindled to plump ovoid cells embedded in a fibrous stroma, with a vaguely plexiform arrangement of small myxoid and adipocytic nodules ("fibroblastic"). CD34 was diffusely positive in all cases tested (21/21), while desmin immunoreactivity was identified in 12 of 21 cases (diffuse = 7, focal = 5). PLAG1 rearrangements were identified in 13 tumors in the entire cohort (59%), including all 5 fibroblastic tumors. RNA sequencing detected eight PLAG1 fusion partners, of which two were known (CHCHD7 and COL3A1) and six were novel (SRSF3, HNRNPC, PCMTD1, YWHAZ, CTDSP2, and PPP2R2A). Twelve cases had follow-up (1-107 months; median 21 months), and no recurrences were reported. Lipoblastomas may occur in older children and adults and may be difficult to recognize due to their predominantly adipocytic or fibrous appearance. Awareness that lipoblastomas may occur in older patients, careful evaluation for foci showing more typical morphologic features, ancillary immunohistochemistry for CD34 and desmin, and molecular genetic studies to identify PLAG1 rearrangements are the keys to recognizing these tumors.

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.

Germline Elongator mutations in Sonic Hedgehog medulloblastoma.

Cancer genomics has revealed many genes and core molecular processes that contribute to human malignancies, but the genetic and molecular bases of many rare cancers remains unclear. Genetic predisposition accounts for 5 to 10% of cancer diagnoses in children1,2, and genetic events that cooperate with known somatic driver events are poorly understood. Pathogenic germline variants in established cancer predisposition genes have been recently identified in 5% of patients with the malignant brain tumour medulloblastoma3. Here, by analysing all protein-coding genes, we identify and replicate rare germline loss-of-function variants across ELP1 in 14% of paediatric patients with the medulloblastoma subgroup Sonic Hedgehog (MBSHH). ELP1 was the most common medulloblastoma predisposition gene and increased the prevalence of genetic predisposition to 40% among paediatric patients with MBSHH. Parent-offspring and pedigree analyses identified two families with a history of paediatric medulloblastoma. ELP1-associated medulloblastomas were restricted to the molecular SHHα subtype4 and characterized by universal biallelic inactivation of ELP1 owing to somatic loss of chromosome arm 9q. Most ELP1-associated medulloblastomas also exhibited somatic alterations in PTCH1, which suggests that germline ELP1 loss-of-function variants predispose individuals to tumour development in combination with constitutive activation of SHH signalling. ELP1 is the largest subunit of the evolutionarily conserved Elongator complex, which catalyses translational elongation through tRNA modifications at the wobble (U34) position5,6. Tumours from patients with ELP1-associated MBSHH were characterized by a destabilized Elongator complex, loss of Elongator-dependent tRNA modifications, codon-dependent translational reprogramming, and induction of the unfolded protein response, consistent with loss of protein homeostasis due to Elongator deficiency in model systems7-9. Thus, genetic predisposition to proteome instability may be a determinant in the pathogenesis of paediatric brain cancers. These results support investigation of the role of protein homeostasis in other cancer types and potential for therapeutic interference.

Loss of glucocorticoid receptor expression mediates in vivo dexamethasone resistance in T-cell acute lymphoblastic leukemia.

Despite decades of clinical use, mechanisms of glucocorticoid resistance are poorly understood. We treated primary murine T lineage acute lymphoblastic leukemias (T-ALLs) with the glucocorticoid dexamethasone (DEX) alone and in combination with the pan-PI3 kinase inhibitor GDC-0941 and observed a robust response to DEX that was modestly enhanced by GDC-0941. Continuous in vivo treatment invariably resulted in outgrowth of drug-resistant clones, ~30% of which showed markedly reduced glucocorticoid receptor (GR) protein expression. A similar proportion of relapsed human T-ALLs also exhibited low GR protein levels. De novo or preexisting mutations in the gene encoding GR (Nr3c1) occurred in relapsed clones derived from multiple independent parental leukemias. CRISPR/Cas9 gene editing confirmed that loss of GR expression confers DEX resistance. Exposing drug-sensitive T-ALLs to DEX in vivo altered transcript levels of multiple genes, and this response was attenuated in relapsed T-ALLs. These data implicate reduced GR protein expression as a frequent cause of glucocorticoid resistance in T-ALL.

Primary bone sarcoma with BCOR internal tandem duplication.

BCOR internal tandem duplications (ITDs) and rearrangements are implicated in the oncogenesis of a subset of undifferentiated sarcomas. To date, BCOR ITD sarcomas have been exclusively found in non-appendicular infantile soft tissues, whereas BCOR-rearranged sarcomas occur in both bones and soft tissues affecting a wider patient age range. Little is known about patient outcome in BCOR ITD sarcomas. We present a BCOR-expressing, primary bone, undifferentiated sarcoma case involving an adolescent male's left tibia that, unexpectedly, harbored a BCOR ITD instead of a BCOR rearrangement. Furthermore, the patient achieved a partial histologic response after receiving a Ewing sarcoma chemotherapy regimen. Our case expands the clinical spectrum of BCOR ITD sarcomas and suggests that childhood and adult BCOR-expressing sarcomas with an undifferentiated histology should be considered for both BCOR rearrangement and ITD screening. Accurate BCOR mutation identification in undifferentiated sarcomas is essential to define their clinical spectrum and to develop effective management strategies.

Enrichment of heterozygous germline RECQL4 loss-of-function variants in pediatric osteosarcoma.

Patients harboring germline pathogenic biallelic variants in genes involved in the recognition and repair of DNA damage are known to have a substantially increased cancer risk. Emerging evidence suggests that individuals harboring heterozygous variants in these same genes may also be at heightened, albeit lesser, risk for cancer. Herein, we sought to determine whether heterozygous variants in RECQL4, the gene encoding an essential DNA helicase that is defective in children with the autosomal recessive cancer-predisposing condition Rothmund-Thomson syndrome (RTS), are associated with increased risk for childhood cancer. To address this question, we interrogated germline sequence data from 4435 pediatric cancer patients at St. Jude Children's Research Hospital and 1127 from the National Cancer Institute Therapeutically Applicable Research to Generate Effective Treatment (TARGET) database and identified 24 (0.43%) who harbored loss-of-function (LOF) RECQL4 variants, including five of 249 (2.0%) with osteosarcoma (OS). These RECQL4 variants were significantly overrepresented in children with OS, the cancer most frequently observed in patients with RTS, as compared to 134,187 noncancer controls in the Genome Aggregation Database (gnomAD v2.1; P = 0.00087, odds ratio [OR] = 7.1, 95% CI, 2.9-17). Nine of the 24 (38%) individuals possessed the same c.1573delT (p.Cys525Alafs) variant located in the highly conserved DNA helicase domain, suggesting that disruption of this domain is central to oncogenesis. Altogether these data expand our understanding of the genetic factors predisposing to childhood cancer and reveal a novel association between heterozygous RECQL4 LOF variants and development of pediatric OS.

Convergent genetic aberrations in murine and human T lineage acute lymphoblastic leukemias.

The lack of predictive preclinical models is a fundamental barrier to translating knowledge about the molecular pathogenesis of cancer into improved therapies. Insertional mutagenesis (IM) in mice is a robust strategy for generating malignancies that recapitulate the extensive inter- and intra-tumoral genetic heterogeneity found in advanced human cancers. While the central role of "driver" viral insertions in IM models that aberrantly increase the expression of proto-oncogenes or disrupt tumor suppressors has been appreciated for many years, the contributions of cooperating somatic mutations and large chromosomal alterations to tumorigenesis are largely unknown. Integrated genomic studies of T lineage acute lymphoblastic leukemias (T-ALLs) generated by IM in wild-type (WT) and Kras mutant mice reveal frequent point mutations and other recurrent non-insertional genetic alterations that also occur in human T-ALL. These somatic mutations are sensitive and specific markers for defining clonal dynamics and identifying candidate resistance mechanisms in leukemias that relapse after an initial therapeutic response. Primary cancers initiated by IM and resistant clones that emerge during in vivo treatment close key gaps in existing preclinical models, and are robust platforms for investigating the efficacy of new therapies and for elucidating how drug exposure shapes tumor evolution and patterns of resistance.

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 BRAF1,2. However, in approximately 50% of cases no genetic driver has been established2. 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 MEK3,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.

PAX5-driven subtypes of B-progenitor acute lymphoblastic leukemia.

Recent genomic studies have identified chromosomal rearrangements defining new subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL), however many cases lack a known initiating genetic alteration. Using integrated genomic analysis of 1,988 childhood and adult cases, we describe a revised taxonomy of B-ALL incorporating 23 subtypes defined by chromosomal rearrangements, sequence mutations or heterogeneous genomic alterations, many of which show marked variation in prevalence according to age. Two subtypes have frequent alterations of the B lymphoid transcription-factor gene PAX5. One, PAX5alt (7.4%), has diverse PAX5 alterations (rearrangements, intragenic amplifications or mutations); a second subtype is defined by PAX5 p.Pro80Arg and biallelic PAX5 alterations. We show that p.Pro80Arg impairs B lymphoid development and promotes the development of B-ALL with biallelic Pax5 alteration in vivo. These results demonstrate the utility of transcriptome sequencing to classify B-ALL and reinforce the central role of PAX5 as a checkpoint in B lymphoid maturation and leukemogenesis.

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.

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.

Diverse and Targetable Kinase Alterations Drive Histiocytic Neoplasms.

UNLABELLED: Histiocytic neoplasms are clonal, hematopoietic disorders characterized by an accumulation of abnormal, monocyte-derived dendritic cells or macrophages in Langerhans cell histiocytosis (LCH) and non-Langerhans cell histiocytosis (non-LCH), respectively. The discovery of BRAF(V600E) mutations in approximately 50% of these patients provided the first molecular therapeutic target in histiocytosis. However, recurrent driving mutations in the majority of patients with BRAF(V600E)-wild-type non-LCH are unknown, and recurrent cooperating mutations in non-MAP kinase pathways are undefined for the histiocytic neoplasms. Through combined whole-exome and transcriptome sequencing, we identified recurrent kinase fusions involving BRAF, ALK, and NTRK1, as well as recurrent, activating MAP2K1 and ARAF mutations in patients with BRAF(V600E)-wild-type non-LCH. In addition to MAP kinase pathway lesions, recurrently altered genes involving diverse cellular pathways were identified. Treatment of patients with MAP2K1- and ARAF-mutated non-LCH using MEK and RAF inhibitors, respectively, resulted in clinical efficacy, demonstrating the importance of detecting and targeting diverse kinase alterations in these disorders. SIGNIFICANCE: We provide the first description of kinase fusions in systemic histiocytic neoplasms and activating ARAF and MAP2K1 mutations in non-Langerhans histiocytic neoplasms. Refractory patients with MAP2K1- and ARAF-mutant histiocytoses had clinical responses to MEK inhibition and sorafenib, respectively, highlighting the importance of comprehensive genomic analysis of these disorders.