Compensatory signals associated with the activation of human GC 5' splice sites

Splicing Dysregulation of Non-Canonical GC-5' Splice Sites of Breast Cancer Susceptibility Genes ATM and PALB2

Background/Objectives: The non-canonical GC-5' splice sites (5'ss) are the most common exception (~1%) to the classical GT/AG splicing rule. They constitute weak 5'ss and can be regulated by splicing factors, so they are especially sensitive to genetic variations inducing the misrecognition of their respective exons. We aimed to investigate the GC-5'ss of the breast/ovarian cancer susceptibility genes, ATM (exon 50), BRIP1 (exon 1), and PALB2 (exon 12), and their dysregulation induced by DNA variants. Methods: Splicing assays of the minigenes, mgATM_49-52, mgBRIP1_1-2, and mgPALB2_5-12, were conducted to study the regulation of the indicated GC-5'ss. Results: A functional map of the splicing regulatory elements (SRE) formed by overlapping exonic microdeletions revealed three essential intervals, ATM c.7335_7344del, PALB2 c.3229_3258del, and c.3293_3322del, which are likely targets for spliceogenic SRE-variants. We then selected 14 ATM and 9 PALB2 variants (Hexplorer score < -40) located at these intervals that were assayed in MCF-7 cells. Nine ATM and three PALB2 variants affected splicing, impairing the recognition of exons 50 and 12, respectively. Therefore, these variants likely disrupt the active SREs involved in the inclusion of both exons in the mature mRNA. DeepCLIP predictions suggested the participation of several splicing factors in exon recognition, including SRSF1, SRSF2, and SRSF7, involved in the recognition of other GC sites. The ATM spliceogenic variants c.7336G>T (p.(Glu2446Ter)) and c.7340T>A (p.(Leu2447Ter)) produced significant amounts of full-length transcripts (55-59%), which include premature termination stop codons, so they would inactivate ATM through both splicing disruption and protein truncation mechanisms. Conclusions: ATM exon 50 and PALB2 exon 12 require specific sequences for efficient recognition by the splicing machinery. The mapping of SRE-rich intervals in minigenes is a valuable approach for the identification of spliceogenic variants that outperforms any prediction software. Indeed, 12 spliceogenic SRE-variants were identified in the critical intervals.

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Inefficient splicing of long non-coding RNAs is associated with higher transcript complexity in human and mouse

Recent reports show that long non-coding RNAs (lncRNAs) have inefficient splicing and fewer alternative splice variants than mRNAs. Here, we have explored the efficiency of lncRNAs and mRNAs in producing various splice variants, given the number of exons in humans and mice. Intriguingly, lncRNAs produce more splice variants per exon, referred to as Transcript Complexity, than mRNAs. Most lncRNA splice variants are the product of the alternative last exon and exon skipping. LncRNAs and mRNAs with higher transcript complexity have shorter intron lengths. Longer exon length and GC/AG at 5'/3' splice sites are associated with higher transcript complexity in lncRNAs. Lastly, our results indicate that inefficient splicing of lncRNAs may facilitate multiple introns splicing and, thus, more spliced products per exon.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

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

RAD51D Aberrant Splicing in Breast Cancer: Identification of Splicing Regulatory Elements and Minigene-Based Evaluation of 53 DNA Variants

RAD51D loss-of-function variants increase lifetime risk of breast and ovarian cancer. Splicing disruption is a frequent pathogenic mechanism associated with variants in susceptibility genes. Herein, we have assessed the splicing and clinical impact of splice-site and exonic splicing enhancer (ESE) variants identified through the study of ~113,000 women of the BRIDGES cohort. A RAD51D minigene with exons 2-9 was constructed in splicing vector pSAD. Eleven BRIDGES splice-site variants (selected by MaxEntScan) were introduced into the minigene by site-directed mutagenesis and tested in MCF-7 cells. The 11 variants disrupted splicing, collectively generating 25 different aberrant transcripts. All variants but one produced negligible levels (<3.4%) of the full-length (FL) transcript. In addition, ESE elements of the alternative exon 3 were mapped by testing four overlapping exonic microdeletions (≥30-bp), revealing an ESE-rich interval (c.202_235del) with critical sequences for exon 3 recognition that might have been affected by germline variants. Next, 26 BRIDGES variants and 16 artificial exon 3 single-nucleotide substitutions were also assayed. Thirty variants impaired splicing with variable amounts (0-65.1%) of the FL transcript, although only c.202G>A demonstrated a complete aberrant splicing pattern without the FL transcript. On the other hand, c.214T>C increased efficiency of exon 3 recognition, so only the FL transcript was detected (100%). In conclusion, 41 RAD51D spliceogenic variants (28 of which were from the BRIDGES cohort) were identified by minigene assays. We show that minigene-based mapping of ESEs is a powerful approach for identifying ESE hotspots and ESE-disrupting variants. Finally, we have classified nine variants as likely pathogenic according to ACMG/AMP-based guidelines, highlighting the complex relationship between splicing alterations and variant interpretation.

Transposon clusters as substrates for aberrant splice-site activation

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5ʹ splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5ʹ splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

Comprehensive genetic diagnosis of patients with Duchenne/Becker muscular dystrophy (DMD/BMD) and pathogenicity analysis of splice site variants in the DMD gene

Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) are caused by mutations in the DMD gene. The aim of this study is to identify pathogenic DMD variants in probands and reduce the risk of recurrence of the disease in affected families. Variations in 100 unrelated DMD/BMD patients were detected by multiplex ligation-dependent probe amplification (MLPA) and next-generation sequencing (NGS). Pathogenic variants in DMD were successfully identified in all cases, and 11 of them were novel. The most common mutations were intragenic deletions (69%), with two hotspots located in the 5' end (exons 2-19) and the central of the DMD gene (exons 45-55), while point mutations were observed in 22% patients. Further, c.1149+1G>A and c.1150-2A>G were confirmed by hybrid minigene splicing assay (HMSA). This two splice site mutations would lead to two aberrant DMD isoforms which give rise to severely truncated protein. Therefore, the clinical use of MLPA, NGS, and HMSA is an effective strategy to identify variants. Importantly, eight embryos were terminated pregnancies according to prenatal diagnosis and a healthy boy was successfully delivered by preimplantation genetic diagnosis (PGD). Early and accurate genetic diagnosis is essential for prenatal diagnosis/PGD to reduce the risk of recurrence of DMD in affected families.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

First estimate of the scale of canonical 5' splice site GT>GC variants capable of generating wild-type transcripts

It has long been known that canonical 5' splice site (5'SS) GT>GC variants may be compatible with normal splicing. However, to date, the actual scale of canonical 5'SSs capable of generating wild-type transcripts in the case of GT>GC substitutions remains unknown. Herein, combining data derived from a meta-analysis of 45 human disease-causing 5'SS GT>GC variants and a cell culture-based full-length gene splicing assay of 103 5'SS GT>GC substitutions, we estimate that ~15-18% of canonical GT 5'SSs retain their capacity to generate between 1% and 84% normal transcripts when GT is substituted by GC. We further demonstrate that the canonical 5'SSs in which substitution of GT by GC-generated normal transcripts exhibit stronger complementarity to the 5' end of U1 snRNA than those sites whose substitutions of GT by GC did not lead to the generation of normal transcripts. We also observed a correlation between the generation of wild-type transcripts and a milder than expected clinical phenotype but found that none of the available splicing prediction tools were capable of reliably distinguishing 5'SS GT>GC variants that generated wild-type transcripts from those that did not. Our findings imply that 5'SS GT>GC variants in human disease genes may not invariably be pathogenic.

Mis-splicing in breast cancer: identification of pathogenic BRCA2 variants by systematic minigene assays

Splicing disruption is a common mechanism of gene inactivation associated with germline variants of susceptibility genes. To study the role of BRCA2 mis-splicing in hereditary breast/ovarian cancer (HBOC), we performed a comprehensive analysis of variants from BRCA2 exons 2-9, as well as the initial characterization of the regulatory mechanisms of such exons. A pSAD-based minigene with exons 2-9 was constructed and validated in MCF-7 cells, producing the expected transcript (1016-nt/V1-BRCA2_exons_2-9-V2). DNA variants from mutational databases were analyzed by NNSplice and Human Splicing Finder softwares. To refine ESE-variant prediction, we mapped the regulatory regions through a functional strategy whereby 26 exonic microdeletions were introduced into the minigene and tested in MCF-7 cells. Thus, we identified nine spliceogenic ESE-rich intervals where ESE-variants may be located. Combining bioinformatics and microdeletion assays, 83 variants were selected and genetically engineered in the minigene. Fifty-three changes impaired splicing: 28 variants disrupted the canonical sites, four created new ones, 10 abrogated enhancers, eight created silencers and three caused a double-effect. Notably, nine spliceogenic-ESE variants were located within ESE-containing intervals. Capillary electrophoresis and sequencing revealed more than 23 aberrant transcripts, where exon skipping was the most common event. Interestingly, variant c.67G>A triggered the usage of a noncanonical GC-donor 4-nt upstream. Thirty-six variants that induced severe anomalies (>60% aberrant transcripts) were analyzed according to the ACMG guidelines. Thus, 28 variants were classified as pathogenic, five as likely pathogenic and three as variants of uncertain significance. Interestingly, 13 VUS were reclassified as pathogenic or likely pathogenic variants. In conclusion, a large fraction of BRCA2 variants (∼64%) provoked splicing anomalies lending further support to the high prevalence of this disease-mechanism. The low accuracy of ESE-prediction algorithms may be circumvented by functional ESE-mapping that represents an optimal strategy to identify spliceogenic ESE-variants. Finally, systematic functional assays by minigenes depict a valuable tool for the initial characterization of splicing anomalies and the clinical interpretation of variants. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Eleven percent intact PGM3 in a severely immunodeficient patient with a novel splice-site mutation, a case report

Background

A novel immunodeficiency, frequently accompanied by high serum-IgE, and caused by mutations in the PGM3 gene was described in 2014. To date there are no unique phenotype characteristics for PGM3 deficiency. PGM3 encodes a carbohydrate-modifying enzyme, phosphoglucomutase 3. Null-mutations are quite likely lethal, and to date only missense mutations or small deletions have been reported. Such mutations frequently cause a combination of reduced enzyme activity and protein instability, complicating determination of the enzyme level needed for survival. Here we present the first patient with a homozygous splice-modifying mutation in the PGM3 gene. An A > G substitution at position c.871 + 3 (transcript NM_001199917) is causing a deletion of exon 7 in the majority of PGM3 transcripts. In addition, this case further increases the clinical phenotypes of immunodeficiency caused by PGM3 mutations.

Case presentation

We describe the symptoms of a 3-year-old girl who was severely growth retarded, had vascular malformations, extensive eczema, multiple food-allergies, and was prone to infections. Unlike the majority of reported PGM3 deficient patients she lacked skeletal dysplasia and had normal neurocognitive development. In addition to the high serum-IgE, she displayed altered T cell numbers with reduced naïve CD4⁺ and CD8⁺ T-cells, increased number of activated effector memory CD8⁺ T cells and aberrant T-cell functions.

The patient was homozygous for a new hypomorphic, splice-modifying mutation in the PGM3 gene, causing severely reduced mRNA levels. In the patient’s cells, we observed 5% intact mRNA and approximately 11% of the protein levels seen in healthy controls.

Treatment with allogeneic hematopoietic stem cell therapy was planned, but unfortunately the clinical condition deteriorated with multi-organ failure, which led to her death at 3 years of age.

Conclusions

There is still no specific phenotype identified that distinguishes immunodeficiency caused by PGM3 mutations from other forms of immunodeficiency. The patient described here yields new information on the phenotypic variability among these patients. In addition, since all the synthesized protein is wild-type, it is possible for the first time to estimate the enzyme activity in vivo. The results suggest that1/10 of the normal PGM3 level is sufficient for survival but that it is insufficient for accurate carbohydrate processing.

Electronic supplementary material

The online version of this article (10.1186/s12887-018-1258-9) contains supplementary material, which is available to authorized users.

Characterization of spliceogenic variants located in regions linked to high levels of alternative splicing: BRCA2 c.7976+5G > T as a case study

Many BRCA1 and BRCA2 (BRCA1/2) genetic variants have been studied at mRNA level and linked to hereditary breast and ovarian cancer due to splicing alteration. In silico tools are reliable when assessing variants located in consensus splice sites, but we may identify variants in complex genomic contexts for which bioinformatics is not precise enough. In this study, we characterize BRCA2 c.7976 + 5G > T variant located in intron 17 which has an atypical donor site (GC). This variant was identified in three unrelated Spanish families and we have detected exon 17 skipping as the predominant transcript occurring in carriers. We have also detected several isoforms (Δ16-18, Δ17,18, Δ18, and ▼17q²²⁴ ) at different expression levels among carriers and controls. This study remarks the challenge of interpreting genetic variants when multiple alternative isoforms are present, and that caution must be taken when using in silico tools to identify potential spliceogenic variants located in GC-AG introns.

The cyclin-dependent kinase G group defines a thermo-sensitive alternative splicing circuit modulating the expression of Arabidopsis ATU2AF65A

The ability to adapt growth and development to temperature variations is crucial to generate plant varieties resilient to predicted temperature changes. However, the mechanisms underlying plant response to progressive increases in temperature have just started to be elucidated. Here, we report that the cyclin-dependent kinase G1 (CDKG1) is a central element in a thermo-sensitive mRNA splicing cascade that transduces changes in ambient temperature into differential expression of the fundamental spliceosome component, ATU2AF65A. CDKG1 is alternatively spliced in a temperature-dependent manner. We found that this process is partly dependent on both the cyclin-dependent kinase G2 (CDKG2) and the interacting co-factor CYCLIN L1 (CYCL1), resulting in two distinct messenger RNAs. The relative abundance of both CDKG1 transcripts correlates with ambient temperature and possibly with different expression levels of the associated protein isoforms. Both CDKG1 alternative transcripts are necessary to fully complement the expression of ATU2AF65A across the temperature range. Our data support a previously unidentified temperature-dependent mechanism based on the alternative splicing (AS) of CDKG1 and regulated by CDKG2 and CYCL1. We propose that changes in ambient temperature affect the relative abundance of CDKG1 transcripts, and this in turn translates into differential CDKG1 protein expression coordinating the AS of ATU2AF65A.

Activation of the intronic cryptic 5' splice site depends on its distance to the upstream cassette exon

Splice site selection is a key step that determines the mRNA isoforms generated from a single transcript. The large diversity in splice site sequences emphasizes the plasticity of splice site recognition and selection. In this report, a cell-based reporter system using a SMN1/2 cassette exon was applied to study the roles governing the activation of a cryptic 5'SS from the intron 4 of the CT/CGRP gene. We found that the cryptic site was activated when placed within 124nt downstream the cassette exon, and the level of activation was negatively correlated with its distance from the exon. In addition, activation was not affected by PTB but was eliminated by an insertion extending the exon length. Activated cryptic 5'SSs in intron or exon could override the original alternative 5'SS, obeying the U1 base-pairing rule. These results suggest that the exon length itself could represent a factor in determining the splice site selection.

Functional classification of DNA variants by hybrid minigenes: Identification of 30 spliceogenic variants of BRCA2 exons 17 and 18

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ex14-20) was constructed in the pSAD vector. Fifty-two candidate variants were selected with splicing prediction programs, introduced in MGBR2_ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ex14-20 produced a stable transcript of the expected size (1,806 nucleotides) and structure (V1-[BRCA2_exons_14–20]–V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3’ end of exon 17 (c.7944-7973) and the 5’ end of exon 18 (c.7979-7988, c.7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with >16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c.7806-2A>G, c.7806-1G>A, c.7806-1G>T, c.7806-1_7806-2dup, c.7976+1G>A, c.7977-3_7978del, c.7977-2A>T, c.7977-1G>T, c.7977-1G>C, c.8009C>A, c.8331+1G>T and c.8331+2T>C) or likely pathogenic (c.7806-9T>G, c.7976G>C, c.7976G>A, c.7977-7C>G, c.7985C>G, c.8023A>G, c.8035G>T and c.8331G>A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17–18 at the BRCA Share database. The remaining 8 variants (c.7975A>G, c.7977-6T>G, c.7988A>T, c.7992T>A, c.8007A>G, c.8009C>T, c.8009C>G, and c.8072C>T) induced partial splicing anomalies with important ratios of the full-length transcript (≥70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA.

A significant proportion of disease-causing mutations of inherited disorders impair splicing. Massive sequencing projects of genetic diseases generate thousands of sequence variations that require functional and clinical interpretations. We have shown that splicing reporter minigenes of the breast cancer genes BRCA1 and BRCA2 are useful tools to functionally test DNA variants. In this work, we have constructed a 7-exon BRCA2 minigene (exons 14 to 20) where we mapped critical splicing regulatory sequences and tested 52 selected variants of exons 17 and 18 detected in breast cancer patients. We finely located three DNA segments on both exons that presumably contain splicing enhancer sequences. We observed that a total of 30 variants of any type disrupted the splicing patterns and, given the severity of their outcomes, we classified 20 of them as pathogenic or likely pathogenic. We also showed that a wide range of splicing elements were affected including canonical and novel 5’ and 3’ splice sites, the polypyrimidine tract and enhancer and silencer sequences. We concluded that splicing aberrations are frequent in Hereditary Breast and Ovarian Cancer and that minigenes are valuable tools to functionally classify DNA variants of any human disease gene under the splicing viewpoint.

Antisense Oligonucleotides Modulating Activation of a Nonsense-Mediated RNA Decay Switch Exon in the ATM Gene

ATM (ataxia-telangiectasia, mutated) is an important cancer susceptibility gene that encodes a key apical kinase in the DNA damage response pathway. ATM mutations in the germ line result in ataxia-telangiectasia (A-T), a rare genetic syndrome associated with hypersensitivity to double-strand DNA breaks and predisposition to lymphoid malignancies. ATM expression is limited by a tightly regulated nonsense-mediated RNA decay (NMD) switch exon (termed NSE) located in intron 28. In this study, we identify antisense oligonucleotides that modulate NSE inclusion in mature transcripts by systematically targeting the entire 3.1-kb-long intron. Their identification was assisted by a segmental deletion analysis of transposed elements, revealing NSE repression upon removal of a distant antisense Alu and NSE activation upon elimination of a long terminal repeat transposon MER51A. Efficient NSE repression was achieved by delivering optimized splice-switching oligonucleotides to embryonic and lymphoblastoid cells using chitosan-based nanoparticles. Together, these results provide a basis for possible sequence-specific radiosensitization of cancer cells, highlight the power of intronic antisense oligonucleotides to modify gene expression, and demonstrate transposon-mediated regulation of NSEs.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Molecular Genetic Characterization of Thyroid Dyshormonogenesis in a French Bulldog

Background

A case of congenital hypothyroidism with goiter (CHG) in a juvenile French bulldog was identified and hypothesized to be caused by dyshormonogenesis of genetic etiology.

Objectives

To describe case management, unusual phenotypic aspects, and a CHG‐causing mutation in a French bulldog.

Animals

Thyroid tissue and blood from a CHG‐affected French bulldog and 4 normal control dogs and buccal brush samples of 125 French bulldogs were studied.

Methods

Standard clinical assessment and laboratory tests were applied. Thyroid peroxidase (TPO) iodide oxidation activity was measured in vitro, and TPO protein was assessed on Western blots. Thyroid peroxidase exons and flanking splice sites were amplified from genomic DNA and sequenced. Thyroid peroxidase cDNA was amplified from thyroid RNA and sequenced.

Results

At 9 months of age, the affected dog had signs of cretinism, but near‐normal skeletal maturation. The enlarged thyroid glands exhibited noninflammatory fibrosis and aberrant follicular organization. Thyroid peroxidase activity and immunocrossreactive protein were undetectable. There was a T>C mutation of the intron 12 splice donor consensus that caused abnormally spliced mRNA, consistent with absent TPO function. The mutant allele was not observed in 125 clinically normal French bulldogs.

Conclusions

Presumptive CHG in a French bulldog with unusual clinical presentation is described. Genetic etiology was confirmed by identifying the underlying TPO mutation.

Splice-correction strategies for treatment of X-linked agammaglobulinemia

X-linked agammaglobulinemia (XLA) is a primary immunodeficiency disease caused by mutations in the gene coding for Bruton’s tyrosine kinase (BTK). Deficiency of BTK leads to a developmental block in B cell differentiation; hence, the patients essentially lack antibody-producing plasma cells and are susceptible to various infections. A substantial portion of the mutations in BTK results in splicing defects, consequently preventing the formation of protein-coding mRNA. Antisense oligonucleotides (ASOs) are therapeutic compounds that have the ability to modulate pre-mRNA splicing and alter gene expression. The potential of ASOs has been exploited for a few severe diseases, both in pre-clinical and clinical studies. Recently, advances have also been made in using ASOs as a personalized therapy for XLA. Splice-correction of BTK has been shown to be feasible for different mutations in vitro, and a recent proof-of-concept study demonstrated the feasibility of correcting splicing and restoring BTK both ex vivo and in vivo in a humanized bacterial artificial chromosome (BAC)-transgenic mouse model. This review summarizes the advances in splice correction, as a personalized medicine for XLA, and outlines the promises and challenges of using this technology as a curative long-term treatment option.

Splice-correcting oligonucleotides restore BTK function in X-linked agammaglobulinemia model

X-linked agammaglobulinemia (XLA) is an inherited immunodeficiency that results from mutations within the gene encoding Bruton's tyrosine kinase (BTK). Many XLA-associated mutations affect splicing of BTK pre-mRNA and severely impair B cell development. Here, we assessed the potential of antisense, splice-correcting oligonucleotides (SCOs) targeting mutated BTK transcripts for treating XLA. Both the SCO structural design and chemical properties were optimized using 2'-O-methyl, locked nucleic acid, or phosphorodiamidate morpholino backbones. In order to have access to an animal model of XLA, we engineered a transgenic mouse that harbors a BAC with an authentic, mutated, splice-defective human BTK gene. BTK transgenic mice were bred onto a Btk knockout background to avoid interference of the orthologous mouse protein. Using this model, we determined that BTK-specific SCOs are able to correct aberrantly spliced BTK in B lymphocytes, including pro-B cells. Correction of BTK mRNA restored expression of functional protein, as shown both by enhanced lymphocyte survival and reestablished BTK activation upon B cell receptor stimulation. Furthermore, SCO treatment corrected splicing and restored BTK expression in primary cells from patients with XLA. Together, our data demonstrate that SCOs can restore BTK function and that BTK-targeting SCOs have potential as personalized medicine in patients with XLA.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

Insulin receptor substrate-1 associates with small nucleolar RNA which contributes to ribosome biogenesis

Insulin receptor substrates (IRSs) are well known to play crucial roles in mediating intracellular signals of insulin-like growth factors (IGFs)/insulin. Previously, we showed that IRS-1 forms high molecular mass complexes containing RNAs. To identify RNAs in IRS-1 complexes, we performed ultraviolet (UV) cross-linking and immunoprecipitation analysis using HEK293 cells expressing FLAG-IRS-1 and FLAG-IRS-2. We detected the radioactive signals in the immunoprecipitates of FLAG-IRS-1 proportional to the UV irradiation, but not in the immunoprecipitates of FLAG-IRS-2, suggesting the direct contact of RNAs with IRS-1. RNAs cross-linked to IRS-1 were then amplified by RT-PCR, followed by sequence analysis. We isolated sequence tags attributed to 25 messenger RNAs and 8 non-coding RNAs, including small nucleolar RNAs (snoRNAs). We focused on the interaction of IRS-1 with U96A snoRNA (U96A) and its host Rack1 (receptor for activated C kinase 1) pre-mRNA. We confirmed the interaction of IRS-1 with U96A, and with RACK1 pre-mRNA by immunoprecipitation with IRS-1 followed by Northern blotting or RT-PCR analyses. Mature U96A in IRS-1(-/-) mouse embryonic fibroblasts was quantitatively less than WT. We also found that a part of nuclear IRS-1 is localized in the Cajal body, a nuclear subcompartment where snoRNA mature. The unanticipated function of IRS-1 in snoRNA biogenesis highlights the potential of RNA-associated IRS-1 complex to open a new line of investigation to dissect the novel mechanisms regulating IGFs/insulin-mediated biological events.

CYCLIN-DEPENDENT KINASE G1 is associated with the spliceosome to regulate CALLOSE SYNTHASE5 splicing and pollen wall formation in Arabidopsis

Arabidopsis thaliana CYCLIN-DEPEDENT KINASE G1 (CDKG1) belongs to the family of cyclin-dependent protein kinases that were originally characterized as cell cycle regulators in eukaryotes. Here, we report that CDKG1 regulates pre-mRNA splicing of CALLOSE SYNTHASE5 (CalS5) and, therefore, pollen wall formation. The knockout mutant cdkg1 exhibits reduced male fertility with impaired callose synthesis and abnormal pollen wall formation. The sixth intron in CalS5 pre-mRNA, a rare type of intron with a GC 5' splice site, is abnormally spliced in cdkg1. RNA immunoprecipitation analysis suggests that CDKG1 is associated with this intron. CDKG1 contains N-terminal Ser/Arg (RS) motifs and interacts with splicing factor Arginine/Serine-Rich Zinc Knuckle-Containing Protein33 (RSZ33) through its RS region to regulate proper splicing. CDKG1 and RS-containing Zinc Finger Protein22 (SRZ22), a splicing factor interacting with RSZ33 and U1 small nuclear ribonucleoprotein particle (snRNP) component U1-70k, colocalize in nuclear speckles and reside in the same complex. We propose that CDKG1 is recruited to U1 snRNP through RSZ33 to facilitate the splicing of the sixth intron of CalS5.

Clinical and genetic analysis of patients with vitamin D-dependent rickets type 1A

CONTEXT

Vitamin D-dependent rickets type 1A (VDDR-IA, OMIM 264700) is a rare autosomal recessive disorder and is caused by mutations in the CYP27B1 gene.

OBJECTIVES

We aim to investigate CYP27B1 mutation in seven patients from four separate families and characterize the genotype-phenotype correlation.

METHODS

The entire coding region of the CYP27B1 gene was sequenced, and genotype-phenotype correlation among patients was assessed.

RESULTS

Sequencing analysis identified biallelic CYP27B1 mutations in all patients and monoallelic mutations in their parents. One patient from the first family was compound heterozygous for c.1166G>A (p.Arg389His) and a novel nonsense mutation c.1079 C>A (p.Ser360*). Two patients from the second family were homozygous for a novel splice donor site mutation in intron 1 (c.195 + 2 T>G), causing partial retention of the intron and a shift in the reading frame. Both novel mutations lead to the complete loss of vitamin D1α-hydroxylase activity. Four patients from families 3 and 4 were homozygous for a previously reported duplication mutation in exon 8 (1319-1325dupCCCACCC, Phe443Profs*24). Interestingly, one patient who was presented with severe hypocalcaemia and seizures at 4 months of age as a result of Phe443Profs*24 has improved spontaneously since 11 years of age and does not need regular treatment. Her laboratory tests showed normal serum calcium and 1,25(OH)(2) D after refusing to take medication for 12 months.

CONCLUSIONS

There is a good genotype-phenotype correlation in VDDR-IA. However, some patients may recover from the loss of CYP27B1 function, probably due to 1α-hydroxylase activity exerted by a non-CYP27B1 enzyme.

Complexities of 5'splice site definition: implications in clinical analyses

In higher eukaryotes, the 5' splice site (5'ss) is initially recognized through an RNA-RNA interaction by U1 small nuclear ribonucleoprotein (U1 snRNP). This event represents one of the key steps in initial spliceosomal assembly and many disease-associated mutations in humans often disrupt this process. Beside base pair complementarity, 5'ss recognition can also be modified by additional factors such as RNA secondary structures or the specific binding of other nuclear proteins. In this work, we have focused on investigating a few examples of changes detected within the 5'ss in patients, that would not be immediately considered "disease causing mutations". We show that the splicing outcome of very similar mutations can be very different due to variations in trans-acting factor(s) interactions and specific context influences. Using several NF1 donor sites and SELEX approaches as experimental models, we have examined the binding properties of particular sequence motifs such as GGGU found in donor sites, and how the sequence context can change their interaction with hnRNPs such as H/F and A1/A2. Our results clearly show that even minor differences in local nucleotide context can differentially affect the binding ability of these factors to the GGGU core. Finally, using a previously identified mutation in KCNH2 that resulted in intron retention we show how very similar 5'ss mutations found in patients can have a very different splicing outcome due to the neighbouring sequence context, thus highlighting the general need to approach splicing problems with suitable experimental approaches.