ESEfinder: A web resource to identify exonic splicing enhancers

A Translational Approach to Spinal Neurofibromatosis: Clinical and Molecular Insights from a Wide Italian Cohort

Spinal neurofibromatosis (SNF), a phenotypic subclass of neurofibromatosis 1 (NF1), is characterized by bilateral neurofibromas involving all spinal roots. In order to deepen the understanding of SNF’s clinical and genetic features, we identified 81 patients with SNF, 55 from unrelated families, and 26 belonging to 19 families with at least 1 member affected by SNF, and 106 NF1 patients aged >30 years without spinal tumors. A comprehensive NF1 mutation screening was performed using NGS panels, including NF1 and several RAS pathway genes. The main features of the SNF subjects were a higher number of internal neurofibromas (p < 0.001), nerve root swelling (p < 0.001), and subcutaneous neurofibromas (p = 0.03), while hyperpigmentation signs were significantly less frequent compared with the classical NF1-affected cohorts (p = 0.012). Fifteen patients underwent neurosurgical intervention. The histological findings revealed neurofibromas in 13 patients and ganglioneuromas in 2 patients. Phenotypic variability within SNF families was observed. The proportion of missense mutations was higher in the SNF cases than in the classical NF1 group (21.40% vs. 7.5%, p = 0.007), conferring an odds ratio (OR) of 3.34 (CI = 1.33−10.78). Two unrelated familial SNF cases harbored in trans double NF1 mutations that seemed to have a subclinical worsening effect on the clinical phenotype. Our study, with the largest series of SNF patients reported to date, better defines the clinical and genetic features of SNF, which could improve the management and genetic counseling of NF1.

Novel biallelic mutations in TMEM126B cause splicing defects and lead to Leigh-like syndrome with severe complex I deficiency

Leigh syndrome (LS)/Leigh-like syndrome (LLS) is one of the most common mitochondrial disease subtypes, caused by mutations in either the nuclear or mitochondrial genomes. Here, we identified a novel intronic mutation (c.82-2 A > G) and a novel exonic insertion mutation (c.290dupT) in TMEM126B from a Chinese patient with clinical manifestations of LLS. In silico predictions, minigene splicing assays and patients’ RNA analyses determined that the c.82-2 A > G mutation resulted in complete exon 2 skipping, and the c.290dupT mutation provoked partial and complete exon 3 skipping, leading to translational frameshifts and premature termination. Functional analysis revealed the impaired mitochondrial function in patient-derived lymphocytes due to severe complex I content and assembly defect. Altogether, this is the first report of LLS in a patient carrying mutations in TMEM126B. Our data uncovers the functional effect and the molecular mechanism of the pathogenic variants c.82-2 A > G and c.290dupT, which expands the gene mutation spectrum of LLS and clinical spectrum caused by TMEM126B mutations, and thus help to clinical diagnosis of TMEM126B mutation‐related mitochondrial diseases.

Correction of the Splicing Defect Caused by a Recurrent Variant in ABCA4 (c.769-784C>T) That Underlies Stargardt Disease

Stargardt disease is an inherited retinal disease caused by biallelic mutations in the ABCA4 gene, many of which affect ABCA4 splicing. In this study, nine antisense oligonucleotides (AONs) were designed to correct pseudoexon (PE) inclusion caused by a recurrent deep-intronic variant in ABCA4 (c.769-784C>T). First, the ability of AONs to skip the PE from the final ABCA4 mRNA transcript was assessed in two cellular models carrying the c.769-784C>T variant: a midigene assay using HEK293T cells and patient-derived fibroblasts. Based on the splicing-correcting ability of each individual AON, the three most efficacious AONs targeting independent regions of the PE were selected for a final assessment in photoreceptor precursor cells (PPCs). The final analysis in the PPC model confirmed high efficacy of AON2, -5, and -7 in promoting PE exclusion. Among the three AONs, AON2 is chosen as the lead candidate for further optimization, hereby showcasing the high potential of AONs to correct aberrant splicing events driven by deep-intronic variants.

Using single molecule Molecular Inversion Probes as a cost-effective, high-throughput sequencing approach to target all genes and loci associated with macular diseases

Macular degenerations (MDs) are a subgroup of retinal disorders characterized by central vision loss. Knowledge is still lacking on the extent of genetic and nongenetic factors influencing inherited MD (iMD) and age-related MD (AMD) expression. Single molecule Molecular Inversion Probes (smMIPs) have proven effective in sequencing the ABCA4 gene in patients with Stargardt disease to identify associated coding and noncoding variation, however many MD patients still remain genetically unexplained. We hypothesized that the missing heritability of MDs may be revealed by smMIPs-based sequencing of all MD-associated genes and risk factors. Using 17,394 smMIPs, we sequenced the coding regions of 105 iMD and AMD-associated genes and noncoding or regulatory loci, known pseudo-exons, and the mitochondrial genome in two test cohorts that were previously screened for variants in ABCA4. Following detailed sequencing analysis of 110 probands, a diagnostic yield of 38% was observed. This established an ''MD-smMIPs panel," enabling a genotype-first approach in a high-throughput and cost-effective manner, whilst achieving uniform and high coverage across targets. Further analysis will identify known and novel variants in MD-associated genes to offer an accurate clinical diagnosis to patients. Furthermore, this will reveal new genetic associations for MD and potential genetic overlaps between iMD and AMD.

SRSF6 balances mitochondrial-driven innate immune outcomes through alternative splicing of BAX

To mount a protective response to infection while preventing hyperinflammation, gene expression in innate immune cells must be tightly regulated. Despite the importance of pre-mRNA splicing in shaping the proteome, its role in balancing immune outcomes remains understudied. Transcriptomic analysis of murine macrophage cell lines identified Serine/Arginine Rich Splicing factor 6 (SRSF6) as a gatekeeper of mitochondrial homeostasis. SRSF6-dependent orchestration of mitochondrial health is directed in large part by alternative splicing of the pro-apoptosis pore-forming protein BAX. Loss of SRSF6 promotes accumulation of BAX-κ, a variant that sensitizes macrophages to undergo cell death and triggers upregulation of interferon stimulated genes through cGAS sensing of cytosolic mitochondrial DNA. Upon pathogen sensing, macrophages regulate SRSF6 expression to control the liberation of immunogenic mtDNA and adjust the threshold for entry into programmed cell death. This work defines BAX alternative splicing by SRSF6 as a critical node not only in mitochondrial homeostasis but also in the macrophage's response to pathogens.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

Background

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

Results

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in “cryptic” splicing regions beyond canonical donor or acceptor splice sites.

Conclusion

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

A Mild Case of Autosomal Recessive Osteopetrosis Masquerading as the Dominant Form Involving Homozygous Deep Intronic Variations in the CLCN7 Gene

Osteopetrosis is a heterogeneous group of rare hereditary diseases characterized by increased bone mass of poor quality. Autosomal-dominant osteopetrosis type II (ADOII) is most often caused by mutation of the CLCN7 gene leading to impaired bone resorption. Autosomal recessive osteopetrosis (ARO) is a more severe form and is frequently accompanied by additional morbidities. We report an adult male presenting with classical clinical and radiological features of ADOII. Genetic analyses showed no amino-acid-converting mutation in CLCN7 but an apparent haploinsufficiency and suppression of CLCN7 mRNA levels in peripheral blood mononuclear cells. Next generation sequencing revealed low-frequency intronic homozygous variations in CLCN7, suggesting recessive inheritance. In silico analysis of an intronic duplication c.595-120_595-86dup revealed additional binding sites for Serine- and Arginine-rich Splicing Factors (SRSF), which is predicted to impair CLCN7 expression. Quantitative backscattered electron imaging and histomorphometric analyses revealed bone tissue and material abnormalities. Giant osteoclasts were present and additionally to lamellar bone, and abundant woven bone and mineralized cartilage were observed, together with increased frequency and thickness of cement lines. Bone mineralization density distribution (BMDD) analysis revealed markedly increased average mineral content of the dense bone (CaMean T-score + 10.1) and frequency of bone with highest mineral content (CaHigh T-score + 19.6), suggesting continued mineral accumulation and lack of bone remodelling. Osteocyte lacunae sections (OLS) characteristics were unremarkable except for an unusually circular shape. Together, our findings suggest that the reduced expression of CLCN7 mRNA in osteoclasts, and possibly also osteocytes, causes poorly remodelled bone with abnormal bone matrix with high mineral content. This together with the lack of adequate bone repair mechanisms makes the material brittle and prone to fracture. While the skeletal phenotype and medical history were suggestive of ADOII, genetic analysis revealed that this is a possible mild case of ARO due to deep intronic mutation.

HvGBSSI mutation at the splicing receptor site affected RNA splicing and decreased amylose content in barley

Granule-bound starch synthase I (HvGBSSI) is encoded by the barley waxy (Wx-1) gene and is the sole enzyme in the synthesis of amylose. Here, a Wx-1 mutant was identified from an ethyl methane sulfonate (EMS)-mutagenized barley population. There were two single-base mutations G1086A and A2424G in Wx-1 in the mutant (M2-1105). The G1086A mutation is located at the 3' splicing receptor (AG) site of the fourth intron, resulting in an abnormal RNA splicing. The A2424G mutation was a synonymous mutation in the ninth intron. The pre-mRNA of Wx-1 was incorrectly spliced and transcribed into two abnormal transcripts. The type I transcript had a 6 bp deletion in the 5' of fifth exon, leading to a translated HvGBSSI protein lacking two amino acids with a decreased starch-binding capacity. In the type II transcript, the fourth intron was incorrectly cleaved and retained, resulting in the premature termination of the barley Wx-1 gene. The mutations in the Wx-1 decreased the enzymatic activity of the HvGBSSI enzyme and resulted in a decreased level in amylose content. This work sheds light on a new Wx-1 gene inaction mechanism.

lncRNA260 siRNA Accelerates M2 Macrophage Polarization and Alleviates Oxidative Stress via Inhibiting IL28RA Gene Alternative Splicing

The macrophage transformation of inflammatory M1 to anti-inflammatory M2 could be promoted by activating PI3K/AKT signaling pathway. In our previous study, it was found that downregulation of lncRNA260 could ameliorate hypoxic cardiomyocyte injury by regulating IL28RA through the activation of PI3K/AKT signaling pathways. It was suggested that lncRNA260 siRNA could promote the macrophages toward M2 polarization by regulating IL28RA. In this study, lncRNA260 siRNA was used to observe its effect on the polarization of murine bone marrow-derived macrophages (BMDM) and investigate its related mechanisms. lncRNA 260 specific siRNA were designed and synthesized which were transfected into murine BMDM with liposomes. The experiment was divided into three groups: Hypoxia group, Hypoxia+lncRNA 260-specific siRNA transfection group, and Normoxia group. The CD206-APC/CD11b-FITC or CD206-FITC/CD107b (Mac-3) double positive proportions were used to compare the M2 polarization proportions in the hypoxia process by using the immunofluorescence staining method. The p-AKT, Arg 1, PI3KCG, IL28RAV1, and IL28RAV2 protein expression changes were observed by using the western blot method. Compared with the Normoxia group, the M2 proportions were significantly decreased in the Hypoxia group (P < 0.05). Compared with the hypoxia group, the M2 proportions were significantly increased in the Hypoxia+lncRNA260 siRNA transfection group (P < 0.05). In the Hypoxia group, the ratios of Arg 1/β-Actin, p-AKT/β-Actin, PI3KCG/β-Actin, and IL28RAV1/β-Actin were significantly lower than those in the Normoxia group (P < 0.05). After transfection with lncRNA260 siRNA, the ratios of Arg1/β-Actin, p-AKT/β-Actin, PI3KCG/β-Actin, and IL28RAV1/β-Actin were significantly higher than those in the Hypoxia group (P < 0.05). Compared with the Normoxia group, the IL28RAV2/β-Actin in the Hypoxia group was significantly increased (P < 0.05). After transfection with lncRNA260 siRNA, the ratio of IL28RAV2/β-Actin was significantly decreased than that in the Hypoxia group (P < 0.05). lncRNA260 siRNA could promote the M2 polarization of the hypoxia macrophages by reducing the IL28RAV2 alternative splicing variant, which might be related to the activation of the JAK-STAT and PI3K/AKT signaling pathways. It will provide a new strategy for the anti-inflammation, antioxidative stress therapy, and cardiac remodeling after AMI.

Paediatric BCOR ‐associated sarcomas with a novel long spliced internal tandem duplication of BCOR exon 15

Clear cell sarcoma of the kidney (CCSK) and primitive myxoid mesenchymal tumour of infancy (PMMTI) are paediatric sarcomas that most commonly harbour internal tandem duplications (ITDs) of exon 15 of the BCOR gene, in the range of 87–114 base pairs (bp). Some cases, instead, have BCOR‐CCNB3 or YWHAE‐NUTM2 gene fusions. About 10% of cases lack any of these genetic alterations when tested by standard methods. Two cases of CCSK and one PMMTI lacking the aforementioned mutations were analysed using Archer FusionPlex technology. Two related BCOR exon 15 RNA transcripts with ITDs of lengths 388 and 96 bp were detected in each case; only the 388 bp transcript was identified when genomic DNA was sequenced. In silico analysis of this transcript revealed acceptor and donor splice sites indicating that, at the RNA level, the 388‐bp transcript was likely spliced to form the 96‐bp transcript. The results were confirmed by Sanger sequencing using primers targeting the ITD breakpoint. This novel and unusually long ITD segment is difficult to identify by DNA sequencing using typical primer design strategies flanking entire duplicated segments because it exceeds the typical read lengths of most sequencing platforms as well as the usual fragment lengths obtained from formalin‐fixed paraffin‐embedded material. As diagnosis of CCSK and PMMTI may be challenging by morphology and immunohistochemistry alone, it is important to identify mutations in these cases. Knowledge of this novel BCOR ITD is important in relation to primer design for detection by sequencing, and using RNA versus DNA for sequencing.

Allele-specific alternative splicing of Drosophila Ribosomal protein S21 suppresses a lethal mutation in the Phosphorylated adaptor for RNA export (Phax) gene

Genetic disruptions to the biogenesis of spliceosomal small-nuclear ribonucleoproteins in Drosophila cause wide-spread alternative splicing changes, including changes to the splicing of pre-mRNA for Ribosomal protein S21 (RpS21). Using a transposon mutant for the Phosphorylated adaptor for RNA export (Phax) gene, we demonstrate that changes in the splicing of RpS21 transcripts have a strong influence on the developmental progression of Phax^SH/SH mutants. Different alleles of the Drosophila RpS21 gene are circulating in common laboratory strains and cell lines. These alleles exhibit differences in RpS21 intron retention and splicing efficiency. Differences in the splicing of RpS21 transcripts account for prior conflicting observations of the phenotypic severity of Phax^SH/SH mutant stocks. The alleles uncover a strong splicing enhancer in RpS21 transcripts that can fully suppress the larval lethality and partially suppress the pupal lethality exhibited by Phax^SH/SH mutant lines. In the absence of the splicing enhancer, the splicing of RpS21 transcripts can be modulated in trans by the SR-rich B52 splicing factor. As Phax^SH/SH mutants exhibit wide-spread splicing changes in transcripts for other genes, findings here establish the importance of a single alternative splicing event, RpS21 splicing or intron retention, to the developmental progression of Drosophila.

Performance evaluation of differential splicing analysis methods and splicing analytics platform construction

A proportion of previously defined benign variants or variants of uncertain significance in humans, which are challenging to identify, may induce an abnormal splicing process. An increasing number of methods have been developed to predict splicing variants, but their performance has not been completely evaluated using independent benchmarks. Here, we manually sourced ∼50 000 positive/negative splicing variants from > 8000 studies and selected the independent splicing variants to evaluate the performance of prediction methods. These methods showed different performances in recognizing splicing variants in donor and acceptor regions, reminiscent of different weight coefficient applications to predict novel splicing variants. Of these methods, 66.67% exhibited higher specificities than sensitivities, suggesting that more moderate cut-off values are necessary to distinguish splicing variants. Moreover, the high correlation and consistent prediction ratio validated the feasibility of integration of the splicing prediction method in identifying splicing variants. We developed a splicing analytics platform called SPCards, which curates splicing variants from publications and predicts splicing scores of variants in genomes. SPCards also offers variant-level and gene-level annotation information, including allele frequency, non-synonymous prediction and comprehensive functional information. SPCards is suitable for high-throughput genetic identification of splicing variants, particularly those located in non-canonical splicing regions.

History of development of the life-saving drug “Nusinersen” in spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive disorder with an incidence of 1/6,000–1/10,000 and is the leading fatal disease among infants. Previously, there was no effective treatment for SMA. The first effective drug, nusinersen, was approved by the US FDA in December 2016, providing hope to SMA patients worldwide. The drug was introduced in the European Union in 2017 and China in 2019 and has so far saved the lives of several patients in most parts of the world. Nusinersen are fixed sequence antisense oligonucleotides with special chemical modifications. The development of nusinersen progressed through major scientific discoveries in medicine, genetics, biology, and other disciplines, wherein several scientists have made substantial contributions. In this article, we will briefly describe the pathogenesis and therapeutic strategies of SMA, summarize the timeline of important scientific findings during the development of nusinersen in a detailed, scientific, and objective manner, and finally discuss the implications of the development of nusinersen for SMA research.

Modeling splicing outcome by combining 5'ss strength and splicing regulatory elements

Correct pre-mRNA processing in higher eukaryotes vastly depends on splice site recognition. Beyond conserved 5'ss and 3'ss motifs, splicing regulatory elements (SREs) play a pivotal role in this recognition process. Here, we present in silico designed sequences with arbitrary a priori prescribed splicing regulatory HEXplorer properties that can be concatenated to arbitrary length without changing their regulatory properties. We experimentally validated in silico predictions in a massively parallel splicing reporter assay on more than 3000 sequences and exemplarily identified some SRE binding proteins. Aiming at a unified 'functional splice site strength' encompassing both U1 snRNA complementarity and impact from neighboring SREs, we developed a novel RNA-seq based 5'ss usage landscape, mapping the competition of pairs of high confidence 5'ss and neighboring exonic GT sites along HBond and HEXplorer score coordinate axes on human fibroblast and endothelium transcriptome datasets. These RNA-seq data served as basis for a logistic 5'ss usage prediction model, which greatly improved discrimination between strong but unused exonic GT sites and annotated highly used 5'ss. Our 5'ss usage landscape offers a unified view on 5'ss and SRE neighborhood impact on splice site recognition, and may contribute to improved mutation assessment in human genetics.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Alternatively spliced CSF3R isoforms in SRSF2 P95H mutated myeloid neoplasms

Alternatively spliced colony stimulating factor 3 receptor (CSF3R) isoforms Class III and Class IV are observed in myelodysplastic syndromes (MDS), but their roles in disease remain unclear. We report that the MDS-associated splicing factor SRSF2 affects the expression of Class III and Class IV isoforms and perturbs granulopoiesis. Add-back of the Class IV isoform in Csf3r-null mouse progenitor cells increased granulocyte progenitors with impaired neutrophil differentiation, while add-back of the Class III produced dysmorphic neutrophils in fewer numbers. These CSF3R isoforms were elevated in patients with myeloid neoplasms harboring SRSF2 mutations. Using in vitro splicing assays, we confirmed increased Class III and Class IV transcripts when SRSF2 P95 mutations were co-expressed with the CSF3R minigene in K562 cells. Since SRSF2 regulates splicing partly by recognizing exonic splicing enhancer (ESE) sequences on pre-mRNA, deletion of either ESE motifs within CSF3R exon 17 decreased Class IV transcript levels without affecting Class III. CD34+ cells expressing SRSF2 P95H showed impaired neutrophil differentiation in response to G-CSF and was accompanied by increased levels of Class IV. Our findings suggest that SRSF2 P95H promotes Class IV splicing by binding to key ESE sequences in CSF3R exon 17, and that SRSF2, when mutated, contributes to dysgranulopoiesis.

Protocol to identify host-viral protein interactions between coagulation-related proteins and their genetic variants with SARS-CoV-2 proteins

Here, we describe a bioinformatics pipeline that evaluates the interactions between coagulation-related proteins and genetic variants with SARS-CoV-2 proteins. This pipeline searches for host proteins that may bind to viral protein and identifies and scores the protein genetic variants to predict the disease pathogenesis in specific subpopulations. Additionally, it is able to find structurally similar motifs and identify potential binding sites within the host-viral protein complexes to unveil viral impact on regulated biological processes and/or host-protein impact on viral invasion or reproduction. For complete details on the use and execution of this protocol, please refer to Holcomb et al. (2021).

Prevalent MLC1 mutation causing autosomal recessive megalencephalic leukoencephalopathy in consanguineous Palestinian families

BACKGROUND

Recessive forms of megalencephalic leukoencephalopathy with subcortical cysts (MLC, OMIM 604004) is a rare early-onset leukodystrophy that presents with macrocephaly, seizures, slowly progressive gross motor deterioration, and MRI evidence of diffuse symmetric white matter swelling and subcortical cysts in the anterior temporal and frontoparietal regions. Later in the disease course, significant spasticity and ataxia develop, which may be accompanied by intellectual deterioration. This disease is caused mostly by biallelic pathogenic variants in the MLC1 gene.

METHODS

In this study, we analysed the clinical and molecular architecture of 6 individuals, belonging to 4 unrelated consanguineous Palestinian families, presenting with consistent MLC features. We sequenced the entire coding and flanking intronic regions of the MLC1 gene.

RESULTS

In all recruited individuals, we detected one recurrent homozygous splice donor mutation NM_015166.4: c.423 + 1G > A. All parents were heterozygous carriers. The mutation abolishes a highly conserved splice site in humans and other species. In silico splice predictors suggested the loss of a canonical splice donor site (CADD score 33.0. SpliceAI: 0.980). The c.423 + 1G > A variant is rare; it was detected in only 4 heterozygous carriers in gnomAD.

CONCLUSION

In this study, we identified a recurrent MLC1 variant (c.423 + 1G > A) as the cause of MLC among a group of Palestinian patients originating from a particular region of the country. Cost-effective studies should be performed to evaluate the implementation of carrier screening in adults originating from this region. Our findings have the potential to contribute to improved genetic diagnosis and carrier testing for individuals within this population and the wider community.

Recognition of functional genetic polymorphism using ESE motif definition: a conservative evolutionary approach to CYP2D6/CYP2C19 gene variants

Although predicting the effects of variants near intron-exon boundaries is relatively straightforward, predicting the functional Exon Splicing Enhancers (ESEs) and the possible effects of variants within ESEs remains a challenge. Considering the essential role of CYP2D6/CYP2C19 genes in drug metabolism, we attempted to identify variants that are most likely to disrupt splicing through their effect on these ESEs. ESEs were predicted in these two genes using ESEfinder 3.0, incorporating a series of filters (increased threshold and evolutionary conservation). Finally, reported mutations were evaluated for their potential to disrupt splicing by affecting these ESEs. Initially, 169 and 243 ESEs were predicted for CYP2C19/CYP2D6, respectively. However, applying the filters, the number of predicted ESEs was reduced to 26 and 19 in CYP2C19/CYP2D6, respectively. Comparing prioritized predicted ESEs with known sequence variants in CYP2C19/CYP2D6 genes highlights 18 variations within conserved ESEs for each gene. We found good agreement in cases where such predictions could be compared to experimental evidence. In total, we prioritized a subset of mutational changes in CYP2C19/CYP2D6 genes that may affect the function of these genes and lead to altered drug responses. Clinical studies and functional analysis for investigating detailed functional consequences of the mentioned mutations and their phenotypic outcomes is mostly recommended.

A novel missense variant in endothelin-2 (EDN2) causes a growth and respiratory lethal syndrome in bovine

The high level of fragmentation of the Spanish Lidia cattle breed, divided into lineages called ‘castas’ and into herds within lineages based on reproductive isolation, increases the risk of homozygosity and the outbreak of recessive genetic defects. Since 2004, an increasing number of calves have been identified in a Lidia herd with signs of severe growth retardation, respiratory alterations and juvenile lethality, which constitutes a novel inherited syndrome in cattle and was subsequently termed growth and respiratory lethal syndrome. We performed a genome‐wide association study on a cohort of 13 affected calves and 24 putative non‐carrier parents, mapping the disease to a wide 6 cM region on bovine chromosome 3 (p < 10⁻⁷). Whole genome re‐sequencing of three affected calves and three putative non‐carrier parents identified a novel missense variant (c.149G>A|p.Cys50Tyr) in exon 2 of the endothelin 2 (EDN2) gene. Bioinformatic analyses of p.Cys50Tyr effects predicted them to be damaging for both the structure and the function of the edn2 protein, and to create a new site of splicing that may also affect the pattern of pre‐mRNA splicing and exon definition. Sanger sequencing of this variant on the rest of the sample set confirmed the segregation pattern obtained with whole genome re‐sequencing. The identification of the causative variant and the development of a diagnostic genetic test enable the efficient design of matings to keep the effective population size as high as possible, as well as providing insights into the first EDN2‐associated hereditary disease in cattle or other species.

Antisense Oligonucleotide Rescue of Deep-Intronic Variants Activating Pseudoexons in the 6-Pyruvoyl-Tetrahydropterin Synthase Gene

We report two new 6-pyruvoyl-tetrahydropterin synthase splicing variants identified through genomic sequencing and transcript analysis in a patient with tetrahydrobiopterin deficiency, presenting with hyperphenylalaninemia and monoamine neurotransmitter deficiency. Variant c.243 + 3A>G causes exon 4 skipping. The deep-intronic c.164-672C>T variant creates a potential 5' splice site that leads to the inclusion of four overlapping pseudoexons, corresponding to exonizations of an antisense short interspersed nuclear element AluSq repeat sequence. Two of the identified pseudoexons have been reported previously, activated by different deep-intronic variants, and were also detected at residual levels in control cells. Interestingly, the predominant pseudoexon is nearly identical to a disease causing activated pseudoexon in the F8 gene, with the same 3' and 5' splice sites. Splice switching antisense oligonucleotides (SSOs) were designed to hybridize with splice sites and/or predicted binding sites for regulatory splice factors. Different SSOs corrected the aberrant pseudoexon inclusion, both in minigenes and in fibroblasts from patients carrying the new variant c.164-672C>T or the previously described c.164-716A>T. With SSO treatment PTPS protein was recovered, illustrating the therapeutic potential of the approach, for patients with different pseudoexon activating variants in the region. In addition, the natural presence of pseudoexons in the wild type context suggests the possibility of applying the antisense strategy in patients with hypomorphic PTS variants with the purpose of upregulating their expression to increase overall protein and activity.

Expanding the NGLY1 deficiency phenotype: Case report of an atypical patient

NGLY1 deficiency is a rare congenital disorder of deglycosylation with a unique constellation of symptoms that include hypo- or alacrima, movement disorder, epilepsy, and severe intellectual disability (OMIM #615273). Here we report a patient with NGLY1 deficiency whose clinical presentation lacks many of the features associated with the disease and has a much milder intellectual disability than had been previously reported, expanding the phenotypic spectrum.

Splicing analyses for variants in MMR genes: best practice recommendations from the European Mismatch Repair Working Group

Over 20% of the DNA mismatch repair (MMR) germline variants in suspected Lynch syndrome patients are classified as variants of uncertain significance (VUS). Well-established functional assays are pivotal for assessing the biological impact of these variants and provide relevant evidence for clinical classification. In our collaborative European Mismatch Repair Working Group (EMMR-WG) we compared three different experimental approaches for evaluating the effect of seven variants on mRNA splicing in MMR genes: (i) RT-PCR of full-length transcripts (FLT), (ii) RT-PCR of targeted transcript sections (TTS), both from patient biological samples and (iii) minigene splicing assays. An overall good concordance was observed between splicing patterns in TTS, FLT and minigene analyses for all variants. The FLT analysis depicted a higher number of different isoforms and mitigated PCR-bias towards shorter isoforms. TTS analyses may miss aberrant isoforms and minigene assays may under/overestimate the severity of certain splicing defects. The interpretation of the experimental findings must be cautious to adequately discriminate abnormal events from physiological complex alternative splicing patterns. A consensus strategy for investigating the impact of MMR variants on splicing was defined. First, RNA should be obtained from patient's cell cultures (such as fresh lymphocyte cultures) incubated with/without a nonsense-mediated decay inhibitor. Second, FLT RT-PCR analysis is recommended to oversee all generated isoforms. Third, TTS analysis and minigene assays are useful independent approaches for verifying and clarifying FLT results. The use of several methodologies is likely to increase the strength of the experimental evidence which contributes to improve variant interpretation.

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.