Comparative analysis identifies exonic splicing regulatory sequences--The complex definition of enhancers and silencers

An interpretable model of pre-mRNA splicing for animal and plant genes

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

CRISPR-dCas13d-based deep screening of proximal and distal splicing-regulatory elements

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is hampered by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically inactive CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identifies not only known SREs but also a previously unknown distal intronic SRE, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Challenges in the Definitive Diagnosis of Niemann-Pick Type C-Leaky Variants and Alternative Transcripts

Niemann-Pick type C (NPC, ORPHA: 646) is a neuro-visceral, psychiatric disease caused predominantly by pathogenic variants in the NPC1 gene or seldom in NPC2. The rarity of the disease, and its wide range of clinical phenotypes and ages of onset, turn the diagnosis into a significant challenge. Other than the detailed clinical history, the typical diagnostic work-up for NPC includes the quantification of pathognomonic metabolites. However, the molecular basis diagnosis is still of utmost importance to fully characterize the disorder. Here, the authors provide an overview of splicing variants in the NPC1 and NPC2 genes and propose a new workflow for NPC diagnosis. Splicing variants cover a significant part of the disease-causing variants in NPC. The authors used cDNA analysis to study the impact of such variants, including the collection of data to classify them as leaky or non-leaky pathogenic variants. However, the presence of naturally occurring spliced transcripts can misdiagnose or mask a pathogenic variant and make the analysis even more difficult. Analysis of the NPC1 cDNA in NPC patients in parallel with controls is vital to assess and detect alternatively spliced forms. Moreover, nonsense-mediated mRNA decay (NMD) analysis plays an essential role in evaluating the naturally occurring transcripts during cDNA analysis and distinguishing them from other pathogenic variants' associated transcripts.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Deep screening of proximal and distal splicing-regulatory elements in a native sequence context

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is impeded by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically dead CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identified not only known SREs, but also a novel distal intronic splicing enhancer, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

Minigene splicing assays reveal new insights into exonic variants of the SLC12A3 gene in Gitelman syndrome

BACKGROUND

Gitelman syndrome (GS) is a type of salt-losing tubular disease, most of which is caused by SLC12A3 gene variants, and missense variants account for the majority. Recently, the phenomenon of exon skipping, in which variants disrupt normal pre-mRNA splicing, has been related to a variety of diseases. Therefore, we hypothesize that a certain proportion of SLC12A3 variants can result in disease via interfering with the normal splicing process.

METHODS

We analyzed 342 previously presumed SLC12A3 missense variants using bioinformatics programs and identified candidate variants that may alter the splicing of pre-mRNA through minigene assays.

RESULTS

Our study revealed that, among ten candidate variants, six variants (c.602G>A, c.602G>T, c.1667C>T, c.1925G>A, c.2548G>C, and c.2549G>C) led to complete or incomplete exon skipping by affecting exonic splicing regulatory elements and/or disturbing canonical splice sites.

CONCLUSION

It is worth mentioning that this is the largest study on pre-mRNA splicing of SLC12A3 exonic variants. In addition, our study emphasizes the importance of detecting splicing function at the mRNA level in GS and indicates that minigene analysis is a valuable tool for splicing functional assays of variants in vitro.

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.

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.

Transgene-design: a web application for the design of mammalian transgenes

Summary

Transgene-design is a web application to help design transgenes for use in mammalian studies. It is predicated on the recent discovery that human intronless transgenes and native retrogenes can be expressed very effectively if the GC content at exonic synonymous sites is high. In addition, as exonic splice enhancers resident in intron containing genes may have different utility in intronless genes, these can be reduced or increased in density. Input can be a native gene or a commercially ‘optimised’ gene. The option to leave in the first intron and to protect or avoid other motifs is also permitted.

Availability and implementation

Transgene-design is based on a ruby for rails platform. The application is available at https://transgene-design.bath.ac.uk. The code is available under GNU General Public License from GitHub (https://github.com/smuehlh/transgenes).

Supplementary information

Supplementary data are available at Bioinformatics online.

Tra2beta-Dependent Regulation of RIO Kinase 3 Splicing During Rift Valley Fever Virus Infection Underscores the Links Between Alternative Splicing and Innate Antiviral Immunity

Rift Valley fever virus (RVFV) is an emerging pathogen that has potential to cause severe disease in humans and domestic livestock. Propagation of RVFV strain MP-12 is negatively impacted by the actions of RIOK3, a protein involved in the cellular immune response to viral infection. During RVFV infection, RIOK3 mRNA is alternatively spliced to produce an isoform that correlates with the inhibition of interferon β signaling. Here, we identify splicing factor TRA2-β (also known as TRA2beta and hTRA2-β) as a key regulator governing the relative abundance of RIOK3 splicing isoforms. Using RT-PCR and minigenes, we determined that TRA2-β interaction with RIOK3 pre-mRNA was necessary for constitutive splicing of RIOK3 mRNA, and conversely, lack of TRA2-β engagement led to increased alternative splicing. Expression of TRA2-β was found to be necessary for RIOK3's antiviral effect against RVFV. Intriguingly, TRA2-β mRNA is also alternatively spliced during RVFV infection, leading to a decrease in cellular TRA2-β protein levels. These results suggest that splicing modulation serves as an immune evasion strategy by RVFV and/or is a cellular mechanism to prevent excessive immune response. Furthermore, the results suggest that TRA2-β can act as a key regulator of additional steps of the innate immune response to viral infection.

Personalized Medicine for Kidney Transplantation: Association of Graft Survival and Acute Transplant Rejection with Genetic Variation in B Cell Activating Factor System Signaling

Kidney transplantation (KT) clinical outcomes are highly variable across patients and would benefit from predictive biomarkers to achieve personalized/precision medicine. The B cell activating factor (BAFF) system signaling plays an essential role in B lymphocytes' homeostasis, and is implicated in activation and survival of B lymphocytes. Single nucleotide polymorphisms (SNPs) in BAFF system genes are therefore strong candidates to identify the genetic mechanisms underpinning variable clinical outcomes in KT. We report here new findings on BAFF system genetic polymorphisms in KT patients in relation to two key phenotypes of clinical interest: graft survival and acute rejection (AR). A total of 168 KT patients, of which 29 suffered AR, participated in this study. The BAFF system polymorphisms in five genes TNFSF13B, TNFSF13, TNFRSF13C, TNFRSF13B, and TNFRSF17 were characterized using TaqMan SNP genotyping. Patients with KT who had an AA genotype in polymorphism rs3803800 of the TNFSF13 gene had a higher risk of suffering AR (p = 0.046; odds ratios = 3.38, 95% CI: 1.02-11.2). Moreover, patients with AA genotype (rs3803800) in the TNFSF13 gene had a significantly lower AR-free time than the GG/GA genotypes (69.2% vs. 85.7%; p = 0.037). Of importance, bioinformatics analysis showed that the polymorphism rs3803800 could alter splicing regulation and affect the proliferation-inducing ligand (APRIL) expression levels. The analysis of graft survival did not show a significant association with the polymorphisms analyzed in this study. In conclusion, the rs3803800 genetic polymorphism from this study of BAFF system genes appears to display importance in AR-free time for KT patients, and thus, warrants further research in independent populations as a putative predictive biomarker of AR. These findings also inform future personalized/precision medicine efforts and functional genomic studies in KT patients.

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

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

Functional analysis of novel genetic variants of NKX2-5 associated with nonsyndromic congenital heart disease

NKX2-5, a master cardiac regulatory transcription factor was the first known genetic cause of congenital heart diseases (CHDs). To further investigate its role in CHD pathogenesis, we performed mutational screening of 285 CHD probands and 200 healthy controls. Five coding sequence variants were identified in six CHD cases (2.1%), including three in the N-terminal region (p.A61G, p.R95L, and p.E131K) and one each in homeodomain (HD) (p.A148E) and tyrosine-rich domain (p.P247A). Variant-p.A148E showed tertiary structure changes and differential DNA binding affinity of mutant compared to wild type. Two N-terminal variants-p.A61G and p.E131K along with HD variant p.A148E demonstrated significantly reduced transcriptional activity of Nppa and Actc1 promoters in dual luciferase promoter assay supported by their reduced expression in qRT-PCR. Nonetheless, variant p.R95L affected the synergy of NKX2-5 with serum response factor and TBX5 leading to significantly decreased Actc1 promoter activity depicting a distinctive role of this region. The aberrant expression of other target genes-Irx4, Mef2c, Bmp10, Myh6, Myh7, and Myocd is also observed in response to NKX2-5 variants, possibly due to the defective gene regulatory network. Severely impaired downstream promoter activities and abnormal expression of target genes due to N-terminal variants supports the emerging role of this region during cardiac-developmental pathways.

F5-Atlanta: A novel mutation in F5 associated with enhanced East Texas splicing and FV-short production

BACKGROUND

Elucidating the molecular pathogenesis underlying East Texas bleeding disorder (ET) led to the discovery of alternatively spliced F5 transcripts harboring large deletions within exon 13. These alternatively spliced transcripts produce a shortened form of coagulation factor V (FV) in which a large portion of its B-domain is deleted. These FV isoforms bind tissue factor pathway inhibitor alpha (TFPIα) with high affinity, prolonging its circulatory half-life and enhancing its anticoagulant effects. While two missense pathogenic variants highlighted this alternative splicing event, similar internally deleted FV proteins are found in healthy controls.

OBJECTIVE

We identified a novel heterozygous 832 base pair deletion within F5 exon 13, termed F5-Atlanta (F5-ATL), in a patient with severe bleeding. Our objective is to investigate the effect of this deletion on F5 and FV expression.

METHODS & RESULTS

Assessment of patient plasma revealed markedly elevated levels of total and free TFPI and a FV isoform similar in size to the FV-short described in ET. Sequencing analyses of cDNA revealed the presence of a transcript alternatively spliced using the ET splice sites, thereby removing the F5-ATL deletion. This alternative splicing pattern was recapitulated by heterologous expression in mammalian cells.

CONCLUSIONS

These findings support a mechanistic model consisting of cis-acting regulatory sequences encoded within F5 exon 13 that control alternative splicing at the ET splice sites and thereby regulate circulating FV-short and TFPIα levels.

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.

NPC1 silent variant induces skipping of exon 11 (p.V562V) and unfolded protein response was found in a specific Niemann-Pick type C patient

Background

Niemann‐Pick type C (NPC, MIM #257220) is a neuro‐visceral disease, caused predominantly by pathogenic variants in the NPC1 gene. Here we studied patients with clinical diagnosis of NPC but inconclusive results regarding the molecular analysis.

Methods

We used a Next‐Generation Sequencing (NGS)‐panel followed by cDNA analysis. Latter, we used massively parallel single‐cell RNA‐seq (MARS‐Seq) to address gene profiling changes and finally the effect of different variants on the protein and cellular levels.

Results

We identified novel variants and cDNA analysis allowed us to establish the functional effect of a silent variant, previously reported as a polymorphism. We demonstrated that this variant induces the skipping of exon 11 leading to a premature stop codon and identified it in NPC patients from two unrelated families. MARS‐Seq analysis showed that a number of upregulated genes were related to the unfolded protein response (UPR) and endoplasmic reticulum (ER) stress in one specific patient. Also, for all analyzed variants, the NPC1 protein was partially retained in the ER.

Conclusion

We showed that the NPC1 silent polymorphism (p.V562V) is a disease‐causing variant in NPC and that the UPR is upregulated in an NPC patient.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

Phylogenetic comparison of 5' splice site determination in central spliceosomal proteins of the U1-70K gene family, in response to developmental cues and stress conditions

Intron-containing genes have the ability to generate multiple transcript isoforms by splicing, thereby greatly expanding the eukaryotic transcriptome and proteome. In eukaryotic cells, precursor mRNA (pre-mRNA) splicing is performed by a mega-macromolecular complex defined as a spliceosome. Among its splicing components, U1 small nuclear ribonucleoprotein (U1 snRNP) is the smallest subcomplex involved in early spliceosome assembly and 5'-splice site recognition. Its central component, named U1-70K, has been extensively characterized in animals and yeast. Very few investigations on U1-70K genes have been conducted in plants, however. To this end, we performed a comprehensive study to systematically identify 115 U1-70K genes from 67 plant species, ranging from algae to angiosperms. Phylogenetic analysis suggested that the expansion of the plant U1-70K gene family was likely to have been driven by whole-genome duplications. Subsequent comparisons of gene structures, protein domains, promoter regions and conserved splicing patterns indicated that plant U1-70Ks are likely to preserve their conserved molecular function across plant lineages and play an important functional role in response to environmental stresses. Furthermore, genetic analysis using T-DNA insertion mutants suggested that Arabidopsis U1-70K may be involved in response to osmotic stress. Our results provide a general overview of this gene family in Viridiplantae and will act as a reference source for future mechanistic studies on this U1 snRNP-specific splicing factor.

A Depletion of Stop Codons in lincRNA is Owing to Transfer of Selective Constraint from Coding Sequences

Although the constraints on a gene’s sequence are often assumed to reflect the functioning of that gene, here we propose transfer selection, a constraint operating on one class of genes transferred to another, mediated by shared binding factors. We show that such transfer can explain an otherwise paradoxical depletion of stop codons in long intergenic noncoding RNAs (lincRNAs). Serine/arginine-rich proteins direct the splicing machinery by binding exonic splice enhancers (ESEs) in immature mRNA. As coding exons cannot contain stop codons in one reading frame, stop codons should be rare within ESEs. We confirm that the stop codon density (SCD) in ESE motifs is low, even accounting for nucleotide biases. Given that serine/arginine-rich proteins binding ESEs also facilitate lincRNA splicing, a low SCD could transfer to lincRNAs. As predicted, multiexon lincRNA exons are depleted in stop codons, a result not explained by open reading frame (ORF) contamination. Consistent with transfer selection, stop codon depletion in lincRNAs is most acute in exonic regions with the highest ESE density, disappears when ESEs are masked, is consistent with stop codon usage skews in ESEs, and is diminished in both single-exon lincRNAs and introns. Owing to low SCD, the maximum lengths of pseudo-ORFs frequently exceed null expectations. This has implications for ORF annotation and the evolution of de novo protein-coding genes from lincRNAs. We conclude that not all constraints operating on genes need be explained by the functioning of the gene but may instead be transferred owing to shared binding factors.

TIMP1 intron 3 retention is a marker of colon cancer progression controlled by hnRNPA1

We previously reported a 40-transcripts signature marking the normal mucosa to colorectal adenocarcinoma transition. Eight of these mRNAs also showed splicing alterations, including a specific intron 3 retention in tissue metalloprotease inhibitor I (TIMP1), which decreased during the early steps of colorectal cancer progression. To decipher the mechanism of intron 3 retention/splicing, we first searched for putative RNA binding protein binding sites onto the TIMP1 sequence. We identified potential serine arginine rich splicing factor 1 (SRSF1) and heterogeneous nuclear RiboNucleoProtein A1 (hnRNPA1) binding sites at the end of intron 3 and the beginning of exon 4, respectively. RNA immunoprecipitation showed that hnRNPA1, but not SRSF1 could bind to the corresponding region in TIMP1 pre-mRNA in live cells. Furthermore, using a TIMP1-based ex vivo minigene approach, together with a plasmon resonance in vitro RNA binding assay, we confirmed that hnRNPA1 could indeed bind to wild type TIMP1 exon 4 pre-mRNA and control TMP1 intron 3 splicing, the interaction being abolished in presence of a mutant sequence that disrupted this site. These results indicated that hnRNPA1, upon binding to TIMP1 exon 4, was a positive regulator of intron 3 splicing. We propose that this TIMP1-intron 3 + transcript belongs to the class of nuclear transcripts with "detained" introns, an abundant molecular class, including in cancer.

Expanding the phenotype and the genotype of Stromme syndrome: A novel variant of the CENPF gene and literature review

This report describes siblings with Stromme syndrome, a rare genetic condition that primarily presents with a triad of intestinal atresia, cranial and ocular malformations, and other organ systems could be involved. This clinical triad was initially named after the first person to describe it in 1993. Here, we report a family with two siblings who presented with unusual intestinal atresia and ocular and CNS abnormalities. The first patient is a 6-year-old-boy with apple peel duodeno-jejunal atresia, unilateral microphthalmia and microcephaly. The second patient, a younger brother, presented with intestinal atresia, corneal opacity and alobar holoprosencephaly and passed away at the age of 3 months. Exome sequencing showed a novel homozygous variant in the CENPF gene, NM_016343.3: c.1195-2 A > G that was detected in both of the affected siblings. This is a report and literature review of CENPF-related ciliopathy, which may result in Stromme syndrome. As this is the fourth report linking the CENPF gene variant with Stromme syndrome and first reported case presented with holoprosencephaly, it will expand the current knowledge on the genotype and the phenotype of Stromme syndrome.

Dido3-dependent SFPQ recruitment maintains efficiency in mammalian alternative splicing

Alternative splicing is facilitated by accessory proteins that guide spliceosome subunits to the primary transcript. Many of these splicing factors recognize the RNA polymerase II tail, but SFPQ is a notable exception even though essential for mammalian RNA processing. This study reveals a novel role for Dido3, one of three Dido gene products, in alternative splicing. Binding of the Dido3 amino terminus to histones and to the polymerase jaw domain was previously reported, and here we show interaction between its carboxy terminus and SFPQ. We generated a mutant that eliminates Dido3 but preserves other Dido gene products, mimicking reduced Dido3 levels in myeloid neoplasms. Dido mutation suppressed SFPQ binding to RNA and increased skipping for a large group of exons. Exons bearing recognition sequences for alternative splicing factors were nonetheless included more efficiently. Reduced SFPQ recruitment may thus account for increased skipping of SFPQ-dependent exons, but could also generate a splicing factor surplus that becomes available to competing splice sites. Taken together, our data indicate that Dido3 is an adaptor that controls SFPQ utilization in RNA splicing. Distributing splicing factor recruitment over parallel pathways provides mammals with a simple mechanism to regulate exon usage while maintaining RNA splicing efficiency.

The Splicing Code Goes Deep

The importance of genomic sequence context in generating transcriptome diversity through RNA splicing is independently unmasked by two studies in this issue (Jaganathan et al., 2019; Baeza-Centurion et al., 2019).

Massive computational identification of somatic variants in exonic splicing enhancers using The Cancer Genome Atlas

Owing to the development of next-generation sequencing (NGS) technologies, a large number of somatic variants have been identified in various types of cancer. However, the functional significance of most somatic variants remains unknown. Somatic variants that occur in exonic splicing enhancer (ESE) regions are thought to prevent serine and arginine-rich (SR) proteins from binding to ESE sequence motifs, which leads to exon skipping. We computationally identified somatic variants in ESEs by compiling numerous open-access datasets from The Cancer Genome Atlas (TCGA). Using somatic variants and RNA-seq data from 9635 patients across 32 TCGA projects, we identified 646 ESE-disrupting variants. The false positive rate of our method, estimated using a permutation test, was approximately 1%. Of these ESE-disrupting variants, approximately 71% were located in the binding motifs of four classical SR proteins. ESE-disrupting variants occurred in proportion to the number of somatic variants, but not necessarily in the specific genes associated with the biological processes of cancer. Existing bioinformatics tools could not predict the pathogenicity of ESE-disrupting variants identified in this study, although these variants could cause exon skipping. We demonstrated that ESE-disrupting nonsense variants tended to escape nonsense-mediated decay surveillance. Using integrated analyses of open access data, we could specifically identify ESE-disrupting variants. We have generated a powerful tool, which can handle datasets without normal samples or raw data, and thus contribute to reducing variants of uncertain significance because our statistical approach only uses the exon-junction read counts from the tumor samples.

Activation of cryptic splice sites in three patients with chronic granulomatous disease

Background

Chronic granulomatous disease (CGD) is a primary immune deficiency caused by mutations in the genes encoding the structural components of the phagocyte NADPH oxidase. As a result, the patients cannot generate sufficient amounts of reactive oxygen species required for killing pathogenic microorganisms.

Methods

We analyzed NADPH oxidase activity and component expression in neutrophils, performed genomic DNA and cDNA analysis, and used mRNA splicing prediction tools to evaluate the impact of mutations.

Results

In two patients with CGD, we had previously found mutations that cause aberrant pre‐mRNA splicing. In one patient an exonic mutation in a cryptic donor splice site caused the deletion of the 3' part of exon 6 from the mRNA of CYBB. This patient suffers from X‐linked CGD. The second patient, with autosomal CGD, has a mutation in the donor splice site of intron 1 of CYBA that activates a cryptic donor splice site downstream in intron 1, causing the insertion of intronic sequences in the mRNA. The third patient, recently analyzed, also with autosomal CGD, has a mutation in intron 4 of CYBA, 15 bp from the acceptor splice site. This mutation weakens a branch site and activates a cryptic acceptor splice site, causing the insertion of 14 intronic nucleotides into the mRNA.

Conclusion

We found three different mutations, one exonic, one in a donor splice site and one intronic, that all caused missplicing of pre‐mRNA. We analyzed these mutations with four different splice prediction programs and found that predictions of splice site strength, splice enhancer and splice silencer protein binding and branch site strength are all essential for correct prediction of pre‐mRNA splicing.

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.

More than a messenger: Alternative splicing as a therapeutic target

Alternative splicing of pre-mRNA is an essential post- and co-transcriptional mechanism of gene expression regulation that produces multiple mature mRNA transcripts from a single gene. Genetic mutations that affect splicing underlie numerous devastating diseases. The complexity of splicing regulation allows for multiple therapeutic approaches to correct disease-associated mis-splicing events. In this review, we first highlight recent findings from therapeutic strategies that have used splice switching antisense oligonucleotides and small molecules that bind directly to RNA. Second, we summarize different genetic and chemical approaches to target components of the spliceosome to correct splicing defects in pathological conditions. Finally, we present an overview of compounds that target kinases and accessory pathways that intersect with the splicing machinery. Advancements in the understanding of disease-specific defects caused by mis-regulation of alternative splicing will certainly increase the development of therapeutic options for the clinic. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Impact of HOTAIR variants on preeclampsia susceptibility based on blood and placenta and in silico analysis

HOX transcript antisense RNA (HOTAIR) as a lncRNA involves in epigenetic regulation of various genes. Several studies have been suggested the effects of HOTAIR polymorphisms on different diseases. The aim of the present study was to evaluate the effect of maternal and placental HOTAIR polymorphisms on risk of preeclampsia (PE). The maternal blood of 203 preeclamptic and 202 nonpreeclamptic pregnant women as well as the placentas of 87 of preeclamptic and 95 nonpreeclamptic pregnant women were genotyped for HOTAIR polymorphisms. There was no association between maternal and placental HOTAIR polymorphisms (rs12826786, rs920778, and rs1899663) and PE risk. However, the maternal rs4759314AG and dominant model genotypes were associated with increased risk of PE. The maternal and placental HOTAIR rs10783618 polymorphism was associated with PE risk in recessive and allelic models. Haplotype analysis showed that, the maternal CTGAT and CCTAT and placental CTGAT haplotypes were significantly higher and maternal CTGAC, TCTAT, and TTGAT and placental CTGAC haplotypes were significantly lower in PE women. In silico analysis revealed that HOTAIR rs1899663 had a main effect on the secondary structure of mRNA, however, HOTAIR rs4759314 variant had potential alteration of splicing. In conclusion, the maternal and placental HOTAIR rs10783618 polymorphism might increase PE susceptibility. © 2019 IUBMB Life, 71(9):1367-1381, 2019.

Interplay between coding and exonic splicing regulatory sequences

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Pathogenicity and selective constraint on variation near splice sites

Mutations that perturb normal pre-mRNA splicing are significant contributors to human disease. We used exome sequencing data from 7833 probands with developmental disorders (DDs) and their unaffected parents, as well as more than 60,000 aggregated exomes from the Exome Aggregation Consortium, to investigate selection around the splice sites and quantify the contribution of splicing mutations to DDs. Patterns of purifying selection, a deficit of variants in highly constrained genes in healthy subjects, and excess de novo mutations in patients highlighted particular positions within and around the consensus splice site of greater functional relevance. By using mutational burden analyses in this large cohort of proband–parent trios, we could estimate in an unbiased manner the relative contributions of mutations at canonical dinucleotides (73%) and flanking noncanonical positions (27%), and calculate the positive predictive value of pathogenicity for different classes of mutations. We identified 18 patients with likely diagnostic de novo mutations in dominant DD-associated genes at noncanonical positions in splice sites. We estimate 35%–40% of pathogenic variants in noncanonical splice site positions are missing from public databases.

Implication of GATA4 synonymous variants in congenital heart disease: A comprehensive in-silico approach

Synonymous variations, previously considered as neutral, are recently shown to have a significant impact on mRNA structure and stability thereby affecting protein expression and function. Their role in disease pathogenesis is also emerging. GATA4 is an important transcription factor involved in cardiac development and a well-known candidate gene associated with congenital heart disease (CHD). In the present study, we sought to conduct molecular screening of GATA4 gene in 285 sporadic and non-syndromic CHD cases. We identified four synonymous (c.27C>A, c.822C>T, c.1233G>A and c.1263C>T) and two intronic variants (g.83217T>G & g.85012T>A) in GATA4. Extensive computational analysis using widely acceptable tools i.e., Mfold, Human Splicing Finder (HSF) and Codon Usage bias was performed with a view to understand their putative downstream effects on GATA4 function. Mfold, a mRNA structure prediction tool showed the alterations of the mRNA structure and stability due to synonymous variants. Similarly, HSF also confidently predicted effect on the cis-acting regulatory elements of splicing due to four synonymous and one donor site intronic variants. Additionally, a significant change in 'Relative Synonymous Codon usage (RSCU) frequencies' and 'log ratio of codon usage frequencies' of variant codon was also noted that might affect the rate of translation. This study establishes that the synonymous variants are possibly associated with disease phenotype in CHD patients. Comprehensive computational analysis, using well-established web based tools, is suggestive of their potential downstream molecular effects on the structure, stability and expression of GATA4 protein.

Human astroviruses: in silico analysis of the untranslated region and putative binding sites of cellular proteins

Human astrovirus (HAstV) constitutes a major cause of acute gastroenteritis in children. The viral 5' and 3' untranslated regions (UTR) have been involved in the regulation of several molecular mechanisms. However, in astrovirues have been less characterized. Here, we analyzed the secondary structures of the 5' and 3' UTR of HAstV, as well as their putative target sites that might be recognized by cellular factors. To our knowledge, this is the first bioinformatic analysis that predicts the HAstV 5' UTR secondary structure. The analysis showed that both the UTR sequence and secondary structure are highly conserved in all HAstVs analyzed, suggesting their regulatory role of viral activities. Notably, the UTRs of HAstVs contain putative binding sites for the serine/arginine-rich factors SRSF2, SRSF5, SRSF6, SRSF3, and the multifunctional hnRNPE2 protein. More importantly, putative binding sites for PTB were localized in single-stranded RNA sequences, while hnRNPE2 sites were localized in double-stranded sequence of the HAstV 5' and 3' UTR structures. These analyses suggest that the combination of SRSF proteins, hnRNPE2 and PTB described here could be involved in the maintenance of the secondary structure of the HAstVs, possibly allowing the recruitment of the replication complex that selects and recruits viral RNA replication templates.

Modeling human RNA spliceosome mutations in the mouse: not all mice were created equal

Myelodysplastic syndromes (MDS) and related myelodysplastic/myeloproliferative neoplasms (MDS/MPNs) are clonal stem cell disorders, primarily affecting patients over 65 years of age. Mapping of the MDS and MDS/MPN genome identified recurrent heterozygous mutations in the RNA splicing machinery, with the SF3B1, SRSF2, and U2AF1 genes being frequently mutated. To better understand how spliceosomal mutations contribute to MDS pathogenesis in vivo, numerous groups have sought to establish conditional murine models of SF3B1, SRSF2, and U2AF1 mutations. The high degree of conservation of hematopoiesis between mice and human and the well-established phenotyping and genetic modification approaches make murine models an effective tool with which to study how a gene mutation contributes to disease pathogenesis. The murine models of spliceosomal mutations described to date recapitulate human MDS or MDS/MPN to varying extents. Reasons for the differences in phenotypes reported between alleles of the same mutation are varied, but the nature of the genetic modification itself and subsequent analysis methods are important to consider. In this review, we summarize recently reported murine models of SF3B1, SRSF2, and U2AF1 mutations, with a particular focus on the genetically engineered modifications underlying the models and the experimental approaches applied.

Exonic splice regulation imposes strong selection at synonymous sites

What proportion of coding sequence nucleotides have roles in splicing, and how strong is the selection that maintains them? Despite a large body of research into exonic splice regulatory signals, these questions have not been answered. This is because, to our knowledge, previous investigations have not explicitly disentangled the frequency of splice regulatory elements from the strength of the evolutionary constraint under which they evolve. Current data are consistent both with a scenario of weak and diffuse constraint, enveloping large swaths of sequence, as well as with well-defined pockets of strong purifying selection. In the former case, natural selection on exonic splice enhancers (ESEs) might primarily act as a slight modifier of codon usage bias. In the latter, mutations that disrupt ESEs are likely to have large fitness and, potentially, clinical effects. To distinguish between these scenarios, we used several different methods to determine the distribution of selection coefficients for new mutations within ESEs. The analyses converged to suggest that ∼15%-20% of fourfold degenerate sites are part of functional ESEs. Most of these sites are under strong evolutionary constraint. Therefore, exonic splice regulation does not simply impose a weak bias that gently nudges coding sequence evolution in a particular direction. Rather, the selection to preserve these motifs is a strong force that severely constrains the evolution of a substantial proportion of coding nucleotides. Thus synonymous mutations that disrupt ESEs should be considered as a potentially common cause of single-locus genetic disorders.

Association of DFNA5, SYK, and NELL1 variants along with HPV infection in oral cancer among the prolonged tobacco-chewers

Southeast Asia, especially India, is well known for the highest use of smokeless tobacco. These products are known to induce oral squamous cell carcinoma. However, not all long-term tobacco-chewers develop oral squamous cell carcinoma. In addition, germline variants play a crucial role in susceptibility, prognosis, development, and progression of the disease. These prompted us to study the genetic susceptibility to oral squamous cell carcinoma among the long-term tobacco-chewers. Here, we presented a retrospective study on prolonged tobacco-chewers of Northeast India to identify the potential protective or risk-associated germline variants in tobacco-related oral squamous cell carcinoma along with HPV infection. Targeted re-sequencing (n = 60) of 170 genetic regions from 75 genes was carried out in Ion-PGM™ and validation (n = 116) of the observed variants was done using Sequenom iPLEX MassARRAY™ platform followed by polymerase chain reaction-based HPV genotyping and p16-immunohistochemistry study. Subsequently, estimation of population structure, different statistical and in silico approaches were undertaken. We identified one nonsense-mediated mRNA decay transcript variant in the DFNA5 region (rs2237306), associated with Benzo(a)pyrene, as a protective factor (odds ratio = 0.33; p = 0.009) and four harmful (odds ratio > 2.5; p < 0.05) intronic variants, rs182361, rs290974, and rs169724 in SYK and rs1670661 in NELL1 region, involved in genetic susceptibility to tobacco- and HPV-mediated oral oncogenesis. Among the oral squamous cell carcinoma patients, 12.6% (11/87) were HPV positive, out of which 45.5% (5/11) were HPV16-infected, 27.3% (3/11) were HPV18-infected, and 27.3% (3/11) had an infection of both subtypes. Multifactor dimensionality reduction analysis showed that the interactions among HPV and NELL1 variant rs1670661 with age and gender augmented the risk of both non-tobacco- and tobacco-related oral squamous cell carcinoma, respectively. These suggest that HPV infection may be one of the important risk factors for oral squamous cell carcinoma in this population. Finally, we newly report a DFNA5 variant probably conferring protection via nonsense-mediated mRNA decay pathway against tobacco-related oral squamous cell carcinoma. Thus, the analytical approach used here can be useful in predicting the population-specific significant variants associated with oral squamous cell carcinoma in any heterogeneous population.

Exon skipping in CYBB mRNA and skewed inactivation of X chromosome cause late-onset chronic granulomatous disease

Chronic granulomatous disease (CGD) is a hereditary immunodeficiency syndrome caused by a defect in the NADPH oxidase complex, which is essential for bactericidal function of phagocytes. Approximately 70% of patients with CGD have a mutation in the CYBB gene on the X chromosome, resulting in defective expression of gp91^phox, one of the membrane-bound subunits of NADPH oxidase. Although most patients with X-linked CGD are males, owing to transmission of this disease as an X-linked recessive trait, there are female patients with X-linked CGD. Here, we report the case of a teenage girl with X-linked CGD associated with a heterozygous mutation in exon 5 of the CYBB gene (c.389G > C; R130P), which causes skipping of exon 5, resulting in a premature stop codon in exon 6 of CYBB. Accurate pro-mRNA splicing for mature mRNA formation is regulated by several splicing mechanisms that are essential for appropriate recognition of exonic sequences. The c.389G > C mutation disrupts exonic-splicing regulator sequences, thereby resulting in the aberrant skipping of exon 5 in the CYBB transcript of the patient. The patient showed an extremely skewed (≥96%) X inactivation pattern of the HUMARA locus; this inactivation is thought to be responsible for the development of CGD not only in neutrophils but also in monocytic, T-cell, and B-cell lineages and in CD34-positive immature hematopoietic cells. Our case and other reports indicate that the onset of X-linked CGD in female patients tends to occur later in life, and that the symptoms tend to be milder as compared to male patients.

The Entamoeba histolytica Syf1 Homolog Is Involved in the Splicing of AG-Dependent and AG-Independent Transcripts

Syf1 is a tetratricopeptide repeat (TPR) protein implicated in transcription elongation, spliceosome conformation, mRNA nuclear-cytoplasmic export and transcription-coupled DNA repair. Recently, we identified the spliceosomal components of the human parasite Entamoeba histolytica, among them is EhSyf. Molecular predictions confirmed that EhSyf contains 15 type 1 TPR tandem α-antiparallel array motifs. Amoeba transformants carrying plasmids overexpressing HA-tagged or EhSyf silencing plasmids were established to monitor the impact of EhSyf on the splicing of several test Entamoeba transcripts. EhSyf Entamoeba transformants efficiently silenced or overexpressed the proteins in the nucleus. The overexpression or absence of EhSyf notably enhanced or blocked splicing of transcripts irrespective of the strength of their 3′ splice site. Finally, the absence of EhSyf negatively affected the transcription of an intron-less transcript. Altogether our data suggest that EhSyf is a bona fide Syf1 ortholog involved in transcription and splicing.

Genetic variants of ALDH2-rs671 and CYP2E1-rs2031920 contributed to risk of hepatocellular carcinoma susceptibility in a Chinese population

Objective

Acetaldehyde dehydrogenase 2 (ALDH2) and cytochrome P450 2E1 (CYP2E1) have been associated with hepatocellular carcinoma (HCC) susceptibility and prognosis. The polymorphisms ALDH2 rs671 and CYP2E1 rs2031920 are reportedly correlated with the prevalence of HCC in other countries. The aim of this study was to investigate associations between ALDH2 and CYP2E1, and HCC susceptibility in a population of Guangxi, southern China, an area with a high incidence of HCC.

Patients and methods

The study cohort included 300 HCC cases, 292 healthy controls for HCC susceptibility analysis, and another 20 HCC cases and 10 healthy controls for ascertainment. Genotyping was performed using the polymerase chain reaction-restriction fragment length polymorphism method.

Results

The study results demonstrated that mutant genotypes of ALDH2 (G/A and A/A) led to significant differences in HCC susceptibility, as compared with the wild genotype (G/G) with the same C1/C1 genotype in non-drinking individuals (adjusted P=0.010, OR=0.20, 95% CI=0.06–0.68). The mutant genotypes of CYP2E1 (C1/C2 and C2/C2) brought about significant differences in HCC susceptibility, as compared with the wild genotype (C1/C1) and the same G/G genotype (adjusted P=0.025, OR=0.42, 95% CI=0.20–0.90). Drinking plays a role in HCC susceptibility in the same G/G genotype individuals (adjusted P=0.004, OR=0.32, 95% CI=0.15–0.69), but had no impact when combined with CYP2E1 for analysis (all P>0.05).

Conclusion

These results suggest that the mutant genotypes of ALDH2 and CYP2E1 may be protective factors for HCC susceptibility in Guangxi province, China.

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

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

Ancient duons may underpin spatial patterning of gene expression in C4 leaves

If the highly efficient C4 photosynthesis pathway could be transferred to crops with the C3 pathway there could be yield gains of up to 50%. It has been proposed that the multiple metabolic and developmental modifications associated with C4 photosynthesis are underpinned by relatively few master regulators that have allowed the evolution of C4 photosynthesis more than 60 times in flowering plants. Here we identify a component of one such regulator that consists of a pair of cis-elements located in coding sequence of multiple genes that are preferentially expressed in bundle sheath cells of C4 leaves. These motifs represent duons as they play a dual role in coding for amino acids as well as controlling the spatial patterning of gene expression associated with the C4 leaf. They act to repress transcription of C4 photosynthesis genes in mesophyll cells. These duons are also present in the C3 model Arabidopsis thaliana, and, in fact, are conserved in all land plants and even some algae that use C3 photosynthesis. C4 photosynthesis therefore appears to have coopted an ancient regulatory code to generate the spatial patterning of gene expression that is a hallmark of C4 photosynthesis. This intragenic transcriptional regulatory sequence could be exploited in the engineering of efficient photosynthesis of crops.

Saturation mutagenesis reveals manifold determinants of exon definition

To illuminate the extent and roles of exonic sequences in the splicing of human RNA transcripts, we conducted saturation mutagenesis of a 51-nt internal exon in a three-exon minigene. All possible single and tandem dinucleotide substitutions were surveyed. Using high-throughput genetics, 5560 minigene molecules were assayed for splicing in human HEK293 cells. Up to 70% of mutations produced substantial (greater than twofold) phenotypes of either increased or decreased splicing. Of all predicted secondary structural elements, only a single 15-nt stem–loop showed a strong correlation with splicing, acting negatively. The in vitro formation of exon-protein complexes between the mutant molecules and proteins associated with spliceosome formation (U2AF35, U2AF65, U1A, and U1-70K) correlated with splicing efficiencies, suggesting exon definition as the step affected by most mutations. The measured relative binding affinities of dozens of human RNA binding protein domains as reported in the CISBP-RNA database were found to correlate either positively or negatively with splicing efficiency, more than could fit on the 51-nt test exon simultaneously. The large number of these functional protein binding correlations point to a dynamic and heterogeneous population of pre-mRNA molecules, each responding to a particular collection of binding proteins.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

Mutation in an exonic splicing enhancer site causing chronic granulomatous disease

In a male patient suffering from X-linked chronic granulomatous disease (CGD) we found a c.389G>T mutation in exon 5 of the CYBB gene. We have analyzed why 95% of the transcripts of this gene lacked exon 5, leading to a frameshift and premature termination codon. The mutation was located in a region comprising three putative exonic splicing enhancer binding sites, for SRSF1, SRFS2 and SRFS6, according to the ESEfinder Tool (http://rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi). With the Analyser Splice Tool we calculated the probability of skipping of exon 5 in CYBB mRNA, and by means of Sroogle the number of putative binding motifs for splicing enhancer and splicing silencer proteins (http://astlab.tau.ac.il/index.php). These analyses clarify why this exon was skipped in the majority of the mRNA. The normally spliced transcript contains an amino acid change p.Arg130Leu. This poorly expressed transcript gives rise to a protein with low expression but presumably normal activity, leading to a respiratory burst activity in the patient's neutrophils of about 15% of normal.

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

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