A deep intronic DLG4 variant resulting in DLG4-related synaptopathy.

The rare autosomal dominant brain disorder DLG4-related synaptopathy is caused by de novo variants in DLG4 (encoding PSD-95), the majority of which are predicted to be protein-truncating. In addition to splice site variants, a number of synonymous and missense DLG4 variants are predicted to exert their effect through altered RNA splicing, although the pathogenicity of these variants is uncertain without functional RNA studies. Here, we describe a young boy with a deep intronic DLG4 variant (c.2105+235C>T) identified using whole genome sequencing. By using reverse-transcription PCR on RNA derived from peripheral blood, we demonstrate that DLG4 mRNA expression is detectable in blood and the deep intronic variant gives rise to two alternative DLG4 transcripts, one of which includes a pseudoexon. Both alternative transcripts are out-of-frame and predicted to result in protein-truncation, thereby establishing the genetic diagnosis for the proband. This adds to the evidence concerning the pathogenic potential of deep intronic variants and underlines the importance of functional studies, even in cases where reported tissue-specific gene expression might suggest otherwise.

Pseudoexon activation by deep intronic variation in GNE myopathy with thrombocytopenia.

INTRODUCTION/AIMS: GNE myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway. Although over 300 GNE variants have been reported, some patients remain undiagnosed with monoallelic pathogenic variants. This study aims to analyze the entire GNE genomic region to identify novel pathogenic variants. METHODS: Patients with clinically compatible GNE myopathy and monoallelic pathogenic variants in the GNE gene were enrolled. The other GNE pathogenic variant was verified using comprehensive methods including exon 2 quantitative polymerase chain reaction and nanopore long-read single-molecule sequencing (LRS). RESULTS: A deep intronic GNE variant, c.862+870C>T, was identified in nine patients from eight unrelated families. This variant generates a cryptic splice site, resulting in the activation of a novel pseudoexon between exons 5 and 6. It results in the insertion of an extra 146 nucleotides into the messengerRNA (mRNA), which is predicted to result in a truncated humanGNE1(hGNE1) protein. Peanut agglutinin（PNA） lectin staining of muscle tissues showed reduced sialylation of mucin O-glycans on sarcolemmal glycoproteins. Notably, a third of patients with the c.862+870C>T variant exhibited thrombocytopenia. A common core haplotype harboring the deep intronic GNE variant was found in all these patients. DISCUSSION: The transcript with pseudoexon activation potentially affects sialic acid biosynthesis via nonsense-mediated mRNA decay, or resulting in a truncated hGNE1 protein, which interferes with normal enzyme function. LRS is expected to be more frequently incorporated in genetic analysis given its efficacy in detecting hard-to-find pathogenic variants.

Exonisation of an intronic L1 element in the dystrophin gene associated with X-linked muscular dystrophy in a Border Collie dog.

X-linked recessive dystrophinopathies are the most common muscular dystrophies (MDs) in humans and dogs. To date, 20 breed-specific MD-associated variants are described in the canine dystrophin gene (DMD), including one associated with dystrophin-deficient MD in the Border Collie mixed breed. Here, we report the diagnosis and follow-up of mild dystrophin-deficient MD in a 5-month-old male Border Collie, associated with a novel DMD variant. Diagnosis was based on neurological examination and laboratory evaluations including creatine kinase activity, electromyography and muscle biopsies with immunofluorescent staining. Inspection of the Sashimi plots of the RNA-seq data from the affected muscle biopsy led to the discovery of a 162-bp L1 pseudoexon in DMD intron 63, introducing a frameshift and a premature stop codon (NM_001003343.1: c.9271_9272insN[162] p.(Ala3091fs*21)). Reduced DMD mRNA levels were detected for both the non-pseudoexon (50× less) and pseudoexon (3× less) containing transcripts in the affected muscle, compared with the level of the non-pseudoexon containing transcript in a control muscle, resulting in very low dystrophin protein levels and the upregulation of utrophin. Because the variant was only found in the affected dog, not in the healthy mother and grandmother, or in 108 unrelated Border Collies from the Belgian population (46 males and 62 females), it was considered a de novo variant. Although the prognosis for dystrophinopathy is generally regarded as poor, the dog stabilised at the age of 6 months and is still clinically stable at the age of 2 years.

Functional Characterization of Splice Variants in the Diagnosis of Albinism.

Albinism is a genetically heterogeneous disease in which 21 genes are known so far. Its inheritance mode is autosomal recessive except for one X-linked form. The molecular analysis of exonic sequences of these genes allows for about a 70% diagnostic rate. About half (15%) of the unsolved cases are heterozygous for one pathogenic or probably pathogenic variant. Assuming that the missing variant may be located in non-coding regions, we performed sequencing for 122 such heterozygous patients of either the whole genome (27 patients) or our NGS panel (95 patients) that includes, in addition to all exons of the 21 genes, the introns and flanking sequences of five genes, TYR, OCA2, SLC45A2, GPR143 and HPS1. Rare variants (MAF < 0.01) in trans to the first variant were tested by RT-PCR and/or minigene assay. Of the 14 variants tested, nine caused either exon skipping or the inclusion of a pseudoexon, allowing for the diagnosis of 11 patients. This represents 9.8% (12/122) supplementary diagnosis for formerly unsolved patients and 75% (12/16) of those in whom the candidate variant was in trans to the first variant. Of note, one missense variant was demonstrated to cause skipping of the exon in which it is located, thus shedding new light on its pathogenic mechanism. Searching for non-coding variants and testing them for an effect on RNA splicing is warranted in order to increase the diagnostic rate.

Deep Intronic ETFDH Variants Represent a Recurrent Pathogenic Event in Multiple Acyl-CoA Dehydrogenase Deficiency.

Multiple acyl-CoA dehydrogenase deficiency (MADD) is a rare inborn error of metabolism affecting fatty acid and amino acid oxidation with an incidence of 1 in 200,000 live births. MADD has three clinical phenotypes: severe neonatal-onset with or without congenital anomalies, and a milder late-onset form. Clinical diagnosis is supported by urinary organic acid and blood acylcarnitine analysis using tandem mass spectrometry in newborn screening programs. MADD is an autosomal recessive trait caused by biallelic mutations in the ETFA, ETFB, and ETFDH genes encoding the alpha and beta subunits of the electron transfer flavoprotein (ETF) and ETF-coenzyme Q oxidoreductase enzymes. Despite significant advancements in sequencing techniques, many patients remain undiagnosed, impacting their access to clinical care and genetic counseling. In this report, we achieved a definitive molecular diagnosis in a newborn by combining whole-genome sequencing (WGS) with RNA sequencing (RNA-seq). Whole-exome sequencing and next-generation gene panels fail to detect variants, possibly affecting splicing, in deep intronic regions. Here, we report a unique deep intronic mutation in intron 1 of the ETFDH gene, c.35-959A>G, in a patient with early-onset lethal MADD, resulting in pseudo-exon inclusion. The identified variant is the third mutation reported in this region, highlighting ETFDH intron 1 vulnerability. It cannot be excluded that these intronic sequence features may be more common in other genes than is currently believed. This study highlights the importance of incorporating RNA analysis into genome-wide testing to reveal the functional consequences of intronic mutations.

Nanopore Deep Sequencing as a Tool to Characterize and Quantify Aberrant Splicing Caused by Variants in Inherited Retinal Dystrophy Genes.

The contribution of splicing variants to molecular diagnostics of inherited diseases is reported to be less than 10%. This figure is likely an underestimation due to several factors including difficulty in predicting the effect of such variants, the need for functional assays, and the inability to detect them (depending on their locations and the sequencing technology used). The aim of this study was to assess the utility of Nanopore sequencing in characterizing and quantifying aberrant splicing events. For this purpose, we selected 19 candidate splicing variants that were identified in patients affected by inherited retinal dystrophies. Several in silico tools were deployed to predict the nature and estimate the magnitude of variant-induced aberrant splicing events. Minigene assay or whole blood-derived cDNA was used to functionally characterize the variants. PCR amplification of minigene-specific cDNA or the target gene in blood cDNA, combined with Nanopore sequencing, was used to identify the resulting transcripts. Thirteen out of nineteen variants caused aberrant splicing events, including cryptic splice site activation, exon skipping, pseudoexon inclusion, or a combination of these. Nanopore sequencing allowed for the identification of full-length transcripts and their precise quantification, which were often in accord with in silico predictions. The method detected reliably low-abundant transcripts, which would not be detected by conventional strategies, such as RT-PCR followed by Sanger sequencing.

An FBN1 deep intronic variant is associated with pseudoexon formation and a variable Marfan phenotype in a five generation family.

Exome sequencing of genes associated with heritable thoracic aortic disease (HTAD) failed to identify a pathogenic variant in a large family with Marfan syndrome (MFS). A genome-wide linkage analysis for thoracic aortic disease identified a peak at 15q21.1, and genome sequencing identified a novel deep intronic FBN1 variant that segregated with thoracic aortic disease in the family (LOD score 2.7) and was predicted to alter splicing. RT-PCR and bulk RNA sequencing of RNA harvested from fibroblasts explanted from the affected proband revealed an insertion of a pseudoexon between exons 13 and 14 of the FBN1 transcript, predicted to lead to nonsense mediated decay (NMD). Treating the fibroblasts with an NMD inhibitor, cycloheximide, greatly improved the detection of the pseudoexon-containing transcript. Family members with the FBN1 variant had later onset aortic events and fewer MFS systemic features than typical for individuals with haploinsufficiency of FBN1. Variable penetrance of the phenotype and negative genetic testing in MFS families should raise the possibility of deep intronic FBN1 variants and the need for additional molecular studies.

Compound heterozygous splicing variants expand the genotypic spectrum of EMC1-related disorders.

EMC1 encodes subunit 1 of the endoplasmic reticulum (ER) membrane protein complex (EMC), a transmembrane domain insertase involved in membrane protein biosynthesis. Variants in EMC1 are described as a cause of global developmental delay, hypotonia, cortical visual impairment, and commonly, cerebral atrophy on MRI scan. We report an individual with severe global developmental delay and progressive cerebellar atrophy in whom exome sequencing identified a heterozygous essential splice-site variant in intron-3 of EMC1 (NM_015047.3:c.287-1G>A). Whole genome sequencing (WGS) identified a deep intronic variant in intron-20 of EMC1 (NM_015047.3:c.2588-771C>G) that was poorly predicted by in silico programs to disrupt pre-mRNA splicing. Reverse Transcription-PCR (RT-PCR) revealed stochastic activation of a pseudo-exon associated with the c.2588-771C>G variant and mis-splicing arising from the c.287-1G>A variant. This case highlights the utility of WGS and RNA studies to identify and assess likely pathogenicity of deep intronic variants and expands the genotypic and phenotypic spectrum of EMC1-related disorders.

Intron mutations and early transcription termination in Duchenne and Becker muscular dystrophy.

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome.

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.