Characteristics and regulatory elements defining constitutive splicing and different modes of alternative splicing in human and mouse

Inverted Alu repeats in loop-out exon skipping across hominoid evolution

Background

Changes in RNA splicing over the course of evolution have profoundly diversified the functional landscape of the human genome. While DNA sequences proximal to intron-exon junctions are known to be critical for RNA splicing, the impact of distal intronic sequences remains underexplored. Emerging evidence suggests that inverted pairs of intronic Alu elements can promote exon skipping by forming RNA stem-loop structures. However, their prevalence and influence throughout evolution remain unknown.

Results

Here, we present a systematic analysis of inverted Alu pairs across the human genome to assess their impact on exon skipping through predicted RNA stem-loop formation and their relevance to hominoid evolution. We found that inverted Alu pairs, particularly pairs of AluY-AluSx1 and AluSz-AluSx, are enriched in the flanking regions of skippable exons genome-wide and are predicted to form stable stem-loop structures. Exons defined by weak 3' acceptor and strong 5' donor splice sites appear especially prone to this skipping mechanism. Through comparative genome analysis across nine primate species, we identified 67,126 hominoid-specific Alu insertions, primarily from AluY and AluS subfamilies, which form inverted pairs enriched across skippable exons in genes of ubiquitination-related pathways. Experimental validation of exon skipping among several hominoid-specific inverted Alu pairs further reinforced their potential evolutionary significance.

Conclusion

This work extends our current knowledge of the roles of RNA secondary structure formed by inverted Alu pairs and details a newly emerging mechanism through which transposable elements have contributed to genomic innovation across hominoid evolution at the transcriptomic level.

ASpdb: an integrative knowledgebase of human protein isoforms from experimental and AI-predicted structures

Alternative splicing is a crucial cellular process in eukaryotes, enabling the generation of multiple protein isoforms with diverse functions from a single gene. To better understand the impact of alternative splicing on protein structures, protein-protein interaction and human diseases, we developed ASpdb (https://biodataai.uth.edu/ASpdb/), a comprehensive database integrating experimentally determined structures and AlphaFold 2-predicted models for human protein isoforms. ASpdb includes over 3400 canonical isoforms, each represented by both experimentally resolved and predicted structures, and >7200 alternative isoforms with AlphaFold 2 predictions. In addition to detailed splicing events, 3D structures, sequence variations and functional annotations, ASpdb uniquely offers comparative analyses and visualization of structural alterations among isoforms. This resource is invaluable for advancing research in alternative splicing, structural biology and disease mechanisms.

Multimodal screen identifies noise-regulatory proteins

Gene-expression noise can influence cell-fate choices across pathology and physiology. However, a crucial question persists: do regulatory proteins or pathways exist that control noise independently of mean expression levels? Our integrative approach, combining single-cell RNA sequencing with proteomics and regulator enrichment analysis, identifies 32 putative noise regulators. SON, a nuclear speckle-associated protein, alters transcriptional noise without changing mean expression levels. Furthermore, SON's noise control can propagate to the protein level. Long-read and total RNA sequencing shows that SON's noise control does not significantly change isoform usage or splicing efficiency. Moreover, SON depletion reduces state switching in pluripotent mouse embryonic stem cells and impacts their fate choice during differentiation. Collectively, we demonstrate a class of proteins that control noise orthogonally to mean expression levels. This work serves as a proof of concept that can identify other functional noise regulators throughout development and disease progression.

A CRISPR-dCas13 RNA-editing tool to study alternative splicing

Alternative splicing allows multiple transcripts to be generated from the same gene to diversify the protein repertoire and gain new functions despite a limited coding genome. It can impact a wide spectrum of biological processes, including disease. However, its significance has long been underestimated due to limitations in dissecting the precise role of each splicing isoform in a physiological context. Furthermore, identifying key regulatory elements to correct deleterious splicing isoforms has proven equally challenging, increasing the difficulty of tackling the role of alternative splicing in cell biology. In this work, we take advantage of dCasRx, a catalytically inactive RNA targeting CRISPR-dCas13 ortholog, to efficiently switch alternative splicing patterns of endogenous transcripts without affecting overall gene expression levels cost-effectively. Additionally, we demonstrate a new application for the dCasRx splice-editing system to identify key regulatory RNA elements of specific splicing events. With this approach, we are expanding the RNA toolkit to better understand the regulatory mechanisms underlying alternative splicing and its physiological impact in various biological processes, including pathological conditions.

Comparative Analysis of the pIgR Gene from the Antarctic Teleost Trematomus bernacchii Reveals Distinctive Features of Cold-Adapted Notothenioidei

The IgM and IgT classes were previously identified and characterized in the Antarctic teleost Trematomus bernacchii, a species belonging to the Perciform suborder Notothenoidei. Herein, we characterized the gene encoding the polymeric immunoglobulin receptor (pIgR) in the same species and compared it to the pIgR of multiple teleost species belonging to five perciform suborders, including 11 Antarctic and 1 non-Antarctic (Cottoperca gobio) notothenioid species, the latter living in the less-cold peri-Antarctic sea. Antarctic pIgR genes displayed particularly long introns marked by sites of transposable elements and transcription factors. Furthermore, analysis of T. bernacchii pIgR cDNA unveiled multiple amino acid substitutions unique to the Antarctic species, all introducing adaptive features, including N-glycosylation sequons. Interestingly, C. gobio shared most features with the other perciforms rather than with the cold-adapted relatives. T. bernacchii pIgR transcripts were predominantly expressed in mucosal tissues, as indicated by q-PCR and in situ hybridization analysis. These results suggest that in cold-adapted species, pIgR preserved its fundamental role in mucosal immune defense, although remarkable gene structure modifications occurred.

Comprehensive Update and Revision of Nomenclature on Complement C6 and C7 Variants

Complement genes encompass a wide array of variants, giving rise to numerous protein isoforms that have often been shown to exhibit clinical significance. Given that these variants have been discovered over a span of 50 y, one challenging consequence is the inconsistency in the terminology used to classify them. This issue is prominently evident in the nomenclature used for complement C6 and C7 variants, for which we observed a great discrepancy between previously published works and variants described in current genome browsers. This report discusses the causes for the discrepancies in C6 and C7 nomenclature and seeks to establish a classification system that would unify existing and future variants. The inconsistency in the methods used to annotate amino acids and the modifications pinpointed in the C6 and C7 primers are some of the factors that contribute greatly to the discrepancy in the nomenclature. Several variants that were classified incorrectly are highlighted in this report, and we showcase first-hand how a unified classification system is important to match previous with current genetic information. Ultimately, we hope that the proposed classification system of nomenclature becomes an incentive for studies on complement variants and their physiological and/or pathological effects.

Alternative splicing patterns of doublesex reveal a missing link between Nix and doublesex in the sex determination cascade of Aedes albopictus

Sexual development in insects is regulated by a complicated hierarchical cascade of sex determination. The primary signals are diverse, whereas the central nexus doublesex (dsx) gene is relatively conserved within the pathway. Aedes (Stegomyia) albopictus is an important vector with an extensive worldwide distribution. We previously reported that Ae. albopictus dsx (Aalbdsx) yields one male- (Aalbdsx^M ) and three female-specific isoforms (Aalbdsx^F1-3 ); however, the spatiotemporal expression profiles and mechanisms regulating sex-specific alternative splicing require further investigation. In this study, we demonstrated that the Aalbdsx^M messenger RNA (mRNA) represents the default pattern when analyzed in human foreskin fibroblasts and HeLa cells. We combined reverse transcription polymerase chain reaction with RNA immunoprecipitation using specific antibodies against tagged Ae. albopictus male-determining factor AalNix and confirmed that AalNix indirectly regulates dsx pre-mRNA and regulates its alternative splicing. During the early embryo stage (0-2 and 4-8 h), maternal dsx^F and default splicing dsx^M were detected in both sexes; the expression of dsx^M then decreased until sufficient AalNix transcripts accumulated in male embryos at 20-24 h. These findings suggest that one or more potential dsx splicing enhancers can shift dsx^M to dsx^F in both sexes; however, the presence of Nix influences the function of this unknown splicing enhancer and ultimately leads to the formation of dsx^M in males. Finally, our results provide important insight into the regulatory mechanism of dsx alternative splicing in the mosquito.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

Pathogenic Intronic Splice-Affecting Variants in MYBPC3 in Three Patients with Hypertrophic Cardiomyopathy

Genetic variants in MYBPC3 are one of the most common causes of hypertrophic cardiomyopathy (HCM). While variants in MYBPC3 affecting canonical splice site dinucleotides are a well-characterised cause of HCM, only recently has work begun to investigate the pathogenicity of more deeply intronic variants. Here, we present three patients with HCM and intronic splice-affecting MYBPC3 variants and analyse the impact of variants on splicing using in vitro minigene assays. We show that the three variants, a novel c.927-8G>A variant and the previously reported c.1624+4A>T and c.3815-10T>G variants, result in MYBPC3 splicing errors. Analysis of blood-derived patient RNA for the c.3815-10T>G variant revealed only wild type spliced product, indicating that mis-spliced transcripts from the mutant allele are degraded. These data indicate that the c.927-8G>A variant of uncertain significance and likely benign c.3815-10T>G should be reclassified as likely pathogenic. Furthermore, we find shortcomings in commonly applied bioinformatics strategies to prioritise variants impacting MYBPC3 splicing and re-emphasise the need for functional assessment of variants of uncertain significance in diagnostic testing.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

Dynamic changes in intron retention are tightly associated with regulation of splicing factors and proliferative activity during B-cell development

Intron retention (IR) has been proposed to modulate the delay between transcription and translation. Here, we provide an exhaustive characterization of IR in differentiated white blood cells from both the myeloid and lymphoid lineage where we observed highest levels of IR in monocytes and B-cells, in addition to previously reported granulocytes. During B-cell differentiation, we found an increase in IR from the bone marrow precursors to cells residing in secondary lymphoid organs. B-cells that undergo affinity maturation to become antibody producing plasma cells steadily decrease retention. In general, we found an inverse relationship between global IR levels and both the proliferative state of cells, and the global levels of expression of splicing factors. IR dynamics during B-cell differentiation appear to be conserved between human and mouse, suggesting that IR plays an important biological role, evolutionary conserved, during blood cell differentiation. By correlating the expression of non-core splicing factors with global IR levels, and analyzing RNA binding protein knockdown and eCLIP data, we identify a few splicing factors likely playing an evolutionary conserved role in IR regulation. Our work provides new insights into the role of IR during hematopoiesis, and on the main factors involved in regulating IR.

CRISPR/Cas9-mediated modulation of splicing efficiency reveals short splicing isoform of Xist RNA is sufficient to induce X-chromosome inactivation

Alternative splicing of mRNA precursors results in multiple protein variants from a single gene and is critical for diverse cellular processes and development. Xist encodes a long noncoding RNA which is a central player to induce X-chromosome inactivation in female mammals and has two major splicing variants: long and short isoforms of Xist RNA. Although a differentiation-specific and a female-specific expression of Xist isoforms have been reported, the functional role of each Xist RNA isoform is largely unexplored. Using CRISPR/Cas9-mediated targeted modification of the 5' splice site in Xist intron 7, we create mutant female ES cell lines which dominantly express the long- or short-splicing isoform of Xist RNA from the inactive X-chromosome (Xi) upon differentiation. Successful execution of CRISPR/Cas-based splicing modulation indicates that our CRISPR/Cas-based targeted modification of splicing sites is a useful approach to study specific isoforms of a transcript generated by alternative splicing. Upon differentiation of splicing-mutant Xist female ES cells, we find that both long and short Xist isoforms can induce X-chromosome inactivation normally during ES cell differentiation, suggesting that the short splicing isoform of Xist RNA is sufficient to induce X-chromosome inactivation.

RNA-Seq approach for accurate characterization of splicing efficiency of yeast introns

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N_EE) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N_EI5 and N_EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N_Total = N_EE + p * N_EI5 + (1-p) * N_EI3, with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N_EI5 is expected to be smaller than N_EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N_EE/N_Total). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated.

FAK alternative splice mRNA variants expression pattern in colorectal cancer

The Focal adhesion kinase (FAK) is a ubiquitous cytoplasmic tyrosine‐kinase promoting tumor progression and metastasis processes by acting in cancer cells and their tumor microenvironment partners. FAK overexpression in primary colon tumors and their metastasis is associated to poor colorectal cancer (CRC) patients’ outcome. Eight FAK mRNA alternative splice variants have been described and contribute to additional level of FAK activity regulation, some of them corresponding to overactivated FAK isoforms. To date, FAK mRNA alternative splice variants expression and implication in CRC processes remain unknown. Here, using different human CRC cells lines displaying differential invasive capacities in an in vivo murine model recapitulating the different steps of CRC development from primary tumors to liver and lung metastasis, we identified three out of the eight mRNA variants (namely FAK⁰, FAK²⁸ and FAK⁶) differentially expressed along the CRC process and the tumor sites. Our results highlight an association between FAK⁰ and FAK⁶ expressions and the metastatic potential of the most aggressive cell lines HT29 and HCT116, suggesting that FAK⁰ and FAK⁶ could represent aggressiveness markers in CRC. Our findings also suggest a more specific role for FAK²⁸ in the interactions between the tumors cells and their microenvironment. In conclusion, targeting FAK⁰, the common form of FAK, might not be a good strategy based on the numerous roles of this kinase in physiological processes. In contrast, FAK⁶ or FAK²⁸ splice variants, or their corresponding protein isoforms, may putatively represent future therapeutic target candidates in the development of CRC primary tumors and metastasis.

What's new?

Overexpression of the focal adhesion kinase (FAK) is associated with poor outcome in patients with colorectal cancer but the role of the eight splice variants of FAK remains unknown. Here the authors correlated FAK splice variant expression in colorectal tumor cell lines with invasiveness in mouse models. FAK⁰ and FAK⁶ splice variant expression was associated with higher aggressiveness and metastatic potential, underscoring that distinct FAK splice variants may represent new targets in the development of drugs against colorectal cancer and associated metastasis.

Landscape of alternative splicing in Capra_hircus

Alternative splicing (AS) is a fundamental regulatory process in all higher eukaryotes. However, AS landscapes for a number of animals, including goats, have not been explored to date. Here, we sequenced 60 samples representing 5 tissues from 4 developmental stages in triplicate using RNA-seq to elucidate the goat AS landscape. In total, 14,521 genes underwent AS (AS genes), accounting for 85.53% of intron-containing genes (16,697). Among these AS genes, 6,342 were differentially expressed in different tissues. Of the AS events identified, retained introns were most prevalent (37.04% of total AS events). Functional enrichment analysis of differential and specific AS genes indicated goat AS mainly involved in organ function and development. Particularly, AS genes identified in leg muscle were associated with the “regulation of skeletal muscle tissue development” GO term. Given genes were associated with this term, four of which (NRG4, IP6K3, AMPD1, and DYSF) might play crucial roles in skeletal muscle development. Further investigation indicated these five genes, harbored 13 ASs, spliced exclusively in leg muscle, likely played a role in goat leg muscle development. These results provide novel insights into goat AS landscapes and a valuable resource for investigation of goat transcriptome complexity and gene regulation.

Alternative splicing links histone modifications to stem cell fate decision

BACKGROUND

Understanding the embryonic stem cell (ESC) fate decision between self-renewal and proper differentiation is important for developmental biology and regenerative medicine. Attention has focused on mechanisms involving histone modifications, alternative pre-messenger RNA splicing, and cell-cycle progression. However, their intricate interrelations and joint contributions to ESC fate decision remain unclear.

RESULTS

We analyze the transcriptomes and epigenomes of human ESC and five types of differentiated cells. We identify thousands of alternatively spliced exons and reveal their development and lineage-dependent characterizations. Several histone modifications show dynamic changes in alternatively spliced exons and three are strongly associated with 52.8% of alternative splicing events upon hESC differentiation. The histone modification-associated alternatively spliced genes predominantly function in G2/M phases and ATM/ATR-mediated DNA damage response pathway for cell differentiation, whereas other alternatively spliced genes are enriched in the G1 phase and pathways for self-renewal. These results imply a potential epigenetic mechanism by which some histone modifications contribute to ESC fate decision through the regulation of alternative splicing in specific pathways and cell-cycle genes. Supported by experimental validations and extended datasets from Roadmap/ENCODE projects, we exemplify this mechanism by a cell-cycle-related transcription factor, PBX1, which regulates the pluripotency regulatory network by binding to NANOG. We suggest that the isoform switch from PBX1a to PBX1b links H3K36me3 to hESC fate determination through the PSIP1/SRSF1 adaptor, which results in the exon skipping of PBX1.

CONCLUSION

We reveal the mechanism by which alternative splicing links histone modifications to stem cell fate decision.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

In silico analysis of the sequence features responsible for alternatively spliced introns in the model green alga Chlamydomonas reinhardtii

Alternatively spliced introns are the ones that are usually spliced but can be occasionally retained in a transcript isoform. They are the most frequently used alternative splice form in plants (~50% of alternative splicing events). Chlamydomonas reinhardtii, a unicellular alga, is a good model to understand alternative splicing (AS) in plants from an evolutionary perspective as it diverged from land plants a billion years ago. Using over 7 million cDNA sequences from both pyrosequencing and Sanger sequencing, we found that a much higher percentage of genes (~20% of multi-exon genes) undergo AS than previously reported (3-5%). We found a full component of SR and SR-like proteins possibly involved in AS. The most prevalent type of AS event (40%) was retention of introns, most of which were supported by multiple cDNA evidence (72%) while only 20% of them have coding capacity. By comparing retained and constitutive introns, we identified sequence features potentially responsible for the retention of introns, in the framework of an "intron definition" model for splicing. We find that retained introns tend to have a weaker 5' splice site, more Gs in their poly-pyrimidine tract and a lesser conservation of nucleotide 'C' at position -3 of the 3' splice site. In addition, the sequence motifs found in the potential branch-point region differed between retained and constitutive introns. Furthermore, the enrichment of G-triplets and C-triplets among the first and last 50 nt of the introns significantly differ between constitutive and retained introns. These could serve as intronic splicing enhancers. All the alternative splice forms can be accessed at http://bioinfolab.miamioh.edu/cgi-bin/PASA_r20140417/cgi-bin/status_report.cgi?db=Chre_AS .

Merging Absolute and Relative Quantitative PCR Data to Quantify STAT3 Splice Variant Transcripts

Human signal transducer and activator of transcription 3 (STAT3) is one of many genes containing a tandem splicing site. Alternative donor splice sites 3 nucleotides apart result in either the inclusion (S) or exclusion (ΔS) of a single residue, Serine-701. Further downstream, splicing at a pair of alternative acceptor splice sites result in transcripts encoding either the 55 terminal residues of the transactivation domain (α) or a truncated transactivation domain with 7 unique residues (β). As outlined in this manuscript, measuring the proportions of STAT3's four spliced transcripts (Sα, Sβ, ΔSα and ΔSβ) was possible using absolute qPCR (quantitative polymerase chain reaction). The protocol therefore distinguishes and measures highly similar splice variants. Absolute qPCR makes use of calibrator plasmids and thus specificity of detection is not compromised for the sake of efficiency. The protocol necessitates primer validation and optimization of cycling parameters. A combination of absolute qPCR and efficiency-dependent relative qPCR of total STAT3 transcripts allowed a description of the fluctuations of STAT3 splice variants' levels in eosinophils treated with cytokines. The protocol also provided evidence of a co-splicing interdependence between the two STAT3 splicing events. The strategy based on a combination of the two qPCR techniques should be readily adaptable to investigation of co-splicing at other tandem splicing sites.

Genome-Wide Analysis of Alternative Splicing Provides Insights into Stress Adaptation of the Pacific Oyster

Alternative splicing (AS) is thought to enhance transcriptome diversity dramatically and play an important role in stress adaptation. While well studied in vertebrates, AS remains poorly understood in invertebrates. Here, we used high-throughput RNA-sequencing data to perform a genome-wide survey of AS in the Pacific oyster (Crassostrea gigas), an economically important mollusk that is cultivated worldwide. This analysis identified 8223 AS events corresponding to 4480 genes in the Pacific oyster, suggesting that about 16 % of oyster multiexonic genes undergo AS. We observed that a majority of the identified AS events were related to skipped exons (37.8 %). Then Gene Ontology analysis was conducted to analyze the function of the genes that undergo AS and the genes that produce more than five AS isoforms. After that, the expression of AS isoforms facing temperature, salinity, and air exposure challenge were examined. To validate our bioinformatic-predicted results and examine whether AS affects stress adaptation, we selected heat-shock protein 60 (HSP60) and HSP90 genes, both of which experience AS, for reverse transcription PCR (RT-PCR). We also performed quantitative real-time PCR (qRT-PCR) to determine the relative expression of each AS isoform among different stress adapted populations. Our study indicates that AS events are likely complex in the Pacific oyster and may be related to stress adaptation. These results will complement the predicted gene database of C. gigas and provide an invaluable resource for future functional genomic studies on molluscs.

Selection preserves Ubiquitin Specific Protease 4 alternative exon skipping in therian mammals

Ubiquitin specific protease 4 (USP4) is a highly networked deubiquitinating enzyme with reported roles in cancer, innate immunity and RNA splicing. In mammals it has two dominant isoforms arising from inclusion or skipping of exon 7 (E₇). We evaluated two plausible mechanisms for the generation of these isoforms: (A) E₇ skipping due to a long upstream intron and (B) E₇ skipping due to inefficient 5′ splice sites (5′SS) and/or branchpoint sites (BPS). We then assessed whether E₇ alternative splicing is maintained by selective pressure or arose from genetic drift. Both transcript variants were generated from a USP4-E₇ minigene construct with short flanking introns, an observation consistent with the second mechanism whereby differential splice signal strengths are the basis of E₇ skipping. Optimization of the downstream 5′SS eliminated E₇ skipping. Experimental validation of the correlation between 5′SS identity and exon skipping in vertebrates pinpointed the +6 site as the key splicing determinant. Therian mammals invariably display a 5′SS configuration favouring alternative splicing and the resulting isoforms have distinct subcellular localizations. We conclude that alternative splicing of mammalian USP4 is under selective maintenance and that long and short USP4 isoforms may target substrates in various cellular compartments.

Correlated Evolution of Nucleotide Positions within Splice Sites in Mammals

Splice sites (SSs)--short nucleotide sequences flanking introns--are under selection for spliceosome binding, and adhere to consensus sequences. However, non-consensus nucleotides, many of which probably reduce SS performance, are frequent. Little is known about the mechanisms maintaining such apparently suboptimal SSs. Here, we study the correlations between strengths of nucleotides occupying different positions of the same SS. Such correlations may arise due to epistatic interactions between positions (i.e., a situation when the fitness effect of a nucleotide in one position depends on the nucleotide in another position), their evolutionary history, or to other reasons. Within both the intronic and the exonic parts of donor SSs, nucleotides that increase (decrease) SS strength tend to co-occur with other nucleotides increasing (respectively, decreasing) it, consistent with positive epistasis. Between the intronic and exonic parts of donor SSs, the correlations of nucleotide strengths tend to be negative, consistent with negative epistasis. In the course of evolution, substitutions at a donor SS tend to decrease the strength of its exonic part, and either increase or do not change the strength of its intronic part. In acceptor SSs, the situation is more complicated; the correlations between adjacent positions appear to be driven mainly by avoidance of the AG dinucleotide which may cause aberrant splicing. In summary, both the content and the evolution of SSs is shaped by a complex network of interdependences between adjacent nucleotides that respond to a range of sometimes conflicting selective constraints.

Promoter-like epigenetic signatures in exons displaying cell type-specific splicing

BACKGROUND

Pre-mRNA splicing occurs mainly co-transcriptionally, and both nucleosome density and histone modifications have been proposed to play a role in splice site recognition and regulation. However, the extent and mechanisms behind this interplay remain poorly understood.

RESULTS

We use transcriptomic and epigenomic data generated by the ENCODE project to investigate the association between chromatin structure and alternative splicing. We find a strong and significant positive association between H3K9ac, H3K27ac, H3K4me3, epigenetic marks characteristic of active promoters, and exon inclusion in a small but well-defined class of exons, representing approximately 4 % of all regulated exons. These exons are systematically maintained at comparatively low levels of inclusion across cell types, but their inclusion is significantly enhanced in particular cell types when in physical proximity to active promoters.

CONCLUSION

Histone modifications and other chromatin features that activate transcription can be co-opted to participate in the regulation of the splicing of exons that are in physical proximity to promoter regions.

Splicing predictions reliably classify different types of alternative splicing

Alternative splicing is a key player in the creation of complex mammalian transcriptomes and its misregulation is associated with many human diseases. Multiple mRNA isoforms are generated from most human genes, a process mediated by the interplay of various RNA signature elements and trans-acting factors that guide spliceosomal assembly and intron removal. Here, we introduce a splicing predictor that evaluates hundreds of RNA features simultaneously to successfully differentiate between exons that are constitutively spliced, exons that undergo alternative 5' or 3' splice-site selection, and alternative cassette-type exons. Surprisingly, the splicing predictor did not feature strong discriminatory contributions from binding sites for known splicing regulators. Rather, the ability of an exon to be involved in one or multiple types of alternative splicing is dictated by its immediate sequence context, mainly driven by the identity of the exon's splice sites, the conservation around them, and its exon/intron architecture. Thus, the splicing behavior of human exons can be reliably predicted based on basic RNA sequence elements.

Mechanism of alternative splicing and its regulation

Alternative splicing of precursor mRNA is an essential mechanism to increase the complexity of gene expression, and it plays an important role in cellular differentiation and organism development. Regulation of alternative splicing is a complicated process in which numerous interacting components are at work, including cis-acting elements and trans-acting factors, and is further guided by the functional coupling between transcription and splicing. Additional molecular features, such as chromatin structure, RNA structure and alternative transcription initiation or alternative transcription termination, collaborate with these basic components to generate the protein diversity due to alternative splicing. All these factors contributing to this one fundamental biological process add up to a mechanism that is critical to the proper functioning of cells. Any corruption of the process may lead to disruption of normal cellular function and the eventuality of disease. Cancer is one of those diseases, where alternative splicing may be the basis for the identification of novel diagnostic and prognostic biomarkers, as well as new strategies for therapy. Thus, an in-depth understanding of alternative splicing regulation has the potential not only to elucidate fundamental biological principles, but to provide solutions for various diseases.

Genome-wide analysis of alternative splicing in Volvox carteri

Background

Alternative splicing is an essential mechanism for increasing transcriptome and proteome diversity in eukaryotes. Particularly in multicellular eukaryotes, this mechanism is involved in the regulation of developmental and physiological processes like growth, differentiation and signal transduction.

Results

Here we report the genome-wide analysis of alternative splicing in the multicellular green alga Volvox carteri. The bioinformatic analysis of 132,038 expressed sequence tags (ESTs) identified 580 alternative splicing events in a total of 426 genes. The predominant type of alternative splicing in Volvox is intron retention (46.5%) followed by alternative 5′ (17.9%) and 3′ (21.9%) splice sites and exon skipping (9.5%). Our analysis shows that in Volvox at least ~2.9% of the intron-containing genes are subject to alternative splicing. Considering the total number of sequenced ESTs, the Volvox genome seems to provide more favorable conditions (e.g., regarding length and GC content of introns) for the occurrence of alternative splicing than the genome of its close unicellular relative Chlamydomonas. Moreover, many randomly chosen alternatively spliced genes of Volvox do not show alternative splicing in Chlamydomonas. Since the Volvox genome contains about the same number of protein-coding genes as the Chlamydomonas genome (~14,500 protein-coding genes), we assumed that alternative splicing may play a key role in generation of genomic diversity, which is required to evolve from a simple one-cell ancestor to a multicellular organism with differentiated cell types (Mol Biol Evol 31:1402-1413, 2014). To confirm the alternative splicing events identified by bioinformatic analysis, several genes with different types of alternatively splicing have been selected followed by experimental verification of the predicted splice variants by RT-PCR.

Conclusions

The results show that our approach for prediction of alternative splicing events in Volvox was accurate and reliable. Moreover, quantitative real-time RT-PCR appears to be useful in Volvox for analyses of relationships between the appearance of specific alternative splicing variants and different kinds of physiological, metabolic and developmental processes as well as responses to environmental changes.

Electronic supplementary material

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

SeeSite: Characterizing Relationships between Splice Junctions and Splicing Enhancers

RNA splicing is a cellular process driven by the interaction between numerous regulatory sequences and binding sites, however, such interactions have been primarily explored by laboratory methods since computational tools largely ignore the relationship between different splicing elements. Current computational methods identify either splice sites or other regulatory sequences, such as enhancers and silencers. We present a novel approach for characterizing co-occurring relationships between splice site motifs and splicing enhancers. Our approach relies on an efficient algorithm for approximately solving Consensus Sequence with Outliers , an NP-complete string clustering problem. In particular, we give an algorithm for this problem that outputs near-optimal solutions in polynomial time. To our knowledge, this is the first formulation and computational attempt for detecting co-occurring sequence elements in RNA sequence data. Further, we demonstrate that SeeSite is capable of showing that certain ESEs are preferentially associated with weaker splice sites, and that there exists a co-occurrence relationship with splice site motifs.

Global remodeling of nucleosome positions in C. elegans

Background

Eukaryotic chromatin architecture is affected by intrinsic histone-DNA sequence preferences, steric exclusion between nucleosome particles, formation of higher-order structures, and in vivo activity of chromatin remodeling enzymes.

Results

To disentangle sequence-dependent nucleosome positioning from the other factors, we have created two high-throughput maps of nucleosomes assembled in vitro on genomic DNA from the nematode worm Caenorhabditis elegans. A comparison of in vitro nucleosome positions with those observed in a mixed-stage, mixed-tissue population of C. elegans cells reveals that in vivo sequence preferences are modified on the genomic scale. Indeed, G/C dinucleotides are predicted to be most favorable for nucleosome formation in vitro but not in vivo. Nucleosome sequence read coverage in vivo is distinctly lower in chromosome arms than in central regions; the observed changes in apparent nucleosome sequence specificity, likely due to genome-wide chromatin remodeler activity, contribute to the formation of these megabase-scale chromatin domains. We also observe that the majority of well-positioned in vivo nucleosomes do not occupy thermodynamically favorable sequences observed in vitro. Finally, we find that exons are intrinsically more amenable to nucleosome formation compared to introns. Nucleosome occupancy of introns and exons consistently increases with G/C content in vitro but not in vivo, in agreement with our observation that G/C dinucleotide enrichment does not strongly promote in vivo nucleosome formation.

Conclusions

Our findings highlight the importance of both sequence specificity and active nucleosome repositioning in creating large-scale chromatin domains, and the antagonistic roles of intrinsic sequence preferences and chromatin remodelers in C. elegans.

Sequence read data has been deposited into Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra; accession number SRA050182). Additional data, software and computational predictions are available on the Nucleosome Explorer website (http://nucleosome.rutgers.edu).

Accurate identification and analysis of human mRNA isoforms using deep long read sequencing

Precise identification of RNA-coding regions and transcriptomes of eukaryotes is a significant problem in biology. Currently, eukaryote transcriptomes are analyzed using deep short-read sequencing experiments of complementary DNAs. The resulting short-reads are then aligned against a genome and annotated junctions to infer biological meaning. Here we use long-read complementary DNA datasets for the analysis of a eukaryotic transcriptome and generate two large datasets in the human K562 and HeLa S3 cell lines. Both data sets comprised at least 4 million reads and had median read lengths greater than 500 bp. We show that annotation-independent alignments of these reads provide partial gene structures that are very much in-line with annotated gene structures, 15% of which have not been obtained in a previous de novo analysis of short reads. For long-noncoding RNAs (i.e., lncRNA) genes, however, we find an increased fraction of novel gene structures among our alignments. Other important aspects of transcriptome analysis, such as the description of cell type-specific splicing, can be performed in an accurate, reliable and completely annotation-free manner, making it ideal for the analysis of transcriptomes of newly sequenced genomes. Furthermore, we demonstrate that long read sequence can be assembled into full-length transcripts with considerable success. Our method is applicable to all long read sequencing technologies.

Removal of retained introns regulates translation in the rapidly developing gametophyte of Marsilea vestita

The utilization of stored RNA is a driving force in rapid development. Here, we show that retention and subsequent removal of introns from pre-mRNAs regulate temporal patterns of translation during rapid and posttranscriptionally controlled spermatogenesis of the fern Marsilea vestita. Analysis of RNAseq-derived transcriptomes revealed a large subset of intron-retaining transcripts (IRTs) that encode proteins essential for gamete development. Genomic and IRT sequence comparisons show that other introns have been previously removed from the IRT pre-mRNAs. Fully spliced isoforms appear at distinct times during development in a spliceosome-dependent and transcription-independent manner. RNA interference knockdowns of 17/17 IRTs produced anomalies after the time points when those transcripts would normally be spliced. Intron retention is a functional mechanism for forestalling precocious translation of transcripts in the male gametophyte of M. vestita. These results have broad implications for plant gene regulation, where intron retention is widespread.

Position weight matrix, gibbs sampler, and the associated significance tests in motif characterization and prediction

Position weight matrix (PWM) is not only one of the most widely used bioinformatic methods, but also a key component in more advanced computational algorithms (e.g., Gibbs sampler) for characterizing and discovering motifs in nucleotide or amino acid sequences. However, few generally applicable statistical tests are available for evaluating the significance of site patterns, PWM, and PWM scores (PWMS) of putative motifs. Statistical significance tests of the PWM output, that is, site-specific frequencies, PWM itself, and PWMS, are in disparate sources and have never been collected in a single paper, with the consequence that many implementations of PWM do not include any significance test. Here I review PWM-based methods used in motif characterization and prediction (including a detailed illustration of the Gibbs sampler for de novo motif discovery), present statistical and probabilistic rationales behind statistical significance tests relevant to PWM, and illustrate their application with real data. The multiple comparison problem associated with the test of site-specific frequencies is best handled by false discovery rate methods. The test of PWM, due to the use of pseudocounts, is best done by resampling methods. The test of individual PWMS for each sequence segment should be based on the extreme value distribution.

Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis

Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. As very few genome-wide studies analyzing AS in plants are available, we have performed high-throughput sequencing of a normalized cDNA library which resulted in a high coverage transcriptome map of Arabidopsis. We detect ∼150,000 splice junctions derived mostly from typical plant introns, including an eightfold increase in the number of U12 introns (2069). Around 61% of multiexonic genes are alternatively spliced under normal growth conditions. Moreover, we provide experimental validation of 540 AS transcripts (from 256 genes coding for important regulatory factors) using high-resolution RT-PCR and Sanger sequencing. Intron retention (IR) is the most frequent AS event (∼40%), but many IRs have relatively low read coverage and are less well-represented in assembled transcripts. Additionally, ∼51% of Arabidopsis genes produce AS transcripts which do not involve IR. Therefore, the significance of IR in generating transcript diversity was generally overestimated in previous assessments. IR analysis allowed the identification of a large set of cryptic introns inside annotated coding exons. Importantly, a significant fraction of these cryptic introns are spliced out in frame, indicating a role in protein diversity. Furthermore, we show extensive AS coupled to nonsense-mediated decay in AFC2, encoding a highly conserved LAMMER kinase which phosphorylates splicing factors, thus establishing a complex loop in AS regulation. We provide the most comprehensive analysis of AS to date which will serve as a valuable resource for the plant community to study transcriptome complexity and gene regulation.

Factors affecting splicing strength of yeast genes

Accurate and efficient splicing is of crucial importance for highly-transcribed intron-containing genes (ICGs) in rapidly replicating unicellular eukaryotes such as the budding yeast Saccharomyces cerevisiae. We characterize the 5' and 3' splice sites (ss) by position weight matrix scores (PWMSs), which is the highest for the consensus sequence and the lowest for splice sites differing most from the consensus sequence and used PWMS as a proxy for splicing strength. HAC1, which is known to be spliced by a nonspliceosomal mechanism, has the most negative PWMS for both its 5' ss and 3' ss. Several genes under strong splicing regulation and requiring additional splicing factors for their splicing also have small or negative PWMS values. Splicing strength is higher for highly transcribed ICGs than for lowly transcribed ICGs and higher for transcripts that bind strongly to spliceosomes than those that bind weakly. The 3' splice site features a prominent poly-U tract before the 3'AG. Our results suggest the potential of using PWMS as a screening tool for ICGs that are either spliced by a nonspliceosome mechanism or under strong splicing regulation in yeast and other fungal species.

Characterization of intronic structures and alternative splicing in Phytophthora sojae by comparative analysis of expressed sequence tags and genomic sequences

The oomycetes, a distinct phylogenetic lineage of fungus-like microorganisms, are heterokonts (stramenopiles) belonging to the supergroup Chromalveolata. Although the complete genomic sequences of a number of oomycetes have been reported, little information regarding the introns therein is available. Here, we investigated the introns of Phytophthora sojae, a pathogen that causes soybean root and stem rot, by a comparative analysis of genomic sequences and expressed sequence tags. A total of 4013 introns were identified, of which 96.6% contained canonical splice sites. The P. sojae genome possessed features distinct from other organisms at 5' splice sites, polypyrimidine tracts, branch sites, and 3' splice sites. Diverse repeating sequences, ranging from 2 to 10 nucleotides in length, were found at more than half of the intron-exon boundaries. Furthermore, 122 genes underwent alternative splicing. These data indicate that P. sojae has unique splicing mechanisms, and recognition of those mechanisms may lead to more accurate predictions of the location of introns in P. sojae and even other oomycete species.

Genome-wide analysis of alternative splicing in Chlamydomonas reinhardtii

Background

Genome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs.

Results

Our analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy.

Conclusions

The extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

Genome-wide evidence for selection acting on single amino acid repeats

Low complexity and homopolymer sequences within coding regions are known to evolve rapidly. While their expansion may be deleterious, there is increasing evidence for a functional role associated with these amino acid sequences. Homopolymer sequences are thought to evolve mostly through replication slippage and, therefore, they may be expected to be longer in regions with relaxed selective constraint. Within the coding sequences of eukaryotes, alternatively spliced exons are known to evolve under relaxed constraints in comparison to those exons that are constitutively spliced because they are not included in all of the mature mRNA of a gene. This relaxed exposure to selection leads to faster rates of evolution for alternatively spliced exons in comparison to constitutively spliced exons. Here, we have tested the effect of splicing on the structure (composition, length) of homopolymer sequences in relation to the splicing pattern in which they are found. We observed a significant relationship between alternative splicing and homopolymer sequences with alternatively spliced genes being enriched in number and length of homopolymer sequences. We also observed lower codon diversity and longer homocodons, suggesting a balance between slippage and point mutations linked to the constraints imposed by selection.

Complex selection on 5' splice sites in intron-rich organisms

In contrast to the typically streamlined genomes of prokaryotes, many eukaryotic genomes are riddled with long intergenic regions, spliceosomal introns, and repetitive elements. What explains the persistence of these and other seemingly suboptimal structures? There are three general hypotheses: (1) the structures in question are not actually suboptimal but optimal, being favored by selection, for unknown reasons; (2) the structures are not suboptimal, but of (essentially) equal fitness to "optimal" ones; or (3) the structures are truly suboptimal, but selection is too weak to systematically eliminate them. The 5' splice sites of introns offer a rare opportunity to directly test these hypotheses. Intron-poor species show a clear consensus splice site; most introns begin with the same six nucleotide sequence (typically GTAAGT or GTATGT), indicating efficient selection for this consensus sequence. In contrast, intron-rich species have much less pronounced boundary consensus sequences, and only small minorities of introns in intron-rich species share the same boundary sequence. We studied rates of evolutionary change of 5' splice sites in three groups of closely related intron-rich species--three primates, five Drosophila species, and four Cryptococcus fungi. Surprisingly, the results indicate that changes from consensus-to-variant nucleotides are generally disfavored by selection, but that changes from variant to consensus are neither favored nor disfavored. This evolutionary pattern is consistent with selective differences across introns, for instance, due to compensatory changes at other sites within the gene, which compensate for the otherwise suboptimal consensus-to-variant changes in splice boundaries.

Evolution of alternative splicing regulation: changes in predicted exonic splicing regulators are not associated with changes in alternative splicing levels in primates

Alternative splicing is tightly regulated in a spatio-temporal and quantitative manner. This regulation is achieved by a complex interplay between spliceosomal (trans) factors that bind to different sequence (cis) elements. cis-elements reside in both introns and exons and may either enhance or silence splicing. Differential combinations of cis-elements allows for a huge diversity of overall splicing signals, together comprising a complex ‘splicing code’. Many cis-elements have been identified, and their effects on exon inclusion levels demonstrated in reporter systems. However, the impact of interspecific differences in these elements on the evolution of alternative splicing levels has not yet been investigated at genomic level. Here we study the effect of interspecific differences in predicted exonic splicing regulators (ESRs) on exon inclusion levels in human and chimpanzee. For this purpose, we compiled and studied comprehensive datasets of predicted ESRs, identified by several computational and experimental approaches, as well as microarray data for changes in alternative splicing levels between human and chimpanzee. Surprisingly, we found no association between changes in predicted ESRs and changes in alternative splicing levels. This observation holds across different ESR exon positions, exon lengths, and 5′ splice site strengths. We suggest that this lack of association is mainly due to the great importance of context for ESR functionality: many ESR-like motifs in primates may have little or no effect on splicing, and thus interspecific changes at short-time scales may primarily occur in these effectively neutral ESRs. These results underscore the difficulties of using current computational ESR prediction algorithms to identify truly functionally important motifs, and provide a cautionary tale for studies of the effect of SNPs on splicing in human disease.

Comparative component analysis of exons with different splicing frequencies

Transcriptional isoforms are not just random combinations of exons. What has caused exons to be differentially spliced and whether exons with different splicing frequencies are subjected to divergent regulation by potential elements or splicing signals? Beyond the conventional classification for alternatively spliced exons (ASEs) and constitutively spliced exons (CSEs), we have classified exons from alternatively spliced human genes and their mouse orthologs (12,314 and 5,464, respectively) into four types based on their splicing frequencies. Analysis has indicated that different groups of exons presented divergent compositional and regulatory properties. Interestingly, with the decrease of splicing frequency, exons tend to have greater lengths, higher GC content, and contain more splicing elements and repetitive elements, which seem to imply that the splicing frequency is influenced by such factors. Comparison of non-alternatively spliced (NAS) mouse genes with alternatively spliced human orthologs also suggested that exons with lower splicing frequencies may be newly evolved ones which gained functions with splicing frequencies altered through the evolution. Our findings have revealed for the first time that certain factors may have critical influence on the splicing frequency, suggesting that exons with lower splicing frequencies may originate from old repetitive sequences, with splicing sites altered by mutation, gaining novel functions and become more frequently spliced.

Systematic analysis of alternative promoters correlated with alternative splicing in human genes

Interactions between various events are essential for complex and delicate transcriptional regulation. To delineate the features and potential roles of alternative promoters (APs) correlated with alternative splicing (AS), we have systematically analyzed 9908 putative alternative promoters (PAPs) from 3797 human genes. Our results showed that approximately 65% of AS events are associated with PAPs. Intriguingly, PAPs per human AS gene only averaged 2.6 for our dataset, which was significantly lower than previously reported. This seems to imply that the human genome contains a small pool of appropriable PAPs for AS genes. Exploration of the characteristics of PAPs such as CpG islands, TATA boxes, GC-content, transcription factor binding sites (TFBSs) and repetitive elements suggested that, respectively, 87% and 90% of PAPs of human AS genes are CpG- and TATA box-poor. The GC-content is significantly higher in the downstream of transcription start sites (TSSs) than upstream (58% vs. 53%), and there is a strong negative correlation between the GC-content and the number of PAPs. These suggested that GC-content around the TSSs plays an important role in the regulation of AS. Moreover, different APs contain distinct densities of repetitive elements and TFBSs, indicating that such sequences have an intrinsic role in the divergent regulation of PAPs and AS. Substantial difference was also found between human AS genes in terms of PAP numbers. A close connection between PAPs and AS may play a critical role in the choice of APs and regulation of AS genes. Furthermore, the distribution of AS genes on different human chromosomes also influences the numbers of PAPs and isoforms of AS genes. Our results may provide important clues for further studies on regulatory network of transcription-related events.

Where do introns come from?

The systems for mRNA surveillance, capping, and cleavage/polyadenylation are proposed to play pivotal roles in the physical establishment and distribution of spliceosomal introns along a transcript.

The effect of intron length on exon creation ratios during the evolution of mammalian genomes

Recent studies report that alternatively spliced exons tend to occur in longer introns, which is attributed to the length constraints for splice site pairing for the two major splicing mechanisms, intron definition versus exon definition. Using genome-wide studies of EST and microarray data from human and mouse, we have analyzed the distribution of various subsets of alternatively spliced exons, based on their inclusion level and evolutionary history, versus increasing intron length. Alternative exons may be included in either a major or minor fraction of all transcripts (known as major-form and minor-form exons, respectively). We find that major-form exons are seven- to eightfold more likely to be contained in short introns (<400 nt) than minor-form exons, which occur preferentially in longer introns. Since minor-form exons are more likely to be novel (approximately 75%), this implied that novel exons arise more frequently in longer introns. To test this hypothesis, we used whole genome alignments to classify exons according to their phylogenetic age. We find that older exons, i.e., exons that are conserved in all mammals, predominate at shorter intron lengths, for both major- and minor-form exons. In contrast, exons that arose recently during primate evolution are more prevalent at longer intron lengths (>1000 nt). This suggests that the observed correlation of longer intron lengths with alternatively spliced exons may be at least partly due to biases in the probability of exon creation, which is higher in long introns.

Width of gene expression profile drives alternative splicing

Alternative splicing generates an enormous amount of functional and proteomic diversity in metazoan organisms. This process is probably central to the macromolecular and cellular complexity of higher eukaryotes. While most studies have focused on the molecular mechanism triggering and controlling alternative splicing, as well as on its incidence in different species, its maintenance and evolution within populations has been little investigated. Here, we propose to address these questions by comparing the structural characteristics as well as the functional and transcriptional profiles of genes with monomorphic or polymorphic splicing, referred to as MS and PS genes, respectively. We find that MS and PS genes differ particularly in the number of tissues and cell types where they are expressed.We find a striking deficit of PS genes on the sex chromosomes, particularly on the Y chromosome where it is shown not to be due to the observed lower breadth of expression of genes on that chromosome. The development of a simple model of evolution of cis-regulated alternative splicing leads to predictions in agreement with these observations. It further predicts the conditions for the emergence and the maintenance of cis-regulated alternative splicing, which are both favored by the tissue specific expression of splicing variants. We finally propose that the width of the gene expression profile is an essential factor for the acquisition of new transcript isoforms that could later be maintained by a new form of balancing selection.

The "alternative" choice of constitutive exons throughout evolution

Alternative cassette exons are known to originate from two processes—exonization of intronic sequences and exon shuffling. Herein, we suggest an additional mechanism by which constitutively spliced exons become alternative cassette exons during evolution. We compiled a dataset of orthologous exons from human and mouse that are constitutively spliced in one species but alternatively spliced in the other. Examination of these exons suggests that the common ancestors were constitutively spliced. We show that relaxation of the 5′ splice site during evolution is one of the molecular mechanisms by which exons shift from constitutive to alternative splicing. This shift is associated with the fixation of exonic splicing regulatory sequences (ESRs) that are essential for exon definition and control the inclusion level only after the transition to alternative splicing. The effect of each ESR on splicing and the combinatorial effects between two ESRs are conserved from fish to human. Our results uncover an evolutionary pathway that increases transcriptome diversity by shifting exons from constitutive to alternative splicing.

Alternative splicing is believed to play a major role in the creation of transcriptomic diversification leading to higher order of organismal complexity, especially in mammals. As much as 80% of human genes generate more than one type of mRNA by alternative splicing. Thus, alternative splicing can bridge the low number of protein coding genes (∼24,500) and the total number of proteins generated in the human proteome (∼90,000). The correlation between the higher order of phenotypic diversity and alternative splicing was recently demonstrated and thus the origin of alternative splicing is of great interest. There are currently two models regarding the origin of alternatively spliced exons—exonization of intronic sequences and exon shuffling. According to these two mechanisms, a protein-coding gene was first established and only then a new alternative exon appeared within it or was added to the gene. Our current study provides evidences for a new mechanism indicating that during evolution constitutively spliced exons became alternatively spliced. Large-scale bioinformatic analyses reveal the magnitude of this process and experimental validation systems provide insights into its mechanisms.

A systematic analysis of intronic sequences downstream of 5' splice sites reveals a widespread role for U-rich motifs and TIA1/TIAL1 proteins in alternative splicing regulation

To identify human intronic sequences associated with 5' splice site recognition, we performed a systematic search for motifs enriched in introns downstream of both constitutive and alternative cassette exons. Significant enrichment was observed for U-rich motifs within 100 nucleotides downstream of 5' splice sites of both classes of exons, with the highest enrichment between positions +6 and +30. Exons adjacent to U-rich intronic motifs contain lower frequencies of exonic splicing enhancers and higher frequencies of exonic splicing silencers, compared with exons not followed by U-rich intronic motifs. These findings motivated us to explore the possibility of a widespread role for U-rich motifs in promoting exon inclusion. Since cytotoxic granule-associated RNA binding protein (TIA1) and TIA1-like 1 (TIAL1; also known as TIAR) were previously shown in vitro to bind to U-rich motifs downstream of 5' splice sites, and to facilitate 5' splice site recognition in vitro and in vivo, we investigated whether these factors function more generally in the regulation of splicing of exons followed by U-rich intronic motifs. Simultaneous knockdown of TIA1 and TIAL1 resulted in increased skipping of 36/41 (88%) of alternatively spliced exons associated with U-rich motifs, but did not affect 32/33 (97%) alternatively spliced exons that are not associated with U-rich motifs. The increase in exon skipping correlated with the proximity of the first U-rich motif and the overall "U-richness" of the adjacent intronic region. The majority of the alternative splicing events regulated by TIA1/TIAL1 are conserved in mouse, and the corresponding genes are associated with diverse cellular functions. Based on our results, we estimate that approximately 15% of alternative cassette exons are regulated by TIA1/TIAL1 via U-rich intronic elements.

Searching for splicing motifs

Intron removal during pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. The sequences constituting the splice sites provide insufficient information to distinguish true splice sites from the greater number of false splice sites that populate transcripts. Additional information used for exon recognition resides in a large number of positively or negatively acting elements that lie both within exons and in the adjacent introns. The identification of such sequence motifs has progressed rapidly in recent years, such that extensive lists are now available for exonic splicing enhancers and exonic splicing silencers. These motifs have been identified both by empirical experiments and by computational predictions, the validity of the latter being confirmed by experimental verification. Molecular searches have been carried out either by the selection of sequences that bind to splicing factors, or enhance or silence splicing in vitro or in vivo. Computational methods have focused on sequences of 6 or 8 nucleotides that are over- or under-represented in exons, compared to introns or transcripts that do not undergo splicing. These various methods have sought to provide global definitions of motifs, yet the motifs are distinctive to the method used for identification and display little overlap. Astonishingly, at least three-quarters of a typical mRNA would be comprised of these motifs. A present challenge lies in understanding how the cell integrates this surfeit of information to generate what is usually a binary splicing decision.

Widespread and subtle: alternative splicing at short-distance tandem sites

Alternative splicing at donor or acceptor sites located just a few nucleotides apart is widespread in many species. It results in subtle changes in the transcripts and often in the encoded proteins. Several of these tandem splice events contribute to the repertoire of functionally different proteins, whereas many are neutral or deleterious. Remarkably, some of the functional events are differentially spliced in tissues or developmental stages, whereas others exhibit constant splicing ratios, indicating that function is not always associated with differential splicing. Stochastic splice site selection seems to play a major role in these processes. Here, we review recent progress in understanding functional and evolutionary aspects as well as the mechanism of splicing at short-distance tandem sites.

Alternative approach to a heavy weight problem

Obesity is reaching epidemic proportions in developed countries and represents a significant risk factor for hypertension, heart disease, diabetes, and dyslipidemia. Splicing mutations constitute at least 14% of disease-causing mutations, thus implicating polymorphisms that affect splicing as likely candidates for disease susceptibility. A recent study suggested that genes associated with obesity were significantly enriched for rare nucleotide variants. Here, we examined these variants and revealed that they are located near splice junctions and tend to affect exonic splicing regulatory sequences. We also show that the majority of the exons that harbor these SNPs are constitutively spliced, yet they exhibit weak splice sites, typical to alternatively spliced exons, and are hence suboptimal for recognition by the splicing machinery and prone to become alternatively spliced. Using ex vivo assays, we tested a few representative variants and show that they indeed affect splicing by causing a shift from a constitutive to an alternative pattern, suggesting a possible link between extreme body mass index and abnormal splicing patterns.