Web-accessible application for identifying pathogenic transcripts with RNA-seq: Increased sensitivity in diagnosis of neurodevelopmental disorders.

For neurodevelopmental disorders (NDDs), a molecular diagnosis is key for management, predicting outcome, and counseling. Often, routine DNA-based tests fail to establish a genetic diagnosis in NDDs. Transcriptome analysis (RNA sequencing [RNA-seq]) promises to improve the diagnostic yield but has not been applied to NDDs in routine diagnostics. Here, we explored the diagnostic potential of RNA-seq in 96 individuals including 67 undiagnosed subjects with NDDs. We performed RNA-seq on single individuals' cultured skin fibroblasts, with and without cycloheximide treatment, and used modified OUTRIDER Z scores to detect gene expression outliers and mis-splicing by exonic and intronic outliers. Analysis was performed by a user-friendly web application, and candidate pathogenic transcriptional events were confirmed by secondary assays. We identified intragenic deletions, monoallelic expression, and pseudoexonic insertions but also synonymous and non-synonymous variants with deleterious effects on transcription, increasing the diagnostic yield for NDDs by 13%. We found that cycloheximide treatment and exonic/intronic Z score analysis increased detection and resolution of aberrant splicing. Importantly, in one individual mis-splicing was found in a candidate gene nearly matching the individual's specific phenotype. However, pathogenic splicing occurred in another neuronal-expressed gene and provided a molecular diagnosis, stressing the need to customize RNA-seq. Lastly, our web browser application allowed custom analysis settings that facilitate diagnostic application and ranked pathogenic transcripts as top candidates. Our results demonstrate that RNA-seq is a complementary method in the genomic diagnosis of NDDs and, by providing accessible analysis with improved sensitivity, our transcriptome analysis approach facilitates wider implementation of RNA-seq in routine genome diagnostics.

An 8q24 Gain in Pancreatic Juice Is a Candidate Biomarker for the Detection of Pancreatic Cancer.

Secretin-stimulated pancreatic juice (PJ), collected from the duodenum, presents a valuable biomarker source for the (earlier) detection of pancreatic cancer (PC). Here, we evaluate the feasibility and performance of shallow sequencing to detect copy number variations (CNVs) in cell-free DNA (cfDNA) from PJ for PC detection. First, we confirmed the feasibility of shallow sequencing in PJ (n = 4), matched plasma (n = 3) and tissue samples (n = 4, microarray). Subsequently, shallow sequencing was performed on cfDNA from PJ of 26 cases (25 sporadic PC, 1 high-grade dysplasia) and 19 controls with a hereditary or familial increased risk of PC. 40 of the 45 PJ samples met the quality criteria for cfDNA analysis. Nine individuals had an 8q24 gain (oncogene MYC; 23%; eight cases (33%) and one control (6%), p = 0.04); six had both a 2q gain (STAT1) and 5p loss (CDH10; 15%; four cases (7%) and two controls (13%), p = 0.72). The presence of an 8q24 gain differentiated the cases and controls, with a sensitivity of 33% (95% CI 16-55%) and specificity of 94% (95% CI 70-100%). The presence of either an 8q24 or 2q gain with a 5p loss was related to a sensitivity of 50% (95% CI 29-71%) and specificity of 81% (95% CI 54-96%). Shallow sequencing of PJ is feasible. The presence of an 8q24 gain in PJ shows promise as a biomarker for the detection of PC. Further research is required with a larger sample size and consecutively collected samples in high-risk individuals prior to implementation in a surveillance cohort.

High-yield identification of pathogenic NF1 variants by skin fibroblast transcriptome screening after apparently normal diagnostic DNA testing.

Neurofibromatosis type 1 (NF1) is caused by inactivating mutations in NF1. Due to the size, complexity, and high mutation rate at the NF1 locus, the identification of causative variants can be challenging. To obtain a molecular diagnosis in 15 individuals meeting diagnostic criteria for NF1, we performed transcriptome analysis (RNA-seq) on RNA obtained from cultured skin fibroblasts. In each case, routine molecular DNA diagnostics had failed to identify a disease-causing variant in NF1. A pathogenic variant or abnormal mRNA splicing was identified in 13 cases: 6 deep intronic variants and 2 transposon insertions causing noncanonical splicing, 3 postzygotic changes, 1 branch point mutation and, in 1 case, abnormal splicing for which the responsible DNA change remains to be identified. These findings helped resolve the molecular findings for an additional 17 individuals in multiple families with NF1, demonstrating the utility of skin-fibroblast-based transcriptome analysis for molecular diagnostics. RNA-seq improves mutation detection in NF1 and provides a powerful complementary approach to DNA-based methods. Importantly, our approach is applicable to other genetic disorders, particularly those caused by a wide variety of variants in a limited number of genes and specifically for individuals in whom routine molecular DNA diagnostics did not identify the causative variant.

Integration of metabolomics with genomics: Metabolic gene prioritization using metabolomics data and genomic variant (CADD) scores.

The integration of metabolomics data with sequencing data is a key step towards improving the diagnostic process for finding the disease-causing genetic variant(s) in patients suspected of having an inborn error of metabolism (IEM). The measured metabolite levels could provide additional phenotypical evidence to elucidate the degree of pathogenicity for variants found in genes associated with metabolic processes. We present a computational approach, called Reafect, that calculates for each reaction in a metabolic pathway a score indicating whether that reaction is deficient or not. When calculating this score, Reafect takes multiple factors into account: the magnitude and sign of alterations in the metabolite levels, the reaction distances between metabolites and reactions in the pathway, and the biochemical directionality of the reactions. We applied Reafect to untargeted metabolomics data of 72 patient samples with a known IEM and found that in 81% of the cases the correct deficient enzyme was ranked within the top 5% of all considered enzyme deficiencies. Next, we integrated Reafect with Combined Annotation Dependent Depletion (CADD) scores (a measure for gene variant deleteriousness) and ranked the metabolic genes of 27 IEM patients. We observed that this integrated approach significantly improved the prioritization of the genes containing the disease-causing variant when compared with the two approaches individually. For 15/27 IEM patients the correct affected gene was ranked within the top 0.25% of the set of potentially affected genes. Together, our findings suggest that metabolomics data improves the identification of affected genes in patients suffering from IEM.

Dominant-acting CSF1R variants cause microglial depletion and altered astrocytic phenotype in zebrafish and adult-onset leukodystrophy.

Tissue-resident macrophages of the brain, including microglia, are implicated in the pathogenesis of various CNS disorders and are possible therapeutic targets by their chemical depletion or replenishment by hematopoietic stem cell therapy. Nevertheless, a comprehensive understanding of microglial function and the consequences of microglial depletion in the human brain is lacking. In human disease, heterozygous variants in CSF1R, encoding the Colony-stimulating factor 1 receptor, can lead to adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP) possibly caused by microglial depletion. Here, we investigate the effects of ALSP-causing CSF1R variants on microglia and explore the consequences of microglial depletion in the brain. In intermediate- and late-stage ALSP post-mortem brain, we establish that there is an overall loss of homeostatic microglia and that this is predominantly seen in the white matter. By introducing ALSP-causing missense variants into the zebrafish genomic csf1ra locus, we show that these variants act dominant negatively on the number of microglia in vertebrate brain development. Transcriptomics and proteomics on relatively spared ALSP brain tissue validated a downregulation of microglia-associated genes and revealed elevated astrocytic proteins, possibly suggesting involvement of astrocytes in early pathogenesis. Indeed, neuropathological analysis and in vivo imaging of csf1r zebrafish models showed an astrocytic phenotype associated with enhanced, possibly compensatory, endocytosis. Together, our findings indicate that microglial depletion in zebrafish and human disease, likely as a consequence of dominant-acting pathogenic CSF1R variants, correlates with altered astrocytes. These findings underscore the unique opportunity CSF1R variants provide to gain insight into the roles of microglia in the human brain, and the need to further investigate how microglia, astrocytes, and their interactions contribute to white matter homeostasis.

Non-invasive prenatal diagnosis for translocation carriers-YES please or NO go?

INTRODUCTION: The presence of an unbalanced familial translocation can be reliably assessed in the cytotrophoblast of chorionic villi. However, carriers of a balanced translocation often decline invasive testing. This study aimed to investigate whether an unbalanced translocation can also be diagnosed in cell free DNA by whole-genome non-invasive prenatal screening (NIPS). MATERIAL AND METHODS: Pregnant women carrying a fetus with an unbalanced familial translocation, for whom NIPS as well as microarray data were available, were included in this retrospective assessment. NIPS was performed in the course of the TRIDENT study. RESULTS: In 12 cases, both NIPS and microarray data were available. In 10 of 12 cases the unbalanced translocation was correctly identified by NIPS without prior knowledge on parental translocation. One was missed because the fetal fraction was too low. One was missed because of technical restrictions in calling 16p gains. CONCLUSIONS: This study supports the hypothesis that routine NIPS may be used for prenatal diagnosis of unbalanced inheritance of familial translocations, especially with prior knowledge of the translocation allowing focused examination of the involved chromosomal regions. Our study showed that routine shallow sequencing designed for aneuploidy detection in cell free DNA may be sufficient for higher resolution NIPS, if specialized copy number software is used and if sufficient fetal fraction is present.

Loss of UGP2 in brain leads to a severe epileptic encephalopathy, emphasizing that bi-allelic isoform-specific start-loss mutations of essential genes can cause genetic diseases.

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early-onset, therapy-resistant seizures and developmental delay. Here we report on 22 individuals from 15 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly, visual disturbance and similar minor dysmorphisms. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A > G) in the essential UDP-glucose pyrophosphorylase (UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform, which is predominant in brain. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in neural stem cells, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modeled during pluripotent stem cell differentiation in vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2b in vivo in zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation in UGP2 as a cause of a novel autosomal recessive DEE syndrome. Importantly, it also shows that isoform-specific start-loss mutations causing expression loss of a tissue-relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies.