Diagnostic utility of transcriptome sequencing for rare Mendelian diseases

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

RNA Sequencing in Disease Diagnosis

RNA sequencing (RNA-seq) enables the accurate measurement of multiple transcriptomic phenotypes for modeling the impacts of disease variants. Advances in technologies, experimental protocols, and analysis strategies are rapidly expanding the application of RNA-seq to identify disease biomarkers, tissue- and cell-type-specific impacts, and the spatial localization of disease-associated mechanisms. Ongoing international efforts to construct biobank-scale transcriptomic repositories with matched genomic data across diverse population groups are further increasing the utility of RNA-seq approaches by providing large-scale normative reference resources. The availability of these resources, combined with improved computational analysis pipelines, has enabled the detection of aberrant transcriptomic phenotypes underlying rare diseases. Further expansion of these resources, across both somatic and developmental tissues, is expected to soon provide unprecedented insights to resolve disease origin, mechanism of action, and causal gene contributions, suggesting the continued high utility of RNA-seq in disease diagnosis.

Expanding the genetic and phenotypic landscape of replication factor C complex-related disorders: RFC4 deficiency is linked to a multisystemic disorder

The precise regulation of DNA replication is vital for cellular division and genomic integrity. Central to this process is the replication factor C (RFC) complex, encompassing five subunits, which loads proliferating cell nuclear antigen onto DNA to facilitate the recruitment of replication and repair proteins and enhance DNA polymerase processivity. While RFC1's role in cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS) is known, the contributions of RFC2-5 subunits on human Mendelian disorders is largely unexplored. Our research links bi-allelic variants in RFC4, encoding a core RFC complex subunit, to an undiagnosed disorder characterized by incoordination and muscle weakness, hearing impairment, and decreased body weight. We discovered across nine affected individuals rare, conserved, predicted pathogenic variants in RFC4, all likely to disrupt the C-terminal domain indispensable for RFC complex formation. Analysis of a previously determined cryo-EM structure of RFC bound to proliferating cell nuclear antigen suggested that the variants disrupt interactions within RFC4 and/or destabilize the RFC complex. Cellular studies using RFC4-deficient HeLa cells and primary fibroblasts demonstrated decreased RFC4 protein, compromised stability of the other RFC complex subunits, and perturbed RFC complex formation. Additionally, functional studies of the RFC4 variants affirmed diminished RFC complex formation, and cell cycle studies suggested perturbation of DNA replication and cell cycle progression. Our integrated approach of combining in silico, structural, cellular, and functional analyses establishes compelling evidence that bi-allelic loss-of-function RFC4 variants contribute to the pathogenesis of this multisystemic disorder. These insights broaden our understanding of the RFC complex and its role in human health and disease.

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.

Whole exome sequencing in Serbian patients with hereditary spastic paraplegia

Hereditary spastic paraplegia (HSP) is a group of neurodegenerative diseases with a high genetic and clinical heterogeneity. Numerous HSP patients remain genetically undiagnosed despite screening for known genetic causes of HSP. Therefore, identification of novel variants and genes is needed. Our previous study analyzed 74 adult Serbian HSP patients from 65 families using panel of the 13 most common HSP genes in combination with a copy number variation analysis. Conclusive genetic findings were established in 23 patients from 19 families (29%). In the present study, nine patients from nine families previously negative on the HSP gene panel were selected for the whole exome sequencing (WES). Further, 44 newly diagnosed adult HSP patients from 44 families were sent to WES directly, since many studies showed WES may be used as the first step in HSP diagnosis. WES analysis of cohort 1 revealed a likely genetic cause in five (56%) of nine HSP families, including variants in the ETHE1, ZFYVE26, RNF170, CAPN1, and WASHC5 genes. In cohort 2, possible causative variants were found in seven (16%) of 44 patients (later updated to 27% when other diagnosis were excluded), comprising six different genes: SPAST, SPG11, WASCH5, KIF1A, KIF5A, and ABCD1. These results expand the genetic spectrum of HSP patients in Serbia and the region with implications for molecular genetic diagnosis and future causative therapies. Wide HSP panel can be the first step in diagnosis, alongside with the copy number variation (CNV) analysis, while WES should be performed after.

Expedited Exome Reanalysis Following Deep Phenotyping and Muscle Biopsy in Suspected Mitochondrial Disorder

BACKGROUND

Exome sequencing (ES) is a useful tool in diagnosing suspected mitochondrial disease but can miss pathogenic variants for several reasons. Additional testing, such as muscle biopsy or biochemical testing, can be helpful in exome-negative cases.

METHODS

We report a patient who presented with repeated episodes of lactic acidosis and failure to thrive.

RESULTS

ES and mitochondrial sequencing were initially negative but clinical suspicion for mitochondrial disease remained high. After muscle biopsy showed evidence of mitochondrial dysfunction, the ES was reanalyzed and revealed novel variants in AARS2.

CONCLUSION

This case demonstrates the importance of muscle biopsy and biochemical testing in evaluating patients with a high suspicion of mitochondrial disease, even in the genomics era. Closed-loop communication between molecular genetics laboratories and clinical geneticists is an important step to help establish diagnosis in unsolved cases.

Reanalysis of RNA sequencing data ends diagnostic odyssey and expands the phenotypic spectrum of congenital titinopathy

Although next-generation sequencing has enabled diagnoses for many patients with Mendelian disorders, the majority remain undiagnosed. Here, we present a sibling pair who were clinically diagnosed with Escobar syndrome, however targeted gene testing was negative. Exome sequencing (ES), and later genome sequencing (GS), revealed compound heterozygous TTN variants in both siblings, a maternally inherited frameshift variant [(NM_133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM_133378.4):c.12322G > A; p.(Asp4108Asn)]. This result was considered nondiagnostic due to poor clinical fit and limited pathogenicity evidence for the missense variant of uncertain significance (VUS). Following initial nondiagnostic RNA sequencing (RNAseq) on muscle and further pursuit of other variants detected on the ES/GS, a reanalysis of noncanonical splice sites in the muscle transcriptome identified an out-of-frame exon retraction in TTN, near the known VUS. Interim literature included reports of patients with similar TTN variants who had phenotypic concordance with the siblings, and a diagnosis of a congenital titinopathy was given 4 years after the TTN variants had been initially reported. This report highlights the value of reanalysis of RNAseq with a different approach, expands the phenotypic spectrum of congenital titinopathy and also illustrates how a perceived phenotypic mismatch, and failure to consider known variants, can result in a prolongation of the diagnostic journey.

The expanding diagnostic toolbox for rare genetic diseases

Genomic technologies, such as targeted, exome and short-read genome sequencing approaches, have revolutionized the care of patients with rare genetic diseases. However, more than half of patients remain without a diagnosis. Emerging approaches from research-based settings such as long-read genome sequencing and optical genome mapping hold promise for improving the identification of disease-causal genetic variants. In addition, new omic technologies that measure the transcriptome, epigenome, proteome or metabolome are showing great potential for variant interpretation. As genetic testing options rapidly expand, the clinical community needs to be mindful of their individual strengths and limitations, as well as remaining challenges, to select the appropriate diagnostic test, correctly interpret results and drive innovation to address insufficiencies. If used effectively - through truly integrative multi-omics approaches and data sharing - the resulting large quantities of data from these established and emerging technologies will greatly improve the interpretative power of genetic and genomic diagnostics for rare diseases.

Combined exome and whole transcriptome sequencing identifies a de novo intronic SRCAP variant causing DEHMBA syndrome with severe sleep disorder

Rare heterozygous variants in exons 33-34 of the SRCAP gene are associated with Floating-Harbor syndrome and have a dominant-negative mechanism of action. At variance, heterozygous null alleles falling in other parts of the same gene cause developmental delay, hypotonia, musculoskeletal defects, and behavioral abnormalities (DEHMBA) syndrome. We report an 18-year-old man with DEHMBA syndrome and obstructive sleep apnea, who underwent exome sequencing (ES) and whole transcriptome sequencing (WTS) on peripheral blood. Trio analysis prioritized the de novo heterozygous c.5658+5 G > A variant. WTS promptly demostrated four different abnormal transcripts affecting >40% of the reads, three of which leading to a frameshift. This study demonstrated the efficacy of a combined ES-WTS approach in solving undiagnosed cases. We also speculated that sleep respiratory disorder may be an underdiagnosed complication of DEHMBA syndrome.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

Identification of Synonymous Pathogenic Variants in Monogenic Disorders by Integrating Exome with Transcriptome Sequencing

Exome sequencing is becoming a first-tier clinical diagnostic test for Mendelian diseases, drastically reducing the time and cost of diagnostic odyssey and improving the diagnosis rate. Despite its success, exome sequencing faces practical challenges in assessing the pathogenicity of numerous intronic and synonymous variants, leaving a significant proportion of patients undiagnosed. In this study, a whole-blood transcriptome database was constructed that showed the expression profile of 2981 Online Mendelian Inheritance in Man disease genes in blood samples. Meanwhile, a workflow integrating exome sequencing, blood transcriptome sequencing, and in silico prediction tools to identify and validate splicing-altering intronic or synonymous variants was proposed. Following this pipeline, seven synonymous variants in eight patients were discovered. Of these, the functional evidence of c.981G>A (PIGN), c.1161A>G (ALPL), c.858G>A (ATP6AP2), and c.1011G>T (MTHFR) have not been reported previously. RNA sequencing validation confirmed that these variants induced aberrant splicing, expanding the disease-causing variant spectrum of these genes. Overall, this study shows the feasibility of combining multi-omics data to identify splicing-altering variants, especially the power of RNA sequencing. It also reveals that synonymous variants, which often are overlooked in standard diagnostic approaches, comprise an important portion of unresolved genetic diseases.

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.

Inborn errors of immunity: A field without frontiers

The study of primary immunodeficiencies or inborn errors of immunity continues to drive our knowledge of the function of the human immune system. From the outset, the study of inborn errors has focused on unraveling genetic etiologies and molecular mechanisms. Aided by the continuous growth in genetic diagnostics, the field has moved from the study of an infection dominated phenotype to embrace and unravel diverse manifestations of autoinflammation, autoimmunity, malignancy, and severe allergy in all medical disciplines. It has now moved from the study of ultrarare presentations to producing meaningful impact in conditions as diverse as inflammatory bowel disease, neurological conditions, and hematology. Beyond offering immunogenetic diagnosis, the study of underlying inborn errors of immunity in these conditions points to targeted treatment which can be lifesaving.

RNA sequencing and target long-read sequencing reveal an intronic transposon insertion causing aberrant splicing

More than half of cases with suspected genetic disorders remain unsolved by genetic analysis using short-read sequencing such as exome sequencing (ES) and genome sequencing (GS). RNA sequencing (RNA-seq) and long-read sequencing (LRS) are useful for interpretation of candidate variants and detection of structural variants containing repeat sequences, respectively. Recently, adaptive sampling on nanopore sequencers enables target LRS more easily. Here, we present a Japanese girl with premature chromatid separation (PCS)/mosaic variegated aneuploidy (MVA) syndrome. ES detected a known pathogenic maternal heterozygous variant (c.1402-5A>G) in intron 10 of BUB1B (NM_001211.6), a known responsive gene for PCS/MVA syndrome with autosomal recessive inheritance. Minigene splicing assay revealed that almost all transcripts from the c.1402-5G allele have mis-splicing with 4-bp insertion. GS could not detect another pathogenic variant, while RNA-seq revealed abnormal reads in intron 2. To extensively explore variants in intron 2, we performed adaptive sampling and identified a paternal 3.0 kb insertion. Consensus sequence of 16 reads spanning the insertion showed that the insertion consists of Alu and SVA elements. Realignment of RNA-seq reads to the new reference sequence containing the insertion revealed that 16 reads have 5' splice site within the insertion and 3' splice site at exon 3, demonstrating causal relationship between the insertion and aberrant splicing. In addition, immunoblotting showed severely diminished BUB1B protein level in patient derived cells. These data suggest that detection of transcriptomic abnormalities by RNA-seq can be a clue for identifying pathogenic variants, and determination of insert sequences is one of merits of LRS.

Next-generation sequencing and bioinformatics in rare movement disorders

The ability to sequence entire exomes and genomes has revolutionized molecular testing in rare movement disorders, and genomic sequencing is becoming an integral part of routine diagnostic workflows for these heterogeneous conditions. However, interpretation of the extensive genomic variant information that is being generated presents substantial challenges. In this Perspective, we outline multidimensional strategies for genetic diagnosis in patients with rare movement disorders. We examine bioinformatics tools and computational metrics that have been developed to facilitate accurate prioritization of disease-causing variants. Additionally, we highlight community-driven data-sharing and case-matchmaking platforms, which are designed to foster the discovery of new genotype-phenotype relationships. Finally, we consider how multiomic data integration might optimize diagnostic success by combining genomic, epigenetic, transcriptomic and/or proteomic profiling to enable a more holistic evaluation of variant effects. Together, the approaches that we discuss offer pathways to the improved understanding of the genetic basis of rare movement disorders.

The first case of intellectual disability caused by novel compound heterozygosity for NUDT2 variants

NUDT2 is an enzyme important for maintaining the intracellular level of the diadenosine tetraphosphate (Ap4A). Bi-allelic loss of function variants in NUDT2 has recently been reported as a rare cause of intellectual disability (ID). Herein, we describe a Chinese girl with ID, attention deficit hyperactivity disorder (ADHD), and motor delays with abnormal walking posture and difficulty climbing stairs, who bears compound heterozygous variants c.34 C > T (p.R12*) and c.194T > G (p.I65R) in NUDT2. Homozygous variants c.34 C > T (p.R12*) or c.186del (p.A63Qfs*3) in NUDT2 were previously reported to cause ID. This is the first patient with ID due to compound heterozygous variants in NUDT2 and p.I65R is a novel missense variant. This study enriched the genotype and phenotype of NUDT2-related ID and supported the critical developmental involvement of NUDT2.

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.

Basis with RNA-Seq and WGCNA to explore the effect of Frankincense essential oil on dextran sodium sulfate-induced ulcerative colitis through MAPK/NF-κB signaling

PURPOSE

Frankincense has been shown in studies to have healing benefits for people with ulcerative colitis (UC). However, its underlying mechanisms have not been fully investigated. The objective of this study was to explore the potential molecular mechanisms of Frankincense essential oil (FREO) in improving dextran sodium sulfate (DSS)-induced UC from multiple perspectives.

METHODS

The FREO components were analyzed by GC-MS, and the interactions between the key active components and the mechanism of FREO were determined based on RNA-seq, "quantity-effect" weighting coefficient network pharmacology, WGCNA and pharmacodynamic experiments. The protection of FREO against DSS-induced UC mice was assessed by behavioral and pathological changes through mice. The expression of pro-inflammatory cytokines was measured using enzyme-linked immunosorbent assay. The expression of MAPK and NF-κB-related proteins by the Western Blotting and immunohistochemistry method.

RESULTS

Treatment with FREO significantly improved the symptoms of weight loss, diarrhea, stool blood, and colon shortening in UC mice. Reduced intestinal mucosal damage and the degree of inflammatory cell infiltration in the colon. Decreased TNF-α and IL-6 levels in mice's serum and inhibited phosphorylation of ERK, p65 in MAPK and NF-κB signaling.

CONCLUSION

FREO may decrease the inflammatory response to reduce the symptoms of UC by modulating the MAPK/ NF-κB pathway. This may be due to the synergistic interaction of the effective ingredient Hepten-2-yl tiglate, 6-methyl-5-, Isoneocembrene A and P-Cymene. This study provides a promising drug candidate and a new concept for the treatment of UC.

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.

The Diagnostic Landscape of Adult Neurogenetic Disorders

Neurogenetic diseases affect individuals across the lifespan, but accurate diagnosis remains elusive for many patients. Adults with neurogenetic disorders often undergo a long diagnostic odyssey, with multiple specialist evaluations and countless investigations without a satisfactory diagnostic outcome. Reasons for these diagnostic challenges include: (1) clinical features of neurogenetic syndromes are diverse and under-recognized, particularly those of adult-onset, (2) neurogenetic syndromes may manifest with symptoms that span multiple neurological and medical subspecialties, and (3) a positive family history may not be present or readily apparent. Furthermore, there is a large gap in the understanding of how to apply genetic diagnostic tools in adult patients, as most of the published literature focuses on the pediatric population. Despite these challenges, accurate genetic diagnosis is imperative to provide affected individuals and their families guidance on prognosis, recurrence risk, and, for an increasing number of disorders, offer targeted treatment. Here, we provide a framework for recognizing adult neurogenetic syndromes, describe the current diagnostic approach, and highlight studies using next-generation sequencing in different neurological disease cohorts. We also discuss diagnostic pitfalls, barriers to achieving a definitive diagnosis, and emerging technology that may increase the diagnostic yield of testing.

Uncovering the prominent role of satellite cells in paravertebral muscle development and aging by single-nucleus RNA sequencing

To uncover the role of satellite cells (SCs) in paravertebral muscle development and aging, we constructed a single-nucleus transcriptomic atlas of mouse paravertebral muscle across seven timepoints spanning the embryo (day 16.5) to old (month 24) stages. Eight cell types, including SCs, fast muscle cells, and slow muscle cells, were identified. An energy metabolism-related gene set, TCA CYCLE IN SENESCENCE, was enriched in SCs. Forty-two skeletal muscle disease-related genes were highly expressed in SCs and exhibited similar expression patterns. Among them, Pdha1 was the core gene in the TCA CYCLE IN SENESCENCE; Pgam2, Sod1, and Suclg1 are transcription factors closely associated with skeletal muscle energy metabolism. Transcription factor enrichment analysis of the 42 genes revealed that Myod1 and Mef2a were also highly expressed in SCs, which regulated Pdha1 expression and were associated with skeletal muscle development. These findings hint that energy metabolism may be pivotal in SCs development and aging. Three ligand-receptor pairs of extracellular matrix (ECM)-receptor interactions, Lamc1-Dag1, Lama2-Dag1, and Hspg2-Dag1, may play a vital role in SCs interactions with slow/fast muscle cells and SCs self-renewal. Finally, we built the first database of a skeletal muscle single-cell transcriptome, the Musculoskeletal Cell Atlas (http://www.mskca.tech), which lists 630,040 skeletal muscle cells and provides interactive visualization, a useful resource for revealing skeletal muscle cellular heterogeneity during development and aging. Our study could provide new targets and ideas for developing drugs to inhibit skeletal muscle aging and treat skeletal muscle diseases.

A common benign intronic deletion masking a pathogenic deep intronic PCCB variant - genome sequencing and RNA studies to the rescue

Propionic acidemia (PA) is an autosomal recessive metabolic disorder caused by variants in PCCA or PCCB, both sub-units of the propionyl-CoA carboxylase (PCC) enzyme. PCC is required for the catabolism of certain amino acids and odd-chain fatty acids. In its absence, the accumulated toxic metabolites cause metabolic acidosis, neurologic symptoms, multi-organ dysfunction and possible death. The clinical presentation of PA is highly variable, with typical onset in the neonatal or early infantile period. We encountered two families, whose children were diagnosed with PA. Exome sequencing (ES) failed to identify a pathogenic variant, and we proceeded with genome sequencing (GS), demonstrating homozygosity to a deep intronic PCCB variant. RNA analysis established that this variant creates a pseudoexon with a premature stop codon. The parents are variant carriers, though three of them display pseudo-homozygosity due to a common large benign intronic deletion on the second allele. The parental presumed homozygosity merits special attention, as it masked the causative variant at first, which was resolved only by RNA studies. Arriving at a rapid diagnosis, whether biochemical or genetic, can be crucial in directing lifesaving care, concluding the diagnostic odyssey, and allowing the family prenatal testing in subsequent pregnancies. This study demonstrates the power of integrative genetic studies in reaching a diagnosis, utilizing GS and RNA analysis to overcome ES limitations and define pathogenicity. Importantly, it highlights that intronic deletions should be taken into consideration when analyzing genomic data, so that pseudo-homozygosity would not be misinterpreted as true homozygosity, and pathogenic variants will not be mislabeled as benign.

Precision medicine in rare diseases: What is next?

Molecular diagnostics is a cornerstone of modern precision medicine, broadly understood as tailoring an individual's treatment, follow-up, and care based on molecular data. In rare diseases (RDs), molecular diagnoses reveal valuable information about the cause of symptoms, disease progression, familial risk, and in certain cases, unlock access to targeted therapies. Due to decreasing DNA sequencing costs, genome sequencing (GS) is emerging as the primary method for precision diagnostics in RDs. Several ongoing European initiatives for precision medicine have chosen GS as their method of choice. Recent research supports the role for GS as first-line genetic investigation in individuals with suspected RD, due to its improved diagnostic yield compared to other methods. Moreover, GS can detect a broad range of genetic aberrations including those in noncoding regions, producing comprehensive data that can be periodically reanalyzed for years to come when further evidence emerges. Indeed, targeted drug development and repurposing of medicines can be accelerated as more individuals with RDs receive a molecular diagnosis. Multidisciplinary teams in which clinical specialists collaborate with geneticists, genomics education of professionals and the public, and dialogue with patient advocacy groups are essential elements for the integration of precision medicine into clinical practice worldwide. It is also paramount that large research projects share genetic data and leverage novel technologies to fully diagnose individuals with RDs. In conclusion, GS increases diagnostic yields and is a crucial step toward precision medicine for RDs. Its clinical implementation will enable better patient management, unlock targeted therapies, and guide the development of innovative treatments.

Downstream Assays for Variant Resolution: Epigenetics, RNA Sequnncing, and Metabolomics

As the availability of advanced molecular testing like whole exome and genome sequencing expands, it comes with the added complication of interpreting inconclusive results, including determining the relevance of variants of uncertain significance or failing to find a variant in an otherwise suspected specific genetic disorder. This complication necessitates the use of alternative testing methods to gather more information in support of, or against, a particular genetic diagnosis. Therefore, new genome-wide approaches, including DNA epigenetic testing, RNA sequencing, and metabolomics, are increasingly being used to increase the diagnostic yield when used in conjunction with more conventional genetic tests.

Full-length Isoform Sequencing for Resolving the Molecular Basis of Charcot-Marie-Tooth 2A

Objectives

Transcript sequencing of patient-derived samples has been shown to improve the diagnostic yield for solving cases of suspected Mendelian conditions, yet the added benefit of full-length long-read transcript sequencing is largely unexplored.

Methods

We applied short-read and full-length transcript sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis.

Results

We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts the branch point critical for intron 6 splicing. Full-length long-read isoform complementary DNA (cDNA) sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates 5 distinct altered splicing transcripts. All 5 altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A).

Discussion

This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

mRNA sequencing provides new insights into the pathogenesis of Hirschsprung's disease in mice

PURPOSE

The aim of this study is to use RNA sequencing and RT-qPCR to identify the main susceptibility genes linked to the occurrence and development of Hirschsprung disease in the colonic tissues of EDNRBm1yzcm and wild mice.

METHODS

RNA was extracted from colon tissues of 3 mutant homozygous mice and 3 wild mice. RNA degradation, contamination concentration, and integrity were then measured. The extracted RNA was then sequenced using the Illumina platform. The obtained sequence data are filtered to ensure data quality and compared to the reference genome for further analysis. DESeq2 was used for gene expression analysis of the raw data. In addition, graphene oxide enrichment analysis and RT-qPCR validation were also performed.

RESULTS

This study identified 8354 differentially expressed genes in EDNRBm1yzcm and wild mouse colon tissues by RNA sequencing, including 4346 upregulated genes and 4005 downregulated genes. Correspondingly, the results of RT-qPCR analysis showed good correlation with the transcriptome data. In addition, GO and KEGG enrichment results suggested that there were 8103 terms and 320 pathways in all DEGs. When P < 0.05, 1081 GO terms and 320 KEGG pathways reached a significant level. Finally, through the existing studies and the enrichment results of differentially expressed genes, it was determined that axon guidance and the focal adhesion pathway may be closely related to the occurrence of HSCR.

CONCLUSIONS

This study analyzed and identified the differential genes in colonic tissues between EDNRBm1yzcm mice and wild mice, which provided new insight for further mining the potential pathogenic genes of Hirschsprung's disease.

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.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

RNA Sequencing of Urine-Derived Cells for the Characterization and Diagnosis of Osteogenesis Imperfecta

DNA sequencing is a reliable tool for identifying genetic variants in osteogenesis imperfecta (OI) but cannot always establish pathogenicity, particularly in variants altering splicing. RNA sequencing can provide functional evidence of the effect of a variant on the transcript but requires cells expressing the relevant genes. Here, we used urine-derived cells (UDC) to characterize genetic variants in patients with suspected or confirmed OI and provide evidence on the pathogenicity of variants of uncertain significance (VUS). Urine samples were obtained from 45 children and adolescents; UDC culture was successful in 40 of these participants (age range 4-20 years, 21 females), including 18 participants with OI or suspected OI who had a candidate variant or VUS on DNA sequencing. RNA was extracted from UDC and sequenced on an Illumina NextSeq550 device. Principal component analysis showed that the gene expression profiles of UDC and fibroblasts (based on Genotype Tissue Expression [GTEx] Consortium data) clustered close together and had less variability than those of whole blood cells. Transcript abundance was sufficient for analysis by RNA sequencing (defined as a median gene expression level of ≥10 transcripts per million) for 25 of the 32 bone fragility genes (78%) that were included in our diagnostic DNA sequencing panel. These results were similar to GTEx data for fibroblasts. Abnormal splicing was identified in 7 of the 8 participants with pathogenic or likely pathogenic variants in the splice region or deeper within the intron. Abnormal splicing was also observed in 2 VUS (COL1A1 c.2829+5G>A and COL1A2 c.693+6T>G), but no splice abnormality was observed in 3 other VUS. Abnormal deletions and duplications could also be observed in UDC transcripts. In conclusion, UDC are suitable for RNA transcript analysis in patients with suspected OI and can provide functional evidence for pathogenicity, in particular of variants affecting splicing. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Poison exon annotations improve the yield of clinically relevant variants in genomic diagnostic testing

PURPOSE

Neurodevelopmental disorders (NDDs) often result from rare genetic variation, but genomic testing yield for NDDs remains below 50%, suggesting that clinically relevant variants may be missed by standard analyses. Here, we analyze "poison exons" (PEs), which are evolutionarily conserved alternative exons often absent from standard gene annotations. Variants that alter PE inclusion can lead to loss of function and may be highly penetrant contributors to disease.

METHODS

We curated published RNA sequencing data from developing mouse cortex to define 1937 conserved PE regions potentially relevant to NDDs, and we analyzed variants found by genome sequencing in multiple NDD cohorts.

RESULTS

Across 2999 probands, we found 6 novel clinically relevant variants in PE regions. Five of these variants are in genes that are part of the sodium voltage-gated channel alpha subunit family (SCN1A, SCN2A, and SCN8A), which is associated with epilepsies. One variant is in SNRPB, associated with cerebrocostomandibular syndrome. These variants have moderate to high computational impact assessments, are absent from population variant databases, and in genes with gene-phenotype associations consistent with each probands reported features.

CONCLUSION

With a very minimal increase in variant analysis burden (average of 0.77 variants per proband), annotation of PEs can improve diagnostic yield for NDDs and likely other congenital conditions.

Transcriptomic analysis of paired healthy human skeletal muscles to identify modulators of disease severity in DMD

Muscle damage and fibro-fatty replacement of skeletal muscles is a main pathologic feature of Duchenne muscular dystrophy (DMD) with more proximal muscles affected earlier and more distal affected later in the disease course, suggesting that different skeletal muscle groups possess distinctive characteristics that influence their susceptibility to disease. To explore transcriptomic factors driving differential gene expression and modulating DMD skeletal muscle severity, we characterized the transcriptome of vastus lateralis (VL), a more proximal and susceptible muscle, relative to tibialis anterior (TA), a more distal and protected muscle, in 15 healthy individuals using bulk RNA sequencing to identify gene expression differences that may mediate their relative susceptibility to damage with loss of dystrophin. Matching single nuclei RNA sequencing data was generated for 3 of the healthy individuals, to infer cell composition in the bulk RNA sequencing dataset and to improve mapping of differentially expressed genes to their cell source of expression. A total of 3,410 differentially expressed genes were identified and mapped to cell type using single nuclei RNA sequencing of muscle, including long non-coding RNAs and protein coding genes. There was an enrichment of genes involved in calcium release from the sarcoplasmic reticulum, particularly in the myofibers and these myofiber genes were higher in the VL. There was an enrichment of genes in "Collagen-Containing Extracellular Matrix" expressed by fibroblasts, endothelial, smooth muscle and pericytes, with most genes higher in the TA, as well as genes in "Regulation Of Apoptotic Process" expressed across all cell types. Previously reported genetic modifiers were also enriched within the differentially expressed genes. We also identify 6 genes with differential isoform usage between the VL and TA. Lastly, we integrate our findings with DMD RNA sequencing data from the TA, and identify "Collagen-Containing Extracellular Matrix" and "Negative Regulation Of Apoptotic Process" as differentially expressed between DMD compared to healthy. Collectively, these findings propose novel candidate mechanisms that may mediate differential muscle susceptibility in muscular dystrophies and provide new insight into potential therapeutic targets.

LncRNA NEAT1: A novel regulator associated with the inflammatory response in acute respiratory distress syndrome

BACKGROUND

Acute respiratory distress syndrome (ARDS) is a life-threatening condition with an unfavorable prognosis. As the pathogenesis of ARDS remains unclear, we aimed to identify the core genes associated with ARDS and the mechanisms by which competing endogenous RNAs (ceRNAs) regulate the disease's progression.

METHODS

Three mRNA microarray datasets (GSE17355, GSE48787, and GSE130936), derived from the Gene Expression Omnibus (GEO) database, were selected. Common differentially expressed genes (DEGs) related to acute lung injury (ALI) were identified and subjected to enrichment analysis. Then, hub genes were figured out through the protein-protein interaction (PPI) network and functional analysis, and targeted miRNAs and lncRNAs were predicted. Finally, the ceRNA networks associated with ALI were constructed and validated experimentally.

RESULTS

A total of 155 upregulated and 93 downregulated DEGs were identified in the three datasets. The TNF signaling pathway and IL-17 signaling pathway were the most enriched pathways. Then, eleven DEGs enriched in the IL-17 signaling pathway were selected as the hub genes. Three miRNAs (mmu-mir-155-5p, mmu-mir-21a-5p, and mmu-mir-122-5p), which were located in the lung tissue and predicted to bind the hub genes at the same time, and two lncRNAs (Neat1 and Tug1), which have binding sites for the aforementioned miRNAs, were filtered. With qPCR verification, we identified a ceRNA network composed of NEAT1, miR-21-5p, MMP9, and CXCL5. NEAT1 knockdown promoted the migration and reduced the expression of pro-inflammatory factor and reactive oxygen species (ROS) in lung epithelial cells. We eventually confirmed that NEAT1/miR-21-5p/CXCL5/MMP9 played a pivotal role in regulating the inflammatory response in ALI.

CONCLUSION

The IL-17 signaling pathway is of great importance in the pathogenesis of ARDS. NEAT1/miR-21-5p is involved in the inflammation of ALI by regulating CXCL5 and MMP9.

RNA sequencing reveals a complete picture of a homozygous missense variant in a patient with VPS13D movement disorder: a case report and review of the literature

RNA sequencing (RNA-seq) is a complementary diagnostic tool to exome sequencing (ES), only recently clinically available to undiagnosed patients post-ES, that provides functional information on variants of unknown significance (VUS) by evaluating its effect on RNA transcription. ES became clinically available in the early 2010s and promised an agnostic platform for patients with a neurological disease, especially for those who believed to have a genetic etiology. However, the massive data generated by ES pose challenges in variant interpretation, especially for rare missense, synonymous, and deep intronic variants that may have a splicing effect. Without functional study and/or family segregation analysis, these rare variants would be likely interpreted as VUS which is difficult for clinicians to use in clinical care. Clinicians are able to assess the VUS for phenotypic overlap, but this additional information alone is usually not enough to re-classify a variant. Here, we report a case of a 14-month-old male who presented to clinic with a history of seizures, nystagmus, cerebral palsy, oral aversion, global developmental delay, and poor weight gain requiring gastric tube placement. ES revealed a previously unreported homozygous missense VUS, c.7406A > G p.(Asn2469Ser), in VPS13D. This variant has not been previously reported in genome aggregation database (gnomAD), ClinVar, or in any peer-reviewed published literature. By RNA-seq, we demonstrated that this variant mainly impacts splicing and results in a frameshift and early termination. It is expected to generate either a truncated protein, p.(Val2468fs*19), or no protein from this transcript due to nonsense-mediated mRNA decay leading to VPS13D deficiency. To our knowledge, this is the first case utilizing RNA-seq to further functionally characterize a homozygous novel missense VUS in VPS13D and confirm its impact on splicing. This confirmed pathogenicity gave the diagnosis of VPS13D movement disorder to this patient. Therefore, clinicians should consider utilizing RNA-seq to clarify VUS by evaluating its effect on RNA transcription.

Integrated multi-omics for rapid rare disease diagnosis on a national scale

Critically ill infants and children with rare diseases need equitable access to rapid and accurate diagnosis to direct clinical management. Over 2 years, the Acute Care Genomics program provided whole-genome sequencing to 290 families whose critically ill infants and children were admitted to hospitals throughout Australia with suspected genetic conditions. The average time to result was 2.9 d and diagnostic yield was 47%. We performed additional bioinformatic analyses and transcriptome sequencing in all patients who remained undiagnosed. Long-read sequencing and functional assays, ranging from clinically accredited enzyme analysis to bespoke quantitative proteomics, were deployed in selected cases. This resulted in an additional 19 diagnoses and an overall diagnostic yield of 54%. Diagnostic variants ranged from structural chromosomal abnormalities through to an intronic retrotransposon, disrupting splicing. Critical care management changed in 120 diagnosed patients (77%). This included major impacts, such as informing precision treatments, surgical and transplant decisions and palliation, in 94 patients (60%). Our results provide preliminary evidence of the clinical utility of integrating multi-omic approaches into mainstream diagnostic practice to fully realize the potential of rare disease genomic testing in a timely manner.

High diagnostic potential of short and long read genome sequencing with transcriptome analysis in exome-negative developmental disorders

Exome sequencing (ES) has become the method of choice for diagnosing rare diseases, while the availability of short-read genome sequencing (SR-GS) in a medical setting is increasing. In addition, new sequencing technologies, such as long-read genome sequencing (LR-GS) and transcriptome sequencing, are being increasingly used. However, the contribution of these techniques compared to widely used ES is not well established, particularly in regards to the analysis of non-coding regions. In a pilot study of five probands affected by an undiagnosed neurodevelopmental disorder, we performed trio-based short-read GS and long-read GS as well as case-only peripheral blood transcriptome sequencing. We identified three new genetic diagnoses, none of which affected the coding regions. More specifically, LR-GS identified a balanced inversion in NSD1, highlighting a rare mechanism of Sotos syndrome. SR-GS identified a homozygous deep intronic variant of KLHL7 resulting in a neoexon inclusion, and a de novo mosaic intronic 22-bp deletion in KMT2D, leading to the diagnosis of Perching and Kabuki syndromes, respectively. All three variants had a significant effect on the transcriptome, which showed decreased gene expression, mono-allelic expression and splicing defects, respectively, further validating the effect of these variants. Overall, in undiagnosed patients, the combination of short and long read GS allowed the detection of cryptic variations not or barely detectable by ES, making it a highly sensitive method at the cost of more complex bioinformatics approaches. Transcriptome sequencing is a valuable complement for the functional validation of variations, particularly in the non-coding genome.

A deep intronic TCTN2 variant activating a cryptic exon predicted by SpliceRover in a patient with Joubert syndrome

The recent introduction of genome sequencing in genetic analysis has led to the identification of pathogenic variants located in deep introns. Recently, several new tools have emerged to predict the impact of variants on splicing. Here, we present a Japanese boy of Joubert syndrome with biallelic TCTN2 variants. Exome sequencing identified only a heterozygous maternal nonsense TCTN2 variant (NM_024809.5:c.916C >T, p.(Gln306Ter)). Subsequent genome sequencing identified a deep intronic variant (c.1033+423G>A) inherited from his father. The machine learning algorithms SpliceAI, Squirls, and Pangolin were unable to predict alterations in splicing by the c.1033+423G>A variant. SpliceRover, a tool for splice site prediction using FASTA sequence, was able to detect a cryptic exon which was 85-bp away from the variant and within the inverted Alu sequence while SpliceRover scores for these splice sites showed slight increase (donor) or decrease (acceptor) between the reference and mutant sequences. RNA sequencing and RT-PCR using urinary cells confirmed inclusion of the cryptic exon. The patient showed major symptoms of TCTN2-related disorders such as developmental delay, dysmorphic facial features and polydactyly. He also showed uncommon features such as retinal dystrophy, exotropia, abnormal pattern of respiration, and periventricular heterotopia, confirming these as one of features of TCTN2-related disorders. Our study highlights usefulness of genome sequencing and RNA sequencing using urinary cells for molecular diagnosis of genetic disorders and suggests that database of cryptic splice sites predicted in introns by SpliceRover using the reference sequences can be helpful in extracting candidate variants from large numbers of intronic variants in genome sequencing.

Genetics of mitochondrial diseases: Current approaches for the molecular diagnosis

Mitochondrial diseases are a genetically and phenotypically variable set of monogenic disorders. The main characteristic of mitochondrial diseases is a defective oxidative phosphorylation. Both nuclear and mitochondrial DNA encode the approximately 1500 mitochondrial proteins. Since identification of the first mitochondrial disease gene in 1988 a total of 425 genes have been associated with mitochondrial diseases. Mitochondrial dysfunctions can be caused both by pathogenic variants in the mitochondrial DNA or the nuclear DNA. Hence, besides maternal inheritance, mitochondrial diseases can follow all modes of Mendelian inheritance. The maternal inheritance and tissue specificity distinguish molecular diagnostics of mitochondrial disorders from other rare disorders. With the advances made in the next-generation sequencing technology, whole exome sequencing and even whole-genome sequencing are now the established methods of choice for molecular diagnostics of mitochondrial diseases. They reach a diagnostic rate of more than 50% in clinically suspected mitochondrial disease patients. Moreover, next-generation sequencing is delivering a constantly growing number of novel mitochondrial disease genes. This chapter reviews mitochondrial and nuclear causes of mitochondrial diseases, molecular diagnostic methodologies, and their current challenges and perspectives.

Full-length isoform sequencing for resolving the molecular basis of Charcot-Marie-Tooth 2A

Objectives

Transcript sequencing of patient derived samples has been shown to improve the diagnostic yield for solving cases of likely Mendelian disorders, yet the added benefit of full-length long-read transcript sequencing is largely unexplored.

Methods

We applied short-read and full-length isoform cDNA sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis.

Results

We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts a branch point critical for intron 6 spicing. Full-length long-read isoform cDNA sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates five distinct altered splicing transcripts. All five altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 protein levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A).

Discussion

This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster

EZH1 ( Enhancer of Zeste, homolog 1) , a Polycomb Repressive Complex-2 (PRC2) component, is involved in a myriad of cellular processes through modifying histone 3 lysine27 (H3K27) residues. EZH1 represses transcription of downstream target genes through H3K27 trimethylation (H3K27me3). Genetic mutations in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo variant in EZH1 , p.Ala678Gly, through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and was the best candidate variant found in the exome. Human EZH1 / 2 are homologous to fly Enhancer of zeste E(z) , an essential gene in flies, and the residue (A678 in humans, A691 in Drosophila ) is conserved. To further study this variant, we obtained Drosophila null alleles and generated transgenic flies expressing wild-type (E(z) WT ) and the variant (E(z) A691G ) . The E(z) A691G variant led to hyper H3K27me3 while the E(z) WT did not, suggesting this is as a gain-of-function allele. When expressed under the tubulin promotor in vivo the variant rescued null-lethality similar to wild-type but the E(z) A691G flies exhibit bang sensitivity and shortened lifespan. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila . Biochemically this allele leads to increased H3K27me3 suggesting gain-of-function, but when expressed in adult flies the E(z) A691G has some characteristics of partial loss-of-function which may suggest it is a more complex allele in vivo .

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.

Beyond the exome: what's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order and emerging technologies, such as optical genome mapping and long-read DNA or RNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to a consortium such as GREGoR, which is focused on elucidating the underlying cause of rare unsolved genetic disorders.

Poison exon annotations improve the yield of clinically relevant variants in genomic diagnostic testing

Purpose

Neurodevelopmental disorders (NDDs) often result from rare genetic variation, but genomic testing yield for NDDs remains around 50%, suggesting some clinically relevant rare variants may be missed by standard analyses. Here we analyze "poison exons" (PEs) which, while often absent from standard gene annotations, are alternative exons whose inclusion results in a premature termination codon. Variants that alter PE inclusion can lead to loss-of-function and may be highly penetrant contributors to disease.

Methods

We curated published RNA-seq data from developing mouse cortex to define 1,937 PE regions conserved between humans and mice and potentially relevant to NDDs. We then analyzed variants found by genome sequencing in multiple NDD cohorts.

Results

Across 2,999 probands, we found six clinically relevant variants in PE regions that were previously overlooked. Five of these variants are in genes that are part of the sodium voltage-gated channel alpha subunit family ( SCN1A, SCN2A , and SCN8A ), associated with epilepsies. One variant is in SNRPB , associated with Cerebrocostomandibular Syndrome. These variants have moderate to high computational impact assessments, are absent from population variant databases, and were observed in probands with features consistent with those reported for the associated gene.

Conclusion

With only a minimal increase in variant analysis burden (most probands had zero or one candidate PE variants in a known NDD gene, with an average of 0.77 per proband), annotation of PEs can improve diagnostic yield for NDDs and likely other congenital conditions.

RNA Sequencing as a Diagnostic Tool

This JAMA Insights discusses RNA sequencing, which allows detection of qualitative and quantitative changes in RNA expression across the genome in clinical samples, and how it may improve molecular diagnostic rates achieved by diagnostic exome sequencing and whole-genome sequencing alone.

Stepwise use of genomics and transcriptomics technologies increases diagnostic yield in Mendelian disorders

Purpose: Multi-omics offer worthwhile and increasingly accessible technologies to diagnostic laboratories seeking potential second-tier strategies to help patients with unresolved rare diseases, especially patients clinically diagnosed with a rare OMIM (Online Mendelian Inheritance in Man) disease. However, no consensus exists regarding the optimal diagnostic care pathway to adopt after negative results with standard approaches. Methods: In 15 unsolved individuals clinically diagnosed with recognizable OMIM diseases but with negative or inconclusive first-line genetic results, we explored the utility of a multi-step approach using several novel omics technologies to establish a molecular diagnosis. Inclusion criteria included a clinical autosomal recessive disease diagnosis and single heterozygous pathogenic variant in the gene of interest identified by first-line analysis (60%-9/15) or a clinical diagnosis of an X-linked recessive or autosomal dominant disease with no causative variant identified (40%-6/15). We performed a multi-step analysis involving short-read genome sequencing (srGS) and complementary approaches such as mRNA sequencing (mRNA-seq), long-read genome sequencing (lrG), or optical genome mapping (oGM) selected according to the outcome of the GS analysis. Results: SrGS alone or in combination with additional genomic and/or transcriptomic technologies allowed us to resolve 87% of individuals by identifying single nucleotide variants/indels missed by first-line targeted tests, identifying variants affecting transcription, or structural variants sometimes requiring lrGS or oGM for their characterization. Conclusion: Hypothesis-driven implementation of combined omics technologies is particularly effective in identifying molecular etiologies. In this study, we detail our experience of the implementation of genomics and transcriptomics technologies in a pilot cohort of previously investigated patients with a typical clinical diagnosis without molecular etiology.

Diagnostic potential of the amniotic fluid cells transcriptome in deciphering mendelian disease: a proof-of-concept

RNA sequencing (RNA-seq) is emerging in genetic diagnoses as it provides functional support for the interpretation of variants of uncertain significance. However, the use of amniotic fluid (AF) cells for RNA-seq has not yet been explored. Here, we examined the expression of clinically relevant genes in AF cells (n = 48) compared with whole blood and fibroblasts. The number of well-expressed genes in AF cells was comparable to that in fibroblasts and much higher than that in blood across different disease categories. We found AF cells RNA-seq feasible and beneficial in prenatal diagnosis (n = 4) as transcriptomic data elucidated the molecular consequence leading to the pathogenicity upgrade of variants in CHD7 and COL1A2 and revising the in silico prediction of a variant in MYRF. AF cells RNA-seq could become a reasonable choice for postnatal patients with advantages over fibroblasts and blood as it prevents invasive procedures.

A recurrent single-exon deletion in TBCK might be under-recognized in patients with infantile hypotonia and psychomotor delay

Advanced bioinformatics algorithms allow detection of multiple-exon copy-number variations (CNVs) from exome sequencing (ES) data, while detection of single-exon CNVs remains challenging. A retrospective review of Baylor Genetics' clinical ES patient cohort identified four individuals with homozygous single-exon deletions of TBCK (exon 23, NM_001163435.2), a gene associated with an autosomal recessive neurodevelopmental phenotype. To evaluate the prevalence of this deletion and its contribution to disease, we retrospectively analyzed single nucleotide polymorphism (SNP) array data for 8194 individuals undergoing ES, followed by PCR confirmation and RT-PCR on individuals carrying homozygous or heterozygous exon 23 TBCK deletions. A fifth individual was diagnosed with the TBCK-related disorder due to a heterozygous exon 23 deletion in trans with a c.1860+1G>A (NM_001163435.2) pathogenic variant, and three additional heterozygous carriers were identified. Affected individuals and carriers were from diverse ethnicities including European Caucasian, South Asian, Middle Eastern, Hispanic American and African American, with only one family reporting consanguinity. RT-PCR revealed two out-of-frame transcripts related to the exon 23 deletion. Our results highlight the importance of identifying single-exon deletions in clinical ES, especially for genes carrying recurrent deletions. For patients with early-onset hypotonia and psychomotor delay, this single-exon TBCK deletion might be under-recognized due to technical limitations of ES.

Deep intronic NIPBL de novo mutations and differential diagnoses revealed by whole genome and RNA sequencing in Cornelia de Lange syndrome patients

Cornelia de Lange syndrome (CdLS; MIM# 122470) is a rare developmental disorder. Pathogenic variants in 5 genes explain approximately 50% cases, leaving the other 50% unsolved. We performed whole genome sequencing (WGS) ± RNA sequencing (RNA-seq) in 5 unsolved trios fulfilling the following criteria: (i) clinical diagnosis of classic CdLS, (ii) negative gene panel sequencing from blood and saliva-isolated DNA, (iii) unaffected parents' DNA samples available and (iv) proband's blood-isolated RNA available. A pathogenic de novo mutation (DNM) was observed in a CdLS differential diagnosis gene in 3/5 patients, namely POU3F3, SPEN, and TAF1. In the other two, we identified two distinct deep intronic DNM in NIPBL predicted to create a novel splice site. RT-PCRs and RNA-Seq showed aberrant transcripts leading to the creation of a novel frameshift exon. Our findings suggest the relevance of WGS in unsolved suspected CdLS cases and that deep intronic variants may account for a proportion of them.

Evaluating Genetic Disorders in the Neonate: The Role of Exome Sequencing in the NICU

With recent advances in the technologies used for genetic diagnosis as well as our understanding of the genetic basis of disease, a growing list of options is available for providers when caring for a newborn with features suggesting an underlying genetic etiology. The choice of the most appropriate genetic test for a specific situation includes clinical considerations such as the phenotypic features and type of genetic abnormality suspected, as well as practical considerations such as cost and turnaround time. In this review, we discuss clinical exome sequencing in the context of genetic evaluation of newborns, including technical considerations, variant interpretation, and incidental/secondary findings. Strengths and limitations of exome sequencing are discussed and compared with those of other commonly known tests such as karyotype analysis, fluorescence in situ hybridization, chromosomal microarray, and sequencing panels, along with integration of results from prenatal testing if available. We also review future directions including genome sequencing and other emerging technologies that are starting to be used in clinical settings.

Proteomics in Inherited Metabolic Disorders

Inherited metabolic disorders (IMD) are rare medical conditions caused by genetic defects that interfere with the body's metabolism. The clinical phenotype is highly variable and can present at any age, although it more often manifests in childhood. The number of treatable IMDs has increased in recent years, making early diagnosis and a better understanding of the natural history of the disease more important than ever. In this review, we discuss the main challenges faced in applying proteomics to the study of IMDs, and the key advances achieved in this field using tandem mass spectrometry (MS/MS). This technology enables the analysis of large numbers of proteins in different body fluids (serum, plasma, urine, saliva, tears) with a single analysis of each sample, and can even be applied to dried samples. MS/MS has thus emerged as the tool of choice for proteome characterization and has provided new insights into many diseases and biological systems. In the last 10 years, sequential window acquisition of all theoretical fragmentation spectra mass spectrometry (SWATH-MS) has emerged as an accurate, high-resolution technique for the identification and quantification of proteins differentially expressed between healthy controls and IMD patients. Proteomics is a particularly promising approach to help obtain more information on rare genetic diseases, including identification of biomarkers to aid early diagnosis and better understanding of the underlying pathophysiology to guide the development of new therapies. Here, we summarize new and emerging proteomic technologies and discuss current uses and limitations of this approach to identify and quantify proteins. Moreover, we describe the use of proteomics to identify the mechanisms regulating complex IMD phenotypes; an area of research essential to better understand these rare disorders and many other human diseases.