Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Analysis of the role of Caenorhabditis elegans GC-AG introns in regulated splicing

GC-AG introns represent 0.7% of total human pre-mRNA introns. To study the function of GC-AG introns in splicing regulation, 196 cDNA-confirmed GC-AG introns were identified in Caenorhabditis elegans. These represent 0.6% of the cDNA- confirmed intron data set for this organism. Eleven of these GC-AG introns are involved in alternative splicing. In a comparison of the genomic sequences of homologous genes between C.elegans and Caenorhabditis briggsae for 26 GC-AG introns, the C at the +2 position is conserved in only five of these introns. A system to experimentally test the function of GC-AG introns in alternative splicing was developed. Results from these experiments indicate that the conserved C at the +2 position of the tenth intron of the let-2 gene is essential for developmentally regulated alternative splicing. This C allows the splice donor to function as a very weak splice site that works in balance with an alternative GT splice donor. A weak GT splice donor can functionally replace the GC splice donor and allow for splicing regulation. These results indicate that while the majority of GC-AG introns appear to be constitutively spliced and have no evolutionary constraints to prevent them from being GT-AG introns, a subset of GC-AG introns is involved in alternative splicing and the C at the +2 position of these introns can have an important role in splicing regulation.     

In vitro and in silico analysis reveals an efficient algorithm to predict the splicing consequences of mutations at the 5' splice sites

We have found that two previously reported exonic mutations in the PINK1 and PARK7 genes affect pre-mRNA splicing. To develop an algorithm to predict underestimated splicing consequences of exonic mutations at the 5′ splice site, we constructed and analyzed 31 minigenes carrying exonic splicing mutations and their derivatives. We also examined 189 249 U2-dependent 5′ splice sites of the entire human genome and found that a new variable, the SD-Score, which represents a common logarithm of the frequency of a specific 5′ splice site, efficiently predicts the splicing consequences of these minigenes. We also employed the information contents (R_i) to improve the prediction accuracy. We validated our algorithm by analyzing 32 additional minigenes as well as 179 previously reported splicing mutations. The SD-Score algorithm predicted aberrant splicings in 198 of 204 sites (sensitivity = 97.1%) and normal splicings in 36 of 38 sites (specificity = 94.7%). Simulation of all possible exonic mutations at positions −3, −2 and −1 of the 189 249 sites predicts that 37.8, 88.8 and 96.8% of these mutations would affect pre-mRNA splicing, respectively. We propose that the SD-Score algorithm is a practical tool to predict splicing consequences of mutations affecting the 5′ splice site.     

Regulation of Sex-Specific Selection of fruitless 5′ Splice Sites by transformer and transformer-2

In Drosophila melanogaster, the fruitless (fru) gene controls essentially all aspects of male courtship behavior. It does this through sex-specific alternative splicing of the fru pre-mRNA, leading to the production of male-specific fru mRNAs capable of expressing male-specific fru proteins. Sex-specific fru splicing involves the choice between alternative 5′ splice sites, one used exclusively in males and the other used only in females. Here we report that the Drosophila sex determination genes transformer (tra) and transformer-2 (tra-2) switch fru splicing from the male-specific pattern to the female-specific pattern through activation of the female-specific fru 5′ splice site. Activation of female-specific fru splicing requires cis-acting tra and tra-2 repeat elements that are part of an exonic splicing enhancer located immediately upstream of the female-specific fru 5′ splice site and are recognized by the TRA and TRA-2 proteins in vitro. This fru splicing enhancer is sufficient to promote the activation by tra and tra-2 of both a 5′ splice site and the female-specific doublesex (dsx) 3′ splice site, suggesting that the mechanisms of 5′ splice site activation and 3′ splice site activation may be similar.     

Delineation of the mechanisms of aberrant splicing caused by two unusual intronic mutations in the RSK2 gene involved in Coffin-Lowry syndrome

Coffin–Lowry syndrome (CLS) is caused by mutations in the RSK2 gene encoding a protein kinase of the Ras signalling pathway. We have studied two point mutations which cause aberrant splicing but do not concern the invariant GT or AG nucleotides of splice sites. The first, an A→G transition at position +3 of the 5′ splice site of exon 6, results in vivo and in vitro in exon skipping and premature translation termination. The natural 5′ splice site, although intrinsically weak, is not transactivated under normal conditions. Consequently, replacement of an A/U by a G/U base pairing with U1 snRNA reduces its strength below a critical threshold. The second mutation, an A→G transition 11 nt upstream of exon 5, creates a new AG near the natural 3′ splice site. In vitro this synthetic 3′ AG is used exclusively by the splicing machinery. In vivo this splicing event is also observed, but is underestimated because the resulting RSK2 mRNA contains premature stop codons which trigger the nonsense-mediated decay process. We show that a particular mechanism is involved in the aberrant splicing of exon 5, implying involvement of the natural 3′ AG during the first catalytic step and the new 3′ AG during the second step. Thus, our results explain how these mutations cause severe forms of CLS.     

Application of Gene Network Analysis Techniques Identifies AXIN1/PDIA2 and Endoglin Haplotypes Associated with Bicuspid Aortic Valve

Bicuspid Aortic Valve (BAV) is a highly heritable congenital heart defect. The low frequency of BAV (1% of general population) limits our ability to perform genome-wide association studies. We present the application of four a priori SNP selection techniques, reducing the multiple-testing penalty by restricting analysis to SNPs relevant to BAV in a genome-wide SNP dataset from a cohort of 68 BAV probands and 830 control subjects. Two knowledge-based approaches, CANDID and STRING, were used to systematically identify BAV genes, and their SNPs, from the published literature, microarray expression studies and a genome scan. We additionally tested Functionally Interpolating SNPs (fitSNPs) present on the array; the fourth consisted of SNPs selected by Random Forests, a machine learning approach. These approaches reduced the multiple testing penalty by lowering the fraction of the genome probed to 0.19% of the total, while increasing the likelihood of studying SNPs within relevant BAV genes and pathways. Three loci were identified by CANDID, STRING, and fitSNPS. A haplotype within the AXIN1-PDIA2 locus (p-value of 2.926×10⁻⁰⁶) and a haplotype within the Endoglin gene (p-value of 5.881×10⁻⁰⁴) were found to be strongly associated with BAV. The Random Forests approach identified a SNP on chromosome 3 in association with BAV (p-value 5.061×10⁻⁰⁶). The results presented here support an important role for genetic variants in BAV and provide support for additional studies in well-powered cohorts. Further, these studies demonstrate that leveraging existing expression and genomic data in the context of GWAS studies can identify biologically relevant genes and pathways associated with a congenital heart defect.     

Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer

Background

The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)_nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.

Results

Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.

Conclusion

We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1143) contains supplementary material, which is available to authorized users.     

Analysis of the splicing machinery in fission yeast: a comparison with budding yeast and mammals

Based on genetic and bioinformatic analysis, 80 proteins from the newly sequenced Schizosaccharomyces pombe genome appear to be splicing factors. The fission yeast splicing factors were compared to those of Homo sapiens and Saccharomyces cerevisiae in order to determine the extent of conservation or divergence that has occurred over the billion years of evolution that separate these organisms. Our results indicate that many of the factors present in all three organisms have been well conserved throughout evolution. It is calculated that 38% of the fission yeast splicing factors are more similar to the human proteins than to the budding yeast proteins (>10% more similar or similar over a greater region). Many of the factors in this category are required for recognition of the 3′ splice site. Ten fission yeast splicing factors, including putative regulatory factors, have human homologs, but no apparent budding yeast homologs based on sequence data alone. Many of the budding yeast factors that are absent in fission yeast are associated with the U1 and U4/U6.U5 snRNP. Collectively the data presented in this survey indicate that of the two yeasts, S.pombe contains a splicing machinery more closely reflecting the archetype of a spliceosome.     

Function of a Bovine Papillomavirus Type 1 Exonic Splicing Suppressor Requires a Suboptimal Upstream 3′ Splice Site

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. Previous studies in our laboratory have identified two purine-rich exonic splicing enhancers (ESEs), SE1 and SE2, located between two alternative 3′ splice sites at nucleotide (nt) 3225 and nt 3605. Further analysis of BPV-1 late-pre-mRNA splicing in vitro revealed a 48-nt pyrimidine-rich region immediately downstream of SE1 that inhibits utilization of the nt 3225 3′ splice site. This inhibitory element, which we named an exonic splicing suppressor (ESS), has a U-rich 5′ end, a C-rich central part, and an AG-rich 3′ end (Z. M. Zheng, P. He, and C. C. Baker, J. Virol. 70:4691–4699, 1996). The present study utilized in vitro splicing of both homologous and heterologous pre-mRNAs to further characterize the ESS. The BPV-1 ESS was inserted downstream of the 3′ splice site in the BPV-1 late pre-mRNA, Rous sarcoma virus src pre-mRNA, human immunodeficiency virus tat-rev pre-mRNA, and Drosophila dsx pre-mRNA, all containing a suboptimal 3′ splice site, and in the human β-globin pre-mRNA, which contains a constitutive 3′ splice site. These studies demonstrated that suppression of splicing by the BPV-1 ESS requires an upstream suboptimal 3′ splice site but not an upstream ESE. Furthermore, the ESS functions when located either upstream or downstream of BPV-1 SE1. Mutational analyses demonstrated that the function of the ESS is sequence dependent and that only the C-rich region of the ESS is essential for suppression of splicing in all the pre-mRNAs tested.