Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Somatic SLC35A2 mosaicism correlates with clinical findings in epilepsy brain tissue.

OBJECTIVE: Many genetic studies of intractable epilepsy in pediatric patients primarily focus on inherited, constitutional genetic deficiencies identified in patient blood. Recently, studies have revealed somatic mosaicism associated with epilepsy in which genetic variants are present only in a subset of brain cells. We hypothesize that tissue-specific, somatic mosaicism represents an important genetic etiology in epilepsy and aim to discover somatic alterations in epilepsy-affected brain tissue. METHODS: We have pursued a research study to identify brain somatic mosaicism, using next-generation sequencing (NGS) technologies, in patients with treatment refractory epilepsy who have undergone surgical resection of affected brain tissue. RESULTS: We used an integrated combination of NGS techniques and conventional approaches (radiology, histopathology, and electrophysiology) to comprehensively characterize multiple brain regions from a single patient with intractable epilepsy. We present a 3-year-old male patient with West syndrome and intractable tonic seizures in whom we identified a pathogenic frameshift somatic variant in SLC35A2, present at a range of variant allele fractions (4.2%-19.5%) in 12 different brain tissues detected by targeted sequencing. The proportion of the SLC35A2 variant correlated with severity and location of neurophysiology and neuroimaging abnormalities for each tissue. CONCLUSIONS: Our findings support the importance of tissue-based sequencing and highlight a correlation in our patient between SLC35A2 variant allele fractions and the severity of epileptogenic phenotypes in different brain tissues obtained from a grid-based resection of clinically defined epileptogenic regions.

Samovar: Single-Sample Mosaic Single-Nucleotide Variant Calling with Linked Reads.

Linked-read sequencing enables greatly improves haplotype assembly over standard paired-end analysis. The detection of mosaic single-nucleotide variants benefits from haplotype assembly when the model is informed by the mapping between constituent reads and linked reads. Samovar evaluates haplotype-discordant reads identified through linked-read sequencing, thus enabling phasing and mosaic variant detection across the entire genome. Samovar trains a random forest model to score candidate sites using a dataset that considers read quality, phasing, and linked-read characteristics. Samovar calls mosaic single-nucleotide variants (SNVs) within a single sample with accuracy comparable with what previously required trios or matched tumor/normal pairs and outperforms single-sample mosaic variant callers at minor allele frequency 5%-50% with at least 30X coverage. Samovar finds somatic variants in both tumor and normal whole-genome sequencing from 13 pediatric cancer cases that can be corroborated with high recall with whole exome sequencing. Samovar is available open-source at https://github.com/cdarby/samovar under the MIT license.

Analysis of hedgehog signaling in cerebellar granule cell precursors in a conditional Nsdhl allele demonstrates an essential role for cholesterol in postnatal CNS development.

NSDHL is a 3ß-hydroxysterol dehydrogenase that is involved in the removal of two C-4 methyl groups in one of the later steps of cholesterol biosynthesis. Mutations in the gene encoding the enzyme are responsible for the X-linked, male lethal mouse mutations bare patches and striated, as well as most cases of human CHILD syndrome. Rare, hypomorphic NSDHL mutations are also associated with X-linked intellectual disability in males with CK syndrome. Since hemizygous male mice with Nsdhl mutations die by midgestation, we generated a conditional targeted Nsdhl mutation (Nsdhl(tm1.1Hrm)) to investigate the essential role of cholesterol in the early postnatal CNS. Ablation of Nsdhl in radial glia using GFAP-cre resulted in live-born, normal appearing affected male pups. However, the pups develop overt ataxia by postnatal day 8-10 and die shortly thereafter. Histological abnormalities include progressive loss of cortical and hippocampal neurons, as well as deficits in the proliferation and migration of cerebellar granule precursors and subsequent massive apoptosis of the cerebellar cortex. We replicated the granule cell precursor proliferation defect in vitro and demonstrate that it results from defective signaling by SHH. Furthermore, this defect is almost completely rescued by supplementation of the culture media with exogenous cholesterol, while methylsterol accumulation above the enzymatic block appears to be associated with increased cell death. These data support the absolute requirement for cholesterol synthesis in situ once the blood-brain-barrier forms and cholesterol transport to the fetus is abolished. They further emphasize the complex ramifications of cholesterogenic enzyme deficiency on cellular metabolism.

Contactin 4 as an autism susceptibility locus.

Structural and sequence variation have been described in several members of the contactin (CNTN) and contactin-associated protein (CNTNAP) gene families in association with neurodevelopmental disorders, including autism. Using array comparative genome hybridization (CGH), we identified a maternally inherited ∼535 kb deletion at 3p26.3 encompassing the 5' end of the contactin 4 gene (CNTN4) in a patient with autism. Based on this finding and previous reports implicating genomic rearrangements of CNTN4 in autism spectrum disorders (ASDs) and 3p- microdeletion syndrome, we undertook sequencing of the coding regions of the gene in a local ASD cohort in comparison with a set of controls. Unique missense variants were identified in 4 of 75 unrelated individuals with ASD, as well as in 1 of 107 controls. All of the amino acid substitutions were nonsynonomous, occurred at evolutionarily conserved positions, and were, thus, felt likely to be deleterious. However, these data did not reach statistical significance, nor did the variants segregate with disease within all of the ASD families. Finally, there was no detectable difference in binding of two of the variants to the interacting protein PTPRG in vitro. Thus, additional larger studies will be necessary to determine whether CNTN4 functions as an autism susceptibility locus in combination with other genetic and/or environmental factors.

Significant contributions of the extraembryonic membranes and maternal genotype to the placental pathology in heterozygous Nsdhl deficient female embryos.

Mutations in the gene encoding the cholesterol biosynthetic enzyme NSDHL are associated with the X-linked male-lethal bare patches (Bpa) mouse. Mutant male embryos for several Nsdhl alleles die in midgestation with placental insufficiency. We examined here a possible role of the maternal genotype in such placental pathology. Pre-pregnancy plasma cholesterol levels were similar between wild-type (WT) and Bpa(1H)/+ dams fed a standard, cholesterol-free diet. However, there was a marked decrease in cholesterol levels between embryonic day (E)8.5 and E10.5 for both genotypes. Further, there was a significant lag between E11.5 and E13.5 (P = 0.0011) in the recovery of levels in Bpa(1H)/+ dams to their pre-pregnancy values. To investigate possible effects of the maternal genotype on fetal placentation, we generated transgenic mice that expressed human NSDHL and rescued the male lethality of the Bpa(1H) null allele. We then compared placenta area at E10.5 in WT and Bpa(1H)/+ female embryos where the mutant X chromosome was transmitted from a heterozygous mother or a rescued mutant father. In mutant conceptuses, placental areas were approximately 50% less than WT. Surprisingly, expression of Nsdhl in trophoblast lineages of the placenta and yolk sac endoderm, which occurs only from the maternally inherited allele in a female embryo, had the largest effect on placental area (-0.681 mm(2); P < 0.0001). The maternal genotype had a smaller effect, independent of the fetal genotype (-0.283 mm(2); P = 0.024). These data demonstrate significant effects of the mother and fetal membranes on pregnancy outcome, with possible implications for cholesterol homeostasis during human pregnancy.