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.

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.

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.

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.

Mutations in ANAPC1, Encoding a Scaffold Subunit of the Anaphase-Promoting Complex, Cause Rothmund-Thomson Syndrome Type 1.

Rothmund-Thomson syndrome (RTS) is an autosomal-recessive disorder characterized by poikiloderma, sparse hair, short stature, and skeletal anomalies. Type 2 RTS, which is defined by the presence of bi-allelic mutations in RECQL4, is characterized by increased cancer susceptibility and skeletal anomalies, whereas the genetic basis of RTS type 1, which is associated with juvenile cataracts, is unknown. We studied ten individuals, from seven families, who had RTS type 1 and identified a deep intronic splicing mutation of the ANAPC1 gene, a component of the anaphase-promoting complex/cyclosome (APC/C), in all affected individuals, either in the homozygous state or in trans with another mutation. Fibroblast studies showed that the intronic mutation causes the activation of a 95 bp pseudoexon, leading to mRNAs with premature termination codons and nonsense-mediated decay, decreased ANAPC1 protein levels, and prolongation of interphase. Interestingly, mice that were heterozygous for a knockout mutation have an increased incidence of cataracts. Our results demonstrate that deficiency in the APC/C is a cause of RTS type 1 and suggest a possible link between the APC/C and RECQL4 helicase because both proteins are involved in DNA repair and replication.

Macroglossia in a pig diagnosed as Becker muscular dystrophy due to dystrophin pseudoexon insertion derived from intron 26.

We report a case of Becker muscular dystrophy in a 6-month-old, mixed-breed, castrated male pig detected with macroglossia at a meat inspection center. The pig presented a severely enlarged tongue extending outside its mouth. The tongue was firm and pale with discolored muscles. Histologically, there was severe fibrosis, fatty replacement, and myofiber necrosis, degeneration, and regeneration. Immunofluorescence showed focal and severely weak labeling for dystrophin at the sarcolemma of myocytes in the tongue. Analysis of dystrophin mRNA showed a 62 base pair insertion between exons 26 and 27. The insertion was derived from intron 26. Based on these findings, we diagnosed the case as Becker muscular dystrophy-the first known muscular dystrophy case induced by pseudoexon insertion in animals.

Identification of Spliceogenic Variants beyond Canonical GT-AG Splice Sites in Hereditary Cancer Genes.

Pathogenic/likely pathogenic variants in susceptibility genes that interrupt RNA splicing are a well-documented mechanism of hereditary cancer syndromes development. However, if RNA studies are not performed, most of the variants beyond the canonical GT-AG splice site are characterized as variants of uncertain significance (VUS). To decrease the VUS burden, we have bioinformatically evaluated all novel VUS detected in 732 consecutive patients tested in the routine genetic counseling process. Twelve VUS that were predicted to cause splicing defects were selected for mRNA analysis. Here, we report a functional characterization of 12 variants located beyond the first two intronic nucleotides using RNAseq in APC, ATM, FH, LZTR1, MSH6, PALB2, RAD51C, and TP53 genes. Based on the analysis of mRNA, we have successfully reclassified 50% of investigated variants. 25% of variants were downgraded to likely benign, whereas 25% were upgraded to likely pathogenic leading to improved clinical management of the patient and the family members.

Novel deep intronic mutation in PLA2G6 causing early-onset Parkinson's disease with brain iron accumulation through pseudo-exon activation.

PLA2G6 is the causative gene for a group of autosomal recessive neurodegenerative disorders known as PLA2G6-associated neurodegeneration (PLAN). We present a case with early-onset parkinsonism, ataxia, cognitive decline, cerebellar atrophy, and brain iron accumulation. Sequencing of PLA2G6 coding regions identified only a heterozygous nonsense variant, but mRNA analysis revealed the presence of an aberrant transcript isoform due to a novel deep intronic variant (c.2035-274G > A) leading to activation of an intronic pseudo-exon. These results expand the genotypic spectrum of PLAN, showing the paramount importance of detecting possible pathogenic variants in deep intronic regions in undiagnosed patients.

Identification and Rescue of Splice Defects Caused by Two Neighboring Deep-Intronic ABCA4 Mutations Underlying Stargardt Disease.

Sequence analysis of the coding regions and splice site sequences in inherited retinal diseases is not able to uncover ∼40% of the causal variants. Whole-genome sequencing can identify most of the non-coding variants, but their interpretation is still very challenging, in particular when the relevant gene is expressed in a tissue-specific manner. Deep-intronic variants in ABCA4 have been associated with autosomal-recessive Stargardt disease (STGD1), but the exact pathogenic mechanism is unknown. By generating photoreceptor precursor cells (PPCs) from fibroblasts obtained from individuals with STGD1, we demonstrated that two neighboring deep-intronic ABCA4 variants (c.4539+2001G>A and c.4539+2028C>T) result in a retina-specific 345-nt pseudoexon insertion (predicted protein change: p.Arg1514Leufs∗36), likely due to the creation of exonic enhancers. Administration of antisense oligonucleotides (AONs) targeting the 345-nt pseudoexon can significantly rescue the splicing defect observed in PPCs of two individuals with these mutations. Intriguingly, an AON that is complementary to c.4539+2001G>A rescued the splicing defect only in PPCs derived from an individual with STGD1 with this but not the other mutation, demonstrating the high specificity of AONs. In addition, a single AON molecule rescued splicing defects associated with different neighboring mutations, thereby providing new strategies for the treatment of persons with STGD1. As many genes associated with human genetic conditions are expressed in specific tissues and pre-mRNA splicing may also rely on organ-specific factors, our approach to investigate and treat splicing variants using differentiated cells derived from individuals with STGD1 can be applied to any tissue of interest.

Missing heritability in Bloom syndrome: First report of a deep intronic variant leading to pseudo-exon activation in the BLM gene.

Pathogenic biallelic variants in the BLM/RECQL3 gene cause a rare autosomal recessive disorder called Bloom syndrome (BS). This syndrome is characterized by severe growth delay, immunodeficiency, dermatological manifestations and a predisposition to a wide variety of cancers, often multiple and very early in life. Literature shows that the main mode of BLM inactivation is protein translation termination. We expanded the molecular spectrum of BS by reporting the first deep intronic variant causing intron exonisation. We describe a patient with a clinical phenotype of BS and a strong increase in sister chromatid exchanges (SCE), who was found to be compound heterozygous for a novel nonsense variant c.3379C>T, p.(Gln1127Ter) in exon 18 and a deep intronic variant c.3020-258A>G in intron 15 of the BLM gene. The deep intronic variant creates a high-quality de novo donor splice site, which leads to retention of two intron segments. Both pseudo-exons introduce a premature stop codon into the reading frame and abolish BLM protein expression, confirmed by Western Blot analysis. These findings illustrate the role of non-coding variation in Mendelian disorders and herewith highlight an unmet need in routine testing of Mendelian disorders, being the added value of RNA-based approaches to provide a complete molecular diagnosis.

In silico identification of pseudo-exon activation events in personal genome and transcriptome data.

Causative mutations for human genetic disorders have mainly been identified in exonic regions that code for amino acid sequences. Recently, however, it has been reported that mutations in deep intronic regions can also cause certain human genetic disorders by creating novel splice sites, leading to pseudo-exon activation. To investigate how frequently pseudo-exon activation events occur in normal individuals, we conducted in silico identification of such events using personal genome data and corresponding high-quality transcriptome data. With rather stringent conditions, on average, 2.6 pseudo-exon activation events per individual were identified. More pseudo-exon activation events were found in 5' donor splice sites than in 3' acceptor splice sites. Although pseudo-exon activation events have sporadically been reported as causative mutations in genetic disorders, it is revealed in this study that such events can be observed in normal individuals at a certain frequency. We estimate that human genomes typically contain on average at least 10 pseudo-exon activation events. The actual number should be higher than this, because we used stringent criteria to identify pseudo-exon activation events. This suggests that it is worth considering the possibility of pseudo-exon activation when searching for causative mutations of genetic disorders if candidate mutations are not identified in coding regions or RNA splice sites.

Antisense Oligonucleotide-Based Rescue of Aberrant Splicing Defects Caused by 15 Pathogenic Variants in ABCA4.

The discovery of novel intronic variants in the ABCA4 locus has contributed significantly to solving the missing heritability in Stargardt disease (STGD1). The increasing number of variants affecting pre-mRNA splicing makes ABCA4 a suitable candidate for antisense oligonucleotide (AON)-based splicing modulation therapies. In this study, AON-based splicing modulation was assessed for 15 recently described intronic variants (three near-exon and 12 deep-intronic variants). In total, 26 AONs were designed and tested in vitro using a midigene-based splice system. Overall, partial or complete splicing correction was observed for two variants causing exon elongation and all variants causing pseudoexon inclusion. Together, our results confirm the high potential of AONs for the development of future RNA therapies to correct splicing defects causing STGD1.

Deep-intronic variants in CNGB3 cause achromatopsia by pseudoexon activation.

Our comprehensive cohort of 1100 unrelated achromatopsia (ACHM) patients comprises a considerable number of cases (~5%) harboring only a single pathogenic variant in the major ACHM gene CNGB3. We sequenced the entire CNGB3 locus in 33 of these patients to find a second variant which eventually explained the patients' phenotype. Forty-seven intronic CNGB3 variants were identified in 28 subjects after a filtering step based on frequency and the exclusion of variants found in cis with pathogenic alleles. In a second step, in silico prediction tools were used to filter out those variants with little odds of being deleterious. This left three variants that were analyzed using heterologous splicing assays. Variant c.1663-1205G>A, found in 14 subjects, and variant c.1663-2137C>T, found in two subjects, were indeed shown to exert a splicing defect by causing pseudoexon insertion into the transcript. Subsequent screening of further unsolved CNGB3 subjects identified four additional cases harboring the c.1663-1205G>A variant which makes it the eighth most frequent CNGB3 variant in our cohort. Compound heterozygosity could be validated in ten cases. Our study demonstrates that whole gene sequencing can be a powerful approach to identify the second pathogenic allele in patients apparently harboring only one disease-causing variant.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars.

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.