Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Geneapp: A web application for visualizing alternative splicing for biomedicine

Alternative Splicing (AS) is an essential mechanism for eukaryotes. However, the consequences of deleting a single exon can be dramatic for the organism and can lead to cancer in humans. Additionally, alternative 5' and 3' splice sites, which define the boundaries of exons, also play key roles to human disorders. Therefore, Investigating AS events is crucial for understanding the molecular basis of human diseases and developing therapeutic strategies. Workflow for AS event analysis can be sampling followed by data analysis with bioinformatics to identify the different AS events in the control and case samples, data visualization for curation, and selection of relevant targets for experimental validation. The raw output of the analysis software does not favor the inspection of events by bioinformaticians requiring custom scripts for data visualization. In this work, we propose the Geneapp application with three modules: GeneappScript, GeneappServer, and GeneappExplorer. GeneappScript is a wrapper that assists in identifying AS in samples compared in two different approaches, while GeneappServer integrates data from AS analysis already performed by the user. In GeneappExplorer, the user visualizes the previous dataset by exploring AS events in genes with functional annotation. This targeted screens that Geneapp allows to perform helps in the identification of targets for experimental validation to confirm the hypotheses under study. The Geneapp is freely available for non-commercial use at https://geneapp.net to advance research on AS for bioinformatics.

Efficient, specific, and combinatorial control of endogenous exon splicing with dCasRx-RBM25

Efficient targeted control of splicing is a major goal of functional genomics and therapeutic applications. Guide (g)RNA-directed, deactivated (d)Cas CRISPR enzymes fused to splicing effectors represent a promising strategy due to the flexibility of these systems. However, efficient, specific, and generalizable activation of endogenous exons using this approach has not been previously reported. By screening over 300 dCasRx-splicing factor fusion proteins tethered to splicing reporters, we identify dCasRx-RBM25 as a potent activator of exons. Moreover, dCasRx-RBM25 efficiently activates the splicing of ∼90% of targeted endogenous alternative exons and displays high on-target specificity. Using gRNA arrays for combinatorial targeting, we demonstrate that dCasRx-RBM25 enables multiplexed activation and repression of exons. Using this feature, the targeting of neural-regulated exons in Ptpb1 and Puf60 in embryonic stem cells reveals combinatorial effects on downstream alternative splicing events controlled by these factors. Collectively, our results enable versatile, combinatorial exon-resolution functional assays and splicing-directed therapeutic applications.

Genomic Mosaicism of the Brain: Origin, Impact, and Utility

Genomic mosaicism describes the phenomenon where some but not all cells within a tissue harbor unique genetic mutations. Traditionally, research focused on the impact of genomic mosaicism on clinical phenotype-motivated by its involvement in cancers and overgrowth syndromes. More recently, we increasingly shifted towards the plethora of neutral mosaic variants that can act as recorders of cellular lineage and environmental exposures. Here, we summarize the current state of the field of genomic mosaicism research with a special emphasis on our current understanding of this phenomenon in brain development and homeostasis. Although the field of genomic mosaicism has a rich history, technological advances in the last decade have changed our approaches and greatly improved our knowledge. We will provide current definitions and an overview of contemporary detection approaches for genomic mosaicism. Finally, we will discuss the impact and utility of genomic mosaicism.

A deep intronic variant in MME causes autosomal recessive Charcot-Marie-Tooth neuropathy through aberrant splicing

BACKGROUND

Loss-of-function variants in MME (membrane metalloendopeptidase) are a known cause of recessive Charcot-Marie-Tooth Neuropathy (CMT). A deep intronic variant, MME c.1188+428A>G (NM_000902.5), was identified through whole genome sequencing (WGS) of two Australian families with recessive inheritance of axonal CMT using the seqr platform. MME c.1188+428A>G was detected in a homozygous state in Family 1, and in a compound heterozygous state with a known pathogenic MME variant (c.467del; p.Pro156Leufs*14) in Family 2.

AIMS

We aimed to determine the pathogenicity of the MME c.1188+428A>G variant through segregation and splicing analysis.

METHODS

The splicing impact of the deep intronic MME variant c.1188+428A>G was assessed using an in vitro exon-trapping assay.

RESULTS

The exon-trapping assay demonstrated that the MME c.1188+428A>G variant created a novel splice donor site resulting in the inclusion of an 83 bp pseudoexon between MME exons 12 and 13. The incorporation of the pseudoexon into MME transcript is predicted to lead to a coding frameshift and premature termination codon (PTC) in MME exon 14 (p.Ala397ProfsTer47). This PTC is likely to result in nonsense mediated decay (NMD) of MME transcript leading to a pathogenic loss-of-function.

INTERPRETATION

To our knowledge, this is the first report of a pathogenic deep intronic MME variant causing CMT. This is of significance as deep intronic variants are missed using whole exome sequencing screening methods. Individuals with CMT should be reassessed for deep intronic variants, with splicing impacts being considered in relation to the potential pathogenicity of variants.

DNA methylation analysis in patients with neurodevelopmental disorders improves variant interpretation and reveals complexity

Analysis of genomic DNA methylation by generating epigenetic signature profiles (episignatures) is increasingly being implemented in genetic diagnosis. Here we report our experience using episignature analysis to resolve both uncomplicated and complex cases of neurodevelopmental disorders (NDDs). We analyzed 97 NDDs divided into (1) a validation cohort of 59 patients with likely pathogenic/pathogenic variants characterized by a known episignature and (2) a test cohort of 38 patients harboring variants of unknown significance or unidentified variants. The expected episignature was obtained in most cases with likely pathogenic/pathogenic variants (53/59 [90%]), a revealing exception being the overlapping profile of two SMARCB1 pathogenic variants with ARID1A/B:c.6200, confirmed by the overlapping clinical features. In the test cohort, five cases showed the expected episignature, including (1) novel pathogenic variants in ARID1B and BRWD3; (2) a deletion in ATRX causing MRXFH1 X-linked mental retardation; and (3) confirmed the clinical diagnosis of Cornelia de Lange (CdL) syndrome in mutation-negative CdL patients. Episignatures analysis of the in BAF complex components revealed novel functional protein interactions and common episignatures affecting homologous residues in highly conserved paralogous proteins (SMARCA2 M856V and SMARCA4 M866V). Finally, we also found sex-dependent episignatures in X-linked disorders. Implementation of episignature profiling is still in its early days, but with increasing utilization comes increasing awareness of the capacity of this methodology to help resolve the complex challenges of genetic diagnoses.

Bioinformatics pipeline for the systematic mining genomic and proteomic variation linked to rare diseases: The example of monogenic diabetes

Monogenic diabetes is characterized as a group of diseases caused by rare variants in single genes. Like for other rare diseases, multiple genes have been linked to monogenic diabetes with different measures of pathogenicity, but the information on the genes and variants is not unified among different resources, making it challenging to process them informatically. We have developed an automated pipeline for collecting and harmonizing data on genetic variants linked to monogenic diabetes. Furthermore, we have translated variant genetic sequences into protein sequences accounting for all protein isoforms and their variants. This allows researchers to consolidate information on variant genes and proteins linked to monogenic diabetes and facilitates their study using proteomics or structural biology. Our open and flexible implementation using Jupyter notebooks enables tailoring and modifying the pipeline and its application to other rare diseases.

Alternative splicing of pre-mRNA modulates the immune response in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Little is known about the role of alternative splicing (AS) in regulating gene expression in Mycobacteria-infected individuals in distinct stages of infection. Pre-mRNA AS consists of the removal of introns and the assembly of exons contained in eukaryotic genes. AS events can influence transcript stability or structure with important physiological consequences. Using RNA-Seq data from peripheral blood (PB) and ileocecal valve (ICV) samples collected from Holstein cattle with focal and diffuse paratuberculosis (PTB)-associated histopathological lesions in gut tissues and without lesions (controls), we detected differential AS profiles between the infected and control groups. Four of the identified AS events were experimentally validated by reverse transcription-digital droplet PCR (RT-ddPCR). AS events in several genes correlated with changes in gene expression. In the ICV of animals with diffuse lesions, for instance, alternatively spliced genes correlated with changes in the expression of genes involved in endocytosis, antigen processing and presentation, complement activation, and several inflammatory and autoimmune diseases in humans. Taken together, our results identified common mechanisms of AS involvement in the pathogenesis of PTB and human diseases and shed light on novel diagnostic and therapeutic interventions to control these diseases.

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants

BACKGROUND

Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis.

RESULTS

Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity.

CONCLUSIONS

By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

Whole genome sequencing in clinical practice

Whole genome sequencing (WGS) is becoming the preferred method for molecular genetic diagnosis of rare and unknown diseases and for identification of actionable cancer drivers. Compared to other molecular genetic methods, WGS captures most genomic variation and eliminates the need for sequential genetic testing. Whereas, the laboratory requirements are similar to conventional molecular genetics, the amount of data is large and WGS requires a comprehensive computational and storage infrastructure in order to facilitate data processing within a clinically relevant timeframe. The output of a single WGS analyses is roughly 5 MIO variants and data interpretation involves specialized staff collaborating with the clinical specialists in order to provide standard of care reports. Although the field is continuously refining the standards for variant classification, there are still unresolved issues associated with the clinical application. The review provides an overview of WGS in clinical practice - describing the technology and current applications as well as challenges connected with data processing, interpretation and clinical reporting.

Predicting the impact of rare variants on RNA splicing in CAGI6

Variants which disrupt splicing are a frequent cause of rare disease that have been under-ascertained clinically. Accurate and efficient methods to predict a variant's impact on splicing are needed to interpret the growing number of variants of unknown significance (VUS) identified by exome and genome sequencing. Here, we present the results of the CAGI6 Splicing VUS challenge, which invited predictions of the splicing impact of 56 variants ascertained clinically and functionally validated to determine splicing impact. The performance of 12 prediction methods, along with SpliceAI and CADD, was compared on the 56 functionally validated variants. The maximum accuracy achieved was 82% from two different approaches, one weighting SpliceAI scores by minor allele frequency, and one applying the recently published Splicing Prediction Pipeline (SPiP). SPiP performed optimally in terms of sensitivity, while an ensemble method combining multiple prediction tools and information from databases exceeded all others for specificity. Several challenge methods equalled or exceeded the performance of SpliceAI, with ultimate choice of prediction method likely to depend on experimental or clinical aims. One quarter of the variants were incorrectly predicted by at least 50% of the methods, highlighting the need for further improvements to splicing prediction methods for successful clinical application.

Prioritizing de novo potential non-canonical splicing variants in neurodevelopmental disorders

BACKGROUND

Genomic variants outside of the canonical splicing site (±2) may generate abnormal mRNA splicing, which are defined as non-canonical splicing variants (NCSVs). However, the clinical interpretation of NCSVs in neurodevelopmental disorders (NDDs) is largely unknown.

METHODS

We investigated the contribution of NCSVs to NDDs from 345,787 de novo variants (DNVs) in 47,574 patients with NDDs. We performed functional enrichment and protein-protein interaction analysis to assess the association between genes carrying prioritised NCSVs and NDDs. Minigene was used to validate the impact of NCSVs on mRNA splicing.

FINDINGS

We observed significantly more NCSVs (p = 0.02, odds ratio [OR] = 2.05) among patients with NDD than in controls. Both canonical splicing variants (CSVs) and NCSVs contributed to an equal proportion of patients with NDD (0.76% vs. 0.82%). The candidate genes carrying NCSVs were associated with glutamatergic synapse and chromatin remodelling. Minigene successfully validated 59 of 79 (74.68%) NCSVs that led to abnormal splicing in 40 candidate genes, and 9 of the genes (ARID1B, KAT6B, TCF4, SMARCA2, SHANK3, PDHA1, WDR45, SCN2A, SYNGAP1) harboured recurrent NCSVs with the same variant present in more than two unrelated patients with NDD. Moreover, 36 of 59 (61.02%) NCSVs are novel clinically relevant variants, including 34 unreported and 2 clinically conflicting interpretations or of uncertain significance NCSVs in the ClinVar database.

INTERPRETATION

This study highlights the common pathology and clinical importance of NCSVs in unsolved patients with NDD.

FUNDING

The present study was funded by grants from the National Natural Science Foundation of China, China Postdoctoral Science Foundation, the Hunan Youth Science and Technology Innovation Talent Project, the Provincial Natural Science Foundation of Hunan, The Scientific Research Program of FuRong laboratory, and the Natural Science Project of the University of Anhui Province.

Cohort analysis of novel SPAST variants in SPG4 patients and implementation of in vitro and in vivo studies to identify the pathogenic mechanism caused by splicing mutations

Introduction

Pure hereditary spastic paraplegia (SPG) type 4 (SPG4) is caused by mutations of SPAST gene. This study aimed to analyze SPAST variants in SPG4 patients to highlight the occurrence of splicing mutations and combine functional studies to assess the relevance of these variants in the molecular mechanisms of the disease.

Methods

We performed an NGS panel in 105 patients, in silico analysis for splicing mutations, and in vitro minigene assay.

Results and discussion

The NGS panel was applied to screen 105 patients carrying a clinical phenotype corresponding to upper motor neuron syndrome (UMNS), selectively affecting motor control of lower limbs. Pathogenic mutations in SPAST were identified in 12 patients (11.42%), 5 missense, 3 frameshift, and 4 splicing variants. Then, we focused on the patients carrying splicing variants using a combined approach of in silico and in vitro analysis through minigene assay and RNA, if available. For two splicing variants (i.e., c.1245+1G>A and c.1414-2A>T), functional assays confirm the types of molecular alterations suggested by the in silico analysis (loss of exon 9 and exon 12). In contrast, the splicing variant c.1005-1delG differed from what was predicted (skipping exon 7), and the functional study indicates the loss of frame and formation of a premature stop codon. The present study evidenced the high splice variants in SPG4 patients and indicated the relevance of functional assays added to in silico analysis to decipher the pathogenic mechanism.

Clinical features and mutational analysis of X-linked agammaglobulinemia patients in Malaysia

Background

Bruton's tyrosine kinase (BTK) is a cytoplasmic protein involved in the B cell development. X-linked agammaglobulinemia (XLA) is caused by mutation in the BTK gene, which results in very low or absent B cells. Affected males have markedly reduced immunoglobulin levels, which render them susceptible to recurrent and severe bacterial infections. Methods: Patients suspected with X-linked agammaglobulinemia were enrolled during the period of 2010-2018. Clinical summary, and immunological profiles of these patients were recorded. Peripheral blood samples were collected for monocyte BTK protein expression detection and BTK genetic analysis. The medical records between January 2020 and June 2023 were reviewed to investigate COVID-19 in XLA.

Results

Twenty-two patients (from 16 unrelated families) were molecularly diagnosed as XLA. Genetic testing revealed fifteen distinct mutations, including four splicing mutations, four missense mutations, three nonsense mutations, three short deletions, and one large indel mutation. These mutations scattered throughout the BTK gene and mostly affected the kinase domain. All mutations including five novel mutations were predicted to be pathogenic or deleterious by in silico prediction tools. Genetic testing confirmed that eleven mothers and seven sisters were carriers for the disease, while three mutations were de novo. Flow cytometric analysis showed that thirteen patients had minimal BTK expression (0-15%) while eight patients had reduced BTK expression (16-64%). One patient was not tested for monocyte BTK expression due to insufficient sample. Pneumonia (n=13) was the most common manifestation, while Pseudomonas aeruginosa was the most frequently isolated pathogen from the patients (n=4). Mild or asymptomatic COVID-19 was reported in four patients.

Conclusion

This report provides the first overview of demographic, clinical, immunological and genetic data of XLA in Malaysia. The combination of flow cytometric assessment and BTK genetic analysis provides a definitive diagnosis for XLA patients, especially with atypical clinical presentation. In addition, it may also allow carrier detection and assist in genetic counselling and prenatal diagnosis.

Unsuspected consequences of synonymous and missense variants in OCA2 can be detected in blood cell RNA samples of patients with albinism

Oculocutaneous albinism type 2 (OCA2) is the second most frequent form of albinism and represents about 30% of OCA worldwide. As with all types of OCA, patients present with hypopigmentation of hair and skin, as well as severe visual abnormalities. We focused on a subgroup of 29 patients for whom genetic diagnosis was pending because at least one of their identified variants in or around exon 10 of OCA2 is of uncertain significance (VUS). By minigene assay, we investigated the effect of these VUS on exon 10 skipping and showed that not only intronic but also some synonymous variants can result in enhanced exon skipping. We further found that excessive skipping of exon 10 could be detected directly on blood samples of patients and of their one parent with the causal variant, avoiding invasive skin biopsies. Moreover, we show that variants, which result in lack of detectable OCA2 mRNA can be identified from blood samples as well, as shown for the most common OCA2 pathogenic missense variant c.1327G>A/p.(Val443Ile). In conclusion, blood cell RNA analysis allows testing the potential effect of any OCA2 VUS on transcription products. This should help to elucidate yet unsolved OCA2 patients and improve genetic counseling.

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.

Identification of a Novel de novo Splicing Mutation in Duchenne Muscular Dystrophy Gene in an Iranian Family

Introduction

Duchenne muscular dystrophy (DMD) (NM_004006.3) is one of the most notable neuromuscular disorders of early years. The majority of DMD cases are caused by deletions or duplications in dystrophin, while point mutations are less prevalent in dystrophin abnormalities. It is a common knowledge that the severity of the disease depends on the effect of the mutation on the translational reading frame of the dystrophin mRNA.

Case Report

We studied an 8-year-old boy with relevant clinical presentations for DMD. Deletion/duplication screening was performed by using multiplex ligation-dependent probe amplification, and whole-exome sequencing was conducted in order to identify potential variants. A novel de novo splice site variant was identified in the DMD gene (DMD: c.8548-2A>G). To explore the effect of a novel variant in DMD, various in silico analyses were carried out to investigate the pathogenicity of the causative variant. To study the structure of a DMD protein and information on how the genetic variant impacts splicing site in models of wild-type and mutated DMD, we carried out different computational studies. Sanger sequencing was performed for the purpose of variant confirmation and familial segregation analysis.

Discussion

This novel de novo variant was predicted to have an effect on splicing, which leads to DMD due to its significant impacts on dystrophin functionality. The novel mutation would be expected to disrupt the protein structure.

Classification of PRSS1 variants responsible for chronic pancreatitis: An expert perspective from the Franco-Chinese GREPAN study group

BACKGROUND

PRSS1 was the first reported chronic pancreatitis (CP) gene. The existence of both gain-of-function (GoF) and gain-of-proteotoxicity (GoP) pathological PRSS1 variants, together with the fact that PRSS1 variants have been identified in CP subtypes spanning the range from monogenic to multifactorial, has made the classification of PRSS1 variants very challenging.

METHODS

All currently reported PRSS1 variants (derived primarily from two databases) were manually reviewed with respect to their clinical genetics, functional analysis and population allele frequency. They were classified by variant type and pathological mechanism within the framework of our recently proposed ACMG/AMP guidelines-based seven-category system.

RESULTS

The total number of distinct germline PRSS1 variants included for analysis was 100, comprising 3 copy number variants (CNVs), 12 5' and 3' variants, 19 intronic variants, 5 nonsense variants, 1 frameshift deletion variant, 6 synonymous variants, 1 in-frame duplication, 3 gene conversions and 50 missense variants. Based upon a combination of clinical genetic and functional analysis, population data and in silico analysis, we classified 26 variants (all 3 CNVs, the in-frame duplication, all 3 gene conversions and 19 missense) as "pathogenic", 3 variants (missense) as "likely pathogenic", 5 variants (four missense and one promoter) as "predisposing", 13 variants (all missense) as "unknown significance", 2 variants (missense) as "likely benign", and all remaining 51 variants as "benign".

CONCLUSIONS

We describe an expert classification of the 100 PRSS1 variants reported to date. The results have immediate implications for reclassifying many ClinVar-registered PRSS1 variants as well as providing optimal guidelines/standards for reporting PRSS1 variants.

Impaired AIF-CHCHD4 interaction and mitochondrial calcium overload contribute to auditory neuropathy spectrum disorder in patient-iPSC-derived neurons with AIFM1 variant

Auditory neuropathy spectrum disorder (ANSD) is a hearing impairment caused by dysfunction of inner hair cells, ribbon synapses, spiral ganglion neurons and/or the auditory nerve itself. Approximately 1/7000 newborns have abnormal auditory nerve function, accounting for 10%-14% of cases of permanent hearing loss in children. Although we previously identified the AIFM1 c.1265 G > A variant to be associated with ANSD, the mechanism by which ANSD is associated with AIFM1 is poorly understood. We generated induced pluripotent stem cells (iPSCs) from peripheral blood mononuclear cells (PBMCs) via nucleofection with episomal plasmids. The patient-specific iPSCs were edited via CRISPR/Cas9 technology to generate gene-corrected isogenic iPSCs. These iPSCs were further differentiated into neurons via neural stem cells (NSCs). The pathogenic mechanism was explored in these neurons. In patient cells (PBMCs, iPSCs, and neurons), the AIFM1 c.1265 G > A variant caused a novel splicing variant (c.1267-1305del), resulting in AIF p.R422Q and p.423-435del proteins, which impaired AIF dimerization. Such impaired AIF dimerization then weakened the interaction between AIF and coiled-coil-helix-coiled-coil-helix domain-containing protein 4 (CHCHD4). On the one hand, the mitochondrial import of ETC complex subunits was inhibited, subsequently leading to an increased ADP/ATP ratio and elevated ROS levels. On the other hand, MICU1-MICU2 heterodimerization was impaired, leading to _mCa²⁺ overload. Calpain was activated by _mCa²⁺ and subsequently cleaved AIF for its translocation into the nucleus, ultimately resulting in caspase-independent apoptosis. Interestingly, correction of the AIFM1 variant significantly restored the structure and function of AIF, further improving the physiological state of patient-specific iPSC-derived neurons. This study demonstrates that the AIFM1 variant is one of the molecular bases of ANSD. Mitochondrial dysfunction, especially _mCa²⁺ overload, plays a prominent role in ANSD associated with AIFM1. Our findings help elucidate the mechanism of ANSD and may lead to the provision of novel therapies.

Loss-of-function mutations in MYO15A and OTOF cause non-syndromic hearing loss in two Yemeni families

BACKGROUND

Hearing loss is a rare hereditary deficit that is rather common among consanguineous populations. Autosomal recessive non-syndromic hearing loss is the predominant form of hearing loss worldwide. Although prevalent, hearing loss is extremely heterogeneous and poses a pitfall in terms of diagnosis and screening. Using next-generation sequencing has enabled a rapid increase in the identification rate of genes and variants in heterogeneous conditions, including hearing loss. We aimed to identify the causative variants in two consanguineous Yemeni families affected with hearing loss using targeted next-generation sequencing (clinical exome sequencing). The proband of each family was presented with sensorineural hearing loss as indicated by pure-tone audiometry results.

RESULTS

We explored variants obtained from both families, and our analyses collectively revealed the presence and segregation of two novel loss-of-function variants: a frameshift variant, c.6347delA in MYO15A in Family I, and a splice site variant, c.5292-2A > C, in OTOF in Family II. Sanger sequencing and PCR-RFLP of DNA samples from 130 deaf and 50 control individuals confirmed that neither variant was present in our in-house database. In silico analyses predicted that each variant has a pathogenic effect on the corresponding protein.

CONCLUSIONS

We describe two novel loss-of-function variants in MYO15A and OTOF that cause autosomal recessive non-syndromic hearing loss in Yemeni families. Our findings are consistent with previously reported pathogenic variants in the MYO15A and OTOF genes in Middle Eastern individuals and suggest their implication in hearing loss.

Improving Hereditary Hemorrhagic Telangiectasia Molecular Diagnosis: A Referral Center Experience

BACKGROUND

Hereditary hemorrhagic telangiectasia (HHT) is a rare vascular disease inherited in an autosomal dominant manner. Disease-causing variants in endoglin (ENG) and activin A receptor type II-like 1 (ACVRL1) genes are detected in more than 90% of the patients undergoing molecular testing. The identification of variants of unknown significance is often seen as a challenge in clinical practice that makes family screening and genetic counseling difficult. Here, we show that the implementation of cDNA analysis to assess the effect of splice site variants on mRNA splicing is a powerful tool.

METHODS

Gene panel sequencing of genes associated with HHT and other arteriovenous malformation-related syndromes was performed. To evaluate the effect of the splice site variants, cDNA analysis of ENG and ACVRL1 genes was carried out.

RESULTS

three novel splice site variants were identified in ENG (c.68-2A > T and c.1311+4_1311+8del) and ACVLR1 (c.526-6C > G) genes correspondingly in three individuals with HHT that met ≥ 3 Curaçao criteria. All three variants led to an aberrant splicing inducing exon skipping (ENG:c.68-2A > T and ACVRL1:c.526-6C > G) or intron retention (ENG:c.1311+4_1311+8del) allowing the confirmation of the predicted effect on splicing and the reclassification from unknown significance to pathogenic/likely pathogenic of two of them.

CONCLUSIONS

RNA analysis should be performed to assess and/or confirm the impact of variants on splicing. The molecular diagnosis of HHT patients is crucial to allow family screening and accurate genetic counseling. A multidisciplinary approach including clinicians and geneticists is crucial when dealing with patients with rare diseases.

VIsoQLR: an interactive tool for the detection, quantification and fine-tuning of isoforms in selected genes using long-read sequencing

DNA variants altering the pre-mRNA splicing process represent an underestimated cause of human genetic diseases. Their association with disease traits should be confirmed using functional assays from patient cell lines or alternative models to detect aberrant mRNAs. Long-read sequencing is a suitable technique to identify and quantify mRNA isoforms. Available isoform detection and/or quantification tools are generally designed for the whole transcriptome analysis. However experiments focusing on genes of interest need more precise data fine-tuning and visualization tools.Here we describe VIsoQLR, an interactive analyzer, viewer and editor for the semi-automated identification and quantification of known and novel isoforms using long-read sequencing data. VIsoQLR is tailored to thoroughly analyze mRNA expression in splicing assays of selected genes. Our tool takes sequences aligned to a reference, and for each gene, it defines consensus splice sites and quantifies isoforms. VIsoQLR introduces features to edit the splice sites through dynamic and interactive graphics and tables, allowing accurate manual curation. Known isoforms detected by other methods can also be imported as references for comparison. A benchmark against two other popular transcriptome-based tools shows VIsoQLR accurate performance on both detection and quantification of isoforms. Here, we present VIsoQLR principles and features and its applicability in a case study example using nanopore-based long-read sequencing. VIsoQLR is available at https://github.com/TBLabFJD/VIsoQLR .

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.

Minigene Splicing Assays and Long-Read Sequencing to Unravel Pathogenic Deep-Intronic Variants in PAX6 in Congenital Aniridia

PAX6 haploinsufficiency causes aniridia, a congenital eye disorder that involves the iris, and foveal hypoplasia. Comprehensive screening of the PAX6 locus, including the non-coding regions, by next-generation sequencing revealed four deep-intronic variants with potential effects on pre-RNA splicing. Nevertheless, without a functional analysis, their pathogenicity could not be established. We aimed to decipher their impact on the canonical PAX6 splicing using in vitro minigene splicing assays and nanopore-based long-read sequencing. Two multi-exonic PAX6 constructs were generated, and minigene assays were carried out. An aberrant splicing pattern was observed for two variants in intron 6, c.357+136G>A and c.357+334G>A. In both cases, several exonization events, such as pseudoexon inclusions and partial intronic retention, were observed due to the creation or activation of new/cryptic non-canonical splicing sites, including a shared intronic donor site. In contrast, two variants identified in intron 11, c.1032+170A>T and c.1033-275A>C, seemed not to affect splicing processes. We confirmed the high complexity of alternative splicing of PAX6 exon 6, which also involves unreported cryptic intronic sites. Our study highlights the importance of integrating functional studies into diagnostic algorithms to decipher the potential implication of non-coding variants, usually classified as variants of unknown significance, thus allowing variant reclassification to achieve a conclusive genetic diagnosis.

Computational prediction of human deep intronic variation

BACKGROUND

The adoption of whole-genome sequencing in genetic screens has facilitated the detection of genetic variation in the intronic regions of genes, far from annotated splice sites. However, selecting an appropriate computational tool to discriminate functionally relevant genetic variants from those with no effect is challenging, particularly for deep intronic regions where independent benchmarks are scarce.

RESULTS

In this study, we have provided an overview of the computational methods available and the extent to which they can be used to analyze deep intronic variation. We leveraged diverse datasets to extensively evaluate tool performance across different intronic regions, distinguishing between variants that are expected to disrupt splicing through different molecular mechanisms. Notably, we compared the performance of SpliceAI, a widely used sequence-based deep learning model, with that of more recent methods that extend its original implementation. We observed considerable differences in tool performance depending on the region considered, with variants generating cryptic splice sites being better predicted than those that potentially affect splicing regulatory elements. Finally, we devised a novel quantitative assessment of tool interpretability and found that tools providing mechanistic explanations of their predictions are often correct with respect to the ground - information, but the use of these tools results in decreased predictive power when compared to black box methods.

CONCLUSIONS

Our findings translate into practical recommendations for tool usage and provide a reference framework for applying prediction tools in deep intronic regions, enabling more informed decision-making by practitioners.

IntroVerse: a comprehensive database of introns across human tissues

Dysregulation of RNA splicing contributes to both rare and complex diseases. RNA-sequencing data from human tissues has shown that this process can be inaccurate, resulting in the presence of novel introns detected at low frequency across samples and within an individual. To enable the full spectrum of intron use to be explored, we have developed IntroVerse, which offers an extensive catalogue on the splicing of 332,571 annotated introns and a linked set of 4,679,474 novel junctions covering 32,669 different genes. This dataset has been generated through the analysis of 17,510 human control RNA samples from 54 tissues provided by the Genotype-Tissue Expression Consortium. IntroVerse has two unique features: (i) it provides a complete catalogue of novel junctions and (ii) each novel junction has been assigned to a specific annotated intron. This unique, hierarchical structure offers multiple uses, including the identification of novel transcripts from known genes and their tissue-specific usage, and the assessment of background splicing noise for introns thought to be mis-spliced in disease states. IntroVerse provides a user-friendly web interface and is freely available at https://rytenlab.com/browser/app/introverse.

Unpredicted Aberrant Splicing Products Identified in Postmortem Sudden Cardiac Death Samples

Molecular screening for pathogenic mutations in sudden cardiac death (SCD)-related genes is common practice for SCD cases. However, test results may lead to uncertainty because of the identification of variants of unknown significance (VUS) occurring in up to 70% of total identified variants due to a lack of experimental studies. Genetic variants affecting potential splice site variants are among the most difficult to interpret. The aim of this study was to examine rare intronic variants identified in the exonic flanking sequence to meet two main objectives: first, to validate that canonical intronic variants produce aberrant splicing; second, to determine whether rare intronic variants predicted as VUS may affect the splicing product. To achieve these objectives, 28 heart samples of cases of SCD carrying rare intronic variants were studied. Samples were analyzed using 85 SCD genes in custom panel sequencing. Our results showed that rare intronic variants affecting the most canonical splice sites displayed in 100% of cases that they would affect the splicing product, possibly causing aberrant isoforms. However, 25% of these cases (1/4) showed normal splicing, contradicting the in silico results. On the contrary, in silico results predicted an effect in 0% of cases, and experimental results showed >20% (3/14) unpredicted aberrant splicing. Thus, deep intron variants are likely predicted to not have an effect, which, based on our results, might be an underestimation of their effect and, therefore, of their pathogenicity classification and family members’ follow-up.

Predictive functional, statistical and structural analysis of CSNK2A1 and CSNK2B variants linked to neurodevelopmental diseases

Okur-Chung Neurodevelopmental Syndrome (OCNDS) and Poirier-Bienvenu Neurodevelopmental Syndrome (POBINDS) were recently identified as rare neurodevelopmental disorders. OCNDS and POBINDS are associated with heterozygous mutations in the CSNK2A1 and CSNK2B genes which encode CK2α, a serine/threonine protein kinase, and CK2β, a regulatory protein, respectively, which together can form a tetrameric enzyme called protein kinase CK2. A challenge in OCNDS and POBINDS is to understand the genetic basis of these diseases and the effect of the various CK2⍺ and CK2β mutations. In this study we have collected all variants available to date in CSNK2A1 and CSNK2B, and identified hotspots. We have investigated CK2⍺ and CK2β missense mutations through prediction programs which consider the evolutionary conservation, functionality and structure or these two proteins, compared these results with published experimental data on CK2α and CK2β mutants, and suggested prediction programs that could help predict changes in functionality of CK2α mutants. We also investigated the potential effect of CK2α and CK2β mutations on the 3D structure of the proteins and in their binding to each other. These results indicate that there are functional and structural consequences of mutation of CK2α and CK2β, and provide a rationale for further study of OCNDS and POBINDS-associated mutations. These data contribute to understanding the genetic and functional basis of these diseases, which is needed to identify their underlying mechanisms.

A catalog of the genetic causes of hereditary angioedema in the Canary Islands (Spain)

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

Short amplicon reverse transcription-polymerase chain reaction detects aberrant splicing in genes with low expression in blood missed by ribonucleic acid sequencing analysis for clinical diagnosis

Use of blood RNA sequencing (RNA-seq) as a splicing analysis tool for clinical interpretation of variants of uncertain significance (VUSs) found via whole-genome and exome sequencing can be difficult for genes that have low expression in the blood due to insufficient read count coverage aligned to specific genes of interest. Here, we present a short amplicon reverse transcription-polymerase chain reaction(RT-PCR) for the detection of genes with low blood expression. Short amplicon RT-PCR, is designed to span three exons where an exon harboring a variant is flanked by one upstream and one downstream exon. We tested short amplicon RT-PCRs for genes that have median transcripts per million (TPM) values less than one according to the genotype-tissue expression database. Median TPM values of genes analyzed in this study are SYN1 = 0.8549, COL1A1 = 0.6275, TCF4 = 0.4009, DSP = .2894, TTN = 0.2851, COL5A2 = 0.1036, TERT = 0.04452, NTRK2 = 0.0344, ABCA4 = 0.00744, PRPH = 0, and WT1 = 0. All these genes show insufficient exon-spanning read coverage in our RNA-seq data to allow splicing analysis. We successfully detected all genes tested except PRPH and WT1. Aberrant splicing was detected in SYN1, TCF4, NTRK2, TTN, and TERT VUSs. Therefore, our results show short amplicon RT-PCR is a useful alternative for the analysis of splicing events in genes with low TPM in blood RNA for clinical diagnostics.

New Developments and Possibilities in Reanalysis and Reinterpretation of Whole Exome Sequencing Datasets for Unsolved Rare Diseases Using Machine Learning Approaches

Rare diseases impact the lives of 300 million people in the world. Rapid advances in bioinformatics and genomic technologies have enabled the discovery of causes of 20–30% of rare diseases. However, most rare diseases have remained as unsolved enigmas to date. Newer tools and availability of high throughput sequencing data have enabled the reanalysis of previously undiagnosed patients. In this review, we have systematically compiled the latest developments in the discovery of the genetic causes of rare diseases using machine learning methods. Importantly, we have detailed methods available to reanalyze existing whole exome sequencing data of unsolved rare diseases. We have identified different reanalysis methodologies to solve problems associated with sequence alterations/mutations, variation re-annotation, protein stability, splice isoform malfunctions and oligogenic analysis. In addition, we give an overview of new developments in the field of rare disease research using whole genome sequencing data and other omics.

Targeted RNAseq Improves Clinical Diagnosis of Very Early-Onset Pediatric Immune Dysregulation

Despite increased use of whole exome sequencing (WES) for the clinical analysis of rare disease, overall diagnostic yield for most disorders hovers around 30%. Previous studies of mRNA have succeeded in increasing diagnoses for clearly defined disorders of monogenic inheritance. We asked if targeted RNA sequencing could provide similar benefits for primary immunodeficiencies (PIDs) and very early-onset inflammatory bowel disease (VEOIBD), both of which are difficult to diagnose due to high heterogeneity and variable severity. We performed targeted RNA sequencing of a panel of 260 immune-related genes for a cohort of 13 patients (seven suspected PID cases and six VEOIBD) and analyzed variants, splicing, and exon usage. Exonic variants were identified in seven cases, some of which had been previously prioritized by exome sequencing. For four cases, allele specific expression or lack thereof provided additional insights into possible disease mechanisms. In addition, we identified five instances of aberrant splicing associated with four variants. Three of these variants had been previously classified as benign in ClinVar based on population frequency. Digenic or oligogenic inheritance is suggested for at least two patients. In addition to validating the use of targeted RNA sequencing, our results show that rare disease research will benefit from incorporating contributing genetic factors into the diagnostic approach.

ShapeSorter: a fully probabilistic method for detecting conserved RNA structure features supported by SHAPE evidence

There is an increased interest in the determination of RNA structures in vivo as it is now possible to probe them in a high-throughput manner, e.g. using SHAPE protocols. By now, there exist a range of computational methods that integrate experimental SHAPE-probing evidence into computational RNA secondary structure prediction. The state-of-the-art in this field is currently provided by computational methods that employ the minimum-free energy strategy for prediction RNA secondary structures with SHAPE-probing evidence. These methods, however, rely on the assumption that transcripts in vivo fold into the thermodynamically most stable configuration and ignore evolutionary evidence for conserved RNA structure features. We here present a new computational method, ShapeSorter, that predicts RNA structure features without employing the thermodynamic strategy. Instead, ShapeSorter employs a fully probabilistic framework to identify RNA structure features that are supported by evolutionary and SHAPE-probing evidence. Our method can capture RNA structure heterogeneity, pseudo-knotted RNA structures as well as transient and mutually exclusive RNA structure features. Moreover, it estimates P-values for the predicted RNA structure features which allows for easy filtering and ranking. We investigate the merits of our method in a comprehensive performance benchmarking and conclude that ShapeSorter has a significantly superior performance for predicting base-pairs than the existing state-of-the-art methods.

Toward transcriptomics as a primary tool for rare disease investigation

In the past 5 years transcriptome or RNA-sequencing (RNA-seq) has steadily emerged as a complementary assay for rare disease diagnosis and discovery. In this perspective, we summarize several recent developments and challenges in the use of RNA-seq for rare disease investigation. Using an accessible patient sample, such as blood, skin, or muscle, RNA-seq enables the assay of expressed RNA transcripts. Analysis of RNA-seq allows the identification of aberrant or outlier gene expression and alternative splicing as functional evidence to support rare disease study and diagnosis. Further, many types of variant effects can be profiled beyond coding variants, as the consequences of noncoding variants that impact gene expression and splicing can be directly observed. This is particularly apparent for structural variants that disproportionately underlie outlier gene expression and for splicing variants in which RNA-seq can both measure aberrant canonical splicing and detect deep intronic effects. However, a major potential limitation of RNA-seq in rare disease investigation is the developmental and cell type specificity of gene expression as a pathogenic variant's effect may be limited to a specific spatiotemporal context and access to a patient's tissue sample from the relevant tissue and timing of disease expression may not be possible. We speculate that as advances in computational methods and emerging experimental techniques overcome both developmental and cell type specificity, there will be broadening use of RNA sequencing and multiomics in rare disease diagnosis and delivery of precision health.

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

Antisense-Mediated Down-Regulation of Factor V-Short Splicing in a Liver Cell Line Model

Coagulation factor V (FV) is a liver-derived protein encoded by the F5 gene. Alternative splicing of F5 exon 13 produces a low-abundance splicing isoform, known as FV-short, which binds the anticoagulant protein tissue factor pathway inhibitor (TFPIα) with high affinity, stabilising it in the circulation and potently enhancing its anticoagulant activity. Accordingly, rare F5 gene mutations that up-regulate FV-short splicing are associated with bleeding. In this study we have explored the possibility of decreasing FV-short splicing by antisense-based splicing modulation. To this end, we have designed morpholino antisense oligonucleotides (MAOs) targeting the FV-short-specific donor and acceptor splice sites and tested their efficacy in a liver cell line (HepG2) that naturally expresses full-length FV and FV-short. Cells were treated with 0–20 µM MAO, and full-length FV and FV-short mRNA expression was analysed by RT-(q)PCR. Both MAOs, alone or in combination, decreased the FV-short/full-length FV mRNA ratio down to ~50% of its original value in a specific and dose-dependent manner. This pilot study provides proof-of-principle for the possibility to decrease FV-short expression by antisense-mediated splicing modulation. In turn, this may form the basis for novel therapeutic approaches to bleeding disorders caused by FV-short over-expression and/or elevated TFPIα (activity) levels.