A polar mechanism coordinates different regions of alternative splicing within a single gene

Understanding isoform expression by pairing long-read sequencing with single-cell and spatial transcriptomics

RNA isoform diversity, produced via alternative splicing, and alternative usage of transcription start and poly(A) sites, results in varied transcripts being derived from the same gene. Distinct isoforms can play important biological roles, including by changing the sequences or expression levels of protein products. The first single-cell approaches to RNA sequencing-and later, spatial approaches-which are now widely used for the identification of differentially expressed genes, rely on short reads and offer the ability to transcriptomically compare different cell types but are limited in their ability to measure differential isoform expression. More recently, long-read sequencing methods have been combined with single-cell and spatial technologies in order to characterize isoform expression. In this review, we provide an overview of the emergence of single-cell and spatial long-read sequencing and discuss the challenges associated with the implementation of these technologies and interpretation of these data. We discuss the opportunities they offer for understanding the relationships between the distinct variable elements of transcript molecules and highlight some of the ways in which they have been used to characterize isoforms' roles in development and pathology. Single-nucleus long-read sequencing, a special case of the single-cell approach, is also discussed. We attempt to cover both the limitations of these technologies and their significant potential for expanding our still-limited understanding of the biological roles of RNA isoforms.

ScISOr-ATAC reveals convergent and divergent splicing and chromatin specificities between matched cell types across cortical regions, evolution, and in Alzheimer's Disease

Multimodal measurements have become widespread in genomics, however measuring open chromatin accessibility and splicing simultaneously in frozen brain tissues remains unconquered. Hence, we devised Single-Cell-ISOform-RNA sequencing coupled with the Assay-for-Transposase-Accessible-Chromatin (ScISOr-ATAC). We utilized ScISOr-ATAC to assess whether chromatin and splicing alterations in the brain convergently affect the same cell types or divergently different ones. We applied ScISOr-ATAC to three major conditions: comparing (i) the Rhesus macaque (Macaca mulatta) prefrontal cortex (PFC) and visual cortex (VIS), (ii) cross species divergence of Rhesus macaque versus human PFC, as well as (iii) dysregulation in Alzheimer's disease in human PFC. We found that among cortical-layer biased excitatory neuron subtypes, splicing is highly brain-region specific for L3-5/L6 IT_RORB neurons, moderately specific in L2-3 IT_CUX2.RORB neurons and unspecific in L2-3 IT_CUX2 neurons. In contrast, at the chromatin level, L2-3 IT_CUX2.RORB neurons show the highest brain-region specificity compared to other subtypes. Likewise, when comparing human and macaque PFC, strong evolutionary divergence on one molecular modality does not necessarily imply strong such divergence on another molecular level in the same cell type. Finally, in Alzheimer's disease, oligodendrocytes show convergently high dysregulation in both chromatin and splicing. However, chromatin and splicing dysregulation most strongly affect distinct oligodendrocyte subtypes. Overall, these results indicate that chromatin and splicing can show convergent or divergent results depending on the performed comparison, justifying the need for their concurrent measurement to investigate complex systems. Taken together, ScISOr-ATAC allows for the characterization of single-cell splicing and chromatin patterns and the comparison of sample groups in frozen brain samples.

From words to complete phrases: insight into single-cell isoforms using short and long reads

The profiling of gene expression patterns to glean biological insights from single cells has become commonplace over the last few years. However, this approach overlooks the transcript contents that can differ between individual cells and cell populations. In this review, we describe early work in the field of single-cell short-read sequencing as well as full-length isoforms from single cells. We then describe recent work in single-cell long-read sequencing wherein some transcript elements have been observed to work in tandem. Based on earlier work in bulk tissue, we motivate the study of combination patterns of other RNA variables. Given that we are still blind to some aspects of isoform biology, we suggest possible future avenues such as CRISPR screens which can further illuminate the function of RNA variables in distinct cell populations.

Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue

Single-nuclei RNA sequencing characterizes cell types at the gene level. However, compared to single-cell approaches, many single-nuclei cDNAs are purely intronic, lack barcodes and hinder the study of isoforms. Here we present single-nuclei isoform RNA sequencing (SnISOr-Seq). Using microfluidics, PCR-based artifact removal, target enrichment and long-read sequencing, SnISOr-Seq increased barcoded, exon-spanning long reads 7.5-fold compared to naive long-read single-nuclei sequencing. We applied SnISOr-Seq to adult human frontal cortex and found that exons associated with autism exhibit coordinated and highly cell-type-specific inclusion. We found two distinct combination patterns: those distinguishing neural cell types, enriched in TSS-exon, exon-polyadenylation-site and non-adjacent exon pairs, and those with multiple configurations within one cell type, enriched in adjacent exon pairs. Finally, we observed that human-specific exons are almost as tightly coordinated as conserved exons, implying that coordination can be rapidly established during evolution. SnISOr-Seq enables cell-type-specific long-read isoform analysis in human brain and in any frozen or hard-to-dissociate sample.

Efficient RNA polymerase II pause release requires U2 snRNP function

Transcription by RNA polymerase II (Pol II) is coupled to pre-mRNA splicing, but the underlying mechanisms remain poorly understood. Co-transcriptional splicing requires assembly of a functional spliceosome on nascent pre-mRNA, but whether and how this influences Pol II transcription remains unclear. Here we show that inhibition of pre-mRNA branch site recognition by the spliceosome component U2 snRNP leads to a widespread and strong decrease in new RNA synthesis from human genes. Multiomics analysis reveals that inhibition of U2 snRNP function increases the duration of Pol II pausing in the promoter-proximal region, impairs recruitment of the pause release factor P-TEFb, and reduces Pol II elongation velocity at the beginning of genes. Our results indicate that efficient release of paused Pol II into active transcription elongation requires the formation of functional spliceosomes and that eukaryotic mRNA biogenesis relies on positive feedback from the splicing machinery to the transcription machinery.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

The Features and Regulation of Co-transcriptional Splicing in Arabidopsis

Precursor mRNA (pre-mRNA) splicing is essential for gene expression in most eukaryotic organisms. Previous studies from mammals, Drosophila, and yeast show that the majority of splicing events occurs co-transcriptionally. In plants, however, the features of co-transcriptional splicing (CTS) and its regulation still remain largely unknown. Here, we used chromatin-bound RNA sequencing to study CTS in Arabidopsis thaliana. We found that CTS is widespread in Arabidopsis seedlings, with a large proportion of alternative splicing events determined co-transcriptionally. CTS efficiency correlated with gene expression level, the chromatin landscape and, most surprisingly, the number of introns and exons of individual genes, but is independent of gene length. In combination with enhanced crosslinking and immunoprecipitation sequencing analysis, we further showed that the hnRNP-like proteins RZ-1B and RZ-1C promote efficient CTS globally through direct binding, frequently to exonic sequences. Notably, this general effect of RZ-1B/1C on splicing promotion is mainly observed at the chromatin level, not at the mRNA level. RZ-1C promotes CTS of multiple-exon genes in association with its binding to regions both proximal and distal to the regulated introns. We propose that RZ-1C promotes efficient CTS of genes with multiple exons through cooperative interactions with many exons, introns, and splicing factors. Our work thus reveals important features of CTS in plants and provides methodologies for the investigation of CTS and RNA-binding proteins in plants.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

Getting the Entire Message: Progress in Isoform Sequencing

The advent of second-generation sequencing and its application to RNA sequencing have revolutionized the field of genomics by allowing quantification of gene expression, as well as the definition of transcription start/end sites, exons, splice sites and RNA editing sites. However, due to the sequencing of fragments of cDNAs, these methods have not given a reliable picture of complete RNA isoforms. Third-generation sequencing has filled this gap and allows end-to-end sequencing of entire RNA/cDNA molecules. This approach to transcriptomics has been a "niche" technology for a couple of years but now is becoming mainstream with many different applications. Here, we review the background and progress made to date in this rapidly growing field. We start by reviewing the progressive realization that alternative splicing is omnipresent. We then focus on long-noncoding RNA isoforms and the distinct combination patterns of exons in noncoding and coding genes. We consider the implications of the recent technologies of direct RNA sequencing and single-cell isoform RNA sequencing. Finally, we discuss the parameters that define the success of long-read RNA sequencing experiments and strategies commonly used to make the most of such data.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells

Full-length RNA sequencing (RNA-Seq) has been applied to bulk tissue, cell lines and sorted cells to characterize transcriptomes, but applying this technology to single cells has proven to be difficult, with less than ten single-cell transcriptomes having been analyzed thus far. Although single splicing events have been described for ≤200 single cells with statistical confidence, full-length mRNA analyses for hundreds of cells have not been reported. Single-cell short-read 3' sequencing enables the identification of cellular subtypes, but full-length mRNA isoforms for these cell types cannot be profiled. We developed a method that starts with bulk tissue and identifies single-cell types and their full-length RNA isoforms without fluorescence-activated cell sorting. Using single-cell isoform RNA-Seq (ScISOr-Seq), we identified RNA isoforms in neurons, astrocytes, microglia, and cell subtypes such as Purkinje and Granule cells, and cell-type-specific combination patterns of distant splice sites. We used ScISOr-Seq to improve genome annotation in mouse Gencode version 10 by determining the cell-type-specific expression of 18,173 known and 16,872 novel isoforms.

Long-read sequencing of nascent RNA reveals coupling among RNA processing events

Pre-mRNA splicing is accomplished by the spliceosome, a megadalton complex that assembles de novo on each intron. Because spliceosome assembly and catalysis occur cotranscriptionally, we hypothesized that introns are removed in the order of their transcription in genomes dominated by constitutive splicing. Remarkably little is known about splicing order and the regulatory potential of nascent transcript remodeling by splicing, due to the limitations of existing methods that focus on analysis of mature splicing products (mRNAs) rather than substrates and intermediates. Here, we overcome this obstacle through long-read RNA sequencing of nascent, multi-intron transcripts in the fission yeast Schizosaccharomyces pombe. Most multi-intron transcripts were fully spliced, consistent with rapid cotranscriptional splicing. However, an unexpectedly high proportion of transcripts were either fully spliced or fully unspliced, suggesting that splicing of any given intron is dependent on the splicing status of other introns in the transcript. Supporting this, mild inhibition of splicing by a temperature-sensitive mutation in prp2, the homolog of vertebrate U2AF65, increased the frequency of fully unspliced transcripts. Importantly, fully unspliced transcripts displayed transcriptional read-through at the polyA site and were degraded cotranscriptionally by the nuclear exosome. Finally, we show that cellular mRNA levels were reduced in genes with a high number of unspliced nascent transcripts during caffeine treatment, showing regulatory significance of cotranscriptional splicing. Therefore, overall splicing of individual nascent transcripts, 3′ end formation, and mRNA half-life depend on the splicing status of neighboring introns, suggesting crosstalk among spliceosomes and the polyA cleavage machinery during transcription elongation.

Microfluidic isoform sequencing shows widespread splicing coordination in the human transcriptome

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k–200k partitions of 10–200 molecules at a time, enabling analysis of 10–100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT. Coordination of internal exons is found to be extensive among protein coding genes: 23.5%–59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein–protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Regulation of alternative splicing through coupling with transcription and chromatin structure

Alternative precursor messenger RNA (pre-mRNA) splicing plays a pivotal role in the flow of genetic information from DNA to proteins by expanding the coding capacity of genomes. Regulation of alternative splicing is as important as regulation of transcription to determine cell- and tissue-specific features, normal cell functioning, and responses of eukaryotic cells to external cues. Its importance is confirmed by the evolutionary conservation and diversification of alternative splicing and the fact that its deregulation causes hereditary disease and cancer. This review discusses the multiple layers of cotranscriptional regulation of alternative splicing in which chromatin structure, DNA methylation, histone marks, and nucleosome positioning play a fundamental role in providing a dynamic scaffold for interactions between the splicing and transcription machineries. We focus on evidence for how the kinetics of RNA polymerase II (RNAPII) elongation and the recruitment of splicing factors and adaptor proteins to chromatin components act in coordination to regulate alternative splicing.

Ratios of Four STAT3 Splice Variants in Human Eosinophils and Diffuse Large B Cell Lymphoma Cells

Signal transducer and activator of transcription 3 (STAT3) is a key mediator of leukocyte differentiation and proliferation. The 3' end of STAT3 transcripts is subject to two alternative splicing events. One results in either full-length STAT3α or in STAT3β, which lacks part of the C-terminal transactivation domain. The other is at a tandem donor (5') splice site and results in the codon for Ser-701 being included (S) or excluded (ΔS). Despite the proximity of Ser-701 to the site of activating phosphorylation at Tyr-705, ΔS/S splicing has barely been studied. Sequencing of cDNA from purified eosinophils revealed the presence of four transcripts (S-α, ΔS-α, S-β, and ΔS-β) rather than the three reported in publically available databases from which ΔS-β is missing. To gain insight into regulation of the two alternative splicing events, we developed a quantitative(q) PCR protocol to compare transcript ratios in eosinophils in which STAT3 is upregulated by cytokines, activated B cell diffuse large B cell Lymphoma (DLBCL) cells in which STAT3 is dysregulated, and in germinal center B cell-like DLBCL cells in which it is not. With the exception of one line of activated B cell DLCBL cells, the four variants were found in roughly the same ratios despite differences in total levels of STAT3 transcripts. S-α was the most abundant, followed by S-β. ΔS-α and ΔS-β together comprised 15.6±4.0 % (mean±SD, n=21) of the total. The percentage of STAT3β variants that were ΔS was 1.5-fold greater than of STAT3α variants that were ΔS. Inspection of Illumina’s “BodyMap” RNA-Seq database revealed that the ΔS variant accounts for 10-26 % of STAT3 transcripts across 16 human tissues, with less variation than three other genes with the identical tandem donor splice site sequence. Thus, it seems likely that all cells contain the S-α, ΔS-α, S-β, and ΔS-β variants of STAT3.

Assessing long-distance RNA sequence connectivity via RNA-templated DNA-DNA ligation

Many RNAs, including pre-mRNAs and long non-coding RNAs, can be thousands of nucleotides long and undergo complex post-transcriptional processing. Multiple sites of alternative splicing within a single gene exponentially increase the number of possible spliced isoforms, with most human genes currently estimated to express at least ten. To understand the mechanisms underlying these complex isoform expression patterns, methods are needed that faithfully maintain long-range exon connectivity information in individual RNA molecules. In this study, we describe SeqZip, a methodology that uses RNA-templated DNA–DNA ligation to retain and compress connectivity between distant sequences within single RNA molecules. Using this assay, we test proposed coordination between distant sites of alternative exon utilization in mouse Fn1, and we characterize the extraordinary exon diversity of Drosophila melanogaster Dscam1.

DOI:http://dx.doi.org/10.7554/eLife.03700.001

A flow chart can show how an outcome can be achieved from a particular start point by breaking down an activity into a list of possible steps. Often, a flow chart contains several alternative steps, not all of which are taken every time the flow chart is used. The same can be said of genes, which are biological instructions that often contain many options within their DNA sequences.

Proteins—which perform many roles in cells—are built following the instructions contained in genes. First, the DNA sequence of the gene is copied. This produces a molecule of ribonucleic acid (RNA), which is able to move around the cell to find the machinery that can use the genetic information to make a protein. Genes and their RNA copies contain instructions with more steps—called exons—than are necessary to make a working protein, so extra exons are removed (‘spliced’) from the RNA copies. Different combinations of exons can be removed, so splicing can make different versions of the RNA called isoforms. These allow a single gene to build many different proteins. In fruit flies, for example, the different exons of the gene Dscam1 can be spliced into one of 38,016 unique RNA isoforms.

Current technology only allows researchers to deduce the sequence of RNA molecules by combining sequences recorded from short fragments of the molecule. However, before splicing, RNA molecules tend to be much longer than this, so this restricts our understanding of the RNA isoforms found in cells. Here, Roy et al. devised and tested a new method called SeqZip to solve this problem.

SeqZip uses short fragments of DNA called ligamers that can only stick to the sections of RNA that will remain after the molecule has been spliced. After splicing, the ligamers can be stuck together to make a DNA replica of the spliced RNA. The end product is at least 49 times shorter than the original RNA, so it is easier to sequence. In addition, the combinations of the ligamers in the DNA replica show which exons of a specific gene are kept and which ones are spliced out.

To test the method, Roy et al. studied a mouse gene that has six RNA isoforms. SeqZip reduced the length of the RNA by five times and made it possible to measure how frequently the different isoforms naturally arise. Roy et al. also used SeqZip to work out which isoforms of the Dscam1 gene are used at different stages in the life of fruit fly larvae. SeqZip can provide insights into how complex organisms like flies, mice, and humans have evolved with relatively few—a little over 20,000—genes in their genomes.

DOI:http://dx.doi.org/10.7554/eLife.03700.002

Quantification of pre-mRNA escape rate and synergy in splicing

Splicing reactions generally combine high speed with accuracy. However, some of the pre-mRNAs escape the nucleus with a retained intron. Intron retention can control gene expression and increase proteome diversity. We calculated the escape rate for the yeast PTC7 intron and pre-mRNA. This prediction was facilitated by the observation that splicing is a linear process and by deriving simple algebraic expressions from a model of co- and post-transcriptional splicing and RNA surveillance that determines the rate of the nonsense-mediated decay (NMD) of the pre-mRNAs with the retained intron. The escape rate was consistent with the observed threshold of splicing rate below which the mature mRNA level declined. When an mRNA contains multiple introns, the outcome of splicing becomes more difficult to predict since not only the escape rate of the pre-mRNA has to be considered, but also the possibility that the splicing of each intron is influenced by the others. We showed that the two adjacent introns in the SUS1 mRNA are spliced cooperatively, but this does not counteract the escape of the partially spliced mRNA. These findings will help to infer promoter activity and to predict the behavior of and to control splicing regulatory networks.

Experimental assessment of splicing variants using expression minigenes and comparison with in silico predictions

Assessment of the functional consequences of variants near splice sites is a major challenge in the diagnostic laboratory. To address this issue, we created expression minigenes (EMGs) to determine the RNA and protein products generated by splice site variants (n = 10) implicated in cystic fibrosis (CF). Experimental results were compared with the splicing predictions of eight in silico tools. EMGs containing the full-length Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) coding sequence and flanking intron sequences generated wild-type transcript and fully processed protein in Human Embryonic Kidney (HEK293) and CF bronchial epithelial (CFBE41o-) cells. Quantification of variant induced aberrant mRNA isoforms was concordant using fragment analysis and pyrosequencing. The splicing patterns of c.1585-1G>A and c.2657+5G>A were comparable to those reported in primary cells from individuals bearing these variants. Bioinformatics predictions were consistent with experimental results for 9/10 variants (MES), 8/10 variants (NNSplice), and 7/10 variants (SSAT and Sroogle). Programs that estimate the consequences of mis-splicing predicted 11/16 (HSF and ASSEDA) and 10/16 (Fsplice and SplicePort) experimentally observed mRNA isoforms. EMGs provide a robust experimental approach for clinical interpretation of splice site variants and refinement of in silico tools.

The exon junction complex is required for definition and excision of neighboring introns in Drosophila

The exon junction complex (EJC) plays a central role in controlling RNA fate and aids faithful splicing of pre-mRNAs containing large introns via an unknown mechanism. Brennecke and colleagues show that the core EJC plus the accessory factors RnpS1 and Acinus aid in the definition and efficient splicing of neighboring introns. Interestingly, the most highly affected intron belongs to the piwi locus, which explains the reported transposon desilencing in EJC-depleted Drosophila ovaries. Based on transcriptome-wide analysis, the authors propose that the dependency of splicing on the EJC is exploited to control the temporal order of splicing events.

Splicing of pre-mRNAs results in the deposition of the exon junction complex (EJC) upstream of exon–exon boundaries. The EJC plays crucial post-splicing roles in export, translation, localization, and nonsense-mediated decay of mRNAs. It also aids faithful splicing of pre-mRNAs containing large introns, albeit via an unknown mechanism. Here, we show that the core EJC plus the accessory factors RnpS1 and Acinus aid in definition and efficient splicing of neighboring introns. This requires prior deposition of the EJC in close proximity to either an upstream or downstream splicing event. If present in isolation, EJC-dependent introns are splicing-defective also in wild-type cells. Interestingly, the most affected intron belongs to the piwi locus, which explains the reported transposon desilencing in EJC-depleted Drosophila ovaries. Based on a transcriptome-wide analysis, we propose that the dependency of splicing on the EJC is exploited as a means to control the temporal order of splicing events.

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

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

The multiplicity of alternative splicing decisions in Caenorhabditis elegans is linked to specific intronic regulatory motifs and minisatellites

Background

Alternative splicing diversifies the pool of messenger RNA molecules encoded by individual genes. This diversity is particularly high when multiple splicing decisions cause a combinatorial arrangement of several alternate exons. We know very little on how the multiple decisions occurring during the maturation of single transcripts are coordinated and whether specific sequence elements might be involved.

Results

Here, the Caenorhabditis elegans genome was surveyed in order to identify sequence elements that might play a specific role in the regulation of multiple splicing decisions. The introns flanking alternate exons in transcripts whose maturation involves multiple alternative splicing decisions were compared to those whose maturation involves a single decision. Fifty-eight penta-, hexa-, and hepta-meric elements, clustered in 17 groups, were significantly over-represented in genes subject to multiple alternative splicing decisions. Most of these motifs relate to known splicing regulatory elements and appear to be well conserved in the related species Caenorhabditis briggsae. The usage of specific motifs is not linked to the gene product function, but rather depends on the gene structure, since it is influenced by the distance separating the multiple splicing decision sites. Two of these motifs are part of the CeRep25B minisatellite, which is also over-represented at the vicinity of alternative splicing regions. Most of the remaining motifs are not part of repeated sequence elements, but tend to occur in specific heterologous pairs in genes subject to multiple alternative splicing decisions.

Conclusions

The existence of specific intronic sequence elements linked to multiple alternative splicing decisions is intriguing and suggests that these elements might have some specialized regulatory role during splicing.

Electronic supplementary material

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

TMEM16A alternative splicing coordination in breast cancer

Background

TMEM16A, also known as Anoctamin-1, is a calcium-activated chloride channel gene overexpressed in many tumors. The role of TMEM16A in cancer is not completely understood and no data are available regarding the potential tumorigenic properties of the multiple isoforms generated by alternative splicing (AS).

Methods

We evaluated TMEM16A AS pattern, isoforms distribution and Splicing Coordination (SC), in normal tissues and breast cancers, through a semi-quantitative PCR-assay that amplifies transcripts across three AS exons, 6b, 13 and 15.

Results

In breast cancer, we did not observe an association either to AS of individual exons or to specific TMEM16A isoforms, and induced expression of the most common isoforms present in tumors in the HEK293 Flp-In Tet-ON system had no effect on cellular proliferation and migration. The analysis of splicing coordination, a mechanism that regulates AS of distant exons, showed a preferential association of exon 6b and 15 in several normal tissues and tumors: isoforms that predominantly include exon 6b tend to exclude exon 15 and vice versa. Interestingly, we found an increase in SC in breast tumors compared to matched normal tissues.

Conclusions

As the different TMEM16A isoforms do not affect proliferation or migration and do not associate with tumors, our results suggest that the resulting channel activities are not directly involved in cell growth and motility. Conversely, the observed increase in SC in breast tumors suggests that the maintenance of the regulatory mechanism that coordinates distant alternative spliced exons in multiple genes other than TMEM16A is necessary for cancer cell viability.

Inference of alternative splicing from RNA-Seq data with probabilistic splice graphs

Motivation: Alternative splicing and other processes that allow for different transcripts to be derived from the same gene are significant forces in the eukaryotic cell. RNA-Seq is a promising technology for analyzing alternative transcripts, as it does not require prior knowledge of transcript structures or genome sequences. However, analysis of RNA-Seq data in the presence of genes with large numbers of alternative transcripts is currently challenging due to efficiency, identifiability and representation issues.

Results: We present RNA-Seq models and associated inference algorithms based on the concept of probabilistic splice graphs, which alleviate these issues. We prove that our models are often identifiable and demonstrate that our inference methods for quantification and differential processing detection are efficient and accurate.

Availability: Software implementing our methods is available at http://deweylab.biostat.wisc.edu/psginfer.

Contact:cdewey@biostat.wisc.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Coupling between transcription and alternative splicing

The scenario of alternative splicing regulation is far more complex than the classical picture of a pre-mRNA being processed post-transcriptionally in more than one way. Introns are efficiently removed while transcripts are still being synthesized, supporting the idea of a co-transcriptional regulation of alternative splicing. Evidence of a functional coupling between splicing and transcription has recently emerged as it was observed that properties of one process may affect the outcome of the other. Co-transcriptionality is thought to improve splicing efficiency and kinetics by directing the nascent pre-mRNA into proper spliceosome assembly and favoring splicing factor recruitment. Two models have been proposed to explain the coupling of transcription and alternative splicing: in the recruitment model, promoters and pol II status affect the recruitment to the transcribing gene of splicing factors or bifunctional factors acting on both transcription and splicing; in the kinetic model, differences in the elongation rate of pol II would determine the timing in which splicing sites are presented, and thus the outcome of alternative splicing decisions. In the later model, chromatin structure has emerged as a key regulator. Although definitive evidence for transcriptionally coupled alternative splicing alterations in tumor development or cancer pathogenesis is still missing, many alternative splicing events altered in cancer might be subject to transcription-splicing coupling regulation.

Cancer-Associated Perturbations in Alternative Pre-messenger RNA Splicing

For most of our 25,000 genes, the removal of introns by pre-messenger RNA (pre-mRNA) splicing represents an essential step toward the production of functional messenger RNAs (mRNAs). Alternative splicing of a single pre-mRNA results in the production of different mRNAs. Although complex organisms use alternative splicing to expand protein function and phenotypic diversity, patterns of alternative splicing are often altered in cancer cells. Alternative splicing contributes to tumorigenesis by producing splice isoforms that can stimulate cell proliferation and cell migration or induce resistance to apoptosis and anticancer agents. Cancer-specific changes in splicing profiles can occur through mutations that are affecting splice sites and splicing control elements, and also by alterations in the expression of proteins that control splicing decisions. Recent progress in global approaches that interrogate splicing diversity should help to obtain specific splicing signatures for cancer types. The development of innovative approaches for annotating and reprogramming splicing events will more fully establish the essential contribution of alternative splicing to the biology of cancer and will hopefully provide novel targets and anticancer strategies. Metazoan genes are usually made up of several exons interrupted by introns. The introns are removed from the pre-mRNA by RNA splicing. In conjunction with other maturation steps, such as capping and polyadenylation, the spliced mRNA is then transported to the cytoplasm to be translated into a functional protein. The basic mechanism of splicing requires accurate recognition of each extremity of each intron by the spliceosome. Introns are identified by the binding of U1 snRNP to the 5' splice site and the U2AF65/U2AF35 complex to the 3' splice site. Following these interactions, other proteins and snRNPs are recruited to generate the complete spliceosomal complex needed to excise the intron. While many introns are constitutively removed by the spliceosome, other splice junctions are not used systematically, generating the phenomenon of alternative splicing. Alternative splicing is therefore the process by which a single species of pre-mRNA can be matured to produce different mRNA molecules (Fig. 1). Depending on the number and types of alternative splicing events, a pre-mRNA can generate from two to several thousands different mRNAs leading to the production of a corresponding number of proteins. It is now believed that the expression of at least 70 % of human genes is subjected to alternative splicing, implying an enormous contribution to proteomic diversity, and by extension, to the development and the evolution of complex animals. Defects in splicing have been associated with human diseases (Caceres and Kornblihtt, Trends Genet 18(4):186-93, 2002, Cartegni et al., Nat Rev Genet 3(4):285-98, 2002, Pagani and Baralle, Nat Rev Genet 5(5):389-96, 2004), including cancer (Brinkman, Clin Biochem 37(7):584-94, 2004, Venables, Bioessays 28(4):378-86, 2006, Srebrow and Kornblihtt, J Cell Sci 119(Pt 13):2635-2641, 2006, Revil et al., Bull Cancer 93(9):909-919, 2006, Venables, Transworld Res Network, 2006, Pajares et al., Lancet Oncol 8(4):349-57, 2007, Skotheim and Nees, Int J Biochem Cell Biol 39:1432-1449, 2007). Numerous studies have now confirmed the existence of specific differences in the alternative splicing profiles between normal and cancer tissues. Although there are a few cases where specific mutations are the primary cause for these changes, global alterations in alternative splicing in cancer cells may be primarily derived from changes in the expression of RNA-binding proteins that control splice site selection. Overall, these cancer-specific differences in alternative splicing offer an immense potential to improve the diagnosis and the prognosis of cancer. This review will focus on the functional impact of cancer-associated alternative splicing variants, the molecular determinants that alter the splicing decisions in cancer cells, and future therapeutic strategies.

Regulated pre-mRNA splicing: the ghostwriter of the eukaryotic genome

Intron removal is at the heart of mRNA synthesis. It is mediated by one of the cell's largest complexes, the spliceosome. Yet, the fundamental chemistry involved is simple. In this review we will address how the spliceosome acts in diverse ways to optimize gene expression in order to meet the cell's needs. This is done largely by regulating the splicing of key transcripts encoding products that control gene expression pathways. This widespread role is evident even in the yeast Saccharomyces cerevisiae, where many introns appear to have been lost; yet how this control is being achieved is known only in a few cases. Here we explore the relevant examples and posit hypotheses whereby regulated splicing fine-tunes gene expression pathways to maintain cell homeostasis. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Intragenic alternative splicing coordination is essential for Caenorhabditis elegans slo-1 gene function

Alternative splicing is critical for diversifying eukaryotic proteomes, but the rules governing and coordinating splicing events among multiple alternate splice sites within individual genes are not well understood. We developed a quantitative PCR-based strategy to quantify the expression of the 12 transcripts encoded by the Caenorhabditis elegans slo-1 gene, containing three alternate splice sites. Using conditional probability-based models, we show that splicing events are coordinated across these sites. Further, we identify a point mutation in an intron adjacent to one alternate splice site that disrupts alternative splicing at all three sites. This mutation leads to aberrant synaptic transmission at the neuromuscular junction. In a genomic survey, we found that a UAAAUC element disrupted by this mutation is enriched in introns flanking alternate exons in genes with multiple alternate splice sites. These results establish that proper coordination of intragenic alternative splicing is essential for normal physiology of slo-1 in vivo and identify putative specialized cis-regulatory elements that regulate the coordination of intragenic alternative splicing.

Key features of the two-intron Saccharomyces cerevisiae gene SUS1 contribute to its alternative splicing

Alternative pre-mRNA splicing allows dramatic expansion of the eukaryotic proteome and facilitates cellular response to changes in environmental conditions. The Saccharomyces cerevisiae gene SUS1, which encodes a protein involved in mRNA export and histone H2B deubiquitination, contains two introns; non-canonical sequences in the first intron contribute to its retention, a common form of alternative splicing in plants and fungi. Here we show that the pattern of SUS1 splicing changes in response to environmental change such as temperature elevation, and the retained intron product is subject to nonsense-mediated decay. The activities of different splicing factors determine the pattern of SUS1 splicing, including intron retention and exon skipping. Unexpectedly, removal of the 3′ intron is affected by splicing of the upstream intron, suggesting that cross-exon interactions influence intron removal. Production of different SUS1 isoforms is important for cellular function, as we find that the temperature sensitivity and histone H2B deubiquitination defects observed in sus1Δ cells are only partially suppressed by SUS1 cDNA, but SUS1 that is able to undergo splicing complements these phenotypes. These data illustrate a role for S. cerevisiae alternative splicing in histone modification and cellular function and reveal important mechanisms for splicing of yeast genes containing multiple introns.

Competing RNA secondary structures are required for mutually exclusive splicing of the Dscam exon 6 cluster

Alternative splicing of eukaryotic pre-mRNAs is an important mechanism for generating proteome diversity and regulating gene expression. The Drosophila melanogaster Down Syndrome Cell Adhesion Molecule (Dscam) gene is an extreme example of mutually exclusive splicing. Dscam contains 95 alternatively spliced exons that potentially encode 38,016 distinct mRNA and protein isoforms. We previously identified two sets of conserved sequence elements, the docking site and selector sequences in the Dscam exon 6 cluster, which contains 48 mutually exclusive exons. These elements were proposed to engage in competing RNA secondary structures required for mutually exclusive splicing, though this model has not yet been experimentally tested. Here we describe a new system that allowed us to demonstrate that the docking site and selector sequences are indeed required for exon 6 mutually exclusive splicing and that the strength of these RNA structures determines the frequency of exon 6 inclusion. We also show that the function of the docking site has been conserved for ~500 million years of evolution. This work demonstrates that conserved intronic sequences play a functional role in mutually exclusive splicing of the Dscam exon 6 cluster.

First come, first served revisited: factors affecting the same alternative splicing event have different effects on the relative rates of intron removal

Alternative splicing accounts for much of the complexity in higher eukaryotes. Thus, its regulation must allow for flexibility without hampering either its specificity or its fidelity. The mechanisms involved in alternative splicing regulation, especially those acting through coupling with transcription, have not been deeply studied in in vivo models. Much of our knowledge comes from in vitro approaches, where conditions can be precisely controlled at the expense of losing several levels of regulation present in intact cells. Here we studied the relative order of removal of the introns flanking a model alternative cassette exon. We show that there is a preferential removal of the intron downstream from the cassette exon before the upstream intron has been removed. Most importantly, both cis-acting mutations and trans-acting factors that regulate the model alternative splicing event differentially affect the relative order of removal. However, reduction of transcriptional elongation causing higher inclusion of the cassette exon does not change the order of intron removal, suggesting that the assumption, according to the "first come, first served" model, that slow elongation promotes preferential excision of the upstream intron has to be revised. We propose instead that slow elongation favors commitment to exon inclusion during spliceosome assembly. Our results reveal that measuring the order of intron removal may be a straightforward read-out to discriminate among different mechanisms of alternative splice site selection.

Regulation of TMEM16A chloride channel properties by alternative splicing

Expression of TMEM16A protein is associated with the activity of Ca(2+)-activated Cl(-) channels. TMEM16A primary transcript undergoes alternative splicing. thus resulting in the generation of multiple isoforms. We have determined the pattern of splicing and assessed the functional properties of the corresponding TMEM16A variants. We found three alternative exons, 6b, 13, and 15, coding for segments of 22, 4, and 26 amino acids, respectively, which are differently spliced in human organs. By patch clamp experiments on transfected cells, we found that skipping of exon 6b changes the Ca(2+) sensitivity by nearly 4-fold, resulting in Cl(-) currents requiring lower Ca(2+) concentrations to be activated. At the membrane potential of 80 mV, the apparent half-effective concentration decreases from 350 to 90 nm when the segment corresponding to exon 6b is excluded. Skipping of exon 13 instead strongly reduces the characteristic time-dependent activation observed for Ca(2+)-activated Cl(-) channels at positive membrane potentials. This effect was also obtained by deleting only the second pair of amino acids corresponding to exon 13. Alternative splicing appears as an important mechanism to regulate the voltage and Ca(2+) dependence of the TMEM16A-dependent Cl(-) channels in a tissue-specific manner.

DNA damage regulates alternative splicing through inhibition of RNA polymerase II elongation

DNA damage induces apoptosis and many apoptotic genes are regulated via alternative splicing (AS), but little is known about the control mechanisms. Here we show that ultraviolet irradiation (UV) affects cotranscriptional AS in a p53-independent way, through the hyperphosphorylation of RNA polymerase II carboxy-terminal domain (CTD) and a subsequent inhibition of transcriptional elongation, estimated in vivo and in real time. Phosphomimetic CTD mutants not only display lower elongation but also duplicate the UV effect on AS. Consistently, nonphosphorylatable mutants prevent the UV effect. Apoptosis promoted by UV in cells lacking p53 is prevented when the change in AS of the apoptotic gene bcl-x is reverted, confirming the relevance of this mechanism. Splicing-sensitive microarrays revealed a significant overlap of the subsets of genes that have changed AS with UV and those that have reduced expression, suggesting that transcriptional coupling to AS is a key feature of the DNA-damage response.

New insights into form and function of fibronectin splice variants

The extracellular matrix (ECM) is a highly dynamic structure that not only provides a physical framework for cells within connective tissues, but also imparts instructive signals for development, tissue homeostasis and basic cell functions through its composition and ability to exert mechanical forces. The ECM of tissues is composed of, in addition to proteoglycans and hyaluronic acid, a number of proteins, most of which are generated after alternative splicing of their pre-mRNA. However, the precise function of these protein isoforms is still obscure in most cases. Fibronectin (FN), one of the main components of the ECM, is also one of the best-known examples of a family of proteins generated by alternative splicing, having at least 20 different isoforms in humans. Over the last few years, considerable progress on elucidating the functions of the alternatively spliced FN isoforms has been achieved with the essential development of key engineered mouse strains. Here we summarize the phenotypes of the mouse strains having targeted mutations in the FN gene, which may lead to novel insights linking function of alternatively spliced isoforms of fibronectin to human pathologies.

Multiple and specific mRNA processing targets for the major human hnRNP proteins

Alternative splicing is a key mechanism regulating gene expression, and it is often used to produce antagonistic activities particularly in apoptotic genes. Heterogeneous nuclear ribonucleoparticle (hnRNP) proteins form a family of RNA-binding proteins that coat nascent pre-mRNAs. Many but not all major hnRNP proteins have been shown to participate in splicing control. The range and specificity of hnRNP protein action remain poorly documented, even for those affecting splice site selection. We used RNA interference and a reverse transcription-PCR screening platform to examine the implications of 14 of the major hnRNP proteins in the splicing of 56 alternative splicing events in apoptotic genes. Out of this total of 784 alternative splicing reactions tested in three human cell lines, 31 responded similarly to a knockdown in at least two different cell lines. On the other hand, the impact of other hnRNP knockdowns was cell line specific. The broadest effects were obtained with hnRNP K and C, two proteins whose role in alternative splicing had not previously been firmly established. Different hnRNP proteins affected distinct sets of targets with little overlap even between closely related hnRNP proteins. Overall, our study highlights the potential contribution of all of these major hnRNP proteins in alternative splicing control and shows that the targets for individual hnRNP proteins can vary in different cellular contexts.

Functional importance of different patterns of correlation between adjacent cassette exons in human and mouse

Background

Alternative splicing expands transcriptome diversity and plays an important role in regulation of gene expression. Previous studies focus on the regulation of a single cassette exon, but recent experiments indicate that multiple cassette exons within a gene may interact with each other. This interaction can increase the potential to generate various transcripts and adds an extra layer of complexity to gene regulation. Several cases of exon interaction have been discovered. However, the extent to which the cassette exons coordinate with each other remains unknown.

Results

Based on EST data, we employed a metric of correlation coefficients to describe the interaction between two adjacent cassette exons and then categorized these exon pairs into three different groups by their interaction (correlation) patterns. Sequence analysis demonstrates that strongly-correlated groups are more conserved and contain a higher proportion of pairs with reading frame preservation in a combinatorial manner. Multiple genome comparison further indicates that different groups of correlated pairs have different evolutionary courses: (1) The vast majority of positively-correlated pairs are old, (2) most of the weakly-correlated pairs are relatively young, and (3) negatively-correlated pairs are a mixture of old and young events.

Conclusion

We performed a large-scale analysis of interactions between adjacent cassette exons. Compared with weakly-correlated pairs, the strongly-correlated pairs, including both the positively and negatively correlated ones, show more evidence that they are under delicate splicing control and tend to be functionally important. Additionally, the positively-correlated pairs bear strong resemblance to constitutive exons, which suggests that they may evolve from ancient constitutive exons, while negatively and weakly correlated pairs are more likely to contain newly emerging exons.

A 5' splice site enhances the recruitment of basal transcription initiation factors in vivo

Transcription and pre-mRNA splicing are interdependent events. Although mechanisms governing the effects of transcription on splicing are becoming increasingly clear, the means by which splicing affects transcription remain elusive. Using cell lines stably expressing HIV-1 or beta-globin mRNAs, harboring wild-type or various 5' splice site mutations, we demonstrate a strong positive correlation between splicing efficiency and transcription activity. Interestingly, a 5' splice site can stimulate transcription even in the absence of splicing. Chromatin immunoprecipitation experiments show enhanced promoter docking of transcription initiation factors TFIID, TFIIB, and TFIIH on a gene containing a functional 5' splice site. In addition to their promoter association, the TFIID and TFIIH components, TBP and p89, are specifically recruited to the 5' splice site region. Our data suggest a model in which a promoter-proximal 5' splice site via its U1 snRNA interaction can feed back to stimulate transcription initiation by enhancing pre-initiation complex assembly.

Functional coordination of alternative splicing in the mammalian central nervous system

A microarray analysis provides new evidence suggesting that specific cellular processes in the mammalian CNS are coordinated at the level of alternative splicing, and that a complex splicing code underlies CNS-specific alternative splicing regulation.

Background

Alternative splicing (AS) functions to expand proteomic complexity and plays numerous important roles in gene regulation. However, the extent to which AS coordinates functions in a cell and tissue type specific manner is not known. Moreover, the sequence code that underlies cell and tissue type specific regulation of AS is poorly understood.

Results

Using quantitative AS microarray profiling, we have identified a large number of widely expressed mouse genes that contain single or coordinated pairs of alternative exons that are spliced in a tissue regulated fashion. The majority of these AS events display differential regulation in central nervous system (CNS) tissues. Approximately half of the corresponding genes have neural specific functions and operate in common processes and interconnected pathways. Differential regulation of AS in the CNS tissues correlates strongly with a set of mostly new motifs that are predominantly located in the intron and constitutive exon sequences neighboring CNS-regulated alternative exons. Different subsets of these motifs are correlated with either increased inclusion or increased exclusion of alternative exons in CNS tissues, relative to the other profiled tissues.

Conclusion

Our findings provide new evidence that specific cellular processes in the mammalian CNS are coordinated at the level of AS, and that a complex splicing code underlies CNS specific AS regulation. This code appears to comprise many new motifs, some of which are located in the constitutive exons neighboring regulated alternative exons. These data provide a basis for understanding the molecular mechanisms by which the tissue specific functions of widely expressed genes are coordinated at the level of AS.

Identification of an intronic splicing enhancer essential for the inclusion of FGFR2 exon IIIc

The ligand specificity of fibroblast growth factor receptor 2 (FGFR2) is determined by the alternative splicing of exons 8 (IIIb) or 9 (IIIc). Exon IIIb is included in epithelial cells, whereas exon IIIc is included in mesenchymal cells. Although a number of cis elements and trans factors have been identified that play a role in exon IIIb inclusion in epithelium, little is known about the activation of exon IIIc in mesenchyme. We report here the identification of a splicing enhancer required for IIIc inclusion. This 24-nucleotide (nt) downstream intronic splicing enhancer (DISE) is located within intron 9 immediately downstream of exon IIIc. DISE was able to activate the inclusion of heterologous exons rat FGFR2 IIIb and human beta-globin exon 2 in cell lines from different tissues and species and also in HeLa cell nuclear extracts in vitro. DISE was capable of replacing the intronic activator sequence 1 (IAS1), a known IIIb splicing enhancer and vice versa. This fact, together with the requirement for DISE to be close to the 5'-splice site and the ability of DISE to promote binding of U1 snRNP, suggested that IAS1 and DISE belong to the same class of cis-acting elements.

Culture of human embryonic stem cells and the extracellular matrix microenvironment

Among the major obstacles impeding successful derivation and continuous culture of human embryonic stem cells (hESC) for therapeutic purposes, are the presence of feeder cells and feeder-conditioned media of animal origin. The risk of contamination with xenopathogens makes hESC cultured in this way unsafe for future use in regenerative medicine. A holy grail for investigators in the field will be to establish and maintain new hESC lines in completely feeder-free and serum-free defined conditions. Recently, propagation of hESC has become possible, using mammalian- or human-derived extracellular matrix (ECM) and conditioned medium from feeder cells. In addition, providing a three-dimensional ECM environment can even support the derivation of new hESC. In this review, we examine recent advances in the use and development of substrates suitable for the derivation and maintenance of hESC, and our current understanding of the effects of a three-dimensional ECM milieu on cellular behavior.

Unfolding the mystery of alternative splicing through a unique method of in vivo selection

Alternative splicing of pre-messenger RNA (pre-mRNA) is a fundamental mechanism of gene regulation in higher eukaryotes. In addition to creating protein diversity, alternative splicing provides the safest mode of gene evolution. Of late, more and more forms of alternatively spliced transcripts (mRNAs) are being discovered for key genes. Some of the alternatively spliced transcripts are also associated with major human diseases. This has created a sense of urgency to find the methods by which regulation of alternative splicing of specific exons could be best understood. Here I review a powerful in vivo selection method that uses a combinatorial library of partially random sequences. Several advantages of this method include in vivo analysis of large sequences, identification of unique sequence motifs, determination of relative strength of splice sites and identification of long-distance interactions including role of RNA structures. This unique method could be applied to identify tissue-specific cis-elements. Similarly, the method is suitable to find cis-elements that become active in response to specific treatments of cells. Considering this unbiased method uses in vivo conditions, it has potential to identify critical regulatory elements as therapeutic targets for a growing number of splicing-associated diseases.

Signals, pathways and splicing regulation

Alternative splicing of messenger RNA precursors is an extraordinary source of protein diversity and the regulation of this process is crucial for diverse cellular functions in both physiological and pathological situations. For many years, several signaling pathways have been implicated in alternative splicing regulation. Recent work has begun to unravel the molecular mechanisms by which extracellular stimuli activate signaling cascades that modulate the activity of the splicing machinery and therefore the splicing pattern of many different target messenger RNA precursors. These experiments are revealing unexpected aspects of the mechanism that control splicing and the consequences of the regulated splicing events. We summarize here the current knowledge about signal-induced alternative splicing regulation of Slo, NR1, CD44, CD45 and fibronectin genes, and also discuss the importance of some of these events in determination of cellular fate. Furthermore, we highlight the relevance of signal-induced changes in phosphorylation state and subcellular distribution of splicing factors as a way of regulating the splicing process. Lastly, we explore new and unexpected findings about regulated splicing in anucleated cells.

Multivariate analysis and visualization of splicing correlations in single-gene transcriptomes

BACKGROUND

RNA metabolism, through 'combinatorial splicing', can generate enormous structural diversity in the proteome. Alternative domains may interact, however, with unpredictable phenotypic consequences, necessitating integrated RNA-level regulation of molecular composition. Splicing correlations within transcripts of single genes provide valuable clues to functional relationships among molecular domains as well as genomic targets for higher-order splicing regulation.

RESULTS

We present tools to visualize complex splicing patterns in full-length cDNA libraries. Developmental changes in pair-wise correlations are presented vectorially in 'clock plots' and linkage grids. Higher-order correlations are assessed statistically through Monte Carlo analysis of a log-linear model with an empirical-Bayes estimate of the true probabilities of observed and unobserved splice forms. Log-linear coefficients are visualized in a 'spliceprint,' a signature of splice correlations in the transcriptome. We present two novel metrics: the linkage change index, which measures the directional change in pair-wise correlation with tissue differentiation, and the accuracy index, a very simple goodness-of-fit metric that is more sensitive than the integrated squared error when applied to sparsely populated tables, and unlike chi-square, does not diverge at low variance. Considerable attention is given to sparse contingency tables, which are inherent to single-gene libraries.

CONCLUSION

Patterns of splicing correlations are revealed, which span a broad range of interaction order and change in development. The methods have a broad scope of applicability, beyond the single gene--including, for example, multiple gene interactions in the complete transcriptome.

RNA polymerase II C-terminal domain mediates regulation of alternative splicing by SRp20

Previous studies have linked the C-terminal domain (CTD) of RNA polymerase II (pol II) with cotranscriptional precursor messenger RNA processing, but little is known about the CTD's function in regulating alternative splicing. We have examined this function using alpha-amanitin-resistant pol II CTD mutants and fibronectin reporter minigenes. We found that the CTD is required for the inhibitory action of the serine/arginine-rich (SR) protein SRp20 on the inclusion of a fibronectin cassette exon in the mature mRNA. CTD phosphorylation controls transcription elongation, which is a major contributor to alternative splicing regulation. However, the effect of SRp20 is still observed when transcription elongation is reduced. These results suggest that the CTD promotes exon skipping by recruiting SRp20 and that this contributes independently of elongation to the transcriptional control of alternative splicing.

Model of stimulation-responsive splicing and strategies in identification of immunogenic isoforms of tumor antigens and autoantigens

We recently proposed a novel model of stimulation-responsive splicing for the selection of autoantigens and self-tumor antigens. Our model theorizes that the significantly higher rates of alternative splicing of autoantigen and self-tumor antigen transcripts that occur in response to stimuli could induce extra-thymic expression of untolerized antigen epitopes for elicitation of autoimmune and anti-tumor responses. To facilitate the identification of immunogenic isoforms of antigens, we have developed strategies using improved SEREX in conjunction with database-mining and immunogenic isoform mapping. Identification of immunogenic isoforms of autoantigens and self-tumor antigens is very important for the development of novel therapeutics and diagnostic tools for autoimmune diseases and tumors, such as: (1) autoantigen isoform microarrays for disease diagnosis and prognosis; (2) autoantigen isoform-specific tolerizing vaccines and splicing-redirection therapies, as well as (3) immunogenic antigen isoform-specific immunotherapy for tumors.

Defective splicing, disease and therapy: searching for master checkpoints in exon definition

The number of aberrant splicing processes causing human disease is growing exponentially and many recent studies have uncovered some aspects of the unexpectedly complex network of interactions involved in these dysfunctions. As a consequence, our knowledge of the various cis- and trans-acting factors playing a role on both normal and aberrant splicing pathways has been enhanced greatly. However, the resulting information explosion has also uncovered the fact that many splicing systems are not easy to model. In fact we are still unable, with certainty, to predict the outcome of a given genomic variation. Nonetheless, in the midst of all this complexity some hard won lessons have been learned and in this survey we will focus on the importance of the wide sequence context when trying to understand why apparently similar mutations can give rise to different effects. The examples discussed in this summary will highlight the fine 'balance of power' that is often present between all the various regulatory elements that define exon boundaries. In the final part, we shall then discuss possible therapeutic targets and strategies to rescue genetic defects of complex splicing systems.