Single base-pair substitutions in exon-intron junctions of human genes: nature, distribution, and consequences for mRNA splicing

ABO gene polymorphism and COVID-19 severity: The impact on haematological complications, inflammatory markers, and lung lesions

PURPOSE

The study aimed to investigate the connection between an intronic variant in the ABO gene (rs657152) and the severity of COVID-19 in terms of clinical symptoms, haematological complications, inflammatory markers, and lung lesions.

METHODS

After applying exclusion criteria, the study included 240 patients divided into 3 groups: 88 Outpatients, 84 Ward-hospitalized, and 68 ICU-admitted/failed patients. The tetra-ARMS PCR method was used to genotype ABO polymorphism in the patient. Paraclinical tests of patients at the time of admission (before receiving conventional treatments) included levels of C-reactive protein (CRP) and erythrocyte sedimentation rate (ESR), as well as a complete blood count (CBC). Also, the severity of lung lesions was evaluated based on the results of spiral computed tomography (CT) of the chest during admission.

RESULTS

The statistical analysis using the ANOVA test revealed significant differences in the mean values of allele frequencies (p-value = 0.0020) and genotype proportions (p-value = 0.0017) among clinical groups. The study also found a notable difference in ABO polymorphism across different levels of the inflammatory marker CRP, but not with the ESR levels. Furthermore, the study showed a significant difference in the distribution of lung lesion severity and ABO polymorphism among different clinical groups.

CONCLUSION

To conclude, our findings supported the substantial impact of ABO polymorphism rs657152 on the severity of COVID-19 in Iranian patients, specifically concerning haematological complications, inflammatory markers, and lung lesions. The study underscored the protective effect of the AC genotype and the detrimental impact of the CC genotype on clinical manifestations.

OpenASO: RNA Rescue - designing splice-modulating antisense oligonucleotides through community science

Splice-modulating antisense oligonucleotides (ASOs) are precision RNA-based drugs that are becoming an established modality to treat human disease. Previously, we reported the discovery of ASOs that target a novel, putative intronic RNA structure to rescue splicing of multiple pathogenic variants of F8 exon 16 that cause hemophilia A. However, the conventional approach to discovering splice-modulating ASOs is both laborious and expensive. Here, we describe an alternative paradigm that integrates data-driven RNA structure prediction and community science to discover splice-modulating ASOs. Using a splicing-deficient pathogenic variant of F8 exon 16 as a model, we show that 25% of the top-scoring molecules designed in the Eterna OpenASO challenge have a statistically significant impact on enhancing exon 16 splicing. Additionally, we show that a distinct combination of ASOs designed by Eterna players can additively enhance the inclusion of the splicing-deficient exon 16 variant. Together, our data suggests that crowdsourcing designs from a community of citizen scientists may accelerate the discovery of splice-modulating ASOs with potential to treat human disease.

Protein isoform-centric therapeutics: expanding targets and increasing specificity

Most protein-coding genes produce multiple protein isoforms; however, these isoforms are commonly neglected in drug discovery. The expression of protein isoforms can be specific to a disease, tissue and/or developmental stage, and this specific expression can be harnessed to achieve greater drug specificity than pan-targeting of all gene products and to enable improved treatments for diseases caused by aberrant protein isoform production. In recent years, several protein isoform-centric therapeutics have been developed. Here, we collate these studies and clinical trials to highlight three distinct but overlapping modes of action for protein isoform-centric drugs: isoform switching, isoform introduction or depletion, and modulation of isoform activity. In addition, we discuss how protein isoforms can be used clinically as targets for cell type-specific drug delivery and immunotherapy, diagnostic biomarkers and sources of cancer neoantigens. Collectively, we emphasize the value of a focus on isoforms as a route to discovering drugs with greater specificity and fewer adverse effects. This approach could enable the targeting of proteins for which pan-inhibition of all isoforms is toxic and poorly tolerated.

Validating the splicing effect of rare variants in the SLC26A4 gene using minigene assay

BACKGROUND

The SLC26A4 gene is the second most common cause of hereditary hearing loss in human. The aim of this study was to utilize the minigene assay in order to identify pathogenic variants of SLC26A4 associated with enlarged vestibular aqueduct (EVA) and hearing loss (HL) in two patients.

METHODS

The patients were subjected to multiplex PCR amplification and next-generation sequencing of common deafness genes (including GJB2, SLC26A4, and MT-RNR1), then bioinformatics analysis was performed on the sequencing data to identify candidate pathogenic variants. Minigene experiments were conducted to determine the potential impact of the variants on splicing.

RESULTS

Genetic testing revealed that the first patient carried compound heterozygous variants c.[1149 + 1G > A]; [919-2 A > G] in the SLC26A4 gene, while the second patient carried compound heterozygous variants c.[2089 + 3 A > T]; [919-2 A > G] in the same gene. Minigene experiments demonstrated that both c.1149 + 1G > A and c.2089 + 3 A > T affected mRNA splicing. According to the ACMG guidelines and the recommendations of the ClinGen Hearing Loss Expert Panel for ACMG variant interpretation, these variants were classified as "likely pathogenic".

CONCLUSIONS

This study identified the molecular etiology of hearing loss in two patients with EVA and elucidated the impact of rare variants on splicing, thus contributing to the mutational spectrum of pathogenic variants in the SLC26A4 gene.

A systematic assessment of the impact of rare canonical splice site variants on splicing using functional and in silico methods

Canonical splice site variants (CSSVs) are often presumed to cause loss-of-function (LoF) and are assigned very strong evidence of pathogenicity (according to American College of Medical Genetics/Association for Molecular Pathology criterion PVS1). The exact nature and predictability of splicing effects of unselected rare CSSVs in blood-expressed genes are poorly understood. We identified 168 rare CSSVs in blood-expressed genes in 112 individuals using genome sequencing, and studied their impact on splicing using RNA sequencing (RNA-seq). There was no evidence of a frameshift, nor of reduced expression consistent with nonsense-mediated decay, for 25.6% of CSSVs: 17.9% had wildtype splicing only and normal junction depths, 3.6% resulted in cryptic splice site usage and in-frame insertions or deletions, 3.6% resulted in full exon skipping (in frame), and 0.6% resulted in full intron inclusion (in frame). Blind to these RNA-seq data, we attempted to predict the precise impact of CSSVs by applying in silico tools and the ClinGen Sequence Variant Interpretation Working Group 2018 guidelines for applying PVS1 criterion. The predicted impact on splicing using (1) SpliceAI, (2) MaxEntScan, and (3) AutoPVS1, an automatic classification tool for PVS1 interpretation of null variants that utilizes Ensembl Variant Effect Predictor and MaxEntScan, was concordant with RNA-seq analyses for 65%, 63%, and 61% of CSSVs, respectively. In summary, approximately one in four rare CSSVs did not show evidence for LoF based on analysis of RNA-seq data. Predictions from in silico methods were often discordant with findings from RNA-seq. More caution may be warranted in applying PVS1-level evidence to CSSVs in the absence of functional data.

APOER2 splicing repertoire in Alzheimer's disease: Insights from long-read RNA sequencing

Disrupted alternative splicing plays a determinative role in neurological diseases, either as a direct cause or as a driver in disease susceptibility. Transcriptomic profiling of aged human postmortem brain samples has uncovered hundreds of aberrant mRNA splicing events in Alzheimer's disease (AD) brains, associating dysregulated RNA splicing with disease. We previously identified a complex array of alternative splicing combinations across apolipoprotein E receptor 2 (APOER2), a transmembrane receptor that interacts with both the neuroprotective ligand Reelin and the AD-associated risk factor, APOE. Many of the human APOER2 isoforms, predominantly featuring cassette splicing events within functionally important domains, are critical for the receptor's function and ligand interaction. However, a comprehensive repertoire and the functional implications of APOER2 isoforms under both physiological and AD conditions are not fully understood. Here, we present an in-depth analysis of the splicing landscape of human APOER2 isoforms in normal and AD states. Using single-molecule, long-read sequencing, we profiled the entire APOER2 transcript from the parietal cortex and hippocampus of Braak stage IV AD brain tissues along with age-matched controls and investigated several functional properties of APOER2 isoforms. Our findings reveal diverse patterns of cassette exon skipping for APOER2 isoforms, with some showing region-specific expression and others unique to AD-affected brains. Notably, exon 15 of APOER2, which encodes the glycosylation domain, showed less inclusion in AD compared to control in the parietal cortex of females with an APOE ɛ3/ɛ3 genotype. Also, some of these APOER2 isoforms demonstrated changes in cell surface expression, APOE-mediated receptor processing, and synaptic number. These variations are likely critical in inducing synaptic alterations and may contribute to the neuronal dysfunction underlying AD pathogenesis.

Identification and characterization of a novel intronic splicing mutation in CSF1R-related leukoencephalopathy

AIMS

Colony stimulating factor 1 receptor (CSF1R)-related leukoencephalopathy is a rapidly progressing neurodegenerative disease caused by CSF1R gene mutations. This study aimed to identify and investigate the effect of a novel intronic mutation (c.1754-3C>G) of CSF1R on splicing.

METHODS

A novel intronic mutation was identified using whole-exome sequencing. To investigate the impact of this mutation, we employed various bioinformatics tools to analyze the transcription of the CSF1R gene and the three-dimensional structure of its encoded protein. Furthermore, reverse transcription polymerase chain reaction (RT-PCR) was performed to validate the findings.

RESULTS

A novel mutation (c.1754-3C>G) in CSF1R was identified, which results in exon 13 skipping due to the disruption of the 3' splice site consensus sequence NYAG/G. This exon skipping event was further validated in the peripheral blood of the mutation carrier through RT-PCR and Sanger sequencing. Protein structure prediction indicated a disruption in the tyrosine kinase domain, with the truncated protein showing significant structural alterations.

CONCLUSIONS

Our findings underscore the importance of intronic mis-splicing mutations in the diagnosis and management of CSF1R-related leukoencephalopathy.

A partial form of AIRE deficiency underlies a mild form of autoimmune polyendocrine syndrome type 1

Autoimmune polyendocrine syndrome type 1 (APS-1) is caused by mutations in the autoimmune regulator (AIRE) gene. Most patients present with severe chronic mucocutaneous candidiasis and organ-specific autoimmunity from early childhood, but the clinical picture is highly variable. AIRE is crucial for negative selection of T cells, and scrutiny of different patient mutations has previously highlighted many of its molecular mechanisms. In patients with a milder adult-onset phenotype sharing a mutation in the canonical donor splice site of intron 7 (c.879+1G>A), both the predicted altered splicing pattern with loss of exon 7 (AireEx7-/-) and normal full-length AIRE mRNA were found, indicating leaky rather than abolished mRNA splicing. Analysis of a corresponding mouse model demonstrated that the AireEx7-/- mutant had dramatically impaired transcriptional capacity of tissue-specific antigens in medullary thymic epithelial cells but still retained some ability to induce gene expression compared with the complete loss-of-function AireC313X-/- mutant. Our data illustrate an association between AIRE activity and the severity of autoimmune disease, with implications for more common autoimmune diseases associated with AIRE variants, such as primary adrenal insufficiency, pernicious anemia, type 1 diabetes, and rheumatoid arthritis.

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.

Disruption of mitochondrial and lysosomal functions by human CACNA1C variants expressed in HEK 293 and CHO cells

Objective

To investigate the pathogenesis of three novel de novo CACNA1C variants (p.E411D, p.V622G, and p.A272V) in causing neurodevelopmental disorders and arrhythmia.

Methods

Several molecular experiments were carried out on transfected human embryonic kidney 293 (HEK 293) and Chinese hamster ovary (CHO) cells to explore the effects of p.E411D, p.V622G, and p.A272V variants on electrophysiology, mitochondrial and lysosomal functions. Electrophysiological studies, RT-qPCR, western blot, apoptosis assay, mito-tracker fluorescence intensity, lyso-tracker fluorescence intensity, mitochondrial calcium concentration test, and cell viability assay were performed. Besides, reactive oxygen species (ROS) levels, ATP levels, mitochondrial copy numbers, mitochondrial complex I, II, and cytochrome c functions were measured.

Results

The p.E411D variant was found in a patient with attention deficit-hyperactive disorder (ADHD), and moderate intellectual disability (ID). This mutant demonstrated reduced calcium current density, mRNA, and protein expression, and it was localized in the nucleus, cytoplasm, lysosome, and mitochondria. It exhibited an accelerated apoptosis rate, impaired autophagy, and mitophagy. It also demonstrated compromised mitochondrial cytochrome c oxidase, complex I, and II enzymes, abnormal mitochondrial copy numbers, low ATP levels, abnormal mitochondria fluorescence intensity, impaired mitochondrial fusion and fission, and elevated mitochondrial calcium ions. The p.V622G variant was identified in a patient who presented with West syndrome and moderate global developmental delay. The p.A272V variant was found in a patient who presented with epilepsy and mild ID. Both mutants (p.V622G and p.A272V) exhibited reduced calcium current densities, decreased mRNA and protein expressions, and they were localized in the nucleus, cytoplasm, lysosome, and mitochondria. They exhibited accelerated apoptosis and proliferation rates, impaired autophagy, and mitophagy. They also exhibited abnormal mitochondrial cytochrome c oxidase, complex I and II enzymes, abnormal mitochondrial copy numbers, low ATP, high ROS levels, abnormal mitochondria fluorescence intensity, impaired mitochondrial fusion and fission, as well as elevated mitochondrial calcium ions.

Conclusion

The p.E411D, p.V622G and p.A272V mutations of human CACNA1C reduce the expression level of CACNA1C proteins, and impair mitochondrial and lysosomal functions. These effects induced by CACNA1C variants may contribute to the pathogenesis of CACNA1C-related disorders.

Improving Hereditary Hemorrhagic Telangiectasia Molecular Diagnosis: A Referral Center Experience

BACKGROUND

Hereditary hemorrhagic telangiectasia (HHT) is a rare vascular disease inherited in an autosomal dominant manner. Disease-causing variants in endoglin (ENG) and activin A receptor type II-like 1 (ACVRL1) genes are detected in more than 90% of the patients undergoing molecular testing. The identification of variants of unknown significance is often seen as a challenge in clinical practice that makes family screening and genetic counseling difficult. Here, we show that the implementation of cDNA analysis to assess the effect of splice site variants on mRNA splicing is a powerful tool.

METHODS

Gene panel sequencing of genes associated with HHT and other arteriovenous malformation-related syndromes was performed. To evaluate the effect of the splice site variants, cDNA analysis of ENG and ACVRL1 genes was carried out.

RESULTS

three novel splice site variants were identified in ENG (c.68-2A > T and c.1311+4_1311+8del) and ACVLR1 (c.526-6C > G) genes correspondingly in three individuals with HHT that met ≥ 3 Curaçao criteria. All three variants led to an aberrant splicing inducing exon skipping (ENG:c.68-2A > T and ACVRL1:c.526-6C > G) or intron retention (ENG:c.1311+4_1311+8del) allowing the confirmation of the predicted effect on splicing and the reclassification from unknown significance to pathogenic/likely pathogenic of two of them.

CONCLUSIONS

RNA analysis should be performed to assess and/or confirm the impact of variants on splicing. The molecular diagnosis of HHT patients is crucial to allow family screening and accurate genetic counseling. A multidisciplinary approach including clinicians and geneticists is crucial when dealing with patients with rare diseases.

Polymorphism of glutathione S-transferase P1 of dogs with mammary tumours

Mammary tumours constitute more than half of neoplasms in female dogs from different countries. Genome sequences are associated with cancer susceptibility but there is little information available about genetic polymorphisms of glutathione S-transferase P1 (GSTP1) in canine cancers. The aim of this study was to find single nucleotide polymorphisms (SNPs) in GSTP1 of dogs (Canis lupus familiaris) with mammary tumours compared to healthy dogs and to determine the association between GSTP1 polymorphisms and the occurrence of these tumours. The study population included 36 client-owned female dogs with mammary tumours and 12 healthy female dogs, with no previous diagnosis of cancer. DNA was extracted from blood and amplified by PCR assay. PCR-products were sequenced by Sanger method and analysed manually. The 33 polymorphisms were found in GSTP1: 1 coding SNP (exon 4), 24 non-coding SNPs (9 in exon 1), 7 deletions and 1 insertion. The 17 polymorphisms have been found in introns 1, 4, 5 and 6. The dogs with mammary tumours have significant difference from healthy in SNPs I4 c.1018 + 123 T > C (OR 13.412, 95%CI 1.574-114.267, P = .001), I5 c.1487 + 27 T > C (OR 10.737, 95%CI 1.260-91.477, P = .004), I5 c.1487 + 842 G > C (OR 4.714, 95% CI 1.086-20.472, P = .046) and I6 c.2481 + 50 A > G (OR 12.000, 95% CI 1.409-102.207, P = .002). SNP E5 c.1487 T > C and I5 c.1487 + 829 delG also differed significantly (P = .03) but not to the confidence interval. The study, for the first time, showed a positive association of SNPs in GSTP1 with mammary tumours of dogs, that can possibly be used to predict the occurrence of this pathology.

Systematic analysis of CNGA3 splice variants identifies different mechanisms of aberrant splicing

Achromatopsia is an autosomal recessive cone photoreceptor disease that is frequently caused by pathogenic variants in the CNGA3 gene. Here, we present a systematic functional analysis of 20 CNGA3 splice site variants detected in our large cohort of achromatopsia patients and/or listed in common variant databases. All variants were analyzed by functional splice assays based on the pSPL3 exon trapping vector. We demonstrated that ten variants, both at canonical and non-canonical splice sites, induced aberrant splicing, including intronic nucleotide retention, exonic nucleotide deletion and exon skipping, resulting in 21 different aberrant transcripts. Of these, eleven were predicted to introduce a premature termination codon. The pathogenicity of all variants was assessed based on established guidelines for variant classification. Incorporation of the results of our functional analyses enabled re-classification of 75% of variants previously classified as variants of uncertain significance into either likely benign or likely pathogenic. Our study is the first in which a systematic characterization of putative CNGA3 splice variants has been performed. We demonstrated the utility of pSPL3 based minigene assays in the effective assessment of putative splice variants. Our findings improve the diagnosis of achromatopsia patients, who may thus benefit from future gene-based therapeutic strategies.

Mitochondrial carbonic anhydrase VA and VB: properties and roles in health and disease

Carbonic anhydrase V (CA V), a mitochondrial enzyme, was first isolated from guinea-pig liver and subsequently identified in mice and humans. Later, studies revealed that the mouse genome contains two mitochondrial CA sequences, named Car5A and Car5B. The CA VA enzyme is most highly expressed in the liver, whereas CA VB shows a broad tissue distribution. Car5A knockout mice demonstrated a predominant role for CA VA in ammonia detoxification, whereas the roles of CA VB in ureagenesis and gluconeogenesis were evident only in the absence of CA VA. Previous studies have suggested that CA VA is mainly involved in the provision of HCO₃ ^- for biosynthetic processes. In children, mutations in the CA5A gene led to reduced CA activity, and the enzyme was sensitive to increased temperature. The metabolic profiles of these children showed a reduced supply of HCO₃ ^- to the enzymes that take part in intermediary metabolism: carbamoylphosphate synthetase, pyruvate carboxylase, propionyl-CoA carboxylase and 3-methylcrotonyl-CoA carboxylase. Although the role of CA VB is still poorly understood, a recent study reported that it plays an essential role in human Sertoli cells, which sustain spermatogenesis. Metabolic disease associated with CA VA appears to be more common than other inborn errors of metabolism and responds well to treatment with N-carbamyl-l-glutamate. Therefore, early identification of hyperammonaemia will allow specific treatment with N-carbamyl-l-glutamate and prevent neurological sequelae. Carbonic anhydrase VA deficiency should therefore be considered a treatable condition in the differential diagnosis of hyperammonaemia in neonates and young children.

Case report: A novel heterozygous synonymous variant in deep exon region of NIPBL gene generating a non-canonical splice donor in a patient with cornelia de lange syndrome

Cornelia de Lange syndrome (CdLS) is an autosomal dominant or X-linked genetic disease with significant genetic heterogeneity. Variants of the NIPBL gene are responsible for CdLS in 60% of patients. Herein, we report the case of a patient with CdLS showing distinctive facial features, microcephaly, developmental delay, and growth retardation. Whole exome sequencing was performed for the patient, and a novel de novo heterozygous synonymous variant was identified in the deep region of exon 40 in the NIPBL gene (NM_133433.4: c. 6819G > T, p. Gly2273 = ). The clinical significance of the variant was uncertain according to the ACMG/AMP guidelines; however, based on in silico analysis, it was predicted to alter mRNA splicing. To validate the prediction, a reverse transcriptase-polymerase chain reaction was conducted. The variant activated a cryptic splice donor, generating a short transcript of NIPBL. A loss of 137 bp at the 3' end of NIPBL exon 40 was detected, which potentially altered the open reading frame by inserting multiple premature termination codons. Quantitative real-time PCR analysis showed that the ratio of the transcription level of the full-length transcript to that of the altered short transcript in the patient was 5:1, instead of 1:1. These findings may explain the relatively mild phenotype of the patient, regardless of the loss of function of the truncated protein due to a frameshift in the mRNA. To the best of our knowledge, this study is the first to report a synonymous variant in the deep exon regions of the NIPBL gene responsible for CdLS. The identified variant expands the mutational spectrum of the NIPBL gene. Furthermore, synonymous variations may be pathogenic, which should not be ignored in the clinical and genetic diagnosis of the disease.

Selection of Diagnostically Significant Regions of the SLC26A4 Gene Involved in Hearing Loss

Screening pathogenic variants in the SLC26A4 gene is an important part of molecular genetic testing for hearing loss (HL) since they are one of the common causes of hereditary HL in many populations. However, a large size of the SLC26A4 gene (20 coding exons) predetermines the difficulties of its complete mutational analysis, especially in large samples of patients. In addition, the regional or ethno-specific prevalence of SLC26A4 pathogenic variants has not yet been fully elucidated, except variants c.919-2A>G and c.2168A>G (p.His723Arg), which have been proven to be most common in Asian populations. We explored the distribution of currently known pathogenic and likely pathogenic (PLP) variants across the SLC26A4 gene sequence presented in the Deafness Variation Database for the selection of potential diagnostically important parts of this gene. As a result of this bioinformatic analysis, we found that molecular testing ten SLC26A4 exons (4, 6, 10, 11, 13−17 and 19) with flanking intronic regions can provide a diagnostic rate of 61.9% for all PLP variants in the SLC26A4 gene. The primary sequencing of these SLC26A4 regions may be applied as an initial effective diagnostic testing in samples of patients of unknown ethnicity or as a subsequent step after the targeted testing of already-known ethno- or region-specific pathogenic SLC26A4 variants.

Unpredicted Aberrant Splicing Products Identified in Postmortem Sudden Cardiac Death Samples

Molecular screening for pathogenic mutations in sudden cardiac death (SCD)-related genes is common practice for SCD cases. However, test results may lead to uncertainty because of the identification of variants of unknown significance (VUS) occurring in up to 70% of total identified variants due to a lack of experimental studies. Genetic variants affecting potential splice site variants are among the most difficult to interpret. The aim of this study was to examine rare intronic variants identified in the exonic flanking sequence to meet two main objectives: first, to validate that canonical intronic variants produce aberrant splicing; second, to determine whether rare intronic variants predicted as VUS may affect the splicing product. To achieve these objectives, 28 heart samples of cases of SCD carrying rare intronic variants were studied. Samples were analyzed using 85 SCD genes in custom panel sequencing. Our results showed that rare intronic variants affecting the most canonical splice sites displayed in 100% of cases that they would affect the splicing product, possibly causing aberrant isoforms. However, 25% of these cases (1/4) showed normal splicing, contradicting the in silico results. On the contrary, in silico results predicted an effect in 0% of cases, and experimental results showed >20% (3/14) unpredicted aberrant splicing. Thus, deep intron variants are likely predicted to not have an effect, which, based on our results, might be an underestimation of their effect and, therefore, of their pathogenicity classification and family members’ follow-up.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Effects of 14 F9 synonymous codon variants on hemophilia B expression: Alteration of splicing along with protein expression

There is growing evidence that synonymous codon variants (SCVs) can cause disease through the disruption of different processes of protein production. The aim of the study is to investigate whether the 14 SCVs reported in the F9 variant database were the pathogenic causes of hemophilia B. The impacts of SCVs on splicing and protein expression were detected using a combination of in silico prediction, in vitro minigene splicing assay and cell expression detection. The splicing transcripts were identified and quantified by co-amplification fluorescent PCR. The mechanism of splicing was verified by a modified pU1snRNA and pU7snRNA approach. Aberrant splicing patterns were found in eight SCVs. Five of the 8 SCVs produced almost all aberrant splicing isoforms, which were expected to truncate protein, three of them presented a partial defect on both splicing and protein secretion, the overall effects were consistent with the residual Factor IX activity of the affected cases. Neither the pre-messenger RNA (mRNA) splicing process nor the protein function was impaired in the rest six SCVs. In conclusion, our study firstly revealed the pathogenic mechanism of the 14 F9 SCVs and highlighted the importance of performing mRNA splicing analysis and protein expression studies of SCVs in inherited disorders.

Splicing mutations in the CFTR gene as therapeutic targets

The marketing approval, about ten years ago, of the first disease modulator for patients with cystic fibrosis harboring specific CFTR genotypes (~5% of all patients) brought new hope for their treatment. To date, several therapeutic strategies have been approved and the number of CFTR mutations targeted by therapeutic agents is increasing. Although these drugs do not reverse the existing disease, they help to increase the median life expectancy. However, on the basis of their CFTR genotype, ~10% of patients presently do not qualify for any of the currently available CFTR modulator therapies, particularly patients with splicing mutations (~12% of the reported CFTR mutations). Efforts are currently made to develop therapeutic agents that target disease-causing CFTR variants that affect splicing. This highlights the need to fully identify them by scanning non-coding regions and systematically determine their functional consequences. In this review, we present some examples of CFTR alterations that affect splicing events and the different therapeutic options that are currently developed and tested for splice switching.

Donor Splice Site Variant in SLC9A6 Causes Christianson Syndrome in a Lithuanian Family: A Case Report

Background and Objectives: The pathogenic variants of SLC9A6 are a known cause of a rare, X-linked neurological disorder called Christianson syndrome (CS). The main characteristics of CS are developmental delay, intellectual disability, and neurological findings. This study investigated the genetic basis and explored the molecular changes that led to CS in two male siblings presenting with intellectual disability, epilepsy, behavioural problems, gastrointestinal dysfunction, poor height, and weight gain. Materials and Methods: Next-generation sequencing of a tetrad was applied to identify the DNA changes and Sanger sequencing of proband’s cDNA was used to evaluate the impact of a splice site variant on mRNA structure. Bioinformatical tools were used to investigate SLC9A6 protein structure changes. Results: Sequencing and bioinformatical analysis revealed a novel donor splice site variant (NC_000023.11(NM_001042537.1):c.899 + 1G > A) that leads to a frameshift and a premature stop codon. Protein structure modelling showed that the truncated protein is unlikely to form any functionally relevant SLC9A6 dimers. Conclusions: Molecular and bioinformatical analysis revealed the impact of a novel donor splice site variant in the SLC9A6 gene that leads to truncated and functionally disrupted protein causing the phenotype of CS in the affected individuals.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

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

RNA Analysis and Clinical Characterization of a Novel Splice Variant in the NSD1 Gene Causing Familial Sotos Syndrome

BACKGROUND

Sotos syndrome is an autosomal dominant disorder characterized by overgrowth, macrocephaly, distinctive facial features and learning disabilities. Haploinsufficiency of the nuclear receptor SET domain-containing protein 1 (NSD1) gene located on chromosome 5q35 is the major cause of the syndrome. This syndrome shares characteristics with other overgrowth syndromes, which can complicate the differential diagnosis.

METHODS

Genomic DNA was extracted from peripheral blood samples of members of the same family and targeted exome analysis was performed. In silico study of the variant found by next-generation sequencing was used to predict disruption/creation of splice sites and the identification of potential cryptic splice sites. RNA was extracted from peripheral blood samples of patients and functional analyses were performed to confirm the pathogenicity.

RESULTS

We found a novel c.6463 + 5G>A heterozygous NSD1 gene pathogenic variant in a son and his father. Molecular analyses revealed that part of the intron 22 of NSD1 is retained due to the destruction of the splicing donor site, causing the appearance of a premature stop codon in the NSD1 protein.

CONCLUSIONS

Our findings underline the importance of performing RNA functional assays in order to determine the clinical significance of intronic variants, and contribute to the genetic counseling and clinical management of patients and their relatives. Our work also highlights the relevance of using in silico prediction tools to detect a potential alteration in the splicing process.

Hereditary spastic paraplegia associated with a novel homozygous intronic noncanonical splice site variant in the AP4B1 gene

Pathogenic variants in the AP4B1 gene lead to a rare form of hereditary spastic paraplegia (HSP) known as SPG47. We report on a patient with a clinical suspicion of complicated HSP of the lower limbs with intellectual disability, as well as a novel homozygous noncanonical splice site variant in the AP4B1 gene, in which the effect on splicing was validated by RNA analysis. We sequenced 152 genes associated with HSP using Next-Generation Sequencing (NGS). We isolated total RNA from peripheral blood and generated cDNA using reverse transcription-polymerase chain reaction (RT-PCR). A region of AP4B1 mRNA was amplified by PCR and the fragments obtained were purified from the agarose gel and sequenced. We found a homozygous variant of uncertain significance in the AP4B1 gene NM_006594.4: c.1511-6C>G in the proband. Two different AP4B1 mRNA fragments were obtained in the patient and his carrier parents. The shorter fragment was the predominant fragment in the patient and revealed a deletion with skipping of the AP4B1 exon 10. The patient's longer fragment corresponded to an insertion of the last five nucleotides of AP4B1 intron 9. We confirmed that this variant affects the normal splicing of RNA, sustaining the molecular diagnosis of SPG47 in the patient.

Characterization of a Novel Splicing Variant in Acylglycerol Kinase (AGK) Associated with Fatal Sengers Syndrome

Mitochondrial functional integrity depends on protein and lipid homeostasis in the mitochondrial membranes and disturbances in their accumulation can cause disease. AGK, a mitochondrial acylglycerol kinase, is not only involved in lipid signaling but is also a component of the TIM22 complex in the inner mitochondrial membrane, which mediates the import of a subset of membrane proteins. AGK mutations can alter both phospholipid metabolism and mitochondrial protein biogenesis, contributing to the pathogenesis of Sengers syndrome. We describe the case of an infant carrying a novel homozygous AGK variant, c.518+1G>A, who was born with congenital cataracts, pielic ectasia, critical congenital dilated myocardiopathy, and hyperlactacidemia and died 20 h after birth. Using the patient’s DNA, we performed targeted sequencing of 314 nuclear genes encoding respiratory chain complex subunits and proteins implicated in mitochondrial oxidative phosphorylation (OXPHOS). A decrease of 96-bp in the length of the AGK cDNA sequence was detected. Decreases in the oxygen consumption rate (OCR) and the OCR:ECAR (extracellular acidification rate) ratio in the patient’s fibroblasts indicated reduced electron flow through the respiratory chain, and spectrophotometry revealed decreased activity of OXPHOS complexes I and V. We demonstrate a clear defect in mitochondrial function in the patient’s fibroblasts and describe the possible molecular mechanism underlying the pathogenicity of this novel AGK variant. Experimental validation using in vitro analysis allowed an accurate characterization of the disease-causing variant.

Compound heterozygous c.598_612del and c.1746-20C > G CAPN3 genotype cause autosomal recessive limb-girdle muscular dystrophy-1: a case report

BACKGROUND

Autosomal recessive limb-girdle muscular dystrophy-1 (LGMDR1), also known as calpainopathy, is a genetically heterogeneous disorder characterised by progression of muscle weakness. Homozygous or compound heterozygous variants in the CAPN3 gene are known genetic causes of this condition. The aim of this study was to confirm the molecular consequences of the CAPN3 variant NG_008660.1(NM_000070.3):c.1746-20C > G of an individual with suspected LGMDR1 by extensive complementary DNA (cDNA) analysis.

CASE PRESENTATION

In the present study, we report on a male with proximal muscular weakness in his lower limbs. Compound heterozygous NM_000070.3:c.598_612del and NG_008660.1(NM_000070.3):c.1746-20C > G genotype was detected on the CAPN3 gene by targeted next-generation sequencing (NGS). To confirm the pathogenicity of the variant c.1746-20C > G, we conducted genetic analysis based on Sanger sequencing of the proband's cDNA sample. The results revealed that this splicing variant disrupts the original 3' splice site on intron 13, thus leading to the skipping of the DNA fragment involving exon 14 and possibly exon 15. However, the lack of exon 15 in the CAPN3 isoforms present in a blood sample was explained by cell-specific alternative splicing rather than an aberrant splicing mechanism. In silico the c.1746-20C > G splicing variant consequently resulted in frameshift and formation of a premature termination codon (NP_000061.1:p.(Glu582Aspfs*62)).

CONCLUSIONS

Based on the results of our study and the literature we reviewed, both c.598_612del and c.1746-20C > G variants are pathogenic and together cause LGMDR1. Therefore, extensive mRNA and/or cDNA analysis of splicing variants is critical to understand the pathogenesis of the disease.

Reanalysis of Exome Data Identifies Novel SLC25A46 Variants Associated with Leigh Syndrome

SLC25A46 (solute carrier family 25 member 46) mutations have been linked to various neurological diseases with recessive inheritance, including Leigh syndrome, optic atrophy, and lethal congenital pontocerebellar hypoplasia. SLC25A46 is expressed in the outer membrane of mitochondria, where it plays a critical role in mitochondrial dynamics. A deceased 7-month-old female infant was suspected to have Leigh syndrome. Clinical exome sequencing was non-diagnostic, but research reanalysis of the sequencing data identified two novel variants in SLC25A46: a missense (c.1039C>T, p.Arg347Cys; NM_138773, hg19) and a donor splice region variant (c.283+5G>A) in intron 1. Both variants were predicted to be damaging. Sanger sequencing of cDNA detected a single missense allele in the patient compared to control, and the SLC25A46 transcript levels were also reduced due to the splice region variant. Additionally, Western blot analysis of whole-cell lysate showed a decrease of SLC25A46 expression in proband fibroblasts, relative to control cells. Further, analysis of mitochondrial morphology revealed evidence of increased fragmentation of the mitochondrial network in proband fibroblasts, compared to control cells. Collectively, our findings suggest that these novel variants in SLC24A46, the donor splice one and the missense variant, are the cause of the neurological phenotype in this proband.

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.

Evaluation of Suspected Autosomal Alport Syndrome Synonymous Variants

Background

Alport syndrome is an inherited disorder characterized by progressive renal disease, variable sensorineural hearing loss, and ocular abnormalities. Although many pathogenic variants in COL4A3 and COL4A4 have been identified in patients with autosomal Alport syndrome, synonymous mutations in these genes have rarely been identified.

Methods

We conducted in silico splicing analysis using Human Splicing Finder (HSF) and Alamut to predict splicing domain strength and disruption of the sites. Furthermore, we performed in vitro splicing assays using minigene constructs and mRNA analysis of patient samples to determine the pathogenicity of four synonymous variants detected in four patients with suspected autosomal dominant Alport syndrome (COL4A3 [c.693G>A (p.Val231=)] and COL4A4 [c.1353C>T (p.Gly451=), c.735G>A (p.Pro245=), and c.870G>A (p.Lys290=)]).

Results

Both in vivo and in vitro splicing assays showed exon skipping in two out of the four synonymous variants identified (c.735G>A and c.870G>A in COL4A4). Prediction analysis of wild-type and mutated COL4A4 sequences using HSF and Alamut suggested these two variants may lead to the loss of binding sites for several splicing factors, e.g., in acceptor sites and exonic splicing enhancers. The other two variants did not induce aberrant splicing.

Conclusions

This study highlights the pitfalls of classifying the functional consequences of variants by a simple approach. Certain synonymous variants, although they do not alter the amino acid sequence of the encoded protein, can dramatically affect pre-mRNA splicing, as shown in two of our patients. Our findings indicate that transcript analysis should be carried out to evaluate synonymous variants detected in patients with autosomal dominant Alport syndrome.

Regulation of alternative splicing in response to temperature variation in plants

Plants have evolved numerous molecular strategies to cope with perturbations in environmental temperature, and to adjust growth and physiology to limit the negative effects of extreme temperature. One of the strategies involves alternative splicing of primary transcripts to encode alternative protein products or transcript variants destined for degradation by nonsense-mediated decay. Here, we review how changes in environmental temperature-cold, heat, and moderate alterations in temperature-affect alternative splicing in plants, including crops. We present examples of the mode of action of various temperature-induced splice variants and discuss how these alternative splicing events enable favourable plant responses to altered temperatures. Finally, we point out unanswered questions that should be addressed to fully utilize the endogenous mechanisms in plants to adjust their growth to environmental temperature. We also indicate how this knowledge might be used to enhance crop productivity in the future.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

Contribution of Noncanonical Splice Variants to TTN Truncating Variant Cardiomyopathy

BACKGROUND

Heterozygous TTN truncating variants cause 10% to 20% of idiopathic dilated cardiomyopathy (DCM). Although variants which disrupt canonical splice signals (ie, invariant dinucleotide of the splice donor site, invariant dinucleotide of the splice acceptor site) at exon-intron junctions are readily recognized as TTN truncating variants, the effects of other nearby sequence variations on splicing and their contribution to disease is uncertain.

METHODS

Rare variants of unknown significance located in the splice regions of highly expressed TTN exons from 203 DCM cases, 3329 normal subjects, and clinical variant databases were identified. The effects of these variants on splicing were assessed using an in vitro splice assay.

RESULTS

Splice-altering variants of unknown significance were enriched in DCM cases over controls and present in 2% of DCM patients (P=0.002). Application of this method to clinical variant databases demonstrated 20% of similar variants of unknown significance in TTN splice regions affect splicing. Noncanonical splice-altering variants were most frequently located at position +5 of the donor site (P=4.4×10⁷) and position -3 of the acceptor site (P=0.002). SpliceAI, an emerging in silico prediction tool, had a high positive predictive value (86%-95%) but poor sensitivity (15%-50%) for the detection of splice-altering variants. Alternate exons spliced out of most TTN transcripts frequently lacked the consensus base at +5 donor and -3 acceptor positions.

CONCLUSIONS

Noncanonical splice-altering variants in TTN explain 1-2% of DCM and offer a 10-20% increase in the diagnostic power of TTN sequencing in this disease. These data suggest rules that may improve efforts to detect splice-altering variants in other genes and may explain the low percent splicing observed for many alternate TTN exons.

Mutation of c.244G>T in NR5A1 gene causing 46, XY DSD by affecting RNA splicing

Objective

To identify the pathogenic mechanism of the c.244G>T mutation in NR5A1 gene found in a Chinese patient with 46, XY disorders of sex development (DSD). Subjects and methods: Genomic DNA was extracted from a Chinese 46, XY DSD patient. Targeted next-generation and Sanger sequencing were performed to investigate and validate the gene mutation causing 46, XY DSD, respectively. In silico tools were used to predict the pathogenicity of the variant. Dual luciferase reporter gene assay and minigene splicing reporter assay were used to identify the pathogenicity of the variant.

Results

A novel heterozygous variant, c.244G>T (p.Ala82Ser), in NR5A1 gene was detected in the 46, XY DSD patient. Four of five silico tools predicting pathogenicity of missense variants indicated that the variant was pathogenic. However, in vitro functional study showed that p.Ala82Ser did not affect the transcriptional activity of NR5A1. In silico tools predicting the potential splicing loci revealed that c.244G>T led to aberrant splicing of NR5A1 RNA. Minigene splicing reporter assay confirmed that c.244G>T resulted in the deletion of exon2 or deletion of 19 nucleotides in 3′ end of exon2.

Conclusions

Mutation of c.244G>T in NR5A1 results in 46, XY DSD by inducing abnormal splicing of NR5A1 RNA instead of amino acid substitution of NR5A1.

Fucosidosis in Tunisian patients: mutational analysis and homology-based modeling of FUCA1 enzyme

BACKGROUND

Fucosidosis is an autosomal recessive lysosomal storage disease caused by defective alpha-L-fucosidase (FUCA1) activity, leading to the accumulation of fucose-containing glycolipids and glycoproteins in various tissues. Clinical features include angiokeratoma, progressive psychomotor retardation, neurologic signs, coarse facial features, and dysostosis multiplex.

METHODS

All exons and flanking intron regions of FUCA1 were screened by direct sequencing to identify mutations and polymorphisms in three unrelated families with fucosidosis. Bioinformatics tools were then used to predict the impacts of novel alterations on the structure and function of proteins. Furthermore, the identified mutations were localized onto a 3D structure model using the DeepView Swiss-PdbViewer 4.1 software, which established a function-structure relationship of the FUCA1 proteins.

RESULTS

Four novel mutations were identified in this study. Two patients (P1 and P2) in Families 1 and 2 who had the severe phenotype were homoallelic for the two identified frameshift mutations p.K57Sfs*75 and p.F77Sfs*55, respectively. The affected patient (P3) from Family 3, who had the milder phenotype, was heterozygous for the novel missense mutation p.G332E and the novel splice site mutation c.662+5g>c. We verified that this sequence variation did not correspond to a polymorphism by testing 50 unrelated individuals. Additionally, 16 FUCA1 polymorphisms were identified. The structure prediction analysis showed that the missense mutation p.G332E would probably lead to a significant conformational change, thereby preventing the expression of the FUCA1 protein indeed; the 3D structural model of the FUCA1 protein reveals that the glycine at position 332 is located near a catalytic nucleophilic residue. This makes it likely that the enzymatic function of the protein with p.G332E is severely impaired.

CONCLUSION

These are the first FUCA1 mutations identified in Tunisia that cause the fucosidosis disease. Bioinformatics analysis allowed us to establish an approximate structure-function relationship for the FUCA1 protein, thereby providing better genotype/phenotype correlation knowledge.

Induction of cryptic pre-mRNA splice-switching by antisense oligonucleotides

Antisense oligomers (AOs) are increasingly being used to modulate RNA splicing in live cells, both for research and for the development of therapeutics. While the most common intended effect of these AOs is to induce skipping of whole exons, rare examples are emerging of AOs that induce skipping of only part of an exon, through activation of an internal cryptic splice site. In this report, we examined seven AO-induced cryptic splice sites in six genes. Five of these cryptic splice sites were discovered through our own experiments, and two originated from other published reports. We modelled the predicted effects of AO binding on the secondary structure of each of the RNA targets, and how these alterations would in turn affect the accessibility of the RNA to splice factors. We observed that a common predicted effect of AO binding was disruption of the exon definition signal within the exon's excluded segment.

Novel splicing (c.6529-1G>T) and missense (c.1667G>A) mutations causing factor V deficiency

Congenital factor V deficiency (FVD) is a rare bleeding disorder. In this study, we investigated the genetic basis in an African American patient with factor V activity 3%. Custom sequence capture and targeted next-generation (NGS) sequencing of the F5 gene were undertaken followed by PCR and Sanger sequencing. Two novel variants were identified. In silico analyses correlated clinically with the patient's factor V activity and hemorrhagic tendency. A review of the literature regarding these genomic alterations is presented. We described two novel mutations causing moderate FVD. The first, Chr1:g.169483698C>A with cDNA change (F5):c.6529-1G>T, occurred in a conserved nucleotide at the canonical acceptor splice site of intron 24. The second, Chr1:g.169515775C>T with cDNA change (F5):c.1667G>A, was a missense variant of exon 11, affecting a highly conserved amino acid in the A2 domain. Further research into the mechanisms of F5 mutations leading to FVD and residual factor V expression are needed.

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.

Mutation spectrum of the OPA1 gene in a large cohort of patients with suspected dominant optic atrophy: Identification and classification of 48 novel variants

Autosomal dominant optic atrophy is one of the most common inherited optic neuropathies. This disease is genetically heterogeneous, but most cases are due to pathogenic variants in the OPA1 gene: depending on the population studied, 32–90% of cases harbor pathogenic variants in this gene. The aim of this study was to provide a comprehensive overview of the entire spectrum of likely pathogenic variants in the OPA1 gene in a large cohort of patients. Over a period of 20 years, 755 unrelated probands with a diagnosis of bilateral optic atrophy were referred to our laboratory for molecular genetic investigation. Genetic testing of the OPA1 gene was initially performed by a combined analysis using either single-strand conformation polymorphism or denaturing high performance liquid chromatography followed by Sanger sequencing to validate aberrant bands or melting profiles. The presence of copy number variations was assessed using multiplex ligation-dependent probe amplification. Since 2012, genetic testing was based on next-generation sequencing platforms. Genetic screening of the OPA1 gene revealed putatively pathogenic variants in 278 unrelated probands which represent 36.8% of the entire cohort. A total of 156 unique variants were identified, 78% of which can be considered null alleles. Variant c.2708_2711del/p.(V903Gfs*3) was found to constitute 14% of all disease-causing alleles. Special emphasis was placed on the validation of splice variants either by analyzing cDNA derived from patients´ blood samples or by heterologous splice assays using minigenes. Splicing analysis revealed different aberrant splicing events, including exon skipping, activation of exonic or intronic cryptic splice sites, and the inclusion of pseudoexons. Forty-eight variants that we identified were novel. Nine of them were classified as pathogenic, 34 as likely pathogenic and five as variant of uncertain significance. Our study adds a significant number of novel variants to the mutation spectrum of the OPA1 gene and will thereby facilitate genetic diagnostics of patients with suspected dominant optic atrophy.

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.

MLH1 intronic variants mapping to + 5 position of splice donor sites lead to deleterious effects on RNA splicing

Germline pathogenic variants in the DNA mismatch repair genes (MMR): MLH1, MSH2, MSH6, and PMS2, are causative of Lynch syndrome (LS). However, many of the variants mapping outside the invariant splice site positions (IVS ± 1, IVS ± 2) are classified as variants of unknown significance (VUS). Three such variants (MLH1 c.588+5G>C, c.588+5G>T and c.677+5G>A) were identified in 8 unrelated LS families from Argentina, Brazil and Chile. Herein, we collected clinical information on these families and performed segregation analysis and RNA splicing studies to assess the implication of these VUS in LS etiology. Pedigrees showed a clear pattern of variant co-segregation with colorectal cancer and/or other LS-associated malignancies. Tumors presented deficient expression of MLH1-PMS2 proteins in 7/7 of the LS families, and MSI-high status in 3/3 cases. Moreover, RNA analyses revealed that c.588+5G>C and c.588+5G>T induce skipping of exon 7 whereas c.677+5G>A causes skipping of exon 8. In sum, we report that the combined clinical findings in the families and the molecular studies provided the evidences needed to demonstrate that the three MLH1 variants are causative of LS and to classify c.588+5G>C and c.677+5G>A as class 5 (pathogenic), and c.588+5G>T as class 4 (likely-pathogenic). Our findings underline the importance of performing clinical and family analyses, as well as RNA splicing assays in order to determine the clinical significance of intronic variants, and contribute to the genetic counseling and clinical management of patients and their relatives.

Antisense Oligonucleotide Therapy for Ophthalmic Conditions

Antisense oligonucleotides (AON) are synthetic single-stranded fragments of nucleic acids that bind to a specific complementary messenger RNA (mRNA) sequence and change the final gene product. AON were initially approved for treating cytomegalovirus retinitis and have shown promise in treating Mendelian systemic disease. AON are currently being investigated as a treatment modality for many ophthalmic diseases, including inherited retinal disorders (IRD), inflammatory response and wound healing after glaucoma surgery, and macular degeneration. They provide a possible solution to gene therapy for IRD that are not candidates for adeno-associated virus (AAV) delivery. This chapter outlines the historical background of AON and reviews clinical applications and ongoing clinical trials.

Therapeutic Modulation of RNA Splicing in Malignant and Non-Malignant Disease

RNA splicing is the enzymatic process by which non-protein coding sequences are removed from RNA to produce mature protein-coding mRNA. Splicing is thereby a major mediator of proteome diversity as well as a dynamic regulator of gene expression. Genetic alterations disrupting splicing of individual genes or altering the function of splicing factors contribute to a wide range of human genetic diseases as well as cancer. These observations have resulted in the development of therapies based on oligonucleotides that bind to RNA sequences and modulate splicing for therapeutic benefit. In parallel, small molecules that bind to splicing factors to alter their function or modify RNA processing of individual transcripts are being pursued for monogenic disorders as well as for cancer.

Generation of pathogenic TPP1 mutations in human stem cells as a model for neuronal ceroid lipofuscinosis type 2 disease

Neuronal ceroid lipofuscinosis type 2 (CLN2 disease) is an autosomal recessive neurodegenerative disorder generally with onset at 2 to 4 years of age and characterized by seizures, loss of vision, progressive motor and mental decline, and premature death. CLN2 disease is caused by loss-of-function mutations in the tripeptidyl peptidase 1 (TPP1) gene leading to deficiency in TPP1 enzyme activity. Approximately 60% of patients have one of two pathogenic variants (c.509–1G > C or c.622C > T [p.(Arg208*)]). In order to generate a human stem cell model of CLN2 disease, we used CRISPR/Cas9-mediated knock-in technology to introduce these mutations in a homozygous state into H9 human embryonic stem cells. Heterozygous lines of the c.622C > T (p.(Arg208*)) mutation were also generated, which included a heterozygous mutant with a wild-type allele and different compound heterozygous coding mutants resulting from indels on one allele. We describe the methodology that led to the generation of the lines and provide data on the initial validation and characterization of these CLN2 disease models. Notably, both mutant lines (c.509–1G > C and c.622C > T [p.(Arg208*)]) in the homozygous state were shown to have reduced or absent protein, respectively, and deficiency of TPP1 enzyme activity. These models, which we have made available for wide-spread sharing, will be useful for future studies of molecular and cellular mechanisms underlying CLN2 disease and for therapeutic development.

Activation of cryptic donor splice sites by non-coding and coding PAX6 variants contributes to congenital aniridia

BACKGROUND

The paired-domain transcription factor paired box gene 6 (PAX6) causes a wide spectrum of ocular developmental anomalies, including congenital aniridia, Peters anomaly and microphthalmia. Here, we aimed to functionally assess the involvement of seven potentially non-canonical splicing variants on missplicing of exon 6, which represents the main hotspot region for loss-of-function PAX6 variants.

METHODS

By locus-specific analysis of PAX6 using Sanger and/or targeted next-generation sequencing, we screened a Spanish cohort of 106 patients with PAX6-related diseases. Functional splicing assays were performed by in vitro minigene approaches or directly in RNA from patient-derived lymphocytes cell line, when available.

RESULTS

Five out seven variants, including three synonymous changes, one small exonic deletion and one non-canonical splice variant, showed anomalous splicing patterns yielding partial exon skipping and/or elongation.

CONCLUSION

We describe new spliceogenic mechanisms for PAX6 variants mediated by creating or strengthening five different cryptic donor sites at exon 6. Our work revealed that the activation of cryptic PAX6 splicing sites seems to be a recurrent and underestimated cause of aniridia. Our findings pointed out the importance of functional assessment of apparently silent PAX6 variants to uncover hidden genetic alterations and to improve variant interpretation for genetic counselling in aniridia.

Improved Variant Detection in Clinical Myeloid NGS Testing by Supplementing a Commercial Myeloid NGS Assay with Custom or Extended Data Filtering and Accessory Fragment Analysis

BACKGROUND

Commercial myeloid next-generation sequencing (NGS) panels may facilitate uniform generation of raw data between laboratories. However, different strategies for data filtering and variant annotation may contribute to differences in variant detection and reporting. Here, we present how custom data filtering or the use of Oncomine extended data filtering improve detection of clinically relevant mutations with the Oncomine Myeloid Research Assay.

METHODS

The study included all patient samples (n = 264) analyzed during the first-year, single-site, clinical use of the Ion Torrent Oncomine Myeloid Research Assay. In data analysis, the default analysis filter was supplemented with our own data filtering algorithm in order to detect additional clinically relevant mutations. In addition, we developed a sensitive supplementary test for the ASXL1 c.1934dupG p.Gly646fs mutation by fragment analysis.

RESULTS

Using our custom filter chain, we found 96 different reportable variants that were not detected by the default filter chain. Twenty-six of these were classified as variants of strong or potential clinical significance (tier I/tier II variants), and the custom filtering discovered otherwise undetected tier I/tier II variants in 25 of 132 patients with clinically relevant mutations (19%). The remaining 70 variants not detected by the default filter chain were classified as variants of unknown significance. Among these were several unique variants with possible pathogenic potential judged by bioinformatic predictions. The recently launched Oncomine 5.14 extended filter algorithm detects most but not all of the tier I/tier II variants that were not detected by the default filter. The supplementary fragment analysis for the ASXL1 c.1934dupG p.Gly646fs confidently detected a variant allele frequency of down to 4.8% (SD 0.83%). The assay also detected the ASXL1 c.1900_1922del23 mutation.

CONCLUSION

Detection of clinically relevant variants with the Oncomine Myeloid Research NGS assay can be significantly improved by supplementing the default filter chain with custom data filtering or the recently launched Oncomine 5.14 extended filter algorithm. Our accessory fragment analysis facilitates easy testing for frequent ASXL1 mutations that are poorly or not covered by the NGS assay.

Spliceosomal snRNA Epitranscriptomics

Small nuclear RNAs (snRNAs) are critical components of the spliceosome that catalyze the splicing of pre-mRNA. snRNAs are each complexed with many proteins to form RNA-protein complexes, termed as small nuclear ribonucleoproteins (snRNPs), in the cell nucleus. snRNPs participate in pre-mRNA splicing by recognizing the critical sequence elements present in the introns, thereby forming active spliceosomes. The recognition is achieved primarily by base-pairing interactions (or nucleotide-nucleotide contact) between snRNAs and pre-mRNA. Notably, snRNAs are extensively modified with different RNA modifications, which confer unique properties to the RNAs. Here, we review the current knowledge of the mechanisms and functions of snRNA modifications and their biological relevance in the splicing process.

Identification of Deep-Intronic Splice Mutations in a Large Cohort of Patients With Inherited Retinal Diseases

High throughput sequencing technologies have revolutionized the identification of mutations responsible for a diverse set of Mendelian disorders, including inherited retinal disorders (IRDs). However, the causal mutations remain elusive for a significant proportion of patients. This may be partially due to pathogenic mutations located in non-coding regions, which are largely missed by capture sequencing targeting the coding regions. The advent of whole-genome sequencing (WGS) allows us to systematically detect non-coding variations. However, the interpretation of these variations remains a significant bottleneck. In this study, we investigated the contribution of deep-intronic splice variants to IRDs. WGS was performed for a cohort of 571 IRD patients who lack a confident molecular diagnosis, and potential deep intronic variants that affect proper splicing were identified using SpliceAI. A total of six deleterious deep intronic variants were identified in eight patients. An in vitro minigene system was applied to further validate the effect of these variants on the splicing pattern of the associated genes. The prediction scores assigned to splice-site disruption positively correlated with the impact of mutations on splicing, as those with lower prediction scores demonstrated partial splicing. Through this study, we estimated the contribution of deep-intronic splice mutations to unassigned IRD patients and leveraged in silico and in vitro methods to establish a framework for prioritizing deep intronic variant candidates for mechanistic and functional analyses.

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

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