Azacitidine as epigenetic priming for chemotherapy is safe and well-tolerated in infants with newly diagnosed KMT2A-rearranged acute lymphoblastic leukemia: Children's Oncology Group trial AALL15P1.

Infants less than 1 year old diagnosed with KMT2A-rearranged (KMT2A-r) acute lymphoblastic leukemia (ALL) are at high risk of remission failure, relapse, and death due to leukemia, despite intensive therapies. Infant KMT2A-r ALL blasts are characterized by DNA hypermethylation. Epigenetic priming with DNA methyltransferase inhibitors increases the cytotoxicity of chemotherapy in preclinical studies. The Children's Oncology Group trial AALL15P1 tested the safety and tolerability of five days of azacitidine immediately prior to the start of chemotherapy on day six, in four post-induction chemotherapy courses for infants with newly diagnosed KMT2A-r ALL. The treatment was welltolerated, with only two of 31 evaluable patients (6.5%) experiencing dose-limiting toxicity. Whole genome bisulfite sequencing of peripheral blood mononuclear cells (PBMCs) demonstrated decreased DNA methylation in 87% of samples tested following five days of azacitidine. Event-free survival was similar to prior studies of newly diagnosed infant ALL. Azacitidine is safe and results in decreased DNA methylation of PBMCs in infants with KMT2A-r ALL, but the incorporation of azacitidine to enhance cytotoxicity did not impact survival. Clinicaltrials.gov identifier: NCT02828358.

Genetic heterogeneity and enrichment of variants in DNA-repair genes in ameloblastoma.

OBJECTIVE: Ameloblastomas are a group of relatively common odontogenic tumors that frequently originate from the dental epithelium. These tumors are aggressive in nature and present as slow-growing painless cortical expansion of the jaw. Histologically, the follicular and plexiform subtypes constitute two-thirds of solid/multicystic ameloblastomas. The objective of this study was to understand the genetic architecture of follicular and plexiform ameloblastomas using deep whole-exome sequencing. METHODS: Archived formalin-fixed paraffin-embedded tissue blocks of follicular (n = 4) and plexiform (n = 6) ameloblastomas were retrieved and genomic DNAs were isolated from the tumor tissue dissected from the formalin-fixed paraffin-embedded block. The exomes were enriched using the Integrated DNA Technologies Exome Research Panel (IDT, Coralville, IA) and paired-end sequencing was completed on an Illumina NovaSeq 6000 with an average output of 20 GB of data resulting in a mean coverage of 400×. Variant analysis was completed using custom-developed software: Rapid Understanding of Nucleotide variant Effect Software and variant integration and knowledge interpretation in genomes. RESULTS: Our analyses focused on examining somatic variants (gnomAD minor allele frequency ≤1%) in genes found on an Food and Drug Administration -approved clinical cancer sequencing panel (FoundationOne®CDx). In follicular tumors, variants (>20% of the reads) were identified in BRAF, KMT2D, and ABL1 genes. In plexiform tumors, variants (>20% of the reads) were identified in ALK, BRAF, KRAS, KMT2D, SMO, KMT2A, and BRCA2 genes. Enrichment analysis showed a significant role of DNA repair genes in the development of these tumors. CONCLUSION: The variants identified in follicular and plexiform ameloblastomas were enriched in DNA-repair genes. The observed genetic heterogeneity in these ameloblastomas may contribute to the aggressive nature and recurrence risk of these tumors.

Clinical Validation of Genome Reference Consortium Human Build 38 in a Laboratory Utilizing Next-Generation Sequencing Technologies.

BACKGROUND: Laboratories utilizing next-generation sequencing align sequence data to a standardized human reference genome (HRG). Several updated versions, or builds, have been released since the original HRG in 2001, including the Genome Reference Consortium Human Build 38 (GRCh38) in 2013. However, most clinical laboratories still use GRCh37, which was released in 2009. We report our laboratory's clinical validation of GRCh38. METHODS: Migration to GRCh38 was validated by comparing the coordinates (lifting over) of 9443 internally curated variants from GRCh37 to GRCh38, globally comparing protein coding sequence variants aligned with GRCh37 vs GRCh38 from 917 exomes, assessing genes with known discrepancies, comparing coverage differences, and establishing the analytic sensitivity and specificity of variant detection using Genome in a Bottle data. RESULTS: Eight discrepancies, due to strand swap or reference base, were observed. Three clinically relevant variants had the GRCh37 alternate allele as the reference allele in GRCh38. A comparison of 88 295 calls between builds identified 8 disease-associated genes with sequence differences: ABO, BNC2, KIZ, NEFL, NR2E3, PTPRQ, SHANK2, and SRD5A2. Discrepancies in coding regions in GRCh37 were resolved in GRCh38. CONCLUSIONS: There were a small number of clinically significant changes between the 2 genome builds. GRCh38 provided improved detection of nucleotide changes due to the resolution of discrepancies present in GRCh37. Implementation of GRCh38 results in more accurate and consistent reporting.

Altered subcellular localization of transcription factor TEAD4 regulates first mammalian cell lineage commitment.

In the preimplantation mouse embryo, TEAD4 is critical to establishing the trophectoderm (TE)-specific transcriptional program and segregating TE from the inner cell mass (ICM). However, TEAD4 is expressed in the TE and the ICM. Thus, differential function of TEAD4 rather than expression itself regulates specification of the first two cell lineages. We used ChIP sequencing to define genomewide TEAD4 target genes and asked how transcription of TEAD4 target genes is specifically maintained in the TE. Our analyses revealed an evolutionarily conserved mechanism, in which lack of nuclear localization of TEAD4 impairs the TE-specific transcriptional program in inner blastomeres, thereby allowing their maturation toward the ICM lineage. Restoration of TEAD4 nuclear localization maintains the TE-specific transcriptional program in the inner blastomeres and prevents segregation of the TE and ICM lineages and blastocyst formation. We propose that altered subcellular localization of TEAD4 in blastomeres dictates first mammalian cell fate specification.

Hepatocyte nuclear factor 4 alpha deletion promotes diethylnitrosamine-induced hepatocellular carcinoma in rodents.

Hepatocyte nuclear factor 4 alpha (HNF4α), the master regulator of hepatocyte differentiation, has been recently shown to inhibit hepatocyte proliferation by way of unknown mechanisms. We investigated the mechanisms of HNF4α-induced inhibition of hepatocyte proliferation using a novel tamoxifen (TAM)-inducible, hepatocyte-specific HNF4α knockdown mouse model. Hepatocyte-specific deletion of HNF4α in adult mice resulted in increased hepatocyte proliferation, with a significant increase in liver-to-body-weight ratio. We determined global gene expression changes using Illumina HiSeq-based RNA sequencing, which revealed that a significant number of up-regulated genes following deletion of HNF4α were associated with cancer pathogenesis, cell cycle control, and cell proliferation. The pathway analysis further revealed that c-Myc-regulated gene expression network was highly activated following HNF4α deletion. To determine whether deletion of HNF4α affects cancer pathogenesis, HNF4α knockdown was induced in mice treated with the known hepatic carcinogen diethylnitrosamine (DEN). Deletion of HNF4α significantly increased the number and size of DEN-induced hepatic tumors. Pathological analysis revealed that tumors in HNF4α-deleted mice were well-differentiated hepatocellular carcinoma (HCC) and mixed HCC-cholangiocarcinoma. Analysis of tumors and surrounding normal liver tissue in DEN-treated HNF4α knockout mice showed significant induction in c-Myc expression. Taken together, deletion of HNF4α in adult hepatocytes results in increased hepatocyte proliferation and promotion of DEN-induced hepatic tumors secondary to aberrant c-Myc activation.

Evaluation of Hi-C sequencing for the detection of gene fusions in hematologic and solid pediatric cancer samples.

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate HiC as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens previously positive for clinically significant genomic rearrangements. Archived specimen types tested included viable and nonviable frozen leukemic cells, as well as formalin-fixed paraffin-embedded (FFPE) tumor tissues. Initially, pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to HiC analysis to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases with no known genomic rearrangements based on prior clinical diagnostic testing were evaluated to determine whether HiC could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, 100% concordance was observed between HiC and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study demonstrates the value of HiC sequencing to medical diagnostic testing as it identified several clinically significant rearrangements, including those that might have been missed by current clinical testing workflows. Key points: HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome, facilitating detection of genomic rearrangements.HiC was 100% concordant with clinical diagnostic testing workflows for detecting clinically significant genomic rearrangements in pediatric leukemia and rhabdomyosarcoma specimens.HiC detected clinically significant genomic rearrangements not previously detected by prior clinical cytogenetic and molecular testing.HiC performed well with archived non-viable and viable frozen leukemic cell samples, as well as archived formalin-fixed paraffin-embedded tumor tissue specimens.

Evaluation of Hi-C Sequencing for Detection of Gene Fusions in Hematologic and Solid Tumor Pediatric Cancer Samples.

Hi-C sequencing is a DNA-based next-generation sequencing method that preserves the 3D genome conformation and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate Hi-C as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens. Archived viable and non-viable frozen leukemic cells and formalin-fixed paraffin-embedded (FFPE) tumor specimens were analyzed. Pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to Hi-C to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases without known genomic rearrangements based on prior clinical diagnostic testing was evaluated to determine whether Hi-C could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, we observed 100% concordance between Hi-C and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study provides an institutional proof of principle evaluation of Hi-C sequencing to medical diagnostic testing as it identified several clinically relevant rearrangements, including those that were missed by current clinical testing workflows.

DeepSomatic: Accurate somatic small variant discovery for multiple sequencing technologies.

Somatic variant detection is an integral part of cancer genomics analysis. While most methods have focused on short-read sequencing, long-read technologies now offer potential advantages in terms of repeat mapping and variant phasing. We present DeepSomatic, a deep learning method for detecting somatic SNVs and insertions and deletions (indels) from both short-read and long-read data, with modes for whole-genome and exome sequencing, and able to run on tumor-normal, tumor-only, and with FFPE-prepared samples. To help address the dearth of publicly available training and benchmarking data for somatic variant detection, we generated and make openly available a dataset of five matched tumor-normal cell line pairs sequenced with Illumina, PacBio HiFi, and Oxford Nanopore Technologies, along with benchmark variant sets. Across samples and technologies (short-read and long-read), DeepSomatic consistently outperforms existing callers, particularly for indels.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads.

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

RNA-sequencing quantification of hepatic ontogeny and tissue distribution of mRNAs of phase II enzymes in mice.

Phase II conjugating enzymes play key roles in the metabolism of xenobiotics. In the present study, RNA sequencing was used to elucidate hepatic ontogeny and tissue distribution of mRNA expression of all major known Phase II enzymes, including enzymes involved in glucuronidation, sulfation, glutathione conjugation, acetylation, methylation, and amino acid conjugation, as well as enzymes for the synthesis of Phase II cosubstrates, in male C57BL/6J mice. Livers from male C57BL/6J mice were collected at 12 ages from prenatal to adulthood. Many of these Phase II enzymes were expressed at much higher levels in adult livers than in perinatal livers, such as Ugt1a6b, -2a3, -2b1, -2b5, -2b36, -3a1, and -3a2; Gsta1, -m1, -p1, -p2, and -z1; mGst1; Nat8; Comt; Nnmt; Baat; Ugdh; and Gclc. In contrast, hepatic mRNA expression of a few Phase II enzymes decreased during postnatal liver development, such as mGst2, mGst3, Gclm, and Mat2a. Hepatic expression of certain Phase II enzymes peaked during the adolescent stage, such as Ugt1a1, Sult1a1, Sult1c2, Sult1d1, Sult2as, Sult5a1, Tpmt, Glyat, Ugp2, and Mat1a. In adult mice, the total transcripts for Phase II enzymes were comparable in liver, kidney, and small intestine; however, individual Phase II enzymes displayed marked tissue specificity among the three organs. In conclusion, this study unveils for the first time developmental changes in mRNA abundance of all major known Phase II enzymes in mouse liver, as well as their tissue-specific expression in key drug-metabolizing organs. The age- and tissue-specific expression of Phase II enzymes indicate that the detoxification of xenobiotics is highly regulated by age and cell type.

RNA-sequencing quantification of hepatic ontogeny of phase-I enzymes in mice.

Phase-I drug metabolizing enzymes catalyze reactions of hydrolysis, reduction, and oxidation of drugs and play a critical role in drug metabolism. However, the functions of most phase-I enzymes are not mature at birth, which markedly affects drug metabolism in newborns. Therefore, characterization of the expression profiles of phase-I enzymes and the underlying regulatory mechanisms during liver maturation is needed for better estimation of using drugs in pediatric patients. The mouse is an animal model widely used for studying the mechanisms in the regulation of developmental expression of phase-I genes. Therefore, we applied RNA sequencing to provide a "true quantification" of the mRNA expression of phase-I genes in the mouse liver during development. Liver samples of male C57BL/6 mice at 12 different ages from prenatal to adulthood were used for defining the ontogenic mRNA profiles of phase-I families, including hydrolysis: carboxylesterase (Ces), paraoxonase (Pon), and epoxide hydrolase (Ephx); reduction: aldo-keto reductase (Akr), quinone oxidoreductase (Nqo), and dihydropyrimidine dehydrogenase (Dpyd); and oxidation: alcohol dehydrogenase (Adh), aldehyde dehydrogenase (Aldh), flavin monooxygenases (Fmo), molybdenum hydroxylase (Aox and Xdh), cytochrome P450 (P450), and cytochrome P450 oxidoreductase (Por). Two rapidly increasing stages of total phase-I gene expression after birth reflect functional transition of the liver during development. Diverse expression patterns were identified, and some large gene families contained the mRNA of genes that are enriched at different stages of development. Our study reveals the mRNA abundance of phase-I genes in the mouse liver during development and provides a valuable foundation for mechanistic studies in the future.

Integrated genetic profiling of archival pediatric high-grade glial tumors and reassessment with 2021 WHO classification of paediatric CNS tumours.

Though rare, pediatric high-grade gliomas (pHGG) are a leading cause of cancer-related mortality in children. We wanted to determine whether our currently available clinical laboratory methods could better define diagnosis for pHGG that had been archived at our institution for the past 20 years (1998 to 2017). We investigated 33 formalin-fixed paraffin-embedded pHGG using ThermoFisher Oncoscan SNP microarray with somatic mutation analysis, Sanger sequencing, and whole genome sequencing. These data were correlated with historical histopathological, chromosomal, clinical, and radiological data. Tumors were subsequently classified according to the 2021 WHO Classification of Paediatric CNS Tumours. All 33 tumors were found to have genetic aberrations that placed them within a 2021 WHO subtype and/or provided prognostic information; 6 tumors were upgraded from WHO CNS grade 3 to grade 4. New pHGG genetic features were found including two small cell glioblastomas with H3 G34 mutations not previously described; one tumor with STRN-NTRK2 fusion; and a congenital diffuse leptomeningeal glioneuronal tumor without a chromosomal 1p deletion but with KIAA1549-BRAF fusion. Overall, the combination of laboratory methods yielded key information for tumor classification. Thus, even small studies of these uncommon tumor types may yield new genetic features and possible new subtypes that warrant future investigations.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

RNA sequencing reveals dynamic changes of mRNA abundance of cytochromes P450 and their alternative transcripts during mouse liver development.

Cytochromes P450 (P450s) are a superfamily of enzymes that have critical functions in liver to catalyze the biotransformation of numerous drugs. However, the functions of most P450s are not mature at birth, which can markedly affect the metabolism of drugs in newborns. Therefore, characterization of the developmental profiles and regulatory mechanisms of P450 expression is needed for more rational drug therapy of pediatric patients. An animal model is indispensable for studying the mechanisms of postnatal development of the P450s. Hence we used RNA sequencing (RNA-Seq) to provide a "true quantification" of mRNA expression of all P450s in mouse liver during development. Liver samples of male C57BL/6 mice at 12 different ages from prenatal to adulthood were used. Total mRNAs of the 103 mouse P450s displayed two rapid increasing stages after birth, reflecting critical functional transition of liver during development. Four ontogenic expression patterns were identified among the 71 significantly expressed P450s, which categorized genes into neonatal-, adolescent-, adolescent/adult-, and adult-enriched groups. The 10 most highly expressed subfamilies of mouse P450s in livers of adult mice were CYP2E, -2C, -2D, -3A, -4A, -2F, -2A, -1A, -4F, and -2B, which showed diverse expression profiles during development. The expression patterns of multiple members within a P450 subfamily were often classified to different groups. RNA-Seq also enabled the quantification of known transcript variants of CYP2C44, CYP2C50, CYP2D22, CYP3A25, and CYP26B1 and identification of novel transcripts for CYP2B10, CYP2D26, and CYP3A13. In conclusion, this study reveals the mRNA abundance of all the P450s in mouse liver during development and provides a foundation for mechanistic studies in the future.

Benchmarking challenging small variants with linked and long reads.

Genome in a Bottle benchmarks are widely used to help validate clinical sequencing pipelines and develop variant calling and sequencing methods. Here we use accurate linked and long reads to expand benchmarks in 7 samples to include difficult-to-map regions and segmental duplications that are challenging for short reads. These benchmarks add more than 300,000 SNVs and 50,000 insertions or deletions (indels) and include 16% more exonic variants, many in challenging, clinically relevant genes not covered previously, such as PMS2. For HG002, we include 92% of the autosomal GRCh38 assembly while excluding regions problematic for benchmarking small variants, such as copy number variants, that should not have been in the previous version, which included 85% of GRCh38. It identifies eight times more false negatives in a short read variant call set relative to our previous benchmark. We demonstrate that this benchmark reliably identifies false positives and false negatives across technologies, enabling ongoing methods development.

Curated variation benchmarks for challenging medically relevant autosomal genes.

The repetitive nature and complexity of some medically relevant genes poses a challenge for their accurate analysis in a clinical setting. The Genome in a Bottle Consortium has provided variant benchmark sets, but these exclude nearly 400 medically relevant genes due to their repetitiveness or polymorphic complexity. Here, we characterize 273 of these 395 challenging autosomal genes using a haplotype-resolved whole-genome assembly. This curated benchmark reports over 17,000 single-nucleotide variations, 3,600 insertions and deletions and 200 structural variations each for human genome reference GRCh37 and GRCh38 across HG002. We show that false duplications in either GRCh37 or GRCh38 result in reference-specific, missed variants for short- and long-read technologies in medically relevant genes, including CBS, CRYAA and KCNE1. When masking these false duplications, variant recall can improve from 8% to 100%. Forming benchmarks from a haplotype-resolved whole-genome assembly may become a prototype for future benchmarks covering the whole genome.

Factors Affecting Migration to GRCh38 in Laboratories Performing Clinical Next-Generation Sequencing.

The most recent build of the human reference genome, GRCh38, was released in 2013. However, many laboratories performing next-generation sequencing (NGS) continue to align to GRCh37. Our aim was to assess the number of clinical diagnostic laboratories that have migrated to GRCh38 and discern factors impeding migration for those still using GRCh37. A brief, five-question survey was electronically administered to 71 clinical laboratories offering constitutional NGS-based testing and analyzed categorically. Twenty-eight responses meeting inclusion criteria were collected from 24 academic and four commercial diagnostic laboratories. Most of these (14; 50%) reported volumes of <500 NGS-based tests in 2019. Only two respondents (7%) had already migrated entirely to GRCh38; most laboratories (15; 54%) had no plans to migrate. The two prevailing reasons for not yet migrating were as follows: laboratories did not feel the benefits outweighed the time and monetary costs (14; 50%); and laboratories had insufficient staff to facilitate the migration (12; 43%). These data, although limited, suggest most clinical molecular laboratories are reluctant to migrate to GRCh38, and there appear to be multiple obstacles to overcome before GRCh38 is widely adopted.

A diploid assembly-based benchmark for variants in the major histocompatibility complex.

Most human genomes are characterized by aligning individual reads to the reference genome, but accurate long reads and linked reads now enable us to construct accurate, phased de novo assemblies. We focus on a medically important, highly variable, 5 million base-pair (bp) region where diploid assembly is particularly useful - the Major Histocompatibility Complex (MHC). Here, we develop a human genome benchmark derived from a diploid assembly for the openly-consented Genome in a Bottle sample HG002. We assemble a single contig for each haplotype, align them to the reference, call phased small and structural variants, and define a small variant benchmark for the MHC, covering 94% of the MHC and 22368 variants smaller than 50 bp, 49% more variants than a mapping-based benchmark. This benchmark reliably identifies errors in mapping-based callsets, and enables performance assessment in regions with much denser, complex variation than regions covered by previous benchmarks.