ESEfinder: A web resource to identify exonic splicing enhancers

Unique Ataxia-Oculomotor Apraxia 2 (AOA2) in Israel with Novel Variants, Atypical Late Presentation, and Possible Identification of a Poison Exon

AOA2 is a rare progressive adolescent-onset disease characterised by cerebellar vermis atrophy, peripheral neuropathy and elevated serum alpha-fetoprotein (AFP) caused by pathogenic bi-allelic variants in SETX, encoding senataxin, involved in DNA repair and RNA maturation. Sanger sequencing of genomic DNA, co-segregation and oxidative stress functional studies were performed in Family 1. Trio whole-exome sequencing (WES), followed by SETX RNA and qRT-PCR analysis, were performed in Family 2. Sanger sequencing in Family 1 revealed two novel in-frame SETX deletion and duplication variants in trans (c.7009_7011del; p.Val2337del and c.7369_7371dup; p.His2457dup). Patients had increased induced chromosomal aberrations at baseline and following exposure to higher mitomycin-C concentration and increased sensitivity to oxidative stress at the lower mitomycin-C concentration in cell viability test. Trio WES in Family 2 revealed two novel SETX variants in trans, a nonsense variant (c.568C > T; p.Gln190*), and a deep intronic variant (c.5549-107A > G). Intronic variant analysis and SETX mRNA expression revealed activation of a cryptic exon introducing a premature stop codon (p.Met1850Lysfs*18) and resulting in aberrant splicing, as shown by qRT-PCR analysis, thus leading to higher levels of cryptic exon activation. Along with a second deleterious allele, this variant leads to low levels of SETX mRNA and disease manifestations. Our report expands the phenotypic spectrum of AOA2. Results provide initial support for the hypomorphic nature of the novel in-frame deletion and duplication variants in Family 1. Deep-intronic variant analysis of Family 2 variants potentially reveals a previously undescribed poison exon in the SETX gene, which may contribute to tailored therapy development.

Noncanonical Splice Site and Deep Intronic FRMD7 Variants Activate Cryptic Exons in X-linked Infantile Nystagmus

Purpose

We aim to report noncoding pathogenic variants in patients with FRMD7-related infantile nystagmus (FIN).

Methods

Genome sequencing (n = 2 families) and reanalysis of targeted panel next generation sequencing (n = 2 families) was performed in genetically unsolved cases of suspected FIN. Previous sequence analysis showed no pathogenic coding variants in genes associated with infantile nystagmus. SpliceAI, SpliceRover, and Alamut consensus programs were used to annotate noncoding variants. Minigene splicing assay was performed to confirm aberrant splicing. In silico analysis of exonic splicing enhancer and silencer was also performed.

Results

FRMD7 intronic variants were identified based on genome sequencing and targeted next-generation sequencing analysis. These included c.285-12A>G (pedigree 1), c.284+63T>A (pedigrees 2 and 3), and c. 383-1368A>G (pedigree 4). All variants were absent in gnomAD, and the both c.285-12A>G and c.284+63T>A variants were predicted to enhance new splicing acceptor gains with SpliceAI, SpliceRover, and Alamut consensus approaches. However, the c.383-1368 A>G variant only had a significant impact score on the SpliceRover program. The c.383-1368A>G variant was predicted to promote pseudoexon inclusion by binding of exonic splicing enhancer. Aberrant exonizations were validated through minigene constructs, and all variants were segregated in the families.

Conclusions

Deep learning–based annotation of noncoding variants facilitates the discovery of hidden genetic variations in patients with FIN. This study provides evidence of effectiveness of combined deep learning–based splicing tools to identify hidden pathogenic variants in previously unsolved patients with infantile nystagmus.

Translational Relevance

These results demonstrate robust analysis using two deep learning splicing predictions and in vitro functional study can lead to finding hidden genetic variations in unsolved patients.

Gene-specific nonsense-mediated mRNA decay targeting for cystic fibrosis therapy

Low CFTR mRNA expression due to nonsense-mediated mRNA decay (NMD) is a major hurdle in developing a therapy for cystic fibrosis (CF) caused by the W1282X mutation in the CFTR gene. CFTR-W1282X truncated protein retains partial function, so increasing its levels by inhibiting NMD of its mRNA will likely be beneficial. Because NMD regulates the normal expression of many genes, gene-specific stabilization of CFTR-W1282X mRNA expression is more desirable than general NMD inhibition. Synthetic antisense oligonucleotides (ASOs) designed to prevent binding of exon junction complexes (EJC) downstream of premature termination codons (PTCs) attenuate NMD in a gene-specific manner. We describe cocktails of three ASOs that specifically increase the expression of CFTR-W1282X mRNA and CFTR protein upon delivery into human bronchial epithelial cells. This treatment increases the CFTR-mediated chloride current. These results set the stage for clinical development of an allele-specific therapy for CF caused by the W1282X mutation.

The W1282X nonsense mutation in the CFTR gene causes cystic fibrosis by reducing its mRNA and functional protein levels. Here the authors developed antisense-oligonucleotide cocktails that restore CFTR protein function by gene-specific stabilization of CFTR mRNA.

Innovative in Silico Approaches for Characterization of Genes and Proteins

Bioinformatics is an amalgamation of biology, mathematics and computer science. It is a science which gathers the information from biology in terms of molecules and applies the informatic techniques to the gathered information for understanding and organizing the data in a useful manner. With the help of bioinformatics, the experimental data generated is stored in several databases available online like nucleotide database, protein databases, GENBANK and others. The data stored in these databases is used as reference for experimental evaluation and validation. Till now several online tools have been developed to analyze the genomic, transcriptomic, proteomics, epigenomics and metabolomics data. Some of them include Human Splicing Finder (HSF), Exonic Splicing Enhancer Mutation taster, and others. A number of SNPs are observed in the non-coding, intronic regions and play a role in the regulation of genes, which may or may not directly impose an effect on the protein expression. Many mutations are thought to influence the splicing mechanism by affecting the existing splice sites or creating a new sites. To predict the effect of mutation (SNP) on splicing mechanism/signal, HSF was developed. Thus, the tool is helpful in predicting the effect of mutations on splicing signals and can provide data even for better understanding of the intronic mutations that can be further validated experimentally. Additionally, rapid advancement in proteomics have steered researchers to organize the study of protein structure, function, relationships, and dynamics in space and time. Thus the effective integration of all of these technological interventions will eventually lead to steering up of next-generation systems biology, which will provide valuable biological insights in the field of research, diagnostic, therapeutic and development of personalized medicine.

Synonymous Variants: Necessary Nuance in our Understanding of Cancer Drivers and Treatment Outcomes

Once called "silent mutations" and assumed to have no effect on protein structure and function, synonymous variants are now recognized to be drivers for some cancers. There have been significant advances in our understanding of the numerous mechanisms by which synonymous single nucleotide variants (sSNVs) can affect protein structure and function by affecting pre-mRNA splicing, mRNA expression, stability, folding, miRNA binding, translation kinetics, and co-translational folding. This review highlights the need for considering sSNVs in cancer biology to gain a better understanding of the genetic determinants of human cancers and to improve their diagnosis and treatment. We surveyed the literature for reports of sSNVs in cancer and found numerous studies on the consequences of sSNVs on gene function with supporting in vitro evidence. We also found reports of sSNVs that have statistically significant associations with specific cancer types but for which in vitro studies are lacking to support the reported associations. Additionally, we found reports of germline and somatic sSNVs that were observed in numerous clinical studies and for which in silico analysis predicts possible effects on gene function. We provide a review of these investigations and discuss necessary future studies to elucidate the mechanisms by which sSNVs disrupt protein function and are play a role in tumorigeneses, cancer progression, and treatment efficacy. As splicing dysregulation is one of the most well recognized mechanisms by which sSNVs impact protein function, we also include our own in silico analysis for predicting which sSNVs may disrupt pre-mRNA splicing.

SRSF6 Regulates the Alternative Splicing of the Apoptotic Fas Gene by Targeting a Novel RNA Sequence

Simple Summary

Alternative splicing (AS) produces multiple mRNA isoforms from a gene to make a large number of proteins. Fas (Apo-1/CD95) pre-mRNA, a member of TNF receptor family that mediates apoptosis, can generate pro-apoptotic and anti-apoptotic proteins through AS. Here, we identified SRSF6 as an essential regulator protein in Fas AS. We further located a new functional target sequence of SRSF6 in Fas splicing. In addition, our large-scale RNA-seq analysis using GTEX and TCGA indicated that while SRSF6 expression was correlated with Fas expression in normal tissues, the correlation was disrupted in tumors. Our results suggest a novel regulatory mechanisms of Fas AS.

Abstract

Alternative splicing (AS) is a procedure during gene expression that allows the production of multiple mRNAs from a single gene, leading to a larger number of proteins with various functions. The alternative splicing (AS) of Fas (Apo-1/CD95) pre-mRNA can generate membrane-bound or soluble isoforms with pro-apoptotic and anti-apoptotic functions. SRSF6, a member of the Serine/Arginine-rich protein family, plays essential roles in both constitutive and alternative splicing. Here, we identified SRSF6 as an important regulatory protein in Fas AS. The cassette exon inclusion of Fas was decreased by SRSF6-targeting shRNA treatment, but increased by SRSF6 overexpression. The deletion and substitution mutagenesis of the Fas minigene demonstrated that the UGCCAA sequence in the cassette exon of the Fas gene causes the functional disruption of SRSF6, indicating that these sequences are essential for SRSF6 function in Fas splicing. In addition, biotin-labeled RNA-pulldown and immunoblotting analysis showed that SRSF6 interacted with these RNA sequences. Mutagenesis in the splice-site strength alteration demonstrated that the 5′ splice-site, but not the 3′ splice-site, was required for the SRSF6 regulation of Fas pre-mRNA. In addition, a large-scale RNA-seq analysis using GTEX and TCGA indicated that while SRSF6 expression was correlated with Fas expression in normal tissues, the correlation was disrupted in tumors. Furthermore, high SRSF6 expression was linked to the high expression of pro-apoptotic and immune activation genes. Therefore, we identified a novel RNA target with 5′ splice-site dependence of SRSF6 in Fas pre-mRNA splicing, and a correlation between SRSF6 and Fas expression.

High Concentration or Combined Treatment of Antisense Oligonucleotides for Spinal Muscular Atrophy Perturbed SMN2 Splicing in Patient Fibroblasts

Spinal muscular atrophy (SMA) is caused by survival motor neuron 1 SMN1 deletion. The survival motor neuron 2 (SMN2) encodes the same protein as SMN1 does, but it has a splicing defect of exon 7. Some antisense oligonucleotides (ASOs) have been proven to correct this defect. One of these, nusinersen, is effective in SMA-affected infants, but not as much so in advanced-stage patients. Furthermore, the current regimen may exhibit a ceiling effect. To overcome these problems, high-dose ASOs or combined ASOs have been explored. Here, using SMA fibroblasts, we examined the effects of high-concentration ASOs and of combining two ASOs. Three ASOs were examined: one targeting intronic splicing suppressor site N1 (ISS-N1) in intron 7, and two others targeting the 3′ splice site and 5′ region of exon 8. In our experiments on all ASO types, a low or intermediate concentration (50 or 100 nM) showed better splicing efficiency than a high concentration (200 nM). In addition, a high concentration of each ASO created a cryptic exon in exon 6. When a mixture of two different ASOs (100 nM each) was added to the cells, the cryptic exon was included in the mRNA. In conclusion, ASOs at a high concentration or used in combination may show less splicing correction and cryptic exon creation.

Autosomal dominant polycystic kidney disease (ADPKD) in Tunisia: From molecular genetics to the development of prognostic tools

A high prevalence of genetic kidney disease in Tunisia has been detected, and their study provides very important clinical and genetic information. Autosomal dominant polycystic kidney disease (ADPKD) is one of the main causes of morbidity and mortality associated with the kidneys in Tunisia. We present here clinical and genetic characteristics of a cohort of Tunisian patients with ADPKD. Nineteen Tunisian patients with ADPKD, among 4 familial cases and 11 sporadic cases, and 50 Healthy individuals were included in this cohort. Genetic studies of PKD1/2 were carried on using Sanger sequencing and MLPA. In our study, the mean age at diagnosis was 47 ± 18 years. In addition, 84.21% of cases present a family history of ADPKD. Overall, 57.89% of the affected individuals had HTA and 26.31% patients had hematuria. 15.78 % of the patient has extra-renal cysts i.e. one patient with splenic cysts and two patients had liver cysts. 57.89 % of patients were diagnosed with various extra-renal clinical presentations i.e. myopia, hernia, deafness, intracranial aneurysm, respiratory distress, hyperthyroidism, urinary tract infection and lower back pains. The PKD1 genotype showed earlier onset of ESRD compared to PKD2 genotype (43 vs. 55 years old). Six mutations have been detected in PKD1 gene. Among them, three were novels e.g. c.688 T>G, p.C230G and c.690C>G, p.C230W among exon 5 and c.8522A>G, p.N2841S among exon 23. In addition, thirteen single nucleotides polymorphisms have been reported in PKD1 gene. Among them, eleven previously reported in heterozygous state and two novel single nucleotides polymorphisms in heterozygous and homozygous state and predicted to be probable polymorphisms by computational tools: c.496C>T, p.L166= among the exon 4, and c.10165G>C and p.E3389Gln among the exon 31. Only three single nucleotides polymorphisms previously reported in ADPKD database have been identified in PKD2 gene. The description and analysis of our cohort can help in rapid and reliable diagnosis for early management of patients in Tunisia. Indeed, predictive genetic testing can facilitate donor evaluation and increase living related kidney transplantation.

Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts

Spinal muscular atrophy (SMA) is a severe, debilitating neuromuscular condition characterised by loss of motor neurons and progressive muscle wasting. SMA is caused by a loss of expression of SMN1 that encodes the survival motor neuron (SMN) protein necessary for the survival of motor neurons. Restoration of SMN expression through increased inclusion of SMN2 exon 7 is known to ameliorate symptoms in SMA patients. As a consequence, regulation of pre-mRNA splicing of SMN2 could provide a potential molecular therapy for SMA. In this study, we explored if splice switching antisense oligonucleotides could redirect the splicing repressor hnRNPA1 to the hnRNPA1b isoform and restore SMN expression in fibroblasts from a type I SMA patient. Antisense oligonucleotides (AOs) were designed to promote exon 7b retention in the mature mRNA and induce the hnRNPA1b isoform. RT-PCR and western blot analysis were used to assess and monitor the efficiency of different AO combinations. A combination of AOs targeting multiple silencing motifs in hnRNPA1 pre-mRNA led to robust hnRNPA1b induction, which, in turn, significantly increased expression of full-length SMN (FL-SMN) protein. A combination of PMOs targeting the same motifs also strongly induced hnRNPA1b isoform, but surprisingly SMN2 exon 5 skipping was detected, and the PMO cocktail did not lead to a significant increase in expression of FL-SMN protein. We further performed RNA sequencing to assess the genome-wide effects of hnRNPA1b induction. Some 3244 genes were differentially expressed between the hnRNPA1b-induced and untreated SMA fibroblasts, which are functionally enriched in cell cycle and chromosome segregation processes. RT-PCR analysis demonstrated that expression of the master regulator of these enrichment pathways, MYBL2 and FOXM1B, were reduced in response to PMO treatment. These findings suggested that induction of hnRNPA1b can promote SMN protein expression, but not at sufficient levels to be clinically relevant.

ABCA4 c.859-25A>G, a Frequent Palestinian Founder Mutation Affecting the Intron 7 Branchpoint, Is Associated With Early-Onset Stargardt Disease

Purpose

The effect of noncoding variants is often unknown in the absence of functional assays. Here, we characterized an ABCA4 intron 7 variant, c.859-25A>G, identified in Palestinian probands with Stargardt disease (STGD) or cone-rod dystrophy (CRD). We investigated the effect of this variant on the ABCA4 mRNA and retinal phenotype, and its prevalence in Palestine.

Methods

The ABCA4 gene was sequenced completely or partially in 1998 cases with STGD or CRD. The effect of c.859-25A>G on splicing was investigated in silico using SpliceAI and in vitro using splice assays. Homozygosity mapping was performed for 16 affected individuals homozygous for c.859-25A>G. The clinical phenotype was assessed using functional and structural analyses including visual acuity, full-field electroretinography, and multimodal imaging.

Results

The smMIPs-based ABCA4 sequencing revealed c.859-25A>G in 10 Palestinian probands from Hebron and Jerusalem. SpliceAI predicted a significant effect of this putative branchpoint-inactivating variant on the nearby intron 7 splice acceptor site. Splice assays revealed exon 8 skipping and two partial inclusions of intron 7, each having a deleterious effect. Additional genotyping revealed another 46 affected homozygous or compound heterozygous individuals carrying variant c.859-25A>G. Homozygotes shared a genomic segment of 59.6 to 87.9 kb and showed severe retinal defects on ophthalmoscopic evaluation.

Conclusions

The ABCA4 variant c.859-25A>G disrupts a predicted branchpoint, resulting in protein truncation because of different splice defects, and is associated with early-onset STGD1 when present in homozygosity. This variant was found in 25/525 Palestinian inherited retinal dystrophy probands, representing one of the most frequent inherited retinal disease-causing variants in West-Bank Palestine.

Whole-transcriptome sequencing-based concomitant detection of viral and human genetic determinants of cutaneous lesions

Severe viral infections of the skin can occur in patients with inborn errors of immunity (IEI). We report an all-in-one whole-transcriptome sequencing–based method by RNA-Seq on a single skin biopsy for concomitantly identifying the cutaneous virome and the underlying IEI. Skin biopsies were obtained from healthy and lesional skin from patients with cutaneous infections suspected to be of viral origin. RNA-Seq was utilized as the first-tier strategy for unbiased human genome-wide rare variant detection. Reads unaligned to the human genome were utilized for the exploration of 926 viruses in a viral genome catalog. In 9 families studied, the patients carried pathogenic variants in 6 human IEI genes, including IL2RG, WAS, CIB1, STK4, GATA2, and DOCK8. Gene expression profiling also confirmed pathogenicity of the human variants and permitted genome-wide homozygosity mapping, which assisted in identification of candidate genes in consanguineous families. This automated, online, all-in-one computational pipeline, called VirPy, enables simultaneous detection of the viral triggers and the human genetic variants underlying skin lesions in patients with suspected IEI and viral dermatosis.

Silent mutations reveal therapeutic vulnerability in RAS Q61 cancers

RAS family members are the most frequently mutated oncogenes in human cancers. Although KRAS(G12C)-specific inhibitors show clinical activity in patients with cancer^1-3, there are no direct inhibitors of NRAS, HRAS or non-G12C KRAS variants. Here we uncover the requirement of the silent KRAS^G60G mutation for cells to produce a functional KRAS(Q61K). In the absence of this G60G mutation in KRAS^Q61K, a cryptic splice donor site is formed, promoting alternative splicing and premature protein termination. A G60G silent mutation eliminates the splice donor site, yielding a functional KRAS(Q61K) variant. We detected a concordance of KRAS^Q61K and a G60G/A59A silent mutation in three independent pan-cancer cohorts. The region around RAS Q61 is enriched in exonic splicing enhancer (ESE) motifs and we designed mutant-specific oligonucleotides to interfere with ESE-mediated splicing, rendering the RAS(Q61) protein non-functional in a mutant-selective manner. The induction of aberrant splicing by antisense oligonucleotides demonstrated therapeutic effects in vitro and in vivo. By studying the splicing necessary for a functional KRAS(Q61K), we uncover a mutant-selective treatment strategy for RAS^Q61 cancer and expose a mutant-specific vulnerability, which could potentially be exploited for therapy in other genetic contexts.

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

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

Exon Skipping Through Chimeric Antisense U1 snRNAs to Correct Retinitis Pigmentosa GTPase-Regulator (RPGR) Splice Defect

Inherited retinal dystrophies are caused by mutations in more than 250 genes, each of them carrying several types of mutations that can lead to different clinical phenotypes. Mutations in Retinitis Pigmentosa GTPase-Regulator (RPGR) cause X-linked Retinitis pigmentosa (RP). A nucleotide substitution in intron 9 of RPGR causes the increase of an alternatively spliced isoform of the mature mRNA, bearing exon 9a (E9a). This introduces a stop codon, leading to truncation of the protein. Aiming at restoring impaired gene expression, we developed an antisense RNA-based therapeutic approach for the skipping of RPGR E9a. We designed a set of specific U1 antisense snRNAs (U1_asRNAs) and tested their efficacy in vitro, upon transient cotransfection with RPGR minigene reporter systems in HEK-293T, 661W, and PC-12 cell lines.

Cytoprotective Activity of Polyamines Is Associated with the Alternative Splicing of RAD51A Pre-mRNA in Normal Human CD4+ T Lymphocytes

Physiological polyamines are ubiquitous polycations with pleiotropic biochemical activities, including regulation of gene expression and cell proliferation as well as modulation of cell signaling. They can also decrease DNA damage and promote cell survival. In the present study, we demonstrated that polyamines have cytoprotective effects on normal human CD4+ T lymphocytes but not on cancer Jurkat or K562 cells. Pretreatment of lymphocytes with polyamines resulted in a significant reduction in cells with DNA damage induced by doxorubicin, cisplatin, or irinotecan, leading to an increase in cell survival and viability. The induction of RAD51A expression was in response to DNA damage in both cancer and normal cells. However, in normal cells, putrescin pretreatment resulted in alternative splicing of RAD51A and the switch of the predominant expression from the splice variant with the deletion of exon 4 to the full-length variant. Induction of RAD51A alternative splicing by splice-switching oligonucleotides resulted in a decrease in DNA damage and cell protection against cisplatin-induced apoptosis. The results of this study suggest that the cytoprotective activity of polyamines is associated with the alternative splicing of RAD51A pre-mRNA in normal human CD4+ T lymphocytes. The difference in the sensitivity of normal and cancer cells to polyamines may become the basis for the use of these compounds to protect normal lymphocytes during lymphoblastic chemotherapy.

PDE4B Proposed as a High Myopia Susceptibility Gene in Chinese Population

Myopia is the most common cause of refractive error worldwide. High myopia is a severe type of myopia, which usually accompanies pathological changes in the fundus. To identify high myopia susceptibility genes, DNA-pooling based genome-wide association analysis was used to search for a correlation between single nucleotide polymorphisms and high myopia in a Han Chinese cohort (cases vs. controls in discovery stage: 507 vs. 294; replication stage 1: 991 vs. 1,025; replication stage 2: 1,021 vs. 52,708). Three variants (rs10889602T/G, rs2193015T/C, rs9676191A/C) were identified as being significantly associated with high myopia in the discovery, and replication stage. rs10889602T/G is located at the third intron of phosphodiesterase 4B (PDE4B), whose functional assays were performed by comparing the effects of rs10889602T/T deletion of this risk allele on PDE4B and COL1A1 gene and protein expression levels in the rs10889602T/T^del/del, rs10889602T/T^del/wt, and normal control A549 cell lines. The declines in the PDE4B and COL1A1 gene expression levels were larger in the rs10889602T/T deleted A549 cells than in the normal control A549 cells (one-way ANOVA, p < 0.001). The knockdown of PDE4B by siRNA in human scleral fibroblasts led to downregulation of COL1A1. This correspondence between the declines in rs10889602 of the PDE4B gene, PDE4B knockdown, and COL1A1 protein expression levels suggest that PDE4B may be a novel high myopia susceptibility gene, which regulates myopia progression through controlling scleral collagen I expression levels. More studies are needed to determine if there is a correlation between PDE4B and high myopia in other larger sample sized cohorts.

Genetic Control of Splicing at SIRPG Modulates Risk of Type 1 Diabetes

Signal regulatory protein SIRPγ (CD172G) is expressed on the surface of lymphocytes, where it acts by engaging its ligand, CD47. SIRPG, which encodes SIRPγ, contains a nonsynonymous coding variant, rs6043409, which is significantly associated with risk for type 1 diabetes. SIRPG produces multiple transcript isoforms via alternative splicing, all encoding potentially functional proteins. We show that rs6043409 alters a predicted exonic splicing enhancer, resulting in significant shifts in the distribution of SIRPG transcript isoforms. All of these transcript isoforms produced protein upon transient expression in vitro. However, CRISPR/Cas9 targeting of one of the alternatively spliced exons in SIRPG eliminated all SIRPγ expression in Jurkat T cells. These targeted cells formed fewer cell-cell conjugates with each other than with wild-type Jurkat cells, expressed reduced levels of genes associated with CD47 signaling, and had significantly increased levels of cell-surface CD47. In primary CD4+ and CD8+ T cells, cell-surface SIRPγ levels in response to anti-CD3 stimulation varied quantitatively by rs6043409 genotype. Our results suggest that SIRPG is the most likely causative gene for type 1 diabetes risk in the 20p13 region and highlight the role of alternative splicing in lymphocytes in mediating the genetic risk for autoimmunity.

Oligonucleotide correction of an intronic TIMMDC1 variant in cells of patients with severe neurodegenerative disorder

TIMMDC1 encodes the Translocase of Inner Mitochondrial Membrane Domain-Containing protein 1 (TIMMDC1) subunit of complex I of the electron transport chain responsible for ATP production. We studied a consanguineous family with two affected children, now deceased, who presented with failure to thrive in the early postnatal period, poor feeding, hypotonia, peripheral neuropathy and drug-resistant epilepsy. Genome sequencing data revealed a known, deep intronic pathogenic variant TIMMDC1 c.597-1340A>G, also present in gnomAD (~1/5000 frequency), that enhances aberrant splicing. Using RNA and protein analysis we show almost complete loss of TIMMDC1 protein and compromised mitochondrial complex I function. We have designed and applied two different splice-switching antisense oligonucleotides (SSO) to restore normal TIMMDC1 mRNA processing and protein levels in patients' cells. Quantitative proteomics and real-time metabolic analysis of mitochondrial function on patient fibroblasts treated with SSOs showed restoration of complex I subunit abundance and function. SSO-mediated therapy of this inevitably fatal TIMMDC1 neurologic disorder is an attractive possibility.

Open reading frame correction using splice-switching antisense oligonucleotides for the treatment of cystic fibrosis

CFTR gene mutations that result in the introduction of premature termination codons (PTCs) are common in cystic fibrosis (CF). This mutation type causes a severe form of the disease, likely because of low CFTR messenger RNA (mRNA) expression as a result of nonsense-mediated mRNA decay, as well as the production of a nonfunctional, truncated CFTR protein. Current therapeutics for CF, which target residual protein function, are less effective in patients with these types of mutations due in part to low CFTR protein levels. Splice-switching antisense oligonucleotides (ASOs), designed to induce skipping of exons in order to restore the mRNA open reading frame, have shown therapeutic promise preclinically and clinically for a number of diseases. We hypothesized that ASO-mediated skipping of CFTR exon 23 would recover CFTR activity associated with terminating mutations in the exon, including CFTR p.W1282X, the fifth most common mutation in CF. Here, we show that CFTR lacking the amino acids encoding exon 23 is partially functional and responsive to corrector and modulator drugs currently in clinical use. ASO-induced exon 23 skipping rescued CFTR expression and chloride current in primary human bronchial epithelial cells isolated from a homozygote CFTR-W1282X patient. These results support the use of ASOs in treating CF patients with CFTR class I mutations in exon 23 that result in unstable CFTR mRNA and truncations of the CFTR protein.

Exon-skipping antisense oligonucleotides for cystic fibrosis therapy

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), and the CFTR-W1282X nonsense mutation causes a severe form of CF. Although Trikafta and other CFTR-modulation therapies benefit most CF patients, targeted therapy for patients with the W1282X mutation is lacking. The CFTR-W1282X protein has residual activity but is expressed at a very low level due to nonsense-mediated messenger RNA (mRNA) decay (NMD). NMD-suppression therapy and read-through therapy are actively being researched for CFTR nonsense mutants. NMD suppression could increase the mutant CFTR mRNA, and read-through therapies may increase the levels of full-length CFTR protein. However, these approaches have limitations and potential side effects: because the NMD machinery also regulates the expression of many normal mRNAs, broad inhibition of the pathway is not desirable, and read-through drugs are inefficient partly because the mutant mRNA template is subject to NMD. To bypass these issues, we pursued an exon-skipping antisense oligonucleotide (ASO) strategy to achieve gene-specific NMD evasion. A cocktail of two splice-site-targeting ASOs induced the expression of CFTR mRNA without the premature-termination-codon-containing exon 23 (CFTR-Δex23), which is an in-frame exon. Treatment of human bronchial epithelial cells with this cocktail of ASOs that target the splice sites flanking exon 23 results in efficient skipping of exon 23 and an increase in CFTR-Δex23 protein. The splice-switching ASO cocktail increases the CFTR-mediated chloride current in human bronchial epithelial cells. Our results set the stage for developing an allele-specific therapy for CF caused by the W1282X mutation.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

Reclassification of a TMC1 synonymous substitution as a variant disrupting splicing regulatory elements associated with recessive hearing loss

Alterations of the transmembrane channel-like 1 gene (TMC1) are involved in autosomal recessive and dominant nonsyndromic hearing loss (NSHL). To date, up to 117 causal variants including substitutions, insertions and splice variants have been reported in families from different populations. In a patient suffering from severe prelingual NSHL, we identified, in the homozygous state, the previously considered likely benign synonymous c.627C>T; p.(Leu209=) substitution. We used in silico tools predicting variant-induced alterations of splicing regulatory elements (SREs) and pinpointed this transition as a candidate splice-altering variation. Functional splicing analysis, using a minigene assay, confirmed that the variant altered a critical regulatory sequence which is essential for the exon 11 inclusion in the TMC1 transcripts. This result was reinforced by the analysis of orthologous TMC1 mammalian sequences for which the deleterious effect on the mRNA processing of a native thymidine was always counteracted by the presence of a stronger donor splice site or additional enhancer motifs. This study demonstrates, for the first time, the pathogenicity of the c.627C>T alteration leading to its reclassification as a causal variant impacting SREs and highlights the major importance of exhaustive studies to accurately evaluate the pathogenicity of a variant, regardless of the variation type.

OUP accepted manuscript

Alternative splicing is frequently involved in the diversification of protein function and can also be modulated for therapeutic purposes. Here we develop a predictive model, called Exon ByPASS (predicting Exon skipping Based on Protein amino acid SequenceS), to assess the criticality of exon inclusion based solely on information contained in the amino acid sequence upstream and downstream of the exon junctions. By focusing on protein sequence, Exon ByPASS predicts exon skipping independent of tissue and species in the absence of any intronic information. We validate model predictions using transcriptomic and proteomic data and show that the model can capture exon skipping in different tissues and species. Additionally, we reveal potential therapeutic opportunities by predicting synthetically skippable exons and neo-junctions arising in cancer cells.

Silent variant in F8:c.222G>T (p.Thr74Thr) causes a partial exon skipping in a patient with mild hemophilia A

One of the challenges of genetic testing in patients with hemophilia A is the interpretation of sequence variants. Here we report a silent variant found in exon 2 in the F8 gene in a 47‐year‐old patient with a previous von Willebrand disease (VWD) type 1 diagnosis. Clinically he had mild bleeding symptoms restricted to prolonged bleeding from minor wounds. Sanger sequencing of F8 gene using genomic DNA showed a hemizygous silent variant in exon 2: c.222G>T, p.Thr74Thr. When applying ACMG criteria, the variant was predicted to be “likely benign” in the analyzing software or VUS after curating. Sanger sequencing of the patient's cDNA after nested polymerase chain reaction showed that the patient had both a normal transcript containing exons 1–4 and a defect transcript lacking exon 2. These findings explain the patient's low FVIII:C level and led to the diagnosis of mild hemophilia A instead of VWD type 1. This case illustrates that mRNA work‐up may be needed to clarify a patient's phenotype–genotype.

Genetic association of ARID5B with the risk of colorectal cancer within Jammu and Kashmir, India

Colorectal cancer (CRC), which includes the development of cancer from the colon or rectum, is one of the highly prevalent cancers in the populations of Jammu and Kashmir (J&K) in India. However, case-control genetic association studies on CRC are lacking in this population. Various genome-wide association studies have previously shown that single-nucleotide polymorphisms (SNPs) of the AT-rich interaction domain 5B (ARID5B) gene located on chromosome 10q21.2 contribute substantially to the development of colorectal cancer. The association between ARID5B and CRC risk in north Indian population groups is still unknown. To understand the role of ARID5B SNPs in CRC in the population of J&K, we designed a case-control study to investigate the association of the cancer susceptibility variant rs10740055 of ARID5B with CRC in the population of J&K. The study included 180 cases and 390 healthy controls. Genotyping of the rs10740055 variant was performed by RT-PCR using the TaqMan assay technique. Hardy-Weinberg equilibrium of the variant was assessed using the chi-squared test. The allele- and genotype-specific risks were estimated by odds ratios (ORs) with 95% confidence intervals (CIs). The rs10740055 variant showed a higher risk for colorectal cancer with an OR of 3.35 (1.99-5.65 at 95% CI) and P = 0.000005 corrected for age, gender, ethnicity, BMI, alcohol intake and smoking. Our results indicate that the A allele of rs10740055 imparts risk to the population and also that a larger sample size is needed for further statistical validation. The association of other variants in other ARID family genes should also be tested as their role cannot be ruled out.

Associations of MCM8 rs3761873 and rs16991617 variants with abnormal uterine bleeding induced by copper intrauterine device

AIM

Intrauterine device (IUD) is a commonly used contraceptive method worldwide. Abnormal uterine bleeding (AUB) is one of the most common side effects of Cu-IUDs. Since AUB varies among Cu-IUD users, changes in the bleeding-related genetic factors may contribute to AUB. This study aimed to determine the genetic risk factors of AUB after Cu-IUD insertion.

METHODS

We conducted a case-control study on women who experienced AUB after Cu-IUD insertion (case:control = 62:59). Six candidate variants were genotyped using the Sequenom MassARRAY. Genotype and allele frequencies were analyzed using SHEsisPlus. We performed Pearson's Chi-squared test to analyze categorical data, and ESEfinder to predict the impact on splicing regulation.

RESULTS

MCM8 coding sequence variants: rs3761873-A>C was in Exon 7 and rs16991617 A>G was in Exon 12 of all 19 exons, both of which were significantly different between cases and controls (p_allele  = 0.039 and p_genotype  = 0.092). rs6022 and rs6029 in F5 gene and rs3761873 and rs16991617 in the MCM8 gene showed strong linkage disequilibrium (R²  > 0.8). ESEfinder indicated that the variants of MCM8 may affect the splicing regulation.

CONCLUSIONS

MCM8 rs376187 and rs16991617 were associated with AUB in Cu-IUDs users. MCM8 may play a role in AUB by regulating functions of reproductive organs and primary ovarian insufficiency. Our findings may improve the understanding of the genetic basis of AUB caused by Cu-IUDs.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Cas9-specific immune responses compromise local and systemic AAV CRISPR therapy in multiple dystrophic canine models

Adeno-associated virus (AAV)-mediated CRISPR-Cas9 editing holds promise to treat many diseases. The immune response to bacterial-derived Cas9 has been speculated as a hurdle for AAV-CRISPR therapy. However, immunological consequences of AAV-mediated Cas9 expression have thus far not been thoroughly investigated in large mammals. We evaluate Cas9-specific immune responses in canine models of Duchenne muscular dystrophy (DMD) following intramuscular and intravenous AAV-CRISPR therapy. Treatment results initially in robust dystrophin restoration in affected dogs but also induces muscle inflammation, and Cas9-specific humoral and cytotoxic T-lymphocyte (CTL) responses that are not prevented by the muscle-specific promoter and transient prednisolone immune suppression. In normal dogs, AAV-mediated Cas9 expression induces similar, though milder, immune responses. In contrast, other therapeutic (micro-dystrophin and SERCA2a) and reporter (alkaline phosphatase, AP) vectors result in persistent expression without inducing muscle inflammation. Our results suggest Cas9 immunity may represent a critical barrier for AAV-CRISPR therapy in large mammals.

Six splice site variations, three of them novel, in the ABO gene occurring in nine individuals with ABO subtypes

Background

Nucleotide mutations in the ABO gene may reduce the activity of glycosyltransferase, resulting in lower levels of A or B antigen expression in red blood cells. Six known splice sites have been identified according to the database of red cell immunogenetics and the blood group terminology of the International Society of Blood Transfusion. Here, we describe six distinct splice site variants in individuals with ABO subtypes.

Methods

The ABO phenotype was examined using a conventional serological method. A polymerase chain reaction sequence-based typing method was used to examine the whole coding sequence of the ABO gene. The ABO gene haplotypes were studied using allele-specific primer amplification or cloning technology. In silico analytic tools were used to assess the functional effect of splice site variations.

Results

Six distinct variants in the ABO gene splice sites were identified in nine individuals with ABO subtypes, including c.28 + 1_2delGT, c.28 + 5G > A, c.28 + 5G > C, c.155 + 5G > A, c.204-1G > A and c.374 + 5G > A. c.28 + 1_2delGT was detected in an A_w individual, while c.28 + 5G > A, c.28 + 5G > C, and c.204-1G > A were detected in B_el individuals. c.155 + 5G > A was detected in one B₃ and two AB₃ individuals, whereas c.374 + 5G > A was identified in two A_el individuals. Three novel splice site variants (c.28 + 1_2delGT, c.28 + 5G > A and c.28 + 5G > C) in the ABO gene were discovered, all of which resulted in low antigen expression. In silico analysis revealed that all variants had the potential to alter splice transcripts.

Conclusions

Three novel splice site variations in the ABO gene were identified in Chinese individuals, resulting in decreased A or B antigen expression and the formation of ABO subtypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03141-5.

Functional pre-therapeutic evaluation by genome editing of variants of uncertain significance of essential tumor suppressor genes

BACKGROUND

Targeted therapies in oncology are promising but variants of uncertain significance (VUS) limit their use for clinical management and necessitate functional testing in vitro. Using BRCA1 and BRCA2 variants, which have consequences on PARP inhibitor sensitivity, and POLE variants, potential biomarkers of immunotherapy response, we developed a rapid functional assay based on CRISPR-Cas9 genome editing to determine the functional consequences of these variants having potentially direct implications on patients' access to targeted therapies.

METHODS

We first evaluated the functional impact of 26 BRCA1 and 7 BRCA2 variants by editing and comparing NGS results between the variant of interest and a silent control variant. Ten of these variants had already been classified as benign or pathogenic and were used as controls. Finally, we extended this method to the characterization of POLE VUS.

RESULTS

For the 23 variants that were unclassified or for which conflicting interpretations had been reported, 15 were classified as functionally normal and 6 as functionally abnormal. Another two variants were found to have intermediate consequences, both with potential impacts on splicing. We then compared these scores to the patients' responses to PARP inhibitors when possible. Finally, to prove the application of our method to the classification of variants from other tumor suppressor genes, we exemplified with three POLE VUS. Among them, two were classified with an intermediate functional impact and one was functionally abnormal. Eventually, four POLE variants previously classified in databases were also evaluated. However, we found evidence of a discordance with the classification, variant p.Leu424Val being found here functionally normal.

CONCLUSIONS

Our new rapid functional assay can be used to characterize the functional implication of BRCA1 and BRCA2 variants, giving patients whose variants were evaluated as functionally abnormal access to PARP inhibitor treatment. Retrospective analysis of patients' responses to PARP inhibitors, where accessible, was consistent with our functional score evaluation and confirmed the accuracy of our protocol. This method could potentially be extended to the classification of VUS from all essential tumor suppressor genes and can be performed within a timeframe compatible with clinical applications, thereby having a direct theranostic impact.

Identification of novel deep intronic PAH gene variants in patients diagnosed with phenylketonuria

Phenylketonuria (PKU) is caused by phenylalanine hydroxylase (PAH) gene variants. Previously, 94.21% of variants were identified using Sanger sequencing and multiplex ligation-dependent probe amplification. To investigate the remaining variants, we performed whole-genome sequencing for four patients with PKU and unknown genotypes to identify deep intronic or structural variants. We identified three novel heterozygous variants (c.706+368T>C, c.1065+241C>A, and c.1199+502A>T) in a deep PAH gene intron. We detected a c.1199+502A>T variant in 60% (6/10) of PKU patients with genetically undetermined PKU. In silico predictions indicated that the three deep variants may impact splice site selection and result in the inclusion of a pseudo-exon. A c.1199+502A>T PAH minigene and reverse transcription PCR (RT-PCR) on blood RNA from a PKU patient with biallelic variants c.1199+502A>T and c.1199G>A confirmed that the c.1199+502A>T variant may strengthen the predicted branch point and leads to the inclusion of a 25-nt pseudo-exon in the PAH mRNA. Reverse transcription polymerase chain reaction (RT-PCR) on the minigene revealed that c.706+368T>C may create an SRSF2 (SC35) binding site via a 313-nt pseudo-exon, whereas c.1065+241C>A may produce an 81-nt pseudo-exon that strengthens the predicted SRSF1 (SF2/ASF) binding site. These results augment current knowledge of PAH genotypes and show that deep intronic analysis of PAH can genetically diagnose PKU.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.

SLC4A11 mutations causative of congenital hereditary endothelial dystrophy (CHED) progressing to Harboyan syndrome in consanguineous Pakistani families

BACKGROUND

Autosomal recessive corneal hereditary endothelial dystrophy (CHED) is a rare congenital disorder of cornea. Mutations in SLC4A11 gene are associated with CHED phenotype. CHED is also an early feature of Harboyan syndrome. The aim of the present study was to identify genetic mutations in the SLC4A11 gene in CHED cases belonging to inbred Pakistani families. Furthermore, all homozygous mutation carriers were investigated for hearing deficit.

METHODS AND RESULTS

This study included consanguineous CHED families presented at Al-Shifa Trust Eye Hospital, Rawalpindi, Pakistan from June 2018 to September 2018. DNA was extracted from blood samples. Direct sequencing of SLC4A11 gene was performed. All identified variants were evaluated by in silico programs i.e., SIFT, PolyPhen-2, and MutationTaster. Pathogenicity of the two identified splice site variants was analyzed by Human Splicing Finder and MaxEnt Scan. Screening of five CHED families revealed a total of three previously un reported (p.Arg128Gly, c.2241-2A > T and c.1898-2A > C in family CHED19, CHED22 and CHED26 respectively) and two already reported homozygous disease causing variants (p.Arg869Cys and p.Val824Met in family CHED24 and CHED25 respectively) as predicted by mutation taster. All of these variants segregated with disease phenotype and were not detected in controls.

CONCLUSION

Affected individuals of the five CHED families screened in this study had the disease due to SLC4A11 mutations and progressing to Harboyan syndrome. Identification of previously unreported mutations aid to heterogeneity of SLC4A11 and CHED pathogenesis as well as helped to provide genetic counseling to affected families.

A novel mouse model of PMS2 founder mutation that causes mismatch repair defect due to aberrant splicing

Hereditary non-polyposis colorectal cancer, now known as Lynch syndrome (LS) is one of the most common cancer predisposition syndromes and is caused by germline pathogenic variants (GPVs) in DNA mismatch repair (MMR) genes. A common founder GPV in PMS2 in the Canadian Inuit population, NM_000535.5: c.2002A>G, leads to a benign missense (p.I668V) but also acts as a de novo splice site that creates a 5 bp deletion resulting in a truncated protein (p.I668*). Individuals homozygous for this GPV are predisposed to atypical constitutional MMR deficiency with a delayed onset of first primary malignancy. We have generated mice with an equivalent germline mutation (Pms2c.1993A>G) and demonstrate that it results in a splicing defect similar to those observed in humans. Homozygous mutant mice are viable like the Pms2 null mice. However, unlike the Pms2 null mice, these mutant mice are fertile, like humans homozygous for this variant. Furthermore, these mice exhibit a significant increase in microsatellite instability and intestinal adenomas on an Apc mutant background. Rectification of the splicing defect in human and murine fibroblasts using antisense morpholinos suggests that this novel mouse model can be valuable in evaluating the efficacy aimed at targeting the splicing defect in PMS2 that is highly prevalent among the Canadian Inuits.

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.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

A multistep approach to the genotype-phenotype analysis of Polish patients with tuberous sclerosis complex

The aim of this study was to evaluate a cost-effective diagnostic strategy for identification of casual variants for tuberous sclerosis complex (TSC) in the Polish population and to correlate the genetic results with selected phenotypic features. Fifty-five patients, aged 3-44 years, with a clinical diagnosis of TSC were enrolled into the study. All patients received a three-step analysis: next generation sequencing screening (NGS), multiplex ligation-dependent probe amplification (MLPA) and deep sequencing. This multistep approach obtained positive results in 51/55 (93%) patients: of the 51 positives TSC1 variants were observed in 16 (31%) and TSC2 variants in 35 (69%); these included 13 novel variants and two patients with mosaicism. Four patients (7%) had no mutation identified (NMI). Among the TSC1 gene variants, there were five nonsense, four frameshift, three large deletions, two missense and two splicing variants. For the TSC2 gene, 11 were missense, eight splicing, six frameshift, four large deletions, two in-frame deletions and four nonsense variants. The patients with TSC2 changes had their clinical diagnosis of TSC at a younger age than those with TSC1 changes (one year vs three years, p = 0.041). The TSC2 group demonstrated a higher number of major symptoms per patient (p = 0.04). Subependymal giant cell astrocytoma with concomitance of other brain lesions was more common in patients with missense mutations in either gene (23% vs 0%, p = 0.02). Such a multistep molecular diagnostic strategy could increase the possibility of detecting causal variants for TSC and may allow detection of mosaicism at low levels. Missense pathogenic variants in TSC1 or TSC2 gene might be associated with a higher risk of brain lesions.

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.

Therapeutic manipulation of IKBKAP mis-splicing with a small molecule to cure familial dysautonomia

Approximately half of genetic disease-associated mutations cause aberrant splicing. However, a widely applicable therapeutic strategy to splicing diseases is yet to be developed. Here, we analyze the mechanism whereby IKBKAP-familial dysautonomia (FD) exon 20 inclusion is specifically promoted by a small molecule splice modulator, RECTAS, even though IKBKAP-FD exon 20 has a suboptimal 5' splice site due to the IVS20 + 6 T > C mutation. Knockdown experiments reveal that exon 20 inclusion is suppressed in the absence of serine/arginine-rich splicing factor 6 (SRSF6) binding to an intronic splicing enhancer in intron 20. We show that RECTAS directly interacts with CDC-like kinases (CLKs) and enhances SRSF6 phosphorylation. Consistently, exon 20 splicing is bidirectionally manipulated by targeting cellular CLK activity with RECTAS versus CLK inhibitors. The therapeutic potential of RECTAS is validated in multiple FD disease models. Our study indicates that small synthetic molecules affecting phosphorylation state of SRSFs is available as a new therapeutic modality for mechanism-oriented precision medicine of splicing diseases.

IntSplice2: Prediction of the Splicing Effects of Intronic Single-Nucleotide Variants Using LightGBM Modeling

Prediction of the effect of a single-nucleotide variant (SNV) in an intronic region on aberrant pre-mRNA splicing is challenging except for an SNV affecting the canonical GU/AG splice sites (ss). To predict pathogenicity of SNVs at intronic positions −50 (Int-50) to −3 (Int-3) close to the 3’ ss, we developed light gradient boosting machine (LightGBM)-based IntSplice2 models using pathogenic SNVs in the human gene mutation database (HGMD) and ClinVar and common SNVs in dbSNP with 0.01 ≤ minor allelic frequency (MAF) < 0.50. The LightGBM models were generated using features representing splicing cis-elements. The average recall/sensitivity and specificity of IntSplice2 by fivefold cross-validation (CV) of the training dataset were 0.764 and 0.884, respectively. The recall/sensitivity of IntSplice2 was lower than the average recall/sensitivity of 0.800 of IntSplice that we previously made with support vector machine (SVM) modeling for the same intronic positions. In contrast, the specificity of IntSplice2 was higher than the average specificity of 0.849 of IntSplice. For benchmarking (BM) of IntSplice2 with IntSplice, we made a test dataset that was not used to train IntSplice. After excluding the test dataset from the training dataset, we generated IntSplice2-BM and compared it with IntSplice using the test dataset. IntSplice2-BM was superior to IntSplice in all of the seven statistical measures of accuracy, precision, recall/sensitivity, specificity, F1 score, negative predictive value (NPV), and matthews correlation coefficient (MCC). We made the IntSplice2 web service at https://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice2.

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.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Antisense oligonucleotide-based drug development for Cystic Fibrosis patients carrying the 3849+10 kb C-to-T splicing mutation

Background:

Antisense oligonucleotide (ASO)-based drugs for splicing modulation were recently approved for various genetic diseases with unmet need. Here we aimed to develop an ASO-based splicing modulation therapy for Cystic Fibrosis (CF) patients carrying the 3849 + 10 kb C-to-T splicing mutation in the CFTR gene.

Methods:

We have screened, in FRT cells expressing the 3849 + 10 kb C-to-T splicing mutation, ~30 2ʹ-O-Methyl-modified phosphorothioate ASOs, targeted to prevent the recognition and inclusion of a cryptic exon generated due to the mutation. The effect of highly potent ASO candidates on the splicing pattern, protein maturation and CFTR function was further analyzed in well differentiated primary human nasal and bronchial epithelial cells, derived from patients carrying at least one 3849 + 10 kb C-to-T allele.

Results:

A highly potent lead ASO, efficiently delivered by free uptake, was able to significantly increase the level of correctly spliced mRNA and completely restore the CFTR function to wild type levels in cells from a homozygote patient. This ASO led to CFTR function with an average of 43% of wild type levels in cells from various heterozygote patients. Optimized efficiency of the lead ASO was further obtained with 2ʹ-Methoxy Ethyl modification (2ʹMOE).

Conclusion:

The highly efficient splicing modulation and functional correction, achieved by free uptake of the selected lead ASO in various patients, demonstrate the ASO therapeutic potential benefit for CF patients carrying splicing mutations and is aimed to serve as the basis for our current clinical development.

A world-wide analysis of reduced sensitivity to DMI fungicides in the banana pathogen Pseudocercospora fijiensis

BACKGROUND

Pseudocercospora fijiensis is the causal agent of the black leaf streak disease (BLSD) of banana. Bananas are important global export commodities and a major staple food. Their susceptibility to BLSD pushes disease management towards excessive fungicide use, largely relying on multisite inhibitors and sterol demethylation inhibitors (DMIs). These fungicides are ubiquitous in plant disease control, targeting the CYP51 enzyme. We examined sensitivity to DMIs in P. fijiensis field isolates collected from various major banana production zones in Colombia, Costa Rica, Dominican Republic, Ecuador, the Philippines, Guadalupe, Martinique and Cameroon and determined the underlying genetic reasons for the observed phenotypes.

RESULTS

We observed a continuous range of sensitivity towards the DMI fungicides difenoconazole, epoxiconazole and propiconazole with clear cross-sensitivity. Sequence analyses of PfCYP51 in 266 isolates showed 28 independent amino acid substitutions, nine of which correlated with reduced sensitivity to DMIs. In addition to the mutations, we observed up to six insertions in the Pfcyp51 promoter. Such promoter insertions contain repeated elements with a palindromic core and correlate with the enhanced expression of Pfcyp51 and hence with reduced DMI sensitivity. Wild-type isolates from unsprayed bananas fields did not contain any promoter insertions.

CONCLUSION

The presented data significantly contribute to understanding of the evolution and global distribution of DMI resistance mechanisms in P. fijiensis field populations and facilitate the prediction of different DMI efficacy. The overall reduced DMI sensitivity calls for the deployment of a wider range of solutions for sustainable control of this major banana disease. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.