Functional assays of non-canonical splice-site variants in inherited retinal dystrophies genes

Inherited retinal dystrophies are a group of disorders characterized by the progressive degeneration of photoreceptors leading to loss of the visual function and eventually to legal blindness. Although next generation sequencing (NGS) has revolutionized the molecular diagnosis of these diseases, the pathogenicity of some mutations casts doubts. After the screening of 208 patients with a panel of 117 genes, we obtained 383 variants that were analysed in silico with bioinformatic prediction programs. Based on the results of these tools, we selected 15 variants for their functional assessment. Therefore, we carried out minigene assays to unveil whether they could affect the splicing of the corresponding gene. As a whole, seven variants were found to induce aberrant splicing in the following genes: BEST1, CACNA2D4, PRCD, RIMS1, FSCN2, MERTK and MAK. This study shows the efficacy of a workflow, based on the association of the Minimum Allele Frequency, family co-segregation, in silico predictions and in vitro assays to determine the effect of potential splice site variants identified by DNA-based NGS. These findings improve the molecular diagnosis of inherited retinal dystrophies and will allow some patients to benefit from the upcoming gene-based therapeutic strategies.

A +3 variant at a donor splice site leads to a skipping of the MYH11 exon 32, a recurrent RNA defect causing Heritable Thoracic Aortic Aneurysm and Dissection and/or Patent Ductus Arteriosus

Background

Pathogenic variants in MYH11 are associated with either heritable thoracic aortic aneurysm and dissection (HTAAD), patent ductus arteriosus (PDA) syndrome, or megacystis‐microcolon‐intestinal hypoperistalsis syndrome (MMIHS).

Methods and Results

We report a family referred for molecular diagnosis with HTAAD/PDA phenotype in which we found a variant at a non‐conserved position of the 5’ donor splice site of intron 32 of MYH11 potentially altering splicing (NM_002474.3:c.4578+3A>C). Although its cosegregation with disease was observed, it remained of unknown significance. Later, aortic surgery in the proband gave us the opportunity to perform a transcript analysis. This showed a skipping of the exon 32, an RNA defect previously reported to be translated to an in‐frame loss of 71 amino acids and a dominant‐negative effect in the smooth muscle myosin rod. This RNA defect is also reported in 3 other HTAAD/PDA pedigrees.

Conclusion

This report confirms that among rare variants in MYH11, skipping of exon 32 is recurrent. This finding is of particular interest to establish complex genotype–phenotype correlations where some alleles are associated with autosomal dominant HTAAD/PDA, while others result in recessive or dominant visceral myopathies.

A Reporter Based Cellular Assay for Monitoring Splicing Efficiency

During gene expression, the vital step of pre-mRNA splicing involves accurate recognition of splice sites and efficient assembly of spliceosomal complexes to join exons and remove introns prior to cytoplasmic export of the mature mRNA. Splicing efficiency can be altered by the presence of mutations at splice sites, the influence of trans-acting splicing factors, or the activity of therapeutics. Here, we describe the protocol for a cellular assay that can be applied for monitoring the splicing efficiency of any given exon. The assay uses an adaptable plasmid encoded 3-exon/2-intron minigene reporter, which can be expressed in mammalian cells by transient transfection. Post-transfection, total cellular RNA is isolated, and the efficiency of exon splicing in the reporter mRNA is determined by either primer extension or semi-quantitative reverse transcriptase-polymerase chain reaction (RT-PCR). We describe how the impact of disease associated 5' splice-site mutations can be determined by introducing them in the reporter; and how the suppression of these mutations can be achieved by co-transfection with U1 small nuclear RNA (snRNA) construct carrying compensatory mutations in its 5' region that basepairs with the 5'-splice sites at exon-intron junctions in pre-mRNAs. Thus, the reporter can be used for the design of therapeutic U1 particles to improve recognition of mutant 5' splice-sites. Insertion of cis-acting regulatory sites, such as splicing enhancer or silencer sequences, into the reporter can also be used to examine the role of U1 snRNP in regulation mediated by a specific alternative splicing factor. Finally, reporter expressing cells can be incubated with small molecules to determine the effect of potential therapeutics on constitutive pre-mRNA splicing or on exons carrying mutant 5' splice sites. Overall, the reporter assay can be applied to monitor splicing efficiency in a variety of conditions to study fundamental splicing mechanisms and splicing-associated diseases.

Association of Messenger RNA Level With Phenotype in Patients With Choroideremia: Potential Implications for Gene Therapy Dose

Importance

Gene therapy is a promising treatment for choroideremia, an X-linked retinal degeneration. The required minimum level of gene expression to ameliorate degeneration rate is unknown. This can be interrogated by exploring the association between messenger RNA (mRNA) levels and phenotype in mildly affected patients with choroideremia.

Objective

To analyze CHM mRNA splicing outcomes in 2 unrelated patients with the same c.940+3delA CHM splice site variant identified as mildly affected from a previous study of patients with choroideremia.

Design, Setting, and Participants

In this retrospective observational case series, 2 patients with c.940+3delA CHM variants treated at a single tertiary referral center were studied. In addition, a third patient with a c.940+2T>A variant that disrupts the canonical dinucleotide sequence at the same donor site served as a positive control. Data were collected from October 2013 to July 2018.

Main Outcomes and Measures

Central area of residual fundus autofluorescence was used as a biomarker for disease progression. CHM transcript splicing was assessed by both end point and quantitative polymerase chain reaction. Rab escort protein 1 (REP1) expression was assessed by immunoblot.

Results

The 2 mildly affected patients with c.940+3delA variants had large areas of residual autofluorescence for their age and longer degeneration half-lives compared with the previous cohort of patients with choroideremia. The control patient with a c.940+2T>A variant had a residual autofluorescence area within the range expected for his age. Both patients with the c.940+3delA variant expressed residual levels of full-length CHM mRNA transcripts relative to the predominant truncated transcript (mean [SEM] residual level: patient 1, 2.3% [0.3]; patient 2, 4.7% [0.2]), equivalent to approximately less than 1% of the level of full-length CHM expressed in nonaffected individuals. Full-length CHM expression was undetectable in the control patient. REP1 expression was less than the threshold for detection both in patients 1 and 2 and the control patient compared with wild-type controls.

Conclusions and Relevance

These results demonstrate the first genotype-phenotype association in choroideremia. A +3 deletion in intron 7 is sufficient to cause choroideremia in a milder form. If replicated with gene therapy, these findings would suggest that relatively low expression (less than 1%) of the wild-type levels of mRNA would be sufficient to slow disease progression.

The position of nonsense mutations can predict the phenotype severity: A survey on the DMD gene

A nonsense mutation adds a premature stop signal that hinders any further translation of a protein-coding gene, usually resulting in a null allele. To investigate the possible exceptions, we used the DMD gene as an ideal model. First, because dystrophin absence causes Duchenne muscular dystrophy (DMD), while its reduction causes Becker muscular dystrophy (BMD). Second, the DMD gene is X-linked and there is no second allele that can interfere in males. Third, databases are accumulating reports on many mutations and phenotypic data. Finally, because DMD mutations may have important therapeutic implications. For our study, we analyzed large databases (LOVD, HGMD and ClinVar) and literature and revised critically all data, together with data from our internal patients. We totally collected 2593 patients. Positioning these mutations along the dystrophin transcript, we observed a nonrandom distribution of BMD-associated mutations within selected exons and concluded that the position can be predictive of the phenotype. Nonsense mutations always cause DMD when occurring at any point in fifty-one exons. In the remaining exons, we found milder BMD cases due to early 5’ nonsense mutations, if reinitiation can occur, or due to late 3’ nonsense when the shortened product retains functionality. In the central part of the gene, all mutations in some in-frame exons, such as in exons 25, 31, 37 and 38 cause BMD, while mutations in exons 30, 32, 34 and 36 cause DMD. This may have important implication in predicting the natural history and the efficacy of therapeutic use of drug-stimulated translational readthrough of premature termination codons, also considering the action of internal natural rescuers. More in general, our survey confirm that a nonsense mutation should be not necessarily classified as a null allele and this should be considered in genetic counselling.

Combined genetic approaches yield a 48% diagnostic rate in a large cohort of French hearing-impaired patients

Hearing loss is the most common sensory disorder and because of its high genetic heterogeneity, implementation of Massively Parallel Sequencing (MPS) in diagnostic laboratories is greatly improving the possibilities of offering optimal care to patients. We present the results of a two-year period of molecular diagnosis that included 207 French families referred for non-syndromic hearing loss. Our multi-step strategy involved (i) DFNB1 locus analysis, (ii) MPS of 74 genes, and (iii) additional approaches including Copy Number Variations, in silico analyses, minigene studies coupled when appropriate with complete gene sequencing, and a specific assay for STRC. This comprehensive screening yielded an overall diagnostic rate of 48%, equally distributed between DFNB1 (24%) and the other genes (24%). Pathogenic genotypes were identified in 19 different genes, with a high prevalence of GJB2, STRC, MYO15A, OTOF, TMC1, MYO7A and USH2A. Involvement of an Usher gene was reported in 16% of the genotyped cohort. Four de novo variants were identified. This study highlights the need to develop several molecular approaches for efficient molecular diagnosis of hearing loss, as this is crucial for genetic counselling, audiological rehabilitation and the detection of syndromic forms.

Intron retention resulting from a silent mutation in the VWF gene that structurally influences the 5' splice site

Disease-associated silent mutations are considered to affect the accurate pre-messenger RNA (mRNA) splicing either by influencing regulatory elements, leading to exon skipping, or by creating a new cryptic splice site. This study describes a new molecular pathological mechanism by which a silent mutation inhibits splicing and leads to intron retention. We identified a heterozygous silent mutation, c.7464C>T, in exon 44 of the von Willebrand factor (VWF) gene in a family with type 1 von Willebrand disease. In vivo and ex vivo transcript analysis revealed an aberrantly spliced transcript, with intron 44 retained in the mRNA, implying disruption of the first catalytic step of splicing at the 5' splice site (5'ss). The abnormal transcript with the retained intronic region coded a truncated protein that lacked the carboxy-terminal end of the VWF protein. Confocal immunofluorescence characterizations of blood outgrowth endothelial cells derived from the patient confirmed the presence of the truncated protein by demonstrating accumulation of VWF in the endoplasmic reticulum. In silico pre-mRNA secondary and tertiary structure analysis revealed that this substitution, despite its distal position from the 5'ss (85 bp downstream), induces cis alterations in pre-mRNA structure that result in the formation of a stable hairpin at the 5'ss. This hairpin sequesters the 5'ss residues involved in U1 small nuclear RNA interactions, thereby inhibiting excision of the pre-mRNA intronic region. This study is the first to show the allosteric-like/far-reaching effect of an exonic variation on pre-mRNA splicing that is mediated by structural changes in the pre-mRNA.

Whole USH2A Gene Sequencing Identifies Several New Deep Intronic Mutations

Deep intronic mutations leading to pseudoexon (PE) insertions are underestimated and most of these splicing alterations have been identified by transcript analysis, for instance, the first deep intronic mutation in USH2A, the gene most frequently involved in Usher syndrome type II (USH2). Unfortunately, analyzing USH2A transcripts is challenging and for 1.8%-19% of USH2 individuals carrying a single USH2A recessive mutation, a second mutation is yet to be identified. We have developed and validated a DNA next-generation sequencing approach to identify deep intronic variants in USH2A and evaluated their consequences on splicing. Three distinct novel deep intronic mutations have been identified. All were predicted to affect splicing and resulted in the insertion of PEs, as shown by minigene assays. We present a new and attractive strategy to identify deep intronic mutations, when RNA analyses are not possible. Moreover, the bioinformatics pipeline developed is independent of the gene size, implying the possible application of this approach to any disease-linked gene. Finally, an antisense morpholino oligonucleotide tested in vitro for its ability to restore splicing caused by the c.9959-4159A>G mutation provided high inhibition rates, which are indicative of its potential for molecular therapy.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

FUBP1: a new protagonist in splicing regulation of the DMD gene

We investigated the molecular mechanisms for in-frame skipping of DMD exon 39 caused by the nonsense c.5480T>A mutation in a patient with Becker muscular dystrophy. RNase-assisted pull down assay coupled with mass spectrometry revealed that the mutant RNA probe specifically recruits hnRNPA1, hnRNPA2/B1 and DAZAP1. Functional studies in a human myoblast cell line transfected with DMD minigenes confirmed the splicing inhibitory activity of hnRNPA1 and hnRNPA2/B1, and showed that DAZAP1, also known to activate splicing, acts negatively in the context of the mutated exon 39. Furthermore, we uncovered that recognition of endogenous DMD exon 39 in muscle cells is promoted by FUSE binding protein 1 (FUBP1), a multifunctional DNA- and RNA-binding protein whose role in splicing is largely unknown. By serial deletion and mutagenesis studies in minigenes, we delineated a functional intronic splicing enhancer (ISE) in intron 38. FUBP1 recruitment to the RNA sequence containing the ISE was established by RNA pull down and RNA EMSA, and further confirmed by RNA-ChIP on endogenous DMD pre-mRNA. This study provides new insights about the splicing regulation of DMD exon 39, highlighting the emerging role of FUBP1 in splicing and describing the first ISE for constitutive exon inclusion in the mature DMD transcript.

Stem-loop 4 of U1 snRNA is essential for splicing and interacts with the U2 snRNP-specific SF3A1 protein during spliceosome assembly

The pairing of 5′ and 3′ splice sites across an intron is a critical step in spliceosome formation and its regulation. Sharma et al. report a new interaction between stem–loop 4 (SL4) of the U1 snRNA, which recognizes the 5′ splice, and a component of the U2 snRNP complex, which assembles across the intron at the 3′ splice site. U1-SL4 interacts with the SF3A1 protein of the U2 snRNP, and this interaction occurs within prespliceosomal complexes assembled on the pre-mRNA.

The pairing of 5′ and 3′ splice sites across an intron is a critical step in spliceosome formation and its regulation. Interactions that bring the two splice sites together during spliceosome assembly must occur with a high degree of specificity and fidelity to allow expression of functional mRNAs and make particular alternative splicing choices. Here, we report a new interaction between stem–loop 4 (SL4) of the U1 snRNA, which recognizes the 5′ splice site, and a component of the U2 small nuclear ribonucleoprotein particle (snRNP) complex, which assembles across the intron at the 3′ splice site. Using a U1 snRNP complementation assay, we found that SL4 is essential for splicing in vivo. The addition of free U1-SL4 to a splicing reaction in vitro inhibits splicing and blocks complex assembly prior to formation of the prespliceosomal A complex, indicating a requirement for a SL4 contact in spliceosome assembly. To characterize the interactions of this RNA structure, we used a combination of stable isotope labeling by amino acids in cell culture (SILAC), biotin/Neutravidin affinity pull-down, and mass spectrometry. We show that U1-SL4 interacts with the SF3A1 protein of the U2 snRNP. We found that this interaction between the U1 snRNA and SF3A1 occurs within prespliceosomal complexes assembled on the pre-mRNA. Thus, SL4 of the U1 snRNA is important for splicing, and its interaction with SF3A1 mediates contact between the 5′ and 3′ splice site complexes within the assembling spliceosome.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Experience of targeted Usher exome sequencing as a clinical test

We show that massively parallel targeted sequencing of 19 genes provides a new and reliable strategy for molecular diagnosis of Usher syndrome (USH) and nonsyndromic deafness, particularly appropriate for these disorders characterized by a high clinical and genetic heterogeneity and a complex structure of several of the genes involved. A series of 71 patients including Usher patients previously screened by Sanger sequencing plus newly referred patients was studied. Ninety-eight percent of the variants previously identified by Sanger sequencing were found by next-generation sequencing (NGS). NGS proved to be efficient as it offers analysis of all relevant genes which is laborious to reach with Sanger sequencing. Among the 13 newly referred Usher patients, both mutations in the same gene were identified in 77% of cases (10 patients) and one candidate pathogenic variant in two additional patients. This work can be considered as pilot for implementing NGS for genetically heterogeneous diseases in clinical service.

A classification model relative to splicing for variants of unknown clinical significance: application to the CFTR gene

Molecular diagnosis of cystic fibrosis and cystic fibrosis transmembrane regulator (CFTR)-related disorders led to the worldwide identification of nearly 1,900 sequence variations in the CFTR gene that consist mainly of private point mutations and small insertions/deletions. Establishing their effect on the function of the encoded protein and therefore their involvement in the disease is still challenging and directly impacts genetic counseling. In this context, we built a decision tree following the international guidelines for the classification of variants of unknown clinical significance (VUCS) in the CFTR gene specifically focused on their consequences on splicing. We applied general and specific criteria, including comprehensive review of literature and databases, familial genetics data, and thorough in silico studies. This model was tested on 15 intronic and exonic VUCS identified in our cohort. Six variants were classified as probably nonpathogenic considering their impact on splicing and eight as probably pathogenic, which include two apparent missense mutations. We assessed the validity of our method by performing minigenes studies and confirmed that 93% (14/15) were correctly classified. We provide in this study a high-performance method that can play a full role in interpreting the results of molecular diagnosis in emergency context, when functional studies are not achievable.

Overexpression of adapted U1snRNA in patients' cells to correct a 5' splice site mutation in propionic acidemia

Splicing defects account for 16% of the mutant alleles in the PCCA and PCCB genes, encoding both subunits of the propionyl-CoA carboxylase (PCC) enzyme, defective in propionic acidemia, one of the most frequent organic acidemias causing variable neurological impairment. Most of the splicing mutations identified affect the conserved 3' splice (3' ss) or 5' splice (5' ss) sites, the latter predictably through lowering the strength of base pairing with U1snRNA. Among the 5' ss mutations we have focused on the c.1209+3A>G (IVS13+3A>G) mutation in the PCCA gene, identified in four patients (three homozygous and one heterozygous) of common geographical origin and causing exon 13 skipping. To study the potential of splicing modulation to restore PCC function, we analyzed the effect of transient transfections in patients' cells with modified U1snRNA adapted to compensate the mutant change and other mismatches at different positions of the 5' ss. Using this strategy normal transcript could be efficiently recovered with the concomitant disappearance of the aberrant exon skipping transcript, as observed after standard RT-PCR and sequence analysis or using fluorescent primers and semiquantitative RT-PCR. Different efficiencies with up to 100% exon inclusion were observed depending on the transfection conditions and specifically on the adapted U1snRNA used, confirming previously reported dependencies between nucleotides at the 5' ss for its correct recognition by the spliceosome. The reversal of the splicing defect did not result in a significant increase in enzyme activity, suggesting other factors must be taken into account for the application of overexpression of splice factors such as U1 as therapeutic strategy for splice defects.

Ex vivo splicing assays of mutations at noncanonical positions of splice sites in USHER genes

Molecular diagnosis in Usher syndrome type 1 and 2 patients led to the identification of 21 sequence variations located in noncanonical positions of splice sites in MYO7A, CDH23, USH1C, and USH2A genes. To establish experimentally the splicing pattern of these substitutions, whose impact on splicing is not always predictable by available softwares, ex vivo splicing assays were performed. The branch-point mapping strategy was also used to investigate further a putative branch-point mutation in USH2A intron 43. Aberrant splicing was demonstrated for 16 of the 21 (76.2%) tested sequence variations. The mutations resulted more frequently in activation of a nearby cryptic splice site or use of a de novo splice site than exon skipping (37.5%). This study allowed the reclassification as splicing mutations of one silent (c.7872G>A (p.Glu2624Glu) in CDH23) and four missense mutations (c.2993G>A (p.Arg998Lys) in USH2A, c.592G>A (p.Ala198Thr), c.3503G>C [p.Arg1168Pro], c.5944G>A (p.Gly1982Arg) in MYO7A), whereas it provided clues about a role in structure/function in four other cases: c.802G>A (p.Gly268Arg), c.653T>A (p.Val218Glu) (USH2A), and c.397C>T (p.His133Tyr), c.3502C>T (p.Arg1168Trp) (MYO7A). Our data provide insights into the contribution of splicing mutations in Usher genes and illustrate the need to define accurately their splicing outcome for diagnostic purposes.