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.

Analysis of urinary oligosaccharide excretion patterns by UHPLC/HRAM mass spectrometry for screening of lysosomal storage disorders.

Oligosaccharidoses, sphingolipidoses and mucolipidoses are lysosomal storage disorders (LSDs) in which defective breakdown of glycan-side chains of glycosylated proteins and glycolipids leads to the accumulation of incompletely degraded oligosaccharides within lysosomes. In metabolic laboratories, these disorders are commonly diagnosed by thin-layer chromatography (TLC) but more recently also mass spectrometry-based approaches have been published. To expand the possibilities to screen for these diseases, we developed an ultra-high-performance liquid chromatography (UHPLC) with a high-resolution accurate mass (HRAM) mass spectrometry (MS) screening platform, together with an open-source iterative bioinformatics pipeline. This pipeline generates comprehensive biomarker profiles and allows for extensive quality control (QC) monitoring. Using this platform, we were able to identify α-mannosidosis, ß-mannosidosis, α-N-acetylgalactosaminidase deficiency, sialidosis, galactosialidosis, fucosidosis, aspartylglucosaminuria, GM1 gangliosidosis, GM2 gangliosidosis (M. Sandhoff) and mucolipidosis II/III in patient samples. Aberrant urinary oligosaccharide excretions were also detected for other disorders, including NGLY1 congenital disorder of deglycosylation, sialic acid storage disease, MPS type IV B and GSD II (Pompe disease). For the latter disorder, we identified heptahexose (Hex7), as a potential urinary biomarker, in addition to glucose tetrasaccharide (Glc4), for the diagnosis and monitoring of young onset cases of Pompe disease. Occasionally, so-called "neonate" biomarker profiles were observed in young patients, which were probably due to nutrition. Our UHPLC/HRAM-MS screening platform can easily be adopted in biochemical laboratories and allows for simple and robust screening and straightforward interpretation of the screening results to detect disorders in which aberrant oligosaccharides accumulate.

Uveal Melanoma Patients Have a Distinct Metabolic Phenotype in Peripheral Blood.

Uveal melanomas (UM) are detected earlier. Consequently, tumors are smaller, allowing for novel eye-preserving treatments. This reduces tumor tissue available for genomic profiling. Additionally, these small tumors can be hard to differentiate from nevi, creating the need for minimally invasive detection and prognostication. Metabolites show promise as minimally invasive detection by resembling the biological phenotype. In this pilot study, we determined metabolite patterns in the peripheral blood of UM patients (n = 113) and controls (n = 46) using untargeted metabolomics. Using a random forest classifier (RFC) and leave-one-out cross-validation, we confirmed discriminatory metabolite patterns in UM patients compared to controls with an area under the curve of the receiver operating characteristic of 0.99 in both positive and negative ion modes. The RFC and leave-one-out cross-validation did not reveal discriminatory metabolite patterns in high-risk versus low-risk of metastasizing in UM patients. Ten-time repeated analyses of the RFC and LOOCV using 50% randomly distributed samples showed similar results for UM patients versus controls and prognostic groups. Pathway analysis using annotated metabolites indicated dysregulation of several processes associated with malignancies. Consequently, minimally invasive metabolomics could potentially allow for screening as it distinguishes metabolite patterns that are putatively associated with oncogenic processes in the peripheral blood plasma of UM patients from controls at the time of diagnosis.

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.

Transcriptomics of Epstein-Barr virus aids to the classification of T-cell evasion in nasopharyngeal carcinoma.

Epstein-Barr virus (EBV) contributes to oncogenesis and immune evasion in nasopharyngeal carcinoma (NPC). At present, an aggregated, higher-level view on the impact of EBV genes toward the immune microenvironment of NPC is lacking. To this end, we have interrogated tumor-derived RNA sequences of 106 treatment-naive NPC patients for 98 EBV transcripts, and captured the presence of 10 different immune cell populations as well as 23 different modes of T-cell evasion. We discovered 3 clusters of EBV genes that each associate with distinct immunophenotypes of NPC. Cluster 1 associated with gene sets related to immune cell recruitment, such as those encoding for chemoattractants and their receptors. Cluster 2 associated with antigen processing and presentation, such as interferon-related genes, whereas cluster 3 associated with presence of M1-like macrophages, absence of CD4+ T cells, and oncogenic pathways, such as the nuclear factor kappa light-chain enhancer of activated B-cell pathway. We discuss these 3 EBV clusters regarding their potential for stratification for T-cell immunity in NPC together with the next steps needed to validate such therapeutic value.

Benchmarking Outlier Detection Methods for Detecting IEM Patients in Untargeted Metabolomics Data.

Untargeted metabolomics (UM) is increasingly being deployed as a strategy for screening patients that are suspected of having an inborn error of metabolism (IEM). In this study, we examined the potential of existing outlier detection methods to detect IEM patient profiles. We benchmarked 30 different outlier detection methods when applied to three untargeted metabolomics datasets. Our results show great differences in IEM detection performances across the various methods. The methods DeepSVDD and R-graph performed most consistently across the three metabolomics datasets. For datasets with a more balanced number of samples-to-features ratio, we found that AE reconstruction error, Mahalanobis and PCA reconstruction error also performed well. Furthermore, we demonstrated the importance of a PCA transform prior to applying an outlier detection method since we observed that this increases the performance of several outlier detection methods. For only one of the three metabolomics datasets, we observed clinically satisfying performances for some outlier detection methods, where we were able to detect 90% of the IEM patient samples while detecting no false positives. These results suggest that outlier detection methods have the potential to aid the clinical investigator in routine screening for IEM using untargeted metabolomics data, but also show that further improvements are needed to ensure clinically satisfying performances.

Using Out-of-Batch Reference Populations to Improve Untargeted Metabolomics for Screening Inborn Errors of Metabolism.

Untargeted metabolomics is an emerging technology in the laboratory diagnosis of inborn errors of metabolism (IEM). Analysis of a large number of reference samples is crucial for correcting variations in metabolite concentrations that result from factors, such as diet, age, and gender in order to judge whether metabolite levels are abnormal. However, a large number of reference samples requires the use of out-of-batch samples, which is hampered by the semi-quantitative nature of untargeted metabolomics data, i.e., technical variations between batches. Methods to merge and accurately normalize data from multiple batches are urgently needed. Based on six metrics, we compared the existing normalization methods on their ability to reduce the batch effects from nine independently processed batches. Many of those showed marginal performances, which motivated us to develop Metchalizer, a normalization method that uses 10 stable isotope-labeled internal standards and a mixed effect model. In addition, we propose a regression model with age and sex as covariates fitted on reference samples that were obtained from all nine batches. Metchalizer applied on log-transformed data showed the most promising performance on batch effect removal, as well as in the detection of 195 known biomarkers across 49 IEM patient samples and performed at least similar to an approach utilizing 15 within-batch reference samples. Furthermore, our regression model indicates that 6.5-37% of the considered features showed significant age-dependent variations. Our comprehensive comparison of normalization methods showed that our Log-Metchalizer approach enables the use out-of-batch reference samples to establish clinically-relevant reference values for metabolite concentrations. These findings open the possibilities to use large scale out-of-batch reference samples in a clinical setting, increasing the throughput and detection accuracy.

Untargeted Metabolomics-Based Screening Method for Inborn Errors of Metabolism using Semi-Automatic Sample Preparation with an UHPLC- Orbitrap-MS Platform.

Routine diagnostic screening of inborn errors of metabolism (IEM) is currently performed by different targeted analyses of known biomarkers. This approach is time-consuming, targets a limited number of biomarkers and will not identify new biomarkers. Untargeted metabolomics generates a global metabolic phenotype and has the potential to overcome these issues. We describe a novel, single platform, untargeted metabolomics method for screening IEM, combining semi-automatic sample preparation with pentafluorophenylpropyl phase (PFPP)-based UHPLC- Orbitrap-MS. We evaluated analytical performance and diagnostic capability of the method by analysing plasma samples of 260 controls and 53 patients with 33 distinct IEM. Analytical reproducibility was excellent, with peak area variation coefficients below 20% for the majority of the metabolites. We illustrate that PFPP-based chromatography enhances identification of isomeric compounds. Ranked z-score plots of metabolites annotated in IEM samples were reviewed by two laboratory specialists experienced in biochemical genetics, resulting in the correct diagnosis in 90% of cases. Thus, our untargeted metabolomics platform is robust and differentiates metabolite patterns of different IEMs from those of controls. We envision that the current approach to diagnose IEM, using numerous tests, will eventually be replaced by untargeted metabolomics methods, which also have the potential to discover novel biomarkers and assist in interpretation of genetic data.