Pediatric Chordoma: A Tale of Two Genomes.

Little is known regarding the genomic alterations in chordoma, with the exception of loss of SMARCB1, a core member of the SWI/SNF complex, in poorly differentiated chordomas. A TBXT duplication and rs2305089 polymorphism, located at 6q27, are known genetic susceptibility loci. A comprehensive genomic analysis of the nuclear and mitochondrial genomes in pediatric chordoma has not yet been reported. In this study, we performed whole exome and mitochondrial DNA (mtDNA) genome sequencing on 29 chordomas from 23 pediatric patients. Findings were compared with that from whole genome sequencing datasets of 80 adult skull base chordoma patients. In the pediatric chordoma cohort, 81% percent of the somatic mtDNA mutations were observed in NADH complex genes, which is significantly enriched compared to the rest of the mtDNA genes (p=0.001). In adult chordomas, mtDNA mutations were also enriched in the NADH complex genes (p<0.0001). Furthermore, a progressive increase in heteroplasmy of non-synonymous mtDNA mutations was noted in patients with multiple tumors (p=0.0007). In the nuclear genome, rare likely germline in-frame indels in ARID1B, a member of the SWI/SNF complex located at 6q25.3, were observed in five pediatric patients (22%) and four patients in the adult cohort (5%). The frequency of rare ARID1B indels in the pediatric cohort is significantly higher than that of the adult cohort (p=0.0236, Fisher's exact test), but they were both significantly higher than that in the ethnicity-matched populations (p<5.9e-07 and p<0.0001174, respectively). Implications: germline ARID1B indels and mtDNA aberrations appear important for chordoma genesis, especially in pediatric chordoma.

Comparative Clinicopathologic and Genomic Analysis of Hepatocellular Neoplasm, Not Otherwise Specified, and Hepatoblastoma.

Accurate diagnosis and treatment of hepatocellular neoplasm, not otherwise specified (HCN-NOS), poses significant challenges. Our study aimed to investigate the clinicopathologic and genomic similarities and differences between HCN-NOS and hepatoblastoma (HB) to guide diagnostic and treatment strategies. The clinicopathologic characteristics of 16 patients with HCN-NOS and 23 patients with HB were compared. Molecular studies, including the OncoKids DNA- and RNA-based next-generation sequencing panel, chromosomal microarray, and targeted Sanger sequencing analyses of CTNNB1 and TERT promoters, were employed. We found that patients with HCN-NOS were older (P < .001) and more frequently classified as high risk (P < .01), yet they showed no significant differences in alpha fetoprotein levels or survival outcomes compared with those with HB. HCN-NOS and HB had a comparable frequency of sequence variants, with CTNNB1 mutations being predominant in both groups. Notably, TERT promoter mutations (37.5%) and rare clinically significant variants (BRAF, NRAS, and KMT2D) were exclusive to HCN-NOS. HCN-NOS demonstrated a higher prevalence of gains in 1q, encompassing the MDM4 locus (17/17 vs 11/24; P < .001), as well as loss/loss of heterozygosity (LOH) of 1p (11/17 vs 6/24; P < .05) and chromosome 11 (7/17 vs 1/24; P < .01) when compared with HB. Furthermore, the recurrent loss/LOH of chromosomes 3, 4p, 9, 15q, and Y was only observed in HCN-NOS. However, no significant differences were noted in gains of chromosomes 2, 8, and 20, or loss/LOH of 4q and 11p between the 2 groups. Notably, no clinically significant gene fusions were detected in either group. In conclusion, our study reveals that HCN-NOS exhibits high-risk clinicopathologic features and greater structural complexity compared with HB. However, patients with HCN-NOS exhibit comparable alpha fetoprotein levels at diagnosis, CTNNB1 mutation rates, and survival outcomes when subjected to aggressive treatment, as compared with those with HB. These findings have the potential to enhance diagnostic accuracy and inform more effective treatments for HCN-NOS.

An Exome Capture-Based RNA-Sequencing Assay for Genome-Wide Identification and Prioritization of Clinically Important Fusions in Pediatric Tumors.

This study reports the development of an exome capture-based RNA-sequencing assay to detect recurring and novel fusions in hematologic, solid, and central nervous system tumors. The assay used Twist Comprehensive Exome capture with either fresh or formalin-fixed samples and a bioinformatic platform that provides fusion detection, prioritization, and downstream curation. A minimum of 50 million uniquely mapped reads, a consensus read alignment/fusion calling approach using four callers (Arriba, FusionCatcher, STAR-Fusion, and Dragen), and custom software were used to integrate, annotate, and rank the candidate fusion calls. In an evaluation of 50 samples, the number of calls varied substantially by caller, from a mean of 24.8 with STAR-Fusion to 259.6 with FusionCatcher; only 1.1% of calls were made by all four callers. Therefore a filtering and ranking algorithm was developed based on multiple criteria, including number of supporting reads, calling consensus, genes involved, and cross-reference against databases of known cancer-associated or likely false-positive fusions. This approach was highly effective in pinpointing known clinically relevant fusions, ranking them first in 47 of 50 samples (94%). Detection of pathogenic gene fusions in three diagnostically challenging cases highlights the importance of a genome-wide and nontargeted method for fusion detection in pediatric cancer.

Transcriptome Sequencing Allows Comprehensive Genomic Characterization of Pediatric B-Acute Lymphoblastic Leukemia in an Academic Clinical Laboratory.

Studies have shown the power of transcriptome sequencing [RNA sequencing (RNA-Seq)] in identifying known and novel oncogenic drivers and molecular subtypes of B-acute lymphoblastic leukemia (B-ALL). The current study investigated whether the clinically validated RNA-Seq assay, coupled with a custom analysis pipeline, could be used for a comprehensive B-ALL classification. Following comprehensive clinical testing, RNA-Seq was performed on 76 retrospective B-ALL cases, 28 of which had known and 48 had undetermined subtype. Subtypes were accurately identified in all 28 known cases, and in 38 of 48 unknown cases (79%). The subtypes of the unknown cases included the following: PAX5alt (n = 12), DUX4-rearranged (n = 6), Philadelphia chromosome-like (n = 5), low hyperdiploid (n = 4), ETV6::RUNX1-like (n = 3), MEF2D-rearranged (n = 2), PAX5 P80R (n = 2), ZEB2/CEBP (n = 1), NUTM1-rearranged (n = 1), ZNF384-rearranged (n = 1), and TCF3::PBX1 (n = 1). In 15 of 38 cases (39%), classification based on expression profile was corroborated by detection of subtype-defining oncogenic drivers missed by clinical testing. RNA-Seq analysis also detected large copy number abnormalities, oncogenic hot-spot sequence variants, and intragenic IKZF1 deletions. This pilot study confirms the feasibility of implementing an RNA-Seq workflow for clinical diagnosis of molecular subtypes in pediatric B-ALL, reinforcing that RNA-Seq represents a promising global genomic assay for this heterogeneous leukemia.

Low-pass whole-genome and targeted sequencing of cell-free DNA from cerebrospinal fluid in pediatric patients with central nervous system tumors.

Background: Central nervous system tumors are the most common pediatric solid tumors and the most frequent cause of cancer-related morbidity in childhood. Significant advances in understanding the molecular features of these tumors have facilitated the development of liquid biopsy assays that may aid in diagnosis and monitoring response to therapy. In this report, we describe our comprehensive liquid biopsy platform for detection of genome-wide copy number aberrations, sequence variants, and gene fusions using cerebrospinal fluid (CSF) from pediatric patients with brain, spinal cord, and peripheral nervous system tumors. Methods: Cell-free DNA was isolated from the CSF from 55 patients, including 47 patients with tumors and 8 controls. Results: Abnormalities in cell-free DNA were detected in 24 (51%) patients including 11 with copy number alterations, 9 with sequence variants, and 7 with KIAA1549::BRAF fusions. Positive findings were obtained in patients spanning histologic subtypes, tumor grades, and anatomic locations. Conclusions: This study demonstrates the feasibility of employing this platform in routine clinical care in upfront diagnostic and monitoring settings. Future studies are required to determine the utility of this approach for assessing response to therapy and long-term surveillance.

Combined low-pass whole genome and targeted sequencing in liquid biopsies for pediatric solid tumors.

We designed a liquid biopsy (LB) platform employing low-pass whole genome sequencing (LP-WGS) and targeted sequencing of cell-free (cf) DNA from plasma to detect genome-wide copy number alterations (CNAs) and gene fusions in pediatric solid tumors. A total of 143 plasma samples were analyzed from 19 controls and 73 patients, including 44 bone or soft-tissue sarcomas and 12 renal, 10 germ cell, five hepatic, and two thyroid tumors. cfDNA was isolated from plasma collected at diagnosis, during and after therapy, and/or at relapse. Twenty-six of 37 (70%) patients enrolled at diagnosis without prior therapy (radiation, surgery, or chemotherapy) had circulating tumor DNA (ctDNA), based on the detection of CNAs from LP-WGS, including 18 of 27 (67%) patients with localized disease and eight of 10 (80%) patients with metastatic disease. None of the controls had detectable somatic CNAs. There was a high concordance of CNAs identified by LP-WGS to CNAs detected by chromosomal microarray analysis in the matching tumors. Mutations identified in tumor samples with our next-generation sequencing (NGS) panel, OncoKids®, were also detected by LP-WGS of ctDNA in 14 of 26 plasma samples. Finally, we developed a hybridization-based capture panel to target EWSR1 and FOXO1 fusions from patients with Ewing sarcoma or alveolar rhabdomyosarcoma (ARMS), respectively. Fusions were detected in the plasma from 10 of 12 patients with Ewing sarcoma and in two of two patients with ARMS. Combined, these data demonstrate the clinical applicability of our LB platform to evaluate pediatric patients with a variety of solid tumors.

Potential methylation-regulated genes and pathways in hepatocellular neoplasm, not otherwise specified.

Background and Aims: The molecular basis of hepatocellular neoplasm, not otherwise specified (HCN-NOS) is unknown. We aimed to identify gene expression patterns, potential methylation-regulated genes and pathways that characterize the tumor, and its possible relationship to hepatoblastoma and hepatocellular carcinoma (HCC). Approach & Results: Parallel genome-wide profiling of gene expression (RNAseq) and DNA methylation (EPIC850) was performed on 4 pairs of pre-treatment HCN-NOS tumors and adjacent non-tumor controls. 2530 significantly differentially expressed genes (DEGs) were identified between tumors and controls. Many of these DEGs were associated with hepatoblastoma and/or HCC. Analysis Match in Ingenuity Pathway Analysis determined that the gene expression profile of HCN-NOS was unique but significantly similar to that of both hepatoblastoma and HCC. A total of 27,195 CpG sites (CpGs) were significantly differentially methylated (DM) between tumors and controls, with a global hypomethylation pattern and predominant CpG island hypermethylation in promotor regions. Aberrant DNA methylation predominated in Developmental Process and Molecular Function Regulator pathways. Embryonic stem cell pathways were significantly enriched. In total, 1055 aberrantly methylated (at CpGs) and differentially expressed genes were identified, including 25 upstream regulators and sixty-one potential CpG island methylation-regulated genes. Eight methylation-regulated genes (TCF3, MYBL2, SRC, HMGA2, PPARGC1A, SLC22A1, COL2A1 and MYCN) had highly consistent gene expression patterns and prognostic value in patients with HCC, based on comparison to publicly available datasets. Conclusions: HCN-NOS has a unique, stem-cell like gene expression and DNA methylation profile related to both hepatoblastoma and HCC but distinct therefrom. Further, 8 methylation-regulated genes associated with prognosis in HCC were identified.

Potential of Aqueous Humor as a Liquid Biopsy for Uveal Melanoma.

Tumor biopsy can identify prognostic biomarkers for metastatic uveal melanoma (UM), however aqueous humor (AH) liquid biopsy may serve as an adjunct. This study investigated whether the AH of UM eyes has sufficient circulating tumor DNA (ctDNA) to perform genetic analysis. This is a case series of 37 AH samples, taken before or after radiation, and one tumor wash sample, from 12 choroidal and 8 ciliary body (CB) melanoma eyes. AH was analyzed for nucleic acid concentrations. AH DNA and one tumor wash sample underwent shallow whole-genome sequencing followed by Illumina sequencing to detect somatic copy number alterations (SCNAs). Four post-radiation AH underwent targeted sequencing of BAP1 and GNAQ genes. Post-radiation AH had significantly higher DNA and miRNA concentrations than paired pre-radiation samples. Highly recurrent UM SCNAs were identified in 0/11 post-radiation choroidal and 6/8 post-radiation CB AH. SCNAs were highly concordant in a CB post-radiation AH with its matched tumor (r = 0.978). BAP1 or GNAQ variants were detected in 3/4 post-radiation AH samples. AH is a source of ctDNA in UM eyes, particularly in post-radiation CB eyes. For the first time, UM SCNAs and mutations were identified in AH-derived ctDNA. Suggesting that AH can serve as a liquid biopsy for UM.

Non-Invasive Profiling of Advanced Prostate Cancer via Multi-Parametric Liquid Biopsy and Radiomic Analysis.

Integrating liquid biopsies of circulating tumor cells (CTCs) and cell-free DNA (cfDNA) with other minimally invasive measures may yield more comprehensive disease profiles. We evaluated the feasibility of concurrent cellular and molecular analysis of CTCs and cfDNA combined with radiomic analysis of CT scans from patients with metastatic castration-resistant PC (mCRPC). CTCs from 22 patients were enumerated, stained for PC-relevant markers, and clustered based on morphometric and immunofluorescent features using machine learning. DNA from single CTCs, matched cfDNA, and buffy coats was sequenced using a targeted amplicon cancer hotspot panel. Radiomic analysis was performed on bone metastases identified on CT scans from the same patients. CTCs were detected in 77% of patients and clustered reproducibly. cfDNA sequencing had high sensitivity (98.8%) for germline variants compared to WBC. Shared and unique somatic variants in PC-related genes were detected in cfDNA in 45% of patients (MAF > 0.1%) and in CTCs in 92% of patients (MAF > 10%). Radiomic analysis identified a signature that strongly correlated with CTC count and plasma cfDNA level. Integration of cellular, molecular, and radiomic data in a multi-parametric approach is feasible, yielding complementary profiles that may enable more comprehensive non-invasive disease modeling and prediction.

The spectrum of mitochondrial DNA (mtDNA) mutations in pediatric CNS tumors.

BACKGROUND: We previously established the landscape of mitochondrial DNA (mtDNA) mutations in 23 subtypes of pediatric malignancies, characterized mtDNA mutation profiles among these subtypes, and provided statistically significant evidence for a contributory role of mtDNA mutations to pediatric malignancies. METHODS: To further delineate the spectrum of mtDNA mutations in pediatric central nervous system (CNS) tumors, we analyzed 545 tumor-normal paired whole-genome sequencing datasets from the Children's Brain Tumor Tissue Consortium. RESULTS: Germline mtDNA variants were used to determine the haplogroup, and maternal ancestry, which was not significantly different among tumor types. Among 166 (30.5%) tumors we detected 220 somatic mtDNA mutations, primarily missense mutations (36.8%), as well as 22 loss-of-function mutations. Different pediatric CNS tumor subtypes had distinct mtDNA mutation profiles. The number of mtDNA mutations per tumor ranged from 0.20 (dysembryoplastic neuroepithelial tumor [DNET]) to 0.75 (meningiomas). The average heteroplasmy was 10.7%, ranging from 4.6% in atypical teratoid/rhabdoid tumor (AT/RT) to 26% in diffuse intrinsic pontine glioma. High-grade gliomas had a significant higher number of mtDNA mutations per sample than low-grade gliomas (0.6 vs 0.27) (P = .004), with almost twice as many missense mtDNA mutations per sample (0.24 vs 0.11), and higher average heteroplasmy levels (16% vs 10%). Recurrent mtDNA mutations may represent hotspots which may serve as biologic markers of disease. CONCLUSIONS: Our findings demonstrate varying contributions of mtDNA mutations in different subtypes of CNS tumors. Sequencing the mtDNA genome may ultimately be used to characterize CNS tumors at diagnosis and monitor disease progression.

Increased viral variants in children and young adults with impaired humoral immunity and persistent SARS-CoV-2 infection: A consecutive case series.

BACKGROUND: There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses. METHODS: We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape. FINDINGS: We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape. INTERPRETATION: Our results highlight the potential need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered.

Persistent SARS-CoV-2 infection and increasing viral variants in children and young adults with impaired humoral immunity.

Background: There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses. Methods: We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape. Findings: We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape. Interpretation: Our results highlight the need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered. Funding: The work was partially funded by The Saban Research Institute at Children's Hospital Los Angeles intramural support for COVID-19 Directed Research (X.G. and J.D.B.), the Johns Hopkins Center of Excellence in Influenza Research and Surveillance HHSN272201400007C (A.P.), NIH/NIAID R01AI127877 (S.D.B.), NIH/NIAID R01AI130398 (S.D.B.), NIH 1U54CA260517 (S.D.B.), an endowment to S.D.B. from the Crown Family Foundation, an Early Postdoc.Mobility Fellowship Stipend to O.F.W. from the Swiss National Science Foundation (SNSF), and a Coulter COVID-19 Rapid Response Award to S.D.B. L.G. is a SHARE Research Fellow in Pediatric Hematology-Oncology.

Detection of mitochondrial DNA variants at low level heteroplasmy in pediatric CNS and extra-CNS solid tumors with three different enrichment methods.

The mitochondrial genome is small, 16.5 kb, and yet complex to study due to an abundance of mitochondria in any given cell or tissue. Mitochondrial DNA (mtDNA) mutations have been previously described in cancer, many of which were detected at low heteroplasmy. In this study we enriched the mitochondrial genome in primary pediatric tumors for detection of mtDNA variants. We completed mtDNA enrichment using REPLI-g, Agilent SureSelect, and long-range polymerase chain reaction (LRPCR) followed by next generation sequencing (NGS) on Illumina platforms. Primary tumor and germline genomic DNA from a variety of pediatric central nervous system (CNS) and extra-CNS solid tumors were analyzed by the three different methods. Although all three methods performed equally well for detecting variants at high heteroplasmy or homoplasmy, only LRPCR and SureSelect-based enrichment methods provided consistent results for variants that were present at less than five percent heteroplasmy. We then applied both LRPCR and SureSelect to three successive samples from a patient with multiply-recurrent gliofibroma and detected a low-level novel mutation as well as a change in heteroplasmy levels of a synonymous variant that was correlated with progression of disease. IMPLICATION: This study demonstrates that LRPCR and SureSelect enrichment, but not REPLI-g, followed by NGS are accurate methods for studying the mtDNA variations at low heteroplasmy, which may be applied to studying mtDNA mutations in cancer.

OncoKids: A Comprehensive Next-Generation Sequencing Panel for Pediatric Malignancies.

The OncoKids panel is an amplification-based next-generation sequencing assay designed to detect diagnostic, prognostic, and therapeutic markers across the spectrum of pediatric malignancies, including leukemias, sarcomas, brain tumors, and embryonal tumors. This panel uses low input amounts of DNA (20 ng) and RNA (20 ng) and is compatible with formalin-fixed, paraffin-embedded and frozen tissue, bone marrow, and peripheral blood. The DNA content of this panel covers the full coding regions of 44 cancer predisposition loci, tumor suppressor genes, and oncogenes; hotspots for mutations in 82 genes; and amplification events in 24 genes. The RNA content includes 1421 targeted gene fusions. We describe the validation of this panel by using a large cohort of 192 unique clinical samples that included a wide range of tumor types and alterations. Robust performance was observed for analytical sensitivity, reproducibility, and limit of detection studies. The results from this study support the use of OncoKids for routine clinical testing of a wide variety of pediatric malignancies.

A neoplasm with FIP1L1-PDGFRA fusion presenting as pediatric T-cell lymphoblastic leukemia/lymphoma without eosinophilia.

The 2016 World Health Organization (2016 WHO) classification of hematopoietic malignancies classifies neoplasms with a fusion between the FIP1L1 and PDGFRA genes in 4q12 into a group called "myeloid and lymphoid neoplasms with eosinophilia and abnormalities of PDGFRA, PDGFRB or FGFR1 or with PCM1-JAK2". Neoplasms characterized by this fusion are pluripotent stem cell disorders that can show both myeloid and lymphoid differentiation. They typically occur in adult patients and most are characterized by eosinophilia. We describe identification of a FIP1L1-PDGFRA fusion in a 13-year-old boy who presented with T-lymphoblastic leukemia/lymphoma without eosinophilia. Detection of FIP1L1-PDGFRA driven neoplasms at diagnosis is usually critical for proper treatment, since almost all reported cases responded to tyrosine kinase inhibitors. However, our patient's leukemia was refractory to standard chemotherapy, and did not show a meaningful response to tyrosine kinase inhibitor therapy. Testing for a FIP1L1-PDGFRA rearrangement is at present limited to patients with idiopathic hypereosinophilia, and we hypothesize that this abnormality may be under-diagnosed in children with acute leukemias.

A Rapid and Sensitive Next-Generation Sequencing Method to Detect RB1 Mutations Improves Care for Retinoblastoma Patients and Their Families.

Retinoblastoma is a childhood eye malignancy that can lead to the loss of vision, eye(s), and sometimes life. The tumors are initiated by inactivating mutations in both alleles of the tumor-suppressor gene, RB1, or, rarely, by MYCN amplification. Timely identification of a germline RB1 mutation in blood samples or either somatic RB1 mutation or MYCN amplification in tumors is important for effective care and management of retinoblastoma patients and their families. However, current procedures to thoroughly test RB1 mutations are complicated and lengthy. Herein, we report a next-generation sequencing-based method capable of detecting point mutations, small indels, and large deletions or duplications across the entire RB1 gene and amplification of MYCN gene on a single platform. From DNA extraction to clinical interpretation requires only 3 days, enabling early molecular diagnosis of retinoblastoma and optimal treatment outcomes. This method can also detect low-level mosaic mutations in blood samples that can be missed by routine Sanger sequencing. In addition, it can differentiate between RB1 mutation- and MYCN amplification-driven retinoblastomas. This rapid, comprehensive, and sensitive method for detecting RB1 mutations and MYCN amplification can readily identify RB1 mutation carriers and thus improve the management and genetic counseling for retinoblastoma patients and their families.