Violating the splicing rules: TG dinucleotides function as alternative 3' splice sites in U2-dependent introns

A-Z of Epigenetic Readers: Targeting Alternative Splicing and Histone Modification Variants in Cancer

Epigenetic 'reader' proteins, which have evolved to interact with specific chromatin modifications, play pivotal roles in gene regulation. There is growing interest in the alternative splicing mechanisms that affect the functionality of such epigenetic readers in cancer etiology. The current review considers how deregulation of epigenetic processes and alternative splicing events contribute to pathophysiology. An A-Z guide of epigenetic readers is provided, delineating the antagonistic 'yin-yang' roles of full-length versus spliced isoforms, where this is known from the literature. The examples discussed underscore the key contributions of epigenetic readers in transcriptional regulation, early development, and cancer. Clinical implications are considered, offering insights into precision oncology and targeted therapies focused on epigenetic readers that have undergone alternative splicing events during disease pathogenesis. This review underscores the fundamental importance of alternative splicing events in the context of epigenetic readers while emphasizing the critical need for improved understanding of functional diversity, regulatory mechanisms, and future therapeutic potential.

Global analysis of binding sites of U2AF1 and ZRSR2 reveals RNA elements required for mutually exclusive splicing by the U2- and U12-type spliceosome

Recurring mutations in genes encoding 3' splice-site recognition proteins, U2AF1 and ZRSR2 are associated with human cancers. Here, we determined binding sites of the proteins to reveal that U2-type and U12-type splice sites are recognized by U2AF1 and ZRSR2, respectively. However, some sites are spliced by both the U2-type and U12-type spliceosomes, indicating that well-conserved consensus motifs in some U12-type introns could be recognized by the U2-type spliceosome. Nucleotides flanking splice sites of U12-type introns are different from those flanking U2-type introns. Remarkably, the AG dinucleotide at the positions -1 and -2 of 5' splice sites of U12-type introns with GT-AG termini is not present. AG next to 5' splice site introduced by a single nucleotide substitution at the -2 position could convert a U12-type splice site to a U2-type site. The class switch of introns by a single mutation and the bias against G at the -1 position of U12-type 5' splice site support the notion that the identities of nucleotides in exonic regions adjacent to splice sites are fine-tuned to avoid recognition by the U2-type spliceosome. These findings may shed light on the mechanism of selectivity in U12-type intron splicing and the mutations that affect splicing.

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

Comprehensive database and evolutionary dynamics of U12-type introns

During nuclear maturation of most eukaryotic pre-messenger RNAs and long non-coding RNAs, introns are removed through the process of RNA splicing. Different classes of introns are excised by the U2-type or the U12-type spliceosomes, large complexes of small nuclear ribonucleoprotein particles and associated proteins. We created intronIC, a program for assigning intron class to all introns in a given genome, and used it on 24 eukaryotic genomes to create the Intron Annotation and Orthology Database (IAOD). We then used the data in the IAOD to revisit several hypotheses concerning the evolution of the two classes of spliceosomal introns, finding support for the class conversion model explaining the low abundance of U12-type introns in modern genomes.

Functionally Overlapping Variants Control Tuberculosis Susceptibility in Collaborative Cross Mice

The variable outcome of Mycobacterium tuberculosis infection observed in natural populations is difficult to model in genetically homogeneous small-animal models. The newly developed Collaborative Cross (CC) represents a reproducible panel of genetically diverse mice that display a broad range of phenotypic responses to infection. We explored the genetic basis of this variation, focusing on a CC line that is highly susceptible to M. tuberculosis infection. This study identified multiple quantitative trait loci associated with bacterial control and cytokine production, including one that is caused by a novel loss-of-function mutation in the Itgal gene, which is necessary for T cell recruitment to the infected lung. These studies verify the multigenic control of mycobacterial disease in the CC panel, identify genetic loci controlling diverse aspects of pathogenesis, and highlight the utility of the CC resource.

Host genetics plays an important role in determining the outcome of Mycobacterium tuberculosis infection. We previously found that Collaborative Cross (CC) mouse strains differ in their susceptibility to M. tuberculosis and that the CC042/GeniUnc (CC042) strain suffered from a rapidly progressive disease and failed to produce the protective cytokine gamma interferon (IFN-γ) in the lung. Here, we used parallel genetic and immunological approaches to investigate the basis of CC042 mouse susceptibility. Using a population derived from a CC001/Unc (CC001) × CC042 intercross, we mapped four quantitative trait loci (QTL) underlying tuberculosis immunophenotypes (Tip1 to Tip4). These included QTL that were associated with bacterial burden, IFN-γ production following infection, and an IFN-γ-independent mechanism of bacterial control. Further immunological characterization revealed that CC042 animals recruited relatively few antigen-specific T cells to the lung and that these T cells failed to express the integrin alpha L (αL; i.e., CD11a), which contributes to T cell activation and migration. These defects could be explained by a CC042 private variant in the Itgal gene, which encodes CD11a and is found within the Tip2 interval. This 15-bp deletion leads to aberrant mRNA splicing and is predicted to result in a truncated protein product. The Itgal^CC042 genotype was associated with all measured disease traits, indicating that this variant is a major determinant of susceptibility in CC042 mice. The combined effect of functionally distinct Tip variants likely explains the profound susceptibility of CC042 mice and highlights the multigenic nature of tuberculosis control in the Collaborative Cross.

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.

Characterisation of Antigen B Protein Species Present in the Hydatid Cyst Fluid of Echinococcus canadensis G7 Genotype

The larva of cestodes belonging to the Echinococcus granulosus sensu lato (s.l.) complex causes cystic echinococcosis (CE). It is a globally distributed zoonosis with significant economic and public health impact. The most immunogenic and specific Echinococcus-genus antigen for human CE diagnosis is antigen B (AgB), an abundant lipoprotein of the hydatid cyst fluid (HF). The AgB protein moiety (apolipoprotein) is encoded by five genes (AgB1-AgB5), which generate mature 8 kDa proteins (AgB8/1-AgB8/5). These genes seem to be differentially expressed among Echinococcus species. Since AgB immunogenicity lies on its protein moiety, differences in AgB expression within E. granulosus s.l. complex might have diagnostic and epidemiological relevance for discriminating the contribution of distinct species to human CE. Interestingly, AgB2 was proposed as a pseudogene in E. canadensis, which is the second most common cause of human CE, but proteomic studies for verifying it have not been performed yet. Herein, we analysed the protein and lipid composition of AgB obtained from fertile HF of swine origin (E. canadensis G7 genotype). AgB apolipoproteins were identified and quantified using mass spectrometry tools. Results showed that AgB8/1 was the major protein component, representing 71% of total AgB apolipoproteins, followed by AgB8/4 (15.5%), AgB8/3 (13.2%) and AgB8/5 (0.3%). AgB8/2 was not detected. As a methodological control, a parallel analysis detected all AgB apolipoproteins in bovine fertile HF (G1/3/5 genotypes). Overall, E. canadensis AgB comprised mostly AgB8/1 together with a heterogeneous mixture of lipids, and AgB8/2 was not detected despite using high sensitivity proteomic techniques. This endorses genomic data supporting that AgB2 behaves as a pseudogene in G7 genotype. Since recombinant AgB8/2 has been found to be diagnostically valuable for human CE, our findings indicate that its use as antigen in immunoassays could contribute to false negative results in areas where E. canadensis circulates. Furthermore, the presence of anti-AgB8/2 antibodies in serum may represent a useful parameter to rule out E. canadensis infection when human CE is diagnosed.

Cystic echinococcosis (CE), a worldwide-spread zoonosis, affects livestock mammals and humans with significant economic and public health impact. It is caused by the infection with the larva of cestodes belonging to Echinococcus granulosus complex, a series of parasite species with preference for different hosts. Among them, Echinococcus canadensis larva uses mainly camels, goats and pigs as hosts. Species/genotypes belonging to E. canadensis are considered the second most common cause of human CE, but its contribution may be underestimated since causes asymptomatic or more benign infections than other E. granulosus complex species. The most relevant antigen for CE diagnosis is a lipoprotein called antigen B (AgB). AgB antigenicity is linked to its protein moiety that is encoded by several genes. One of these genes, AgB2, seems to be differentially expressed within E. granulosus complex. Using high sensitivity proteomic tools we analysed the composition of AgB obtained from E. canadensis larva, detecting the protein products of all AgB genes, except AgB2. Since AgB subunits have been widely used as antigens in immunoassays for human CE diagnosis, our results indicate that using AgB2 protein product in these assays may lead to false-negative results, particularly in geographical areas where E. canadensis species/genotypes circulate.

Splice Sites Seldom Slide: Intron Evolution in Oomycetes

We examine exon junctions near apparent amino acid insertions and deletions in alignments of orthologous protein-coding genes. In 1,917 ortholog families across nine oomycete genomes, 10–20% of introns are near an alignment gap, indicating at first sight that splice-site displacements are frequent. We designed a robust algorithmic procedure for the delineation of intron-containing homologous regions, and combined it with a parsimony-based reconstruction of intron loss, gain, and splice-site shift events on a phylogeny. The reconstruction implies that 12% of introns underwent an acceptor-site shift, and 10% underwent a donor-site shift. In order to offset gene annotation problems, we amended the procedure with the reannotation of intron boundaries using alignment evidence. The corresponding reconstruction involves much fewer intron gain and splice-site shift events. The frequency of acceptor- and donor-side shifts drops to 4% and 3%, respectively, which are not much different from what one would expect by random codon insertions and deletions. In other words, gaps near exon junctions are mostly artifacts of gene annotation rather than evidence of sliding intron boundaries. Our study underscores the importance of using well-supported gene structure annotations in comparative studies. When transcription evidence is not available, we propose a robust ancestral reconstruction procedure that corrects misannotated intron boundaries using sequence alignments. The results corroborate the view that boundary shifts and complete intron sliding are only accidental in eukaryotic genome evolution and have a negligible impact on protein diversity.

It's a bit over, is that ok? The subtle surplus from tandem alternative splicing

Tandem alternative splice sites (TASS) form a defined class of alternative splicing and give rise to mRNA insertion/deletion variants with only small size differences. Previous work has confirmed evolutionary conservation of TASS elements while many cases show only low tissue specificity of isoform ratios. We pinpoint stochasticity and noise as important methodological issues for the dissection of TASS isoform patterns. Resolving such uncertainties, a recent report showed regulation in a cell culture system, with shifts of alternative splicing isoform ratios dependent on cell density. This novel type of regulation affects not only multiple TASS isoforms, but also other alternative splicing classes, in a concerted manner. Here, we discuss how specific regulatory network architectures may be realized through the novel regulation type and highlight the role of differential isoform functions as a key step in order to better understand the functional role of TASS.

Gene Model Annotations for Drosophila melanogaster: The Rule-Benders

In the context of the FlyBase annotated gene models in Drosophila melanogaster, we describe the many exceptional cases we have curated from the literature or identified in the course of FlyBase analysis. These range from atypical but common examples such as dicistronic and polycistronic transcripts, noncanonical splices, trans-spliced transcripts, noncanonical translation starts, and stop-codon readthroughs, to single exceptional cases such as ribosomal frameshifting and HAC1-type intron processing. In FlyBase, exceptional genes and transcripts are flagged with Sequence Ontology terms and/or standardized comments. Because some of the rule-benders create problems for handlers of high-throughput data, we discuss plans for flagging these cases in bulk data downloads.

A novel laccase from basidiomycete Cerrena sp.: Cloning, heterologous expression, and characterization

A novel laccase gene Lac1 and its cDNA were cloned from a white-rot fungus Cerrena sp. and characterized. The 1554-bp cDNA of Lac1 encoded a mature protein with 497 amino acids, preceded by a signal peptide of 20 amino acids. An unconventional intron splice site and incomplete splicing variants of Lac1 were observed. Lac1 was heterologously expressed in the yeast host Pichia pastoris, and a maximal laccase activity of 6.3UmL(-1) in the fermentation broth was achieved after fermentation for 9 days. The recombinant protein rLac1 was purified, and its enzymatic properties and functional characteristics were investigated. When ABTS was used as the substrate, the enzyme was most active at pH 3.5 and 55°C, and stable at pH 4-10 and 20-60°C. The Km and kcat values of rLac1 toward ABTS were 28.9 μM and 332.4s(-1), respectively. Furthermore, rLac1 was tolerant to common metal ions up to 100mM concentration and capable of decolorizing structurally different dyes in the absence of a redox mediator. Hence, Lac1 may be useful for industrial applications, such as dye decolorization and bioremediation.

A comprehensive survey of non-canonical splice sites in the human transcriptome

We uncovered the diversity of non-canonical splice sites at the human transcriptome using deep transcriptome profiling. We mapped a total of 3.7 billion human RNA-seq reads and developed a set of stringent filters to avoid false non-canonical splice site detections. We identified 184 splice sites with non-canonical dinucleotides and U2/U12-like consensus sequences. We selected 10 of the herein identified U2/U12-like non-canonical splice site events and successfully validated 9 of them via reverse transcriptase-polymerase chain reaction and Sanger sequencing. Analyses of the 184 U2/U12-like non-canonical splice sites indicate that 51% of them are not annotated in GENCODE. In addition, 28% of them are conserved in mouse and 76% are involved in alternative splicing events, some of them with tissue-specific alternative splicing patterns. Interestingly, our analysis identified some U2/U12-like non-canonical splice sites that are converted into canonical splice sites by RNA A-to-I editing. Moreover, the U2/U12-like non-canonical splice sites have a differential distribution of splicing regulatory sequences, which may contribute to their recognition and regulation. Our analysis provides a high-confidence group of U2/U12-like non-canonical splice sites, which exhibit distinctive features among the total human splice sites.

The complex intron landscape and massive intron invasion in a picoeukaryote provides insights into intron evolution

Genes in pieces and spliceosomal introns are a landmark of eukaryotes, with intron invasion usually assumed to have happened early on in evolution. Here, we analyze the intron landscape of Micromonas, a unicellular green alga in the Mamiellophyceae lineage, demonstrating the coexistence of several classes of introns and the occurrence of recent massive intron invasion. This study focuses on two strains, CCMP1545 and RCC299, and their related individuals from ocean samplings, showing that they not only harbor different classes of introns depending on their location in the genome, as for other Mamiellophyceae, but also uniquely carry several classes of repeat introns. These introns, dubbed introner elements (IEs), are found at novel positions in genes and have conserved sequences, contrary to canonical introns. This IE invasion has a huge impact on the genome, doubling the number of introns in the CCMP1545 strain. We hypothesize that each IE class originated from a single ancestral IE that has been colonizing the genome after strain divergence by inserting copies of itself into genes by intron transposition, likely involving reverse splicing. Along with similar cases recently observed in other organisms, our observations in Micromonas strains shed a new light on the evolution of introns, suggesting that intron gain is more widespread than previously thought.

Physiological state co-regulates thousands of mammalian mRNA splicing events at tandem splice sites and alternative exons

Thousands of tandem alternative splice sites (TASS) give rise to mRNA insertion/deletion variants with small size differences. Recent work has concentrated on the question of biological relevance in general, and the physiological regulation of TASS in particular. We have quantitatively studied 11 representative TASS cases in comparison to one mutually exclusive exon case and two cassette exons (CEs) using a panel of human and mouse tissues, as well as cultured cell lines. Tissues show small but significant differences in TASS isoform ratios, with a variance 4- to 20-fold lower than seen for CEs. Remarkably, in cultured cells, all studied alternative splicing (AS) cases showed a cell-density-dependent shift of isoform ratios with similar time series profiles. A respective genome-wide co-regulation of TASS splicing was shown by next-generation mRNA sequencing data. Moreover, data from human and mouse organs indicate that this co-regulation of TASS occurs in vivo, with brain showing the strongest difference to other organs. Together, the results indicate a physiological AS regulation mechanism that functions almost independently from the splice site context and sequence.

Chromosome 10q26 locus and age-related macular degeneration: a progress update

Age-related macular degeneration (AMD) is the leading cause of late-onset central vision loss in developed countries. Both genetic and environmental factors contribute to the onset of AMD. Variation at a locus on chromosome 10q26 has been consistently associated with this disease and represents one of the two strongest genetic effects being identified in AMD. At least three genes are located within the bounds of the locus: pleckstrin homology domain containing family A member 1 (PLEKHA1), age-related maculopathy susceptibility 2 (ARMS2) and high-temperature requirement A serine peptidase 1 (HTRA1), all of which are associated with AMD. Due to the strong linkage disequilibrium (LD) across this region, statistical genetic analysis alone is incapable of distinguishing the effect of an individual gene in the locus. Uncertainty remains, however, in regards to which gene is responsible for the linkage and association of the locus with AMD. Investigating functional consequences of the associated variants and related genes tends to be essential to identifying the biologically responsible gene(s) underlying AMD. This review examines the recent progress and current uncertainty on the genetic and functional analyses of the 10q26 locus in AMD with a focus on ARMS2 and HTRA1. A discussion, which entails the possible multi-faceted approaches for pinpointing the gene(s) in the locus underlying the pathogenesis of AMD, is also included.

Intron retention in mRNA encoding ancillary subunit of insect voltage-gated sodium channel modulates channel expression, gating regulation and drug sensitivity

Insect voltage-gated sodium (Na_v) channels are formed by a well-known pore-forming α-subunit encoded by para-like gene and ancillary subunits related to TipE from the mutation “temperature-induced-paralysis locus E.” The role of these ancillary subunits in the modulation of biophysical and pharmacological properties of Na⁺ currents are not enough documented. The unique neuronal ancillary subunit TipE-homologous protein 1 of Drosophila melanogaster (DmTEH1) strongly enhances the expression of insect Na_v channels when heterologously expressed in Xenopus oocytes. Here we report the cloning and functional expression of two neuronal DmTEH1-homologs of the cockroach, Periplaneta americana, PaTEH1A and PaTEH1B, encoded by a single bicistronic gene. In PaTEH1B, the second exon encoding the last 11-amino-acid residues of PaTEH1A is shifted to 3′UTR by the retention of a 96-bp intron-containing coding-message, thus generating a new C-terminal end. We investigated the gating and pharmacological properties of the Drosophila Na_v channel variant (DmNa_v1-1) co-expressed with DmTEH1, PaTEH1A, PaTEH1B or a truncated mutant PaTEH1Δ(270-280) in Xenopus oocytes. PaTEH1B caused a 2.2-fold current density decrease, concomitant with an equivalent α-subunit incorporation decrease in the plasma membrane, compared to PaTEH1A and PaTEH1Δ(270-280). PaTEH1B positively shifted the voltage-dependences of activation and slow inactivation of DmNa_v1-1 channels to more positive potentials compared to PaTEH1A, suggesting that the C-terminal end of both proteins may influence the function of the voltage-sensor and the pore of Na_v channel. Interestingly, our findings showed that the sensitivity of DmNa_v1-1 channels to lidocaine and to the pyrazoline-type insecticide metabolite DCJW depends on associated TEH1-like subunits. In conclusion, our work demonstrates for the first time that density, gating and pharmacological properties of Na_v channels expressed in Xenopus oocytes can be modulated by an intron retention process in the transcription of the neuronal TEH1-like ancillary subunits of P. americana.

Characterization of two closely related α-amylase paralogs in the bark beetle, Ips typographus (L.)

Ips typographus (L.), the eight-spined spruce bark beetle, causes severe damage throughout Eurasian spruce forests and suitable nuclear markers are needed in order to study its population structure on a genetic level. Two closely related genes encoding α-amylase in I. typographus were characterized and named AmyA and AmyB. Both α-amylase paralogs consisted of six exons and five introns. AmyA encodes a polypeptide of 483 amino acids, whereas AmyB has two alternative transcripts encoding polypeptides of 483 and 370 amino acids. The expression levels of both genes were high during larval stage and adulthood. The AmyB transcripts were absent in the pupal stage. A modification of the allozyme staining method allowed us to detect two clusters of bands on the electrophoretic gel that may correspond to the two α-amylase genes. There was a correlation between the lack of AmyB expression in pupa and the absence of the fast migrating isozyme cluster at this stage, suggesting that the faster migrating isoforms are products of the AmyB gene, whereas the slowly migrating bands are derived from the AmyA.

Constant splice-isoform ratios in human lymphoblastoid cells support the concept of a splico-stat

Splicing generates mature transcripts from genes in pieces in eukaryotic cells. Overwhelming evidence has accumulated that alternative routes in splicing are possible for most human and mammalian genes, thereby allowing formation of different transcripts from one gene. No function has been assigned to the majority of identified alternative splice forms, and it has been assumed that they compose inert or tolerated waste from aberrant or noisy splicing. Here we demonstrate that five human transcription units (WT1, NOD2, GNAS, RABL2A, RABL2B) have constant splice-isoform ratios in genetically diverse lymphoblastoid cell lines independent of the type of alternative splicing (exon skipping, alternative donor/acceptor, tandem splice sites) and gene expression level. Even splice events that create premature stop codons and potentially trigger nonsense-mediated mRNA decay are found at constant fractions. The analyzed alternative splicing events were qualitatively but not quantitatively conserved in corresponding chimpanzee cell lines. Additionally, subtle splicing at tandem acceptor splice sites (GNAS, RABL2A/B) was highly constrained and strongly depends on the upstream donor sequence content. These results also demonstrate that unusual and unproductive splice variants are produced in a regulated manner.

TassDB2 - A comprehensive database of subtle alternative splicing events

BACKGROUND

Subtle alternative splicing events involving tandem splice sites separated by a short (2-12 nucleotides) distance are frequent and evolutionarily widespread in eukaryotes, and a major contributor to the complexity of transcriptomes and proteomes. However, these events have been either omitted altogether in databases on alternative splicing, or only the cases of experimentally confirmed alternative splicing have been reported. Thus, a database which covers all confirmed cases of subtle alternative splicing as well as the numerous putative tandem splice sites (which might be confirmed once more transcript data becomes available), and allows to search for tandem splice sites with specific features and download the results, is a valuable resource for targeted experimental studies and large-scale bioinformatics analyses of tandem splice sites. Towards this goal we recently set up TassDB (Tandem Splice Site DataBase, version 1), which stores data about alternative splicing events at tandem splice sites separated by 3 nt in eight species.

DESCRIPTION

We have substantially revised and extended TassDB. The currently available version 2 contains extensive information about tandem splice sites separated by 2-12 nt for the human and mouse transcriptomes including data on the conservation of the tandem motifs in five vertebrates. TassDB2 offers a user-friendly interface to search for specific genes or for genes containing tandem splice sites with specific features as well as the possibility to download result datasets. For example, users can search for cases of alternative splicing where the proportion of EST/mRNA evidence supporting the minor isoform exceeds a specific threshold, or where the difference in splice site scores is specified by the user. The predicted impact of each event on the protein is also reported, along with information about being a putative target for the nonsense-mediated decay (NMD) pathway. Links are provided to the UCSC genome browser and other external resources.

CONCLUSION

TassDB2, available via http://www.tassdb.info, provides comprehensive resources for researchers interested in both targeted experimental studies and large-scale bioinformatics analyses of short distance tandem splice sites.

Comparison of methods for quantification of subtle splice variants

Alternative splicing is capable of generating multiple mRNA variants from a single gene and is hence a key mediator of molecular diversity generated at the transcript level. Consequently, delivering quantitative information on the fractions of splice variants is essential for the understanding of their biological roles. Here we compare techniques for subtle splice variant quantification that are able to resolve length differences as small as one nucleotide: PAGE with ethidium-bromide densitometry, pyrosequencing, and CE-LIF. We give comprehensive descriptions of assay designs and calibration procedures and present an evaluation of these methods in terms of accuracy, reproducibility and applicability. We also examined template concentrations and reverse transcription-coupled PCR conditions as potential cause of biased results as they were observed for extreme low template concentrations and/or PCR amplicons with size differences of 195 nt. As proof of concept, we determine the splice ratios of variants differing by 3 and 12 nt in five human tissues. We demonstrate that CE-LIF is the most precise and also the most labor- and time-efficient method.

Analysis of splicing patterns by pyrosequencing

Several different mRNAs can be produced from a given pre-mRNA by regulated alternative splicing, or as the result of deregulations that may lead to pathological states. Analysing splicing patterns is therefore of importance to describe and understand developmental programs, cellular responses to internal or external cues, or human diseases. We describe here a method, Pyrosequencing Analysis of Splicing Patterns (PASP), that combines RT–PCR and pyrosequencing of PCR products. We demonstrated that: (i) Ratios of two pure RNAs mixed in various proportions were accurately measured by PASP; (ii) PASP can be adapted to virtually any splicing event, including mutually exclusive exons, complex patterns of exon skipping or inclusion, and alternative 3′ terminal exons; (iii) In extracts from different organs, the proportions of RNA isoforms measured by PASP reflected those measured by other methods. The PASP method is therefore reliable for analysing splicing patterns. All steps are done in 96-wells microplates, without gel electrophoresis, opening the way to high-throughput comparisons of RNA from several sources.

Accurate prediction of NAGNAG alternative splicing

Alternative splicing (AS) involving NAGNAG tandem acceptors is an evolutionarily widespread class of AS. Recent predictions of alternative acceptor usage reported better results for acceptors separated by larger distances, than for NAGNAGs. To improve the latter, we aimed at the use of Bayesian networks (BN), and extensive experimental validation of the predictions. Using carefully constructed training and test datasets, a balanced sensitivity and specificity of ≥92% was achieved. A BN trained on the combined dataset was then used to make predictions, and 81% (38/47) of the experimentally tested predictions were verified. Using a BN learned on human data on six other genomes, we show that while the performance for the vertebrate genomes matches that achieved on human data, there is a slight drop for Drosophila and worm. Lastly, using the prediction accuracy according to experimental validation, we estimate the number of yet undiscovered alternative NAGNAGs. State of the art classifiers can produce highly accurate prediction of AS at NAGNAGs, indicating that we have identified the major features of the ‘NAGNAG-splicing code’ within the splice site and its immediate neighborhood. Our results suggest that the mechanism behind NAGNAG AS is simple, stochastic, and conserved among vertebrates and beyond.

Diverse splicing pathways of the membrane IgHM pre-mRNA in a Chondrostean, the Siberian sturgeon

Teleosts and tetrapods have evolved different splice patterns to generate their membrane-bound IgM. In the tetrapod lineage, the first transmembrane exon is spliced to an internal cryptic site located close to the end of the fourth constant exon. Because teleosts lack this site they use the regular 3'-splice site of the CH3 exon instead. We characterized the mum splicing patterns in a Chondrostean, the Siberian sturgeon. We observed a surprising diversity of splice patterns, the TM1 exon being spliced to a cryptic site at the end of CH4, to a cryptic site in CH3 or to the 3'-end of CH1. These different pathways lead to mIGHM transcripts encoding four, two or one complete C-domain(s), respectively. The short variant CH1-TM1 was found only in VH2 positive transcripts, while the two other variants were observed for IgHM transcripts expressing all VH families. These results shed light on the evolution of IgM splicing pathways.

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.

Selection against tandem splice sites affecting structured protein regions

BACKGROUND

Alternative selection of splice sites in tandem donors and acceptors is a major mode of alternative splicing. Here, we analyzed whether in-frame tandem sites leading to subtle mRNA insertions/deletions of 3, 6, or 9 nucleotides are under natural selection.

RESULTS

We found multiple lines of evidence that the human protein coding sequences are under selection against such in-frame tandem splice events, indicating that these events are often deleterious. The strength of selection is not homogeneous within the coding sequence as protein regions that fold into a fixed 3D structure (intrinsically ordered) are under stronger selection, especially against sites with a strong minor splice site. Investigating structures of functional protein domains, we found that tandem acceptors are preferentially located at the domain surface and outside structural elements such as helices and sheets. Using three-species comparisons, we estimate that more than half of all mutations that create NAGNAG acceptors in the coding region have been eliminated by selection.

CONCLUSION

We estimate that ~2,400 introns are under selection against possessing a tandem site.