Cryptic exon activation causes dystrophinopathy in two Chinese families

Data-driven insights to inform splice-altering variant assessment

Disease-causing genetic variants often disrupt mRNA splicing, an intricate process that is incompletely understood. Thus, accurate inference of which genetic variants will affect splicing and what their functional consequences will be is challenging, particularly for variants outside of the essential splice sites. Here, we describe a set of data-driven heuristics that inform the interpretation of human splice-altering variants (SAVs) based on the analysis of annotated exons, experimentally validated SAVs, and the currently understood principles of splicing biology. We defined requisite splicing criteria by examining around 202,000 canonical protein-coding exons and 19,000 experimentally validated splicing branchpoints. This analysis defined the sequence, spacing, and motif strength required for splicing, with 95.9% of the exons examined meeting these criteria. By considering over 12,000 experimentally validated variants from the SpliceVarDB, we defined a set of heuristics that inform the evaluation of putative SAVs. To ensure the applicability of each heuristic, only those supported by at least 10 experimentally validated variants were considered. This allowed us to establish a measure of spliceogenicity: the proportion of variants at a location (or motif site) that affected splicing in a given context. This study makes considerable advances toward bridging the gap between computational predictions and the biological process of splicing, offering an evidence-based approach to identifying SAVs and evaluating their impact. Our splicing heuristics enhance the current framework for genetic variant evaluation with a robust, detailed, and comprehensible analysis by adding valuable context over traditional binary prediction tools.

Minigene Assay as an Effective Molecular Diagnostic Strategy in Determining the Pathogenicity of Noncanonical Splice-Site Variants in FLCN

Primary spontaneous pneumothorax (PSP) or pulmonary cyst is one of the manifestations of Birt-Hogg-Dubé syndrome, which is caused by pathogenic variants in FLCN gene. Genetic testing in patients with PSP identifies a certain number of missense or intronic variants. These variants are usually considered as variants of uncertain significance, whose functional interpretations pose a challenge in clinical genetics. To improve recognition of pathogenic splice-altering variants in FLCN gene, computational tools are used to prioritize potential splice-altering variants and then a hybrid minigene assay is performed to verify the RNA splicing pattern. Herein, variants in FLCN exon 11 and its flanking sequence are focused. Eight variants detected in 11 patients with PSP are evaluated, and six variants are prioritized by in silico tools as potential splice-altering variants of uncertain significance. Four variants (c.1177-5_1177-3delCTC, c.1292_1300+4del, c.1300+4C>T, and c.1300+5G>A) are demonstrated by minigene assay to alter RNA splicing of FLCN, and the last three of them are novel. RT-PCR of patient-derived RNA gives consistent results. Genotype-phenotype correlation analysis in patients with PSP with these variants demonstrates good concordance. Our results underline the importance of RNA analysis, which could provide molecular evidence for pathogenicity of a variant, and provide essential information for the clinical interpretation of variants. Combining the clinical information, a definitive diagnosis could be made.

Generation and characterization of an induced pluripotent stem cell line (FJMUNi001-A) from a patient with Duchenne muscular dystrophy carrying c.4518 + 512 T > A variant in the DMD gene

Duchenne muscular dystrophy (DMD) is an X-linked recessive degenerative disease characterized by progressive weakness of limbs. In this study, we generated an induced pluripotent stem cell line from a DMD patient's dermal fibroblasts with non-integrating Sendai virus. The patient carried a rare c.4518 + 512 T > A variant in the DMD gene. This iPSC line displayed normal iPSC morphology, karyotypes and pluripotency expression markers, and also can be differentiated into the three germ layers.

Novel Intronic Mutations of TBK1 Promote Aberrant Splicing Modes in Amyotrophic Lateral Sclerosis

TANK-binding kinase 1 (TBK1) has been identified as a causative gene of amyotrophic lateral sclerosis (ALS) in the Caucasian population in 2015. Here, we sequenced for TBK1 variants in a cohort of 15 familial ALS (fALS) and 275 sporadic ALS (sALS) of Chinese origin by targeted next-generation sequencing. We identified one likely benign missense variant (p. Ser398Pro), two missense variants of uncertain significance (p. Ile37Leu and p. Tyr677Asn), and two novel heterozygous variants in introns of TBK1, c.1522-3T > G and c.2066 + 4A > G. We performed splicing assays through minigene plasmids and RNA pull-down assay to determine that the two substitutions of nucleotides disrupted the binding of the important splicing regulator hnRNPA1 and promoted aberrant pre-mRNA splicing modes. The c.1522-3T > G variant promoted nearly 50.0% of abnormal transcripts (3 different types of insertions and deletions (indels) in junction of intron 13-exon 14) and the c.2066 + 4A > G variant inhibited about 75.0% inclusion of exon 19, both causing premature stop codon and producing TBK1 protein without CCD2. Immunofluorescence analysis showed that the expression of TBK1 with intronic variants was lower since less TBK1 distribution was observed in HEK293T cells. Both patients carrying TBK1 c.1522-3T > G and c.2066 + 4A > G variants developed a rapidly progressive ALS, with a survival of 31 and 10 months, respectively. The frequency of loss of function (LoF) variants in TBK1 was 0.73% in sALS in our cohort. We emphasize that intronic sequencing and pre-mRNA splicing analysis cannot be ignored to demonstrate the complex mutational spectrum and pathogenesis of ALS.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

First Identification of Rare Exonic and Deep Intronic Splice-Altering Variants in Patients With Beta-Sarcoglycanopathy

BACKGROUND

The precise genetic diagnosis of a sarcoglycanopathy or dystrophinopathy is sometimes extremely challenging, as pathogenic non-coding variants and/or complex structural variants do exist in DMD or sarcoglycan genes. This study aimed to determine the genetic diagnosis of three patients from two unrelated families with a suspected sarcoglycanopathy or dystrophinopathy based on their clinical, radiological, and pathological features, for whom routine genomic detection approaches failed to yield a definite genetic diagnosis.

METHODS

Muscle-derived reverse transcription-polymerase chain reaction analysis and/or TA cloning of DMD, SGCA, SGCB, SGCD, and SGCG mRNA were performed to identify aberrant transcripts. Genomic Sanger sequencing around the aberrant transcripts was performed to detect possible splice-altering variants. Bioinformatic and segregation studies of the detected genomic variants were performed in both families.

RESULTS

In patients F1-II1 and F1-II2, we identified two novel pathogenic compound heterozygous variants in SGCB. One is a deep intronic splice-altering variant (DISV), c.243 + 1558C > T in intron 2 causing the activation of an 87-base pair (bp) pseudoexon, and the other one is a non-canonical splicing site variant, c.243 + 6T > A leading to the partial intron inclusion of 10-bp sequence. A novel DISV, c.243 + 1576C > G causing a 106-bp pseudoexon activation, and a nonsense variant in SGCB were identified in compound heterozygous state in patient F2-II1. Unexpectedly, the predicted nonsense variant, c.334C > T in exon 3, created a new donor splice site in exon 3 that was stronger than the natural one, resulting in a 97-bp deletion of exon 3 (r.333_429del).

CONCLUSION

This is the first identification of rare exonic and DISVs in the SGCB gene.

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.

Splicing Characteristics of Dystrophin Pseudoexons and Identification of a Novel Pathogenic Intronic Variant in the DMD Gene

Pseudoexon (PE) inclusion has been implicated in various dystrophinopathies; however, its splicing characteristics have not been fully investigated. This study aims to analyze the splicing characteristics of dystrophin PEs and compare them with those of dystrophin canonical exons (CEs). Forty-two reported dystrophin PEs were divided into a splice site (ss) group and a splicing regulatory element (SRE) group. Five dystrophin PEs with characteristics of poison exons were identified and categorized as the possible poison exon group. The comparative analysis of each essential splicing signal among different groups of dystrophin PEs and dystrophin CEs revealed that the possible poison exon group had a stronger 3′ ss compared to any other group. As for auxiliary SREs, different groups of dystrophin PEs were found to have a smaller density of diverse types of exonic splicing enhancers and a higher density of several types of exonic splicing silencers compared to dystrophin CEs. In addition, the possible poison exon group had a smaller density of 3′ ss intronic splicing silencers compared to dystrophin CEs. To our knowledge, our findings indicate for the first time that poison exons might exist in DMD (the dystrophin gene) and present with different splicing characteristics than other dystrophin PEs and CEs.