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

Re-analysis of whole genome sequencing ends a diagnostic odyssey: Case report of an RNU4-2 related neurodevelopmental disorder

Despite increasing knowledge of disease-causing genes in human genetics, approximately half of the individuals affected by neurodevelopmental disorders remain genetically undiagnosed. Part of this missing heritability might be caused by genetic variants outside of protein-coding genes, which are not routinely diagnostically investigated. A recent preprint identified de novo variants in the non-coding spliceosomal snRNA gene RNU4-2 as a cause of a frequent novel syndromic neurodevelopmental disorder. Here we mined 164 whole genome sequencing (WGS) trios from individuals with neurodevelopmental or multiple congenital anomaly disorders that received diagnostic genomic investigations at our clinic. We identify a recurrent de novo RNU4-2 variant (NR_003137.2(RNU4-2):n.64_65insT) in a 5-year-old girl with severe global developmental delay, hypotonia, microcephaly, and seizures that likely explains her phenotype, given that extensive previous genetic investigations failed to identify an alternative cause. We present detailed phenotyping of the individual obtained during a 5-year follow-up. This includes photographs showing recognizable facial features for this novel disorder, which might allow prioritizing other currently unexplained affected individuals sharing similar facial features for targeted investigations of RNU4-2. This case illustrates the power of re-analysis to solve previously unexplained cases even when a diagnostic genome remains negative.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Analysis of 3760 hematologic malignancies reveals rare transcriptomic aberrations of driver genes

BACKGROUND

Rare oncogenic driver events, particularly affecting the expression or splicing of driver genes, are suspected to substantially contribute to the large heterogeneity of hematologic malignancies. However, their identification remains challenging.

METHODS

To address this issue, we generated the largest dataset to date of matched whole genome sequencing and total RNA sequencing of hematologic malignancies from 3760 patients spanning 24 disease entities. Taking advantage of our dataset size, we focused on discovering rare regulatory aberrations. Therefore, we called expression and splicing outliers using an extension of the workflow DROP (Detection of RNA Outliers Pipeline) and AbSplice, a variant effect predictor that identifies genetic variants causing aberrant splicing. We next trained a machine learning model integrating these results to prioritize new candidate disease-specific driver genes.

RESULTS

We found a median of seven expression outlier genes, two splicing outlier genes, and two rare splice-affecting variants per sample. Each category showed significant enrichment for already well-characterized driver genes, with odds ratios exceeding three among genes called in more than five samples. On held-out data, our integrative modeling significantly outperformed modeling based solely on genomic data and revealed promising novel candidate driver genes. Remarkably, we found a truncated form of the low density lipoprotein receptor LRP1B transcript to be aberrantly overexpressed in about half of hairy cell leukemia variant (HCL-V) samples and, to a lesser extent, in closely related B-cell neoplasms. This observation, which was confirmed in an independent cohort, suggests LRP1B as a novel marker for a HCL-V subclass and a yet unreported functional role of LRP1B within these rare entities.

CONCLUSIONS

Altogether, our census of expression and splicing outliers for 24 hematologic malignancy entities and the companion computational workflow constitute unique resources to deepen our understanding of rare oncogenic events in hematologic cancers.

The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

Genomes in clinical care

In the era of precision medicine, genome sequencing (GS) has become more affordable and the importance of genomics and multi-omics in clinical care is increasingly being recognized. However, how to scale and effectively implement GS on an institutional level remains a challenge for many. Here, we present Genome First and Ge-Med, two clinical implementation studies focused on identifying the key pillars and processes that are required to make routine GS and predictive genomics a reality in the clinical setting. We describe our experience and lessons learned for a variety of topics including test logistics, patient care processes, data reporting, and infrastructure. Our model of providing clinical care and comprehensive genomic analysis from a single source may be used by other centers with a similar structure to facilitate the implementation of omics-based personalized health concepts in medicine.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases

PURPOSE

Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION).

METHODS

PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation.

RESULTS

A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants.

CONCLUSION

GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network (UDN) and Genomics Research to Elucidate the Genetics of Rare Disease (GREGoR) Consortium. We identified 2,800 genes with build-dependent quantification across six routinely-collected biospecimens, including 1,391 protein-coding genes and 341 known rare disease genes. We further observed multiple genes that only have detectable expression in a subset of genome builds. Finally, we characterized how genome build impacts the detection of outlier transcriptomic events. Combined, we provide a database of genes impacted by build choice, and recommend that transcriptomics-guided analyses and diagnoses are cross-referenced with these data for robustness.

Integrating RNA-Seq into genome sequencing workflow enhances the analysis of structural variants causing neurodevelopmental disorders

BACKGROUND

Molecular diagnosis of neurodevelopmental disorders (NDDs) is mainly based on exome sequencing (ES), with a diagnostic yield of 31% for isolated and 53% for syndromic NDD. As sequencing costs decrease, genome sequencing (GS) is gradually replacing ES for genome-wide molecular testing. As many variants detected by GS only are in deep intronic or non-coding regions, the interpretation of their impact may be difficult. Here, we showed that integrating RNA-Seq into the GS workflow can enhance the analysis of the molecular causes of NDD, especially structural variants (SVs), by providing valuable complementary information such as aberrant splicing, aberrant expression and monoallelic expression.

METHODS

We performed trio-GS on a cohort of 33 individuals with NDD for whom ES was inconclusive. RNA-Seq on skin fibroblasts was then performed in nine individuals for whom GS was inconclusive and optical genome mapping (OGM) was performed in two individuals with an SV of unknown significance.

RESULTS

We identified pathogenic or likely pathogenic variants in 16 individuals (48%) and six variants of uncertain significance. RNA-Seq contributed to the interpretation in three individuals, and OGM helped to characterise two SVs.

CONCLUSION

Our study confirmed that GS significantly improves the diagnostic performance of NDDs. However, most variants detectable by GS alone are structural or located in non-coding regions, which can pose challenges for interpretation. Integration of RNA-Seq data overcame this limitation by confirming the impact of variants at the transcriptional or regulatory level. This result paves the way for new routinely applicable diagnostic protocols.

Solving the unsolved genetic epilepsies: Current and future perspectives

Many patients with epilepsy undergo exome or genome sequencing as part of a diagnostic workup; however, many remain genetically unsolved. There are various factors that account for negative results in exome/genome sequencing for patients with epilepsy: (1) the underlying cause is not genetic; (2) there is a complex polygenic explanation; (3) the illness is monogenic but the causative gene remains to be linked to a human disorder; (4) family segregation with reduced penetrance; (5) somatic mosaicism or the complexity of, for example, a structural rearrangement; or (6) limited knowledge or diagnostic tools that hinder the proper classification of a variant, resulting in its designation as a variant of unknown significance. The objective of this review is to outline some of the diagnostic options that lie beyond the exome/genome, and that might become clinically relevant within the foreseeable future. These options include: (1) re-analysis of older exome/genome data as knowledge increases or symptoms change; (2) looking for somatic mosaicism or long-read sequencing to detect low-complexity repeat variants or specific structural variants missed by traditional exome/genome sequencing; (3) exploration of the non-coding genome including disruption of topologically associated domains, long range non-coding RNA, or other regulatory elements; and finally (4) transcriptomics, DNA methylation signatures, and metabolomics as complementary diagnostic methods that may be used in the assessment of variants of unknown significance. Some of these tools are currently not integrated into standard diagnostic workup. However, it is reasonable to expect that they will become increasingly available and improve current diagnostic capabilities, thereby enabling precision diagnosis in patients who are currently undiagnosed.

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.