Experimental assessment of splicing variants using expression minigenes and comparison with in silico predictions

Identified five variants in CFTR gene that alter RNA splicing by minigene assay

Background

Cystic fibrosis (CF) is a common monogenic multisystem disease caused primarily by variants in the CFTR gene. Emerging evidence suggests that some variants, which are described as missense, synonymous or nonsense variants in the literature or databases, may be deleterious by affecting the pre-mRNA splicing process.

Methods

We analyzed 27 exonic variants in the CFTR gene utilizing bioinformatics tools and identified candidate variants that could lead to splicing changes through minigene assays. Ultimately, we selected eight candidate variants to assess their effects on pre-mRNA splicing. The numbering of DNA variants is based on the complementary DNA (cDNA)sequence of CFTR (Ref Seq NM_000492.4).

Results

This study assessed the impact of CFTR variants on exon splicing by combining predictive bioinformatics tools with minigene assays. Among the eight candidate single nucleotide alterations, five variants (c.488A>T,c.1117G>T, c.1209G>T, c.3239A>G and c.3367G>C) were identified as causing exon skipping.

Conclusion

Our study employed a minigene system, which offers great flexibility for assessing aberrant splicing patterns when patient mRNA samples are not accessible, to investigate the effects of exonic variants on pre-mRNA splicing. Our experimental outcomes highlight the importance of analyzing exonic variations at the mRNA level.

A De Novo Noncoding RARB Variant Associated with Complex Microphthalmia Alters a Putative Regulatory Element

Retinoic acid receptor beta (RARB) is a transcriptional regulator crucial for coordinating retinoic acid- (RA-) mediated morphogenic movements, cell growth, and differentiation during eye development. Loss- or gain-of-function RARB coding variants have been associated with microphthalmia, coloboma, and anterior segment defects. We identified a de novo variant c.157+1895G>A located within a conserved region (CR1) in the first intron of RARB in an individual with complex microphthalmia and significant global developmental delay. Based on the phenotypic overlap, we further investigated the possible effects of the variant on mRNA splicing and/or transcriptional regulation through in silico and functional studies. In silico analysis identified the possibility of alternative splicing, suggested by one out of three (HSF, SpliceAI, and MaxEntScan) splicing prediction programs, and a strong indication of regulatory function based on publicly available DNase hypersensitivity, histone modification, chromatin folding, and ChIP-seq data sets. Consistent with the predictions of SpliceAI and MaxEntScan, in vitro minigene assays showed no effect on RARB mRNA splicing. Evaluation of CR1 for a regulatory role using luciferase reporter assays in human lens epithelial cells demonstrated a significant increase in the activity of the RARB promoter in the presence of wild-type CR1. This activity was further significantly increased in the presence of CR1 carrying the c.157+1895G>A variant, suggesting that the variant may promote RARB overexpression in human cells. Induction of RARB overexpression in human lens epithelial cells resulted in increased cell proliferation and elevated expression of FOXC1, a known downstream target of RA signaling and a transcription factor whose down- and upregulation is associated with ocular phenotypes overlapping the RARB spectrum. These results support a regulatory role for the CR1 element and suggest that the de novo c.157+1895G>A variant affecting this region may alter the proper regulation of RARB and, as a result, its downstream genes, possibly leading to abnormal development.

Protospacer modification improves base editing of a canonical splice site variant and recovery of CFTR function in human airway epithelial cells

Canonical splice site variants affecting the 5' GT and 3' AG nucleotides of introns result in severe missplicing and account for about 10% of disease-causing genomic alterations. Treatment of such variants has proven challenging due to the unstable mRNA or protein isoforms that typically result from disruption of these sites. Here, we investigate CRISPR-Cas9-mediated adenine base editing for such variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. We validate a CFTR expression minigene (EMG) system for testing base editing designs for two different targets. We then use the EMG system to test non-standard single-guide RNAs with either shortened or lengthened protospacers to correct the most common cystic fibrosis-causing variant in individuals of African descent (c.2988+1G>A). Varying the spacer region length allowed placement of the editing window in a more efficient context and enabled use of alternate protospacer adjacent motifs. Using these modifications, we restored clinically significant levels of CFTR function to human airway epithelial cells from two donors bearing the c.2988+1G>A variant.

Functional restoration of a CFTR splicing mutation through RNA delivery of CRISPR adenine base editor

Cystic fibrosis (CF) is a genetic disease caused by mutations in the CF transmembrane conductance regulator (CFTR) gene. The 2789+5G>A CFTR mutation is a quite frequent defect causing an aberrant splicing and a non-functional CFTR protein. Here we used a CRISPR adenine base editing (ABE) approach to correct the mutation in the absence of DNA double-strand breaks (DSB). To select the strategy, we developed a minigene cellular model reproducing the 2789+5G>A splicing defect. We obtained up to 70% editing in the minigene model by adapting the ABE to the PAM sequence optimal for targeting 2789+5G>A with a SpCas9-NG (NG-ABE). Nonetheless, the on-target base correction was accompanied by secondary (bystander) A-to-G conversions in nearby nucleotides, which affected the wild-type CFTR splicing. To decrease the bystander edits, we used a specific ABE (NG-ABEmax), which was delivered as mRNA. The NG-ABEmax RNA approach was validated in patient-derived rectal organoids and bronchial epithelial cells showing sufficient gene correction to recover the CFTR function. Finally, in-depth sequencing revealed high editing precision genome-wide and allele-specific correction. Here we report the development of a base editing strategy to precisely repair the 2789+5G>A mutation resulting in restoration of the CFTR function, while reducing bystander and off-target activities.

Q1291H-CFTR molecular dynamics simulations and ex vivo theratyping in nasal epithelial models and clinical response to elexacaftor/tezacaftor/ivacaftor in a Q1291H/F508del patient

Background: Cystic fibrosis (CF) is caused by a wide spectrum of mutations in the CF transmembrane conductance regulator (CFTR) gene, with some leading to non-classical clinical presentations. We present an integrated in vivo, in silico and in vitro investigation of an individual with CF carrying the rare Q1291H-CFTR allele and the common F508del allele. At age 56 years, the participant had obstructive lung disease and bronchiectasis, qualifying for Elexacaftor/Tezacaftor/Ivacaftor (ETI) CFTR modulator treatment due to their F508del allele. Q1291H CFTR incurs a splicing defect, producing both a normally spliced but mutant mRNA isoform and a misspliced isoform with a premature termination codon, causing nonsense mediated decay. The effectiveness of ETI in restoring Q1291H-CFTR is largely unknown. Methods: We collected clinical endpoint measurements, including forced expiratory volume in 1 s percent predicted (FEV1pp) and body mass index (BMI), and examined medical history. In silico simulations of the Q1291H-CFTR were compared to Q1291R, G551D, and wild-type (WT)-CFTR. We quantified relative Q1291H CFTR mRNA isoform abundance in patient-derived nasal epithelial cells. Differentiated pseudostratified airway epithelial cell models at air liquid interface were created and ETI treatment impact on CFTR was assessed by electrophysiology assays and Western blot. Results: The participant ceased ETI treatment after 3 months due to adverse events and no improvement in FEV1pp or BMI. In silico simulations of Q1291H-CFTR identified impairment of ATP binding similar to known gating mutants Q1291R and G551D-CFTR. Q1291H and F508del mRNA transcripts composed 32.91% and 67.09% of total mRNA respectively, indicating 50.94% of Q1291H mRNA was misspliced and degraded. Mature Q1291H-CFTR protein expression was reduced (3.18% ± 0.60% of WT/WT) and remained unchanged with ETI. Baseline CFTR activity was minimal (3.45 ± 0.25 μA/cm²) and not enhanced with ETI (5.73 ± 0.48 μA/cm²), aligning with the individual's clinical evaluation as a non-responder to ETI. Conclusion: The combination of in silico simulations and in vitro theratyping in patient-derived cell models can effectively assess CFTR modulator efficacy for individuals with non-classical CF manifestations or rare CFTR mutations, guiding personalized treatment strategies and optimizing clinical outcomes.

Downstream Alternate Start Site Allows N-Terminal Nonsense Variants to Escape NMD and Results in Functional Recovery by Readthrough and Modulator Combination

Genetic variants that introduce premature termination codons (PTCs) have remained difficult to therapeutically target due to lack of protein product. Nonsense mediated mRNA decay (NMD) targets PTC-bearing transcripts to reduce the potentially damaging effects of truncated proteins. Readthrough compounds have been tested on PTC-generating variants in attempt to permit translation through a premature stop. However, readthrough compounds have not proved efficacious in a clinical setting due to lack of stable mRNA. Here, we investigate N-terminal variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which have been shown to escape NMD, potentially through a mechanism of alternative translation initiation at downstream AUG codons. We hypothesized that N-terminal variants in CFTR that evade NMD will produce stable transcript, allowing CFTR function to be restored by a combination of readthrough and protein modulator therapy. We investigate this using two cell line models expressing CFTR-expression minigenes (EMG; HEK293s and CFBEs) and primary human nasal epithelial (NE) cells, and we test readthrough compounds G418 and ELX-02 in combination with CFTR protein modulators. HEK293 cells expressing the variants E60X and L88X generate CFTR-specific core glycosylated products that are consistent with downstream translation initiation. Mutation of downstream methionines at codons 150 and 152 does not result in changes in CFTR protein processing in cells expressing L88X-CFTR-EMG. However, mutation of methionine at 265 results in loss of detectable CFTR protein in cells expressing E60X, L88X, and Y122X CFTR-EMGs, indicating that downstream translation initiation is occurring at the AUG codon at position M265. In HEK293 stable cells harboring L88X, treatment with readthrough compounds alone allows for formation of full-length, but misfolded CFTR protein. Upon addition of protein modulators in combination with readthrough, we observe formation of mature, complex-glycosylated CFTR. In CFBE and NE cells, addition of readthrough ELX-02 and modulator therapy results in substantial recovery of CFTR function. Our work indicates that N-terminal variants generate stable CFTR transcript due to translation initiation at a downstream AUG codon. Thus, individuals with CF bearing 5′ nonsense variants that evade NMD are ideal candidates for treatment with clinically safe readthrough compounds and modulator therapy.

One Size Does Not Fit All: The Past, Present and Future of Cystic Fibrosis Causal Therapies

Cystic fibrosis (CF) is the most common monogenic disorder, caused by mutations in the CF transmembrane conductance regulator (CFTR) gene. Over the last 30 years, tremendous progress has been made in understanding the molecular basis of CF and the development of treatments that target the underlying defects in CF. Currently, a highly effective CFTR modulator treatment (Kalydeco™/Trikafta™) is available for 90% of people with CF. In this review, we will give an extensive overview of past and ongoing efforts in the development of therapies targeting the molecular defects in CF. We will discuss strategies targeting the CFTR protein (i.e., CFTR modulators such as correctors and potentiators), its cellular environment (i.e., proteostasis modulation, stabilization at the plasma membrane), the CFTR mRNA (i.e., amplifiers, nonsense mediated mRNA decay suppressors, translational readthrough inducing drugs) or the CFTR gene (gene therapies). Finally, we will focus on how these efforts can be applied to the 15% of people with CF for whom no causal therapy is available yet.

Modulation of RNA Splicing by Oligonucleotides: Mechanisms of Action and Therapeutic Implications

Dysregulation of RNA splicing causes many diseases and disorders. Several therapeutic approaches have been developed to correct aberrant alternative splicing events for the treatment of cancers and hereditary diseases, including gene therapy and redirecting splicing, using small molecules or splice switching oligonucleotides (SSO). Significant advances in the chemistry and pharmacology of nucleic acid have led to the development of clinically approved SSO drugs for the treatment of spinal muscular dystrophy and Duchenne muscular dystrophy (DMD). In this review, we discuss the mechanisms of SSO action with emphasis on "less common" approaches to modulate alternative splicing, including bipartite and bifunctional SSO, oligonucleotide decoys for splice factors and SSO-mediated mRNA degradation via AS-NMD and NGD pathways. We briefly discuss the current progress and future perspectives of SSO therapy for rare and ultrarare diseases.

An oncogenic splice variant of PDGFRα in adult glioblastoma as a therapeutic target for selective CDK4/6 inhibitors

Understanding human genome alterations is necessary to optimize genome-based cancer therapeutics. However, some newly discovered mutations remain as variants of unknown significance (VUS). Here, the mutation c.1403A > G in exon 10 of the platelet-derived growth factor receptor-alpha (PDGFRA) gene, a VUS found in adult glioblastoma multiforme (GBM), was introduced in human embryonal kidney 293 T (HEK293T) cells using genome editing to investigate its potential oncogenic functions. Genome editing was performed using CRISPR/Cas9; the proliferation, drug sensitivity, and carcinogenic potential of genome-edited cells were investigated. We also investigated the mechanism underlying the observed phenotypes. Three GBM patients carrying the c.1403A > G mutation were studied to validate the in vitro results. The c.1403A > G mutation led to a splice variant (p.K455_N468delinsN) because of the generation of a 3’-acceptor splice site in exon 10. PDGFRA-mutated HEK293T cells exhibited a higher proliferative activity via PDGFRα and the cyclin-dependent kinase (CDK)4/CDK6-cyclin D1 signaling pathway in a ligand-independent manner. They showed higher sensitivity to multi-kinase, receptor tyrosine kinase, and CDK4/CDK6 inhibitors. Of the three GBM patients studied, two harbored the p.K455_N468delinsN splice variant. The splicing mutation c.1403A > G in PDGFRA is oncogenic in nature. Kinase inhibitors targeting PDGFRα and CDK4/CDK6 signaling should be evaluated for treating GBM patients harboring this mutation.

Functional assessment of somatic STK11 variants identified in primary human non-small cell lung cancers

Serine/Threonine Kinase 11 (STK11) encodes an important tumor suppressor that is frequently mutated in lung adenocarcinoma. Clinical studies have shown that mutations in STK11 resulting in loss of function correlate with resistance to anti-PD-1 monoclonal antibody therapy in KRAS-driven non-small cell lung cancer (NSCLC), but the molecular mechanisms responsible remain unclear. Despite this uncertainty, STK11 functional status is emerging as a reliable biomarker for predicting non-response to anti-PD-1 therapy in NSCLC patients. The clinical utility of this biomarker ultimately depends upon accurate classification of STK11 variants. For nonsense variants occurring early in the STK11 coding region, this assessment is straightforward. However, rigorously demonstrating the functional impact of missense variants remains an unmet challenge. Here we present data characterizing four STK11 splice-site variants by analyzing tumor mRNA, and 28 STK11 missense variants using an in vitro kinase assay combined with a cell-based p53-dependent luciferase reporter assay. The variants we report were identified in primary human NSCLC biopsies in collaboration with the University of Vermont Genomic Medicine group. Additionally, we compare our experimental results with data from 22 in silico predictive algorithms. Our work highlights the power, utility and necessity of functional variant assessment and will aid STK11 variant curation, provide a platform to assess novel STK11 variants and help guide anti-PD-1 therapy utilization in KRAS-driven NSCLCs.

Accurate assignment of disease liability to genetic variants using only population data

PURPOSE

The growing size of public variant repositories prompted us to test the accuracy of pathogenicity prediction of DNA variants using population data alone.

METHODS

Under the a priori assumption that the ratio of the prevalence of variants in healthy population vs that in affected populations form 2 distinct distributions (pathogenic and benign), we used a Bayesian method to assign probability to a variant belonging to either distribution.

RESULTS

The approach, termed Bayesian prevalence ratio (BayPR), accurately parsed 300 of 313 expertly curated CFTR variants: 284 of 296 pathogenic/likely pathogenic variants in 1 distribution and 16 of 17 benign/likely benign variants in another. BayPR produced an area under the receiver operating characteristic curve of 0.99 for 103 functionally confirmed missense CFTR variants, which is equal to or exceeds 10 commonly used algorithms (area under the receiver operating characteristic curve range = 0.54-0.99). Application of BayPR to expertly curated variants in 8 genes associated with 7 Mendelian conditions led to the assignment of a disease-causing probability of ≥80% to 1350 of 1374 (98.3%) pathogenic/likely pathogenic variants and of ≤20% to 22 of 23 (95.7%) benign/likely benign variants.

CONCLUSION

Irrespective of the variant type or functional effect, the BayPR approach provides probabilities of pathogenicity for DNA variants responsible for Mendelian disorders using only the variant counts in affected and unaffected population samples.

Identification of novel single-nucleotide variants altering RNA splicing of PKD1 and PKD2

The development of sequencing techniques identified numerous genetic variants, and accurate evaluation of the clinical significance of these variants facilitates the diagnosis of Mendelian diseases. In the present study, 549 rare single- nucleotide variants of uncertain significance were extracted from the ADPKD and ClinVar databases. MaxEntScan scoresplice is an in silico splicing prediction tool that was used to analyze rare PKD1 and PKD2 variants of unknown significance. An in vitro minigene splicing assay was used to verify 37 splicing-altering candidates that were located within seven residues of the splice donor sequence excluding canonical GT dinucleotides or within 21 residues of the acceptor sequence excluding canonical AG dinucleotides of PKD1 and PKD2. We demonstrated that eight PKD1 variants alter RNA splicing and were predicted to be pathogenic.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.

A deep learning approach to identify gene targets of a therapeutic for human splicing disorders

Pre-mRNA splicing is a key controller of human gene expression. Disturbances in splicing due to mutation lead to dysregulated protein expression and contribute to a substantial fraction of human disease. Several classes of splicing modulator compounds (SMCs) have been recently identified and establish that pre-mRNA splicing represents a target for therapy. We describe herein the identification of BPN-15477, a SMC that restores correct splicing of ELP1 exon 20. Using transcriptome sequencing from treated fibroblast cells and a machine learning approach, we identify BPN-15477 responsive sequence signatures. We then leverage this model to discover 155 human disease genes harboring ClinVar mutations predicted to alter pre-mRNA splicing as targets for BPN-15477. Splicing assays confirm successful correction of splicing defects caused by mutations in CFTR, LIPA, MLH1 and MAPT. Subsequent validations in two disease-relevant cellular models demonstrate that BPN-15477 increases functional protein, confirming the clinical potential of our predictions.

On the Corner of Models and Cure: Gene Editing in Cystic Fibrosis

Cystic fibrosis (CF) is a severe genetic disease for which curative treatment is still lacking. Next generation biotechnologies and more efficient cell-based and in vivo disease models are accelerating the development of novel therapies for CF. Gene editing tools, like CRISPR-based systems, can be used to make targeted modifications in the genome, allowing to correct mutations directly in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. Alternatively, with these tools more relevant disease models can be generated, which in turn will be invaluable to evaluate novel gene editing-based therapies for CF. This critical review offers a comprehensive description of currently available tools for genome editing, and the cell and animal models which are available to evaluate them. Next, we will give an extensive overview of proof-of-concept applications of gene editing in the field of CF. Finally, we will touch upon the challenges that need to be addressed before these proof-of-concept studies can be translated towards a therapy for people with CF.

DICER1-associated embryonal rhabdomyosarcoma and adenosarcoma of the gynecologic tract: Pathology, molecular genetics, and indications for molecular testing

Gynecologic sarcomas are uncommon neoplasms, the majority occurring in the uterus. Due to the diverse nature of these, the description of "new" morphological types and the rarity of some of them, pathological diagnosis and treatment is often challenging. Finding genetic alterations specific to, and frequently occurring, in a certain type can aid in the diagnosis. DICER1 is a highly conserved ribonuclease crucial in the biogenesis of microRNAs and mutations in DICER1 (either somatic or germline) have been detected in a wide range of sarcomas including genitourinary embryonal rhabdomyosarcomas (ERMS) and adenosarcomas. Importantly, DICER1-associated sarcomas share morphological features irrespective of the site of origin such that the pathologist can strongly suspect a DICER1 association. A review of the literature shows that almost all gynecologic ERMS reported (outside of the vagina) harbor DICER1 alterations, while approximately 20% of adenosarcomas also do so. These two tumor types exhibit significant morphological overlap and DICER1 tumor testing may be helpful in distinguishing between them, because a negative result makes ERMS unlikely. Given that germline pathogenic DICER1 variants are frequent in uterine (corpus and cervix) ERMS and pathogenic germline variants in this gene cause a hereditary cancer predisposition syndrome (DICER1 syndrome), patients diagnosed with these neoplasms should be referred to medical genetic services. Cooperation between pathologists and geneticists is crucial and will help in improving the diagnosis and management of these uncommon sarcomas.

Organoids as a personalized medicine tool for ultra-rare mutations in cystic fibrosis: The case of S955P and 1717-2A>G

BACKGROUND

For most of the >2000 CFTR gene variants reported, neither the associated disease liability nor the underlying basic defect are known, and yet these are essential for disease prognosis and CFTR-based therapeutics. Here we aimed to characterize two ultra-rare mutations - 1717-2A > G (c.1585-2A > G) and S955P (p.Ser955Pro) - as case studies for personalized medicine.

METHODS

Patient-derived rectal biopsies and intestinal organoids from two individuals with each of these mutations and F508del (p.Phe508del) in the other allele were used to assess CFTR function, response to modulators and RNA splicing pattern. In parallel, we used cellular models to further characterize S955P independently of F508del and to assess its response to CFTR modulators.

RESULTS

Results in both rectal biopsies and intestinal organoids from both patients evidence residual CFTR function. Further characterization shows that 1717-2A > G leads to alternative splicing generating <1% normal CFTR mRNA and that S955P affects CFTR gating. Finally, studies in organoids predict that both patients are responders to VX-770 alone and even more to VX-770 combined with VX-809 or VX-661, although to different levels.

CONCLUSION

This study demonstrates the high potential of personalized medicine through theranostics to extend the label of approved drugs to patients with rare mutations.

Evaluation of both exonic and intronic variants for effects on RNA splicing allows for accurate assessment of the effectiveness of precision therapies

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.

Genetic variants that impact pre-mRNA splicing are a common cause of genetic disease and have varying downstream molecular consequences. As a result, precision therapies that function at the protein level are not always effective for these variants and thus careful assessment is necessary. Here we evaluate RNA-level effects of 52 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene and show that study of splicing and its consequences allows for more accurate assignment of precision therapies.

Seven novel genetic variants in a North Indian cohort with classical homocystinuria

Classical homocystinuria is the most common cause of isolated homocystinuria. The variants of the CBS gene remain unidentified in Indian children with this disorder. Based on the hallmark clinical features, family history, and/or biochemical clues for classical homocystinuria, 16 children below the age of 18 years were evaluated by Sanger sequencing of the coding exons of CBS gene with flanking intronic regions. The common C677T variant of the MTHFR gene was also screened by restriction fragment length polymorphism. Fifteen children were clinically suspected of having classical homocystinuria and one asymptomatic child with positive family history. Only seven children had biochemical features of classical homocystinuria. Sanger sequencing of the CBS gene confirmed 15 different pathogenic or likely pathogenic variants in 14 cases. Of these, seven variants were novel (three frameshift deletions, two nonsense, one missense, one splice site variant) and were predicted to be deleterious by Mutation Taster software. Seven cases were homozygous, another six were compound heterozygous, and one case was single heterozygous in the study. None of the three most frequent mutations reported worldwide viz., I278T, G307S, and IVS 11-2A>C were found in our cohort. No variants were detected in the exons 2, 8, 12, and 14 as compared to reported literature. Eleven out of 15 variants were associated with the conserved catalytic domain of the CBS polypeptide. The MTHFR polymorphism C677T was observed in heterozygous state in six cases. Our study reports the detailed genotype and seven novel variants in the CBS gene, causing classical homocystinuria in Indian children. The genetic analysis will help to offer accurate genetic counseling, prenatal diagnosis, and development of mutation-based novel therapeutic strategies.

A novel splicing pathogenic variant in COL1A1 causing osteogenesis imperfecta (OI) type I in a Chinese family

Background

Osteogenesis imperfecta (OI), a rare autosomal inheritable disorder characterized by bone fragility and skeletal deformity, is caused by pathogenic variants in genes impairing the synthesis and processing of extracellular matrix protein collagen type I. With the use of next‐generation sequencing and panels approaches, an increasing number of OI patients can be confirmed and new pathogenic variants can be discovered. This study sought to identify pathogenic gene variants in a Chinese family with OI I.

Methods

Whole‐exome sequencing was used to identify pathogenic variants in the proband, which is confirmed by Sanger sequencing and cosegregation analysis; MES, HSF, and Spliceman were used to analyze this splicing variant；qRT‐PCR was performed to identify the mRNA expression level of COL1A1 in patient peripheral blood samples; Minigene splicing assay was performed to mimic the splicing process of COL1A1 variants in vitro; Analysis of evolutionary conservation of amino acid residues and structure prediction of the mutant protein.

Results

A novel splicing pathogenic variant (c.3814+1G>T) was identified in this OI family by using whole‐exome sequencing, Sanger sequencing, and cosegregation analysis. Sequencing of RT‐PCR products from the COL1A1 minigene variant reveals a 132‐nucleotide (nt) insertion exists at the junction between exons 48 and exon 49 of the COL1A1 cDNA. Splicing assay indicates that the mutated minigene produces an alternatively spliced transcript which may cause a frameshift resulting in early termination of protein expression. The molecular analysis suggested that the altered amino acid is located at the C‐terminus of type I procollagen.

Conclusion

Our study reveals the pathogenesis of a novel COL1A1 splicing pathogenic variant c.3814+1G>T in a Chinese family with OI I.

Molecular Diagnosis and Genetic Counseling of Cystic Fibrosis and Related Disorders: New Challenges

Identification of the cystic fibrosis transmembrane conductance regulator (CFTR) gene and its numerous variants opened the way to fantastic breakthroughs in diagnosis, research and treatment of cystic fibrosis (CF). The current and future challenges of molecular diagnosis of CF and CFTR-related disorders and of genetic counseling are here reviewed. Technological advances have enabled to make a diagnosis of CF with a sensitivity of 99% by using next generation sequencing in a single step. The detection of heretofore unidentified variants and ethnic-specific variants remains challenging, especially for newborn screening (NBS), CF carrier testing and genotype-guided therapy. Among the criteria for assessing the impact of variants, population genetics data are insufficiently taken into account and the penetrance of CF associated with CFTR variants remains poorly known. The huge diversity of diagnostic and genetic counseling indications for CFTR studies makes assessment of variant disease-liability critical. This is especially discussed in the perspective of wide genome analyses for NBS and CF carrier screening in the general population, as future challenges.

A non-invasive diagnostic assay for rapid detection and characterization of aberrant mRNA-splicing by nonsense mediated decay inhibition

BACKGROUND

Interpretation of genetic variants detected by sequencing of genomic DNA, which may cause splicing defects, regularly requires mRNA analysis. Usually, only bioinformatic testing is provided, because simple and non-invasive assay protocols are lacking. Furthermore, the detection of mis-splicing is often hampered by nonsense mediated mRNA decay (NMD).

METHODS

Starting from a case of Pompe disease with two potential splicing variants an assay for the analysis of splice defects in general was developed. We analyzed the transcripts from the gene of interest by standard methods after short-term culture of the patient's lymphocytes in the presence and absence of a NMD inhibitor. Variant and wild type transcript expression were quantified by allele specific PCR in the patient and both parents and the expression ratio with/without NMD inhibition was calculated for each transcript.

RESULTS

NMD detection in lymphocytes was optimized and evaluated by analyzing a naturally occurring NMD transcript. Several compounds inhibited NMD successfully, including potential therapeutic agents. Sample storage for up to 4 days at room temperature prior to lymphocyte isolation did not affect results. In a proof of concept we identified two candidate variants as severe splicing variants in a patient with Pompe disease, but the strategy can also be used to screen for any mis-spliced transcripts prone to NMD.

CONCLUSIONS

We developed a simple, non-invasive assay for the detection and characterization of potential splicing variants. This is essential, because early and near-term diagnosis and disease classification is required to facilitate therapy in many genetic diseases.

Nonsense-mediated RNA Decay Pathway Inhibition Restores Expression and Function of W1282X CFTR

The recessive genetic disease cystic fibrosis (CF) is caused by loss-of-function mutations in the CFTR (CF transmembrane conductance regulator) gene. Approximately 10% of patients with CF have at least one allele with a nonsense mutation in CFTR. Nonsense mutations generate premature termination codons that can subject mRNA transcripts to rapid degradation through the nonsense-mediated mRNA decay (NMD) pathway. Currently, there are no approved therapies that specifically target nonsense mutations in CFTR. Here, we identified antisense oligonucleotides (ASOs) that target the NMD factor SMG1 to inhibit the NMD pathway, and determined their effects on the W1282X CFTR mutation. First, we developed and validated two in vitro models of the W1282X CFTR mutation. Next, we treated these cells with antisense oligonucleotides to inhibit NMD and measured the effects of these treatments on W1282X expression and function. SMG1-ASO-mediated NMD inhibition upregulated the RNA, protein, and surface-localized protein expression of the truncated W1282X gene product. Additionally, these ASOs increased the CFTR chloride channel function in cells homozygous for the W1282X mutation. Our approach suggests a new therapeutic strategy for patients harboring nonsense mutations and may be beneficial as a single agent in patients with CF and the W1282X mutation.

UGT1A1 Variants c.864+5G>T and c.996+2_996+5del of a Crigler-Najjar Patient Induce Aberrant Splicing in Minigene Assays

A large fraction of DNA variants impairs pre-mRNA splicing in human hereditary disorders. Crigler-Najjar syndrome (CNS) is characterized by a severe unconjugated hyperbilirubinemia caused by variants in the UGT1A1 gene. We previously reported one CNS-type II patient with two splice-site variants in trans (c.864+5G>T and c.996+2_996+5del). According to MaxEntScan, both disrupt their corresponding donor sites (c.864+5G>T: 6.99 → 2.28; c.996+2_996+5del: 5.96 → -11.02), so they were selected for subsequent functional tests. Given the unavailability of patient RNA, we constructed an UGT1A1 splicing-reporter minigene with exons 1-4 to characterize the underlying splicing anomaly. The variant c.996+2_996+5del generated two aberrant transcripts, Δ(E2) (exon 2 skipping/64%) and ▼(E2q135) (intron retention of 135-nt/36%), which lead to the loss of 18 conserved amino-acids and the gain of 45 new ones of a critical functional domain, respectively. The c.864+5G>T variant mainly produced the aberrant transcript Δ(E1q141) (141-nt deletion/70.4%) and the full-length isoform (29.6%). Δ(E1q141) would provoke the loss of 47 amino-acids of the N-terminal domain that encodes for substrate specificity. Thus, the three anomalous transcripts are likely to inactivate UGT1A1. Moreover, this patient is also homozygous for the promoter variant A(TA)7TAA that decreases UGT1A1 expression by 70%, so the full-length transcript produced by c.864+5G>T would be even more reduced (<9%), thus supporting the diagnosis of CNS-type II. Therefore, minigenes represent valuable tools for the functional and clinical classifications of genetic variants.

Potentially lethal cystic fibrosis gene variant in the orangutan

A syndrome of chronic upper and lower airway disease leading to increased morbidity and mortality occurs primarily in captive orangutans. Similarities in symptoms to the inherited human respiratory disease, cystic fibrosis, led us to hypothesize that orangutan respiratory disease is a result of variants in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene. We identified the nonsense variant, c.484A>T (p.Lys162X), in heterozygosity in an unaffected orangutan. Analysis of the pedigree of this orangutan confirmed that both his sire and deceased fetus also harbored the c.484A>T allele. An expression minigene harboring c.484A>T produced no full-length CFTR protein in HEK293 cells. Finally, the c.484A>T CFTR messenger RNA abundance was severely reduced in primary nasal epithelial cells of the orangutan indicating that c.484A>T (p.Lys162X) is potentially lethal. Genetic screening of the captive orangutan population could be used to prevent transmission of this potentially lethal variant, and thus aid in the conservation of this critically endangered species.

Expert consensus guidelines for the genetic diagnosis of Alport syndrome

Recent expert guidelines recommend genetic testing for the diagnosis of Alport syndrome. Here, we describe current best practice and likely future developments. In individuals with suspected Alport syndrome, all three COL4A5, COL4A3 and COL4A4 genes should be examined for pathogenic variants, probably by high throughput-targeted next generation sequencing (NGS) technologies, with a customised panel for simultaneous testing of the three Alport genes. These techniques identify up to 95% of pathogenic COL4A variants. Where causative pathogenic variants cannot be demonstrated, the DNA should be examined for deletions or insertions by re-examining the NGS sequencing data or with multiplex ligation-dependent probe amplification (MLPA). These techniques identify a further 5% of variants, and the remaining few changes include deep intronic splicing variants or cases of somatic mosaicism. Where no pathogenic variants are found, the basis for the clinical diagnosis should be reviewed. Genes in which mutations produce similar clinical features to Alport syndrome (resulting in focal and segmental glomerulosclerosis, complement pathway disorders, MYH9-related disorders, etc.) should be examined. NGS approaches have identified novel combinations of pathogenic variants in Alport syndrome. Two variants, with one in COL4A3 and another in COL4A4, produce a more severe phenotype than an uncomplicated heterozygous change. NGS may also identify further coincidental pathogenic variants in genes for podocyte-expressed proteins that also modify the phenotype. Our understanding of the genetics of Alport syndrome is evolving rapidly, and both genetic and non-genetic factors are likely to contribute to the observed phenotypic variability.

Minigene Splicing Assays Identify 12 Spliceogenic Variants of BRCA2 Exons 14 and 15

A relevant fraction of BRCA2 variants is associated with splicing alterations and with an increased risk of hereditary breast and ovarian cancer (HBOC). In this work, we have carried out a thorough study of variants from BRCA2 exons 14 and 15 reported at mutation databases. A total of 294 variants from exons 14 and 15 and flanking intronic sequences were analyzed with the online splicing tools NNSplice and Human Splicing Finder. Fifty-three out of these 294 variants were selected as candidate splicing variants. All variants but one, were introduced into the minigene MGBR2_ex14-20 (with exons 14–20) by site-directed mutagenesis and assayed in MCF-7 cells. Twelve of the remaining 52 variants (23.1%) impaired splicing at different degrees, yielding from 5 to 100% of aberrant transcripts. Nine variants affected the natural acceptor or donor sites of both exons and three affected putative enhancers or silencers. Fluorescent capillary electrophoresis revealed at least 10 different anomalous transcripts: (E14q5), Δ (E14p10), Δ(E14p246), Δ(E14q256), Δ(E14), Δ(E15p12), Δ(E15p13), Δ(E15p83), Δ(E15) and a 942-nt fragment of unknown structure. All transcripts, except for Δ(E14q256) and Δ(E15p12), are expected to truncate the BRCA2 protein. Nine variants induced severe splicing aberrations with more than 90% of abnormal transcripts. Thus, according to the guidelines of the American College of Medical Genetics and Genomics, eight variants should be classified as pathogenic (c.7008-2A > T, c.7008-1G > A, c.7435+1G > C, c.7436-2A > T, c.7436-2A > G, c.7617+1G > A, c.7617+1G > T, and c.7617+2T > G), one as likely pathogenic (c.7008-3C > G) and three remain as variants of uncertain clinical significance or VUS (c.7177A > G, c.7447A > G and c.7501C > T). In conclusion, functional assays by minigenes constitute a valuable strategy to primarily check the splicing impact of DNA variants and their clinical interpretation. While bioinformatics predictions of splice site variants were accurate, those of enhancer or silencer variants were poor (only 3/23 spliceogenic variants) which showed weak impacts on splicing (∼5–16% of aberrant isoforms). So, the Exonic Splicing Enhancer and Silencer (ESE and ESS, respectively) prediction algorithms require further improvement.

Capitalizing on the heterogeneous effects of CFTR nonsense and frameshift variants to inform therapeutic strategy for cystic fibrosis

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic 'nulls' as some may allow generation of protein that can be targeted to achieve clinical benefit.

Splicing mutations in human genetic disorders: examples, detection, and confirmation

Precise pre-mRNA splicing, essential for appropriate protein translation, depends on the presence of consensus "cis" sequences that define exon-intron boundaries and regulatory sequences recognized by splicing machinery. Point mutations at these consensus sequences can cause improper exon and intron recognition and may result in the formation of an aberrant transcript of the mutated gene. The splicing mutation may occur in both introns and exons and disrupt existing splice sites or splicing regulatory sequences (intronic and exonic splicing silencers and enhancers), create new ones, or activate the cryptic ones. Usually such mutations result in errors during the splicing process and may lead to improper intron removal and thus cause alterations of the open reading frame. Recent research has underlined the abundance and importance of splicing mutations in the etiology of inherited diseases. The application of modern techniques allowed to identify synonymous and nonsynonymous variants as well as deep intronic mutations that affected pre-mRNA splicing. The bioinformatic algorithms can be applied as a tool to assess the possible effect of the identified changes. However, it should be underlined that the results of such tests are only predictive, and the exact effect of the specific mutation should be verified in functional studies. This article summarizes the current knowledge about the "splicing mutations" and methods that help to identify such changes in clinical diagnosis.

Functional assessment of the BMPR2 gene in lymphoblastoid cell lines from Graves' disease patients

In this study, we analysed the possible influence of the c.419-43delT BMPR2 variant in patients with Graves' disease (GD), in a molecular basis, focusing our efforts on possible alterations in the mRNA processing and synthesis. The molecular assessment of this variant in patients with GD would shed light on the association between the BMPR2 gene and the disease. The variant was detected in 18%, 55% and 10% of patients with pulmonary arterial hypertension, GD and in general population, respectively. Patients with GD fold change showed increased BMPR2 expression when matched against the controls, with a mean of 4.21 ± 1.73 (P = 0.001); BMPR2 was overexpressed in the analysed cell cycle stages. Fold change analysis of variant carriers and non-carriers showed slight overexpression and differences between phases, but none of them were statistically significant. BMPR2 expression was confirmed in the lymphoblastoid cell lines (LCLs) with a molecular weight of 115 kD, and no differences between variant carriers and non-carriers were detected. To conclude, the BMPR2 variant c.419-19delT appears in high frequency in patients with GD, and independently of its presence, BMPR2 is overexpressed in the LCLs from the GD patients tested. This increase could be paired with the described decreased expression of transforming growth factor-β1 in thyroid tissue from patients with GD.

Cas9/gRNA targeted excision of cystic fibrosis-causing deep-intronic splicing mutations restores normal splicing of CFTR mRNA

Cystic Fibrosis is an autosomal recessive disorder caused by mutations in the CFTR gene. CRISPR mediated, template-dependent homology-directed gene editing has been used to correct the most common mutation, c.1521_1523delCTT / p.Phe508del (F508del) which affects ~70% of individuals, but the efficiency was relatively low. Here, we describe a high efficiency strategy for editing of three different rare CFTR mutations which together account for about 3% of individuals with Cystic Fibrosis. The mutations cause aberrant splicing of CFTR mRNA due to the creation of cryptic splice signals that result in the formation of pseudoexons containing premature stop codons c.1679+1634A>G (1811+1.6kbA>G) and c.3718-2477C>T (3849+10kbC>T), or an out-of-frame 5' extension to an existing exon c.3140-26A>G (3272-26A>G). We designed pairs of Cas9 guide RNAs to create targeted double-stranded breaks in CFTR either side of each mutation which resulted in high efficiency excision of the target genomic regions via non-homologous end-joining repair. When evaluated in a mini-gene splicing assay, we showed that targeted excision restored normal splicing for all three mutations. This approach could be used to correct aberrant splicing signals or remove disruptive transcription regulatory motifs caused by deep-intronic mutations in a range of other genetic disorders.

Leveraging splice-affecting variant predictors and a minigene validation system to identify Mendelian disease-causing variants among exon-captured variants of uncertain significance

The genetic heterogeneity of Mendelian disorders results in a significant proportion of patients that are unable to be assigned a confident molecular diagnosis after conventional exon sequencing and variant interpretation. Here, we evaluated how many patients with an inherited retinal disease (IRD) have variants of uncertain significance (VUS) that are disrupting splicing in a known IRD gene by means other than affecting the canonical dinucleotide splice site. Three in silico splice-affecting variant predictors were leveraged to annotate and prioritize variants for splicing functional validation. An in vitro minigene system was used to assay each variant's effect on splicing. Starting with 745 IRD patients lacking a confident molecular diagnosis, we validated 23 VUS as splicing variants that likely explain disease in 26 patients. Using our results, we optimized in silico score cutoffs to guide future variant interpretation. Variants that alter base pairs other than the canonical GT-AG dinucleotide are often not considered for their potential effect on RNA splicing but in silico tools and a minigene system can be utilized for the prioritization and validation of such splice-disrupting variants. These variants can be overlooked causes of human disease but can be identified using conventional exon sequencing with proper interpretation guidelines.

Mutations of Pre-mRNA Splicing Regulatory Elements: Are Predictions Moving Forward to Clinical Diagnostics?

For more than three decades, researchers have known that consensus splice sites alone are not sufficient regulatory elements to provide complex splicing regulation. Other regulators, so-called splicing regulatory elements (SREs) are needed. Most importantly, their sequence variants often underlie the development of various human disorders. However, due to their variable location and high degeneracy, these regulatory sequences are also very difficult to recognize and predict. Many different approaches aiming to identify SREs have been tried, often leading to the development of in silico prediction tools. While these tools were initially expected to be helpful to identify splicing-affecting mutations in genetic diagnostics, we are still quite far from meeting this goal. In fact, most of these tools are not able to accurately discern the SRE-affecting pathological variants from those not affecting splicing. Nonetheless, several recent evaluations have given appealing results (namely for EX-SKIP, ESRseq and Hexplorer predictors). In this review, we aim to summarize the history of the different approaches to SRE prediction, and provide additional validation of these tools based on patients' clinical data. Finally, we evaluate their usefulness for diagnostic settings and discuss the challenges that have yet to be met.

Molecular and functional characterization of the BMPR2 gene in Pulmonary Arterial Hypertension

Pulmonary arterial hypertension is a progressive disease that causes the obstruction of precapillary pulmonary arteries and a sustained increase in pulmonary vascular resistance. The aim was to analyze functionally the variants found in the BMPR2 gene and to establish a genotype-phenotype correlation. mRNA expression studies were performed using pSPL3 vector, studies of subcellular localization were performed using pEGFP-N1 vector and luciferase assays were performed using pGL3-Basic vector. We have identified 30 variants in the BMPR2 gene in 27 of 55 patients. In 16 patients we detected pathogenic mutations. Minigene assays revealed that 6 variants (synonymous, missense) result in splicing defect. By immunofluorescence assay, we observed that 4 mutations affect the protein localization. Finally, 4 mutations located in the 5'UTR region showed a decreased transcriptional activity in luciferase assays. Genotype-phenotype correlation, revealed that patients with pathogenic mutations have a more severe phenotype (sPaP p = 0.042, 6MWT p = 0.041), a lower age at diagnosis (p = 0.040) and seemed to have worse response to phosphodiesterase-5-inhibitors (p = 0.010). Our study confirms that in vitro expression analysis is a suitable approach in order to investigate the phenotypic consequences of the nucleotide variants, especially in cases where the involved genes have a pattern of expression in tissues of difficult access.

Functional analysis by minigene assay of putative splicing variants found in Bardet-Biedl syndrome patients

Bardet–Biedl syndrome (BBS) and Alström syndrome (ALMS) are rare diseases belonging to the group of ciliopathies. Although mutational screening studies of BBS/ALMS cohorts have been extensively reported, little is known about the functional effect of those changes. Thus, splicing variants are estimated to represent 15% of disease‐causing mutations, and there is growing evidence that many exonic changes are really splicing variants misclassified. In this study, we aimed to analyse for the first time several variants in BBS2,ARL6/BBS3,BBS4 and ALMS1 genes predicted to produce aberrant splicing by minigene assay. We found discordance between bioinformatics analysis and experimental data when comparing wild‐type and mutant constructs. Remarkably, we identified nonsense variants presumably resistant to nonsense‐mediated decay, even when a premature termination codon would be introduced in the second amino acid (p.(G2*) mutation in ARL6/BBS3 gene). As a whole, we report one of the first functional studies of BBS/ALMS1 variants using minigene assay, trying to elucidate their role in disease. Functional studies of variants identified in BBS and ALMS patients are essential for their proper classification and subsequent genetic counselling and could also be the start point for new therapeutic approaches, currently based only on symptomatic treatment.

Systematic Computational Identification of Variants That Activate Exonic and Intronic Cryptic Splice Sites

We developed a variant-annotation method that combines sequence-based machine-learning classification with a context-dependent algorithm for selecting splice variants. Our approach is distinctive in that it compares the splice potential of a sequence bearing a variant with the splice potential of the reference sequence. After training, classification accurately identified 168 of 180 (93.3%) canonical splice sites of five genes. The combined method, CryptSplice, identified and correctly predicted the effect of 18 of 21 (86%) known splice-altering variants in CFTR, a well-studied gene whose loss-of-function variants cause cystic fibrosis (CF). Among 1,423 unannotated CFTR disease-associated variants, the method identified 32 potential exonic cryptic splice variants, two of which were experimentally evaluated and confirmed. After complete CFTR sequencing, the method found three cryptic intronic splice variants (one known and two experimentally verified) that completed the molecular diagnosis of CF in 6 of 14 individuals. CryptSplice interrogation of sequence data from six individuals with X-linked dyskeratosis congenita caused by an unknown disease-causing variant in DKC1 identified two splice-altering variants that were experimentally verified. To assess the extent to which disease-associated variants might activate cryptic splicing, we selected 458 pathogenic variants and 348 variants of uncertain significance (VUSs) classified as high confidence from ClinVar. Splice-site activation was predicted for 129 (28%) of the pathogenic variants and 75 (22%) of the VUSs. Our findings suggest that cryptic splice-site activation is more common than previously thought and should be routinely considered for all variants within the transcribed regions of genes.

A novel splicing mutation in GALT gene causing Galactosemia in Ecuadorian family

Classic Galactosemia (OMIM 230400) is an autosomal recessive disorder of galactose metabolism caused by mutations in the galactose-1-phosphate uridyl transferase (GALT) gene. This disease caused by the inability to metabolize galactose is potentially life-threatening but its pathophysiology has not been clearly defined. GALT gene presents high allelic heterogeneity and around 336 variations have been identified. Here, we report the case of a patient with Classic Galactosemia who was detected during a neonatal screening in Ecuador. Molecular study revealed a mutation in GALT gene intron 1, c.82+3A>G in homozygous condition, this mutation has not been previously reported. This gene variation was not found in any of the 119 healthy Ecuadorian individuals used as control. Furthermore, the mutation was the only alteration detected in the propositus's GALT after sequencing all exons and introns of this gene. In silico modeling predicted that the mutation was pathogenic.

A synonymous splicing mutation in the SF3B4 gene segregates in a family with highly variable Nager syndrome

Nager syndrome is a rare preaxial acrofacial dysostosis that is caused by heterozygous loss-of-function variants in SF3B4. This gene encodes for a protein required for the assembly of spliceosomal complexes, being a master gene for splicing regulation. The main clinical features of Nager syndrome include facial-mandibular and preaxial limb malformations, with normal cognitive functioning. Most Nager patients are sporadic, but few familial cases with a highly variable phenotype have been reported. In this work, we report a novel synonymous variant within exon 3 of the SF3B4 gene in a family with three members affected by Nager syndrome. No pathogenic variants have been detected in other 24 genes associated with syndromes characterized by mandibulo-facial anomalies. The pathogenicity of the mutation was demonstrated through a hybrid minigene assay, which confirmed an aberrant splicing with the creation of a cryptic splice site, and showed that this allele is hypomorphic. Our findings emphasize the importance to perform functional analyses to assess the possible consequences of synonymous variants and confirmed that hybrid minigenes represent an effective tool to evaluate the effects of variants on splicing, particularly when RNA is not available.

Whole-gene CFTR sequencing combined with digital RT-PCR improves genetic diagnosis of cystic fibrosis

Despite extensive screening, 1-5% of cystic fibrosis (CF) patients lack a definite molecular diagnosis. Next-generation sequencing (NGS) is making affordable genetic testing based on the identification of variants in extended genomic regions. In this frame, we analyzed 23 CF patients and one carrier by whole-gene CFTR resequencing: 4 were previously characterized and served as controls; 17 were cases lacking a complete diagnosis after a full conventional CFTR screening; 3 were consecutive subjects referring to our centers, not previously submitted to any screening. We also included in the custom NGS design the coding portions of the SCNN1A, SCNN1B and SCNN1G genes, encoding the subunits of the sodium channel ENaC, which were found to be mutated in CF-like patients. Besides 2 novel SCNN1B missense mutations, we identified 22 previously-known CFTR mutations, including 2 large deletions (whose breakpoints were precisely mapped), and novel deep-intronic variants, whose role on splicing was excluded by ex-vivo analyses. Finally, for 2 patients, compound heterozygotes for a CFTR mutation and the intron-9c.1210-34TG_[11-12]T₅ allele-known to be associated with decreased CFTR mRNA levels-the molecular diagnosis was implemented by measuring the residual level of wild-type transcript by digital reverse transcription polymerase chain reaction performed on RNA extracted from nasal brushing.

Correction of a Cystic Fibrosis Splicing Mutation by Antisense Oligonucleotides

Cystic fibrosis (CF), the most common life-threatening genetic disease in Caucasians, is caused by ∼2,000 different mutations in the CF transmembrane conductance regulator (CFTR) gene. A significant fraction of these (∼13%) affect pre-mRNA splicing for which novel therapies have been somewhat neglected. We have previously described the effect of the CFTR splicing mutation c.2657+5G>A in IVS16, showing that it originates transcripts lacking exon 16 as well as wild-type transcripts. Here, we tested an RNA-based antisense oligonucleotide (AON) strategy to correct the aberrant splicing caused by this mutation. Two AONs (AON1/2) complementary to the pre-mRNA IVS16 mutant region were designed and their effect on splicing was assessed at the RNA and protein levels, on intracellular protein localization and function. To this end, we used the 2657+5G>A mutant CFTR minigene stably expressed in HEK293 Flp-In cells that express a single copy of the transgene. RNA data from AON1-treated mutant cells show that exon 16 inclusion was almost completely restored (to 95%), also resulting in increased levels of correctly localized CFTR protein at the plasma membrane (PM) and with increased function. A novel two-color CFTR splicing reporter minigene developed here allowed the quantitative monitoring of splicing by automated microscopy localization of CFTR at the PM. The AON strategy is thus a promising therapeutic approach for the specific correction of alternative splicing.

Epidemiology and genetics of cystic fibrosis in Asia: In preparation for the next-generation treatments

Cystic fibrosis (CF) in the Asian population is less frequently reported due to under-diagnosis and lack of centralized CF patient registries. Clinical studies on CF cases from Asia have documented a severe course of the disease. The spectrum of the cystic fibrosis transmembrane conductance regulator (CFTR) variants in this population is quite heterogeneous. In total, 166 variants have been reported on approximately 3700 Asian CF chromosomes. The frequency of F508del among Asians is low compared with Caucasians. Recent in vitro studies have shown promise of small molecule correction and potentiation of 45 different CFTR variants. Of these variants, 16 (including G551D and F508del) have also been observed among Asian CF individuals. We suggest undertaking molecular studies extensively to annotate CFTR variants that will help Asian CF individuals to benefit from the precision medicine gaining momentum in the Western countries.

A rigorous approach for selection of optimal variant sets for carrier screening with demonstration of clinical utility

Carrier screening for certain diseases is recommended by major medical and Ashkenazi Jewish (AJ) societies. Most carrier screening panels test only for common, ethnic-specific variants. However, with formerly isolated ethnic groups becoming increasingly intermixed, this approach is becoming inadequate. Our objective was to develop a rigorous process to curate all variants, for relevant genes, into a database and then apply stringent clinical validity classification criteria to each in order to retain only those with clear evidence for pathogenicity. The resulting variant set, in conjunction with next-generation DNA sequencing (NGS), then affords the capability for an ethnically diverse, comprehensive, highly specific carrier-screening assay. The clinical utility of our approach was demonstrated by screening a pan-ethnic population of 22,864 individuals for Bloom syndrome carrier status using a BLM variant panel comprised of 50 pathogenic variants. In addition to carriers of the common AJ founder variant, we identified 57 carriers of other pathogenic BLM variants. All variants reported had previously been curated and their clinical validity documented, or were of a type that met our stringent, preassigned validity criteria. Thus, it was possible to confidently report an increased number of Bloom’s syndrome carriers compared to traditional, ethnicity-based screening, while not reducing the specificity of the screening due to reporting variants of unknown clinical significance.

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.

Small-scale high-throughput sequencing-based identification of new therapeutic tools in cystic fibrosis

PURPOSE

Although 97-99% of CFTR mutations have been identified, great efforts must be made to detect yet-unidentified mutations.

METHODS

We developed a small-scale next-generation sequencing approach for reliably and quickly scanning the entire gene, including noncoding regions, to identify new mutations. We applied this approach to 18 samples from patients suffering from cystic fibrosis (CF) in whom only one mutation had hitherto been identified.

RESULTS

Using an in-house bioinformatics pipeline, we could rapidly identify a second disease-causing CFTR mutation for 16 of 18 samples. Of them, c.1680-883A>G was found in three unrelated CF patients. Analysis of minigenes and patients' transcripts showed that this mutation results in aberrantly spliced transcripts because of the inclusion of a pseudoexon. It is located only three base pairs from the c.1680-886A>G mutation (1811+1.6kbA>G), the fourth most frequent mutation in southwestern Europe. We next tested the effect of antisense oligonucleotides targeting splice sites on these two mutations on pseudoexon skipping. Oligonucleotide transfection resulted in the restoration of the full-length, in-frame CFTR transcript, demonstrating the effect of antisense oligonucleotide-induced pseudoexon skipping in CF.

CONCLUSION

Our data confirm the importance of analyzing noncoding regions to find unidentified mutations, which is essential to designing targeted therapeutic approaches.

Identification and characterization of aberrant GAA pre-mRNA splicing in pompe disease using a generic approach

Identification of pathogenic variants in monogenic diseases is an important aspect of diagnosis, genetic counseling, and prediction of disease severity. Pathogenic mechanisms involved include changes in gene expression, RNA processing, and protein translation. Variants affecting pre-mRNA splicing are difficult to predict due to the complex mechanism of splicing regulation. A generic approach to systematically detect and characterize effects of sequence variants on splicing would improve current diagnostic practice. Here, it is shown that such approach is feasible by combining flanking exon RT-PCR, sequence analysis of PCR products, and exon-internal quantitative RT-PCR for all coding exons. Application of this approach to one novel and six previously published variants in the acid-alpha glucosidase (GAA) gene causing Pompe disease enabled detection of a total of 11 novel splicing events. Aberrant splicing included cryptic splice-site usage, intron retention, and exon skipping. Importantly, the extent of leaky wild-type splicing correlated with disease onset and severity. These results indicate that this approach enables sensitive detection and in-depth characterization of variants affecting splicing, many of which are still unrecognized or poorly understood. The approach is generic and should be adaptable for application to other monogenic diseases to aid in improved diagnostics.

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.