Alternative splicing caused by RNA secondary structure

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

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.

Intronic RNA: Ad'junk' mediator of post-transcriptional gene regulation

RNA splicing, the process through which intervening segments of noncoding RNA (introns) are excised from pre-mRNAs to allow for the formation of a mature mRNA product, has long been appreciated for its capacity to add complexity to eukaryotic proteomes. However, evidence suggests that the utility of this process extends beyond protein output and provides cells with a dynamic tool for gene regulation. In this review, we aim to highlight the role that intronic RNA plays in mediating specific splicing outcomes in pre-mRNA processing, as well as explore an emerging class of stable intronic sequences that have been observed to act in gene expression control. Building from underlying flexibility in both sequence and structure, intronic RNA provides mechanisms for post-transcriptional gene regulation that are amenable to the tissue and condition specific needs of eukaryotic cells. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

β-Mannosidosis caused by a novel homozygous intragenic inverted duplication in MANBA

β-Mannosidosis is a lysosomal storage disorder characterized by accumulation of disaccharides due to deficiency of the lysosomal enzyme β-mannosidase. The disease is caused by mutations in MANBA and is extremely rare in humans. Although the clinical presentation is heterogeneous, common symptoms include various degrees of developmental delay, behavioral disturbances, hearing loss, and frequent infections. We report a 15-yr-old girl presenting with mild intellectual disability, sensorineural hearing loss, severe behavioral disturbances, dysmorphic traits, and evolving angiokeratomas. Copy-number variation analysis of next-generation sequencing (NGS) data indicated increased coverage in exons 8-11 of MANBA Low β-mannosidase activity (1 µkatal/kg protein, refv 25-40) established the diagnosis of β-mannosidosis. Whole-genome sequencing (WGS) and cDNA analysis revealed a novel homozygous intragenic inverted duplication in MANBA, where a 13.1-kb region between introns 7 and 11 was duplicated and inserted in an inverted orientation, creating a 67-base nonduplicated gap at the insertion point. Both junctions showed microhomology regions. The inverted duplication resulted in exon skipping of exons 8-9 or 8-10. Our report highlights the importance of copy-number variation analysis of data from NGS and in particular the power of WGS in the identification and characterization of copy-number variants.

SpliceV: analysis and publication quality printing of linear and circular RNA splicing, expression and regulation

BACKGROUND

In eukaryotes, most genes code for multiple transcript isoforms that are generated through the complex and tightly regulated process of RNA splicing. Despite arising from identical precursor transcripts, alternatively spliced RNAs can have dramatically different functions. Transcriptome complexity is elevated further by the production of circular RNAs (circRNAs), another class of mature RNA that results from the splicing of a downstream splice donor to an upstream splice acceptor. While there has been a rapid expansion of circRNA catalogs in the last few years through the utilization of next generation sequencing approaches, our understanding of the mechanisms and regulation of circular RNA biogenesis, the impact that circRNA generation has on parental transcript processing, and the functions carried out by circular RNAs remains limited.

RESULTS

Here, we present a visualization and analysis tool, SpliceV, that rapidly plots all relevant forward- and back-splice data, with exon and single nucleotide level coverage information from RNA-seq experiments in a publication quality format. SpliceV also integrates analysis features that assist investigations into splicing regulation and transcript functions through the display of predicted RNA binding protein sites and the configuration of repetitive elements along the primary transcript.

CONCLUSIONS

SpliceV is an easy-to-use splicing visualization tool, compatible with both Python 2.7 and 3+, and distributed under the GNU Public License. The source code is freely available for download at https://github.com/flemingtonlab/SpliceV and can be installed from PyPI using pip.

RNAs in the spliceosome: Insight from cryoEM structures

Pre-mRNA splicing is catalyzed by the spliceosome, a multimegadalton RNA-protein complex. The spliceosome undergoes dramatic compositional and conformational changes through the splicing cycle, forming at least 10 distinct complexes. Recent high-resolution cryoEM structures of various spliceosomal complexes revealed unprecedented details of this large molecular machine. This review highlights insight into the structure and function of the spliceosomal RNA components obtained from these new structures, with a focus on the yeast spliceosome. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Gene Therapy via Trans-Splicing for LMNA-Related Congenital Muscular Dystrophy

We assessed the potential of Lmna-mRNA repair by spliceosome-mediated RNA trans-splicing as a therapeutic approach for LMNA-related congenital muscular dystrophy. This gene therapy strategy leads to reduction of mutated transcript expression for the benefit of corresponding wild-type (WT) transcripts. We developed 5′-RNA pre-trans-splicing molecules containing the first five exons of Lmna and targeting intron 5 of Lmna pre-mRNA. Among nine pre-trans-splicing molecules, differing in the targeted sequence in intron 5 and tested in C2C12 myoblasts, three induced trans-splicing events on endogenous Lmna mRNA and confirmed at protein level. Further analyses performed in primary myotubes derived from an LMNA-related congenital muscular dystrophy (L-CMD) mouse model led to a partial rescue of the mutant phenotype. Finally, we tested this approach in vivo using adeno-associated virus (AAV) delivery in newborn mice and showed that trans-splicing events occurred in WT mice 50 days after AAV delivery, although at a low rate. Altogether, while these results provide the first evidence for reprogramming LMNA mRNA in vitro, strategies to improve the rate of trans-splicing events still need to be developed for efficient application of this therapeutic approach in vivo.

The effects of structure on pre-mRNA processing and stability

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

RNA structure in splicing: An evolutionary perspective

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic β sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

A Phylogenetic Survey on the Structure of the HIV-1 Leader RNA Domain That Encodes the Splice Donor Signal

RNA splicing is a critical step in the human immunodeficiency virus type 1 (HIV-1) replication cycle because it controls the expression of the complex viral proteome. The major 5′ splice site (5′ss) that is positioned in the untranslated leader of the HIV-1 RNA transcript is of particular interest because it is used for the production of the more than 40 differentially spliced subgenomic mRNAs. HIV-1 splicing needs to be balanced tightly to ensure the proper levels of all viral proteins, including the Gag-Pol proteins that are translated from the unspliced RNA. We previously presented evidence that the major 5′ss is regulated by a repressive local RNA structure, the splice donor (SD) hairpin, that masks the 11 nucleotides (nts) of the 5′ss signal for recognition by U1 small nuclear RNA (snRNA) of the spliceosome machinery. A strikingly different multiple-hairpin RNA conformation was recently proposed for this part of the HIV-1 leader RNA. We therefore inspected the sequence of natural HIV-1 isolates in search for support, in the form of base pair (bp) co-variations, for the different RNA conformations.

Viral interactions with components of the splicing machinery

Eukaryotic genes are often interrupted by stretches of sequence with no protein coding potential or obvious function. After transcription, these interrupting sequences must be removed to give rise to the mature messenger RNA. This fundamental process is called RNA splicing and is achieved by complicated machinery made of protein and RNA that assembles around the RNA to be edited. Viruses also use RNA splicing to maximize their coding potential and economize on genetic space, and use clever strategies to manipulate the splicing machinery to their advantage. This article gives an overview of the splicing process and provides examples of viral strategies that make use of various components of the splicing system to promote their replicative cycle. Representative virus families have been selected to illustrate the interaction with various regulatory proteins and ribonucleoproteins. The unifying theme is fine regulation through protein-protein and protein-RNA interactions with the spliceosome components and associated factors to promote or prevent spliceosome assembly on given splice sites, in addition to a strong influence from cis-regulatory sequences on viral transcripts. Because there is an intimate coupling of splicing with the processes that direct mRNA biogenesis, a description of how these viruses couple the regulation of splicing with the retention or stability of mRNAs is also included. It seems that a unique balance of suppression and activation of splicing and nuclear export works optimally for each family of viruses.

The Conservation and Function of RNA Secondary Structure in Plants

RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.

Diverse regulation of 3' splice site usage

The regulation of splice site (SS) usage is important for alternative pre-mRNA splicing and thus proper expression of protein isoforms in cells; its disruption causes diseases. In recent years, an increasing number of novel regulatory elements have been found within or nearby the 3'SS in mammalian genes. The diverse elements recruit a repertoire of trans-acting factors or form secondary structures to regulate 3'SS usage, mostly at the early steps of spliceosome assembly. Their mechanisms of action mainly include: (1) competition between the factors for RNA elements, (2) steric hindrance between the factors, (3) direct interaction between the factors, (4) competition between two splice sites, or (5) local RNA secondary structures or longer range loops, according to the mode of protein/RNA interactions. Beyond the 3'SS, chromatin remodeling/transcription, posttranslational modifications of trans-acting factors and upstream signaling provide further layers of regulation. Evolutionarily, some of the 3'SS elements seem to have emerged in mammalian ancestors. Moreover, other possibilities of regulation such as that by non-coding RNA remain to be explored. It is thus likely that there are more diverse elements/factors and mechanisms that influence the choice of an intron end. The diverse regulation likely contributes to a more complex but refined transcriptome and proteome in mammals.

HIV-1 splicing is controlled by local RNA structure and binding of splicing regulatory proteins at the major 5' splice site

The 5' leader region of the human immunodeficiency virus 1 (HIV-1) RNA genome contains the major 5' splice site (ss) that is used in the production of the many spliced viral RNAs. This splice-donor (SD) region can fold into a stable stem-loop structure and the thermodynamic stability of this RNA hairpin influences splicing efficiency. In addition, splicing may be modulated by binding of splicing regulatory (SR) proteins, in particular SF2/ASF (SRSF1), SC35 (SRSF2), SRp40 (SRSF5) and SRp55 (SRSF6), to sequence elements in the SD region. The role of RNA structure and SR protein binding in splicing control was previously studied by functional analysis of mutant SD sequences. The interpretation of these studies was complicated by the fact that most mutations simultaneously affect both structure and sequence elements. We therefore tried to disentangle the contribution of these two variables by designing more precise SD region mutants with a single effect on either the sequence or the structure. The current analysis indicates that HIV-1 splicing at the major 5'ss is modulated by both the stability of the local RNA structure and the binding of splicing regulatory proteins.

The separation between the 5'-3' ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.

A novel mutation in CLCN1 associated with feline myotonia congenita

Myotonia congenita (MC) is a skeletal muscle channelopathy characterized by inability of the muscle to relax following voluntary contraction. Worldwide population prevalence in humans is 1∶100,000. Studies in mice, dogs, humans and goats confirmed myotonia associated with functional defects in chloride channels and mutations in a skeletal muscle chloride channel (CLCN1). CLCN1 encodes for the most abundant chloride channel in the skeletal muscle cell membrane. Five random bred cats from Winnipeg, Canada with MC were examined. All cats had a protruding tongue, limited range of jaw motion and drooling with prominent neck and proximal limb musculature. All cats had blepharospasm upon palpebral reflex testing and a short-strided gait. Electromyograms demonstrated myotonic discharges at a mean frequency of 300 Hz resembling the sound of a ‘swarm of bees’. Muscle histopathology showed hypertrophy of all fiber types. Direct sequencing of CLCN1 revealed a mutation disrupting a donor splice site downstream of exon 16 in only the affected cats. In vitro translation of the mutated protein predicted a premature truncation and partial lack of the highly conserved CBS1 (cystathionine β-synthase) domain critical for ion transport activity and one dimerization domain pivotal in channel formation. Genetic screening of the Winnipeg random bred population of the cats' origin identified carriers of the mutation. A genetic test for population screening is now available and carrier cats from the feral population can be identified.

HIV-1 splicing at the major splice donor site is restricted by RNA structure

The 5' leader region of the HIV-1 RNA contains the major 5' splice site (ss) that is used in the production of all spliced viral RNAs. This splice-donor (SD) region can fold a stem-loop structure. We demonstrate that whereas stabilization of this SD hairpin reduces splicing efficiency, destabilization increases splicing. Both stabilization and destabilization reduce viral fitness. These results demonstrate that the stability of the SD hairpin can modulate the level of splicing, most likely by controlling the accessibility of the 5'ss for the splicing machinery. The natural stability of the SD hairpin restricts splicing and this stability seems to be fine-tuned to reach the optimal balance between unspliced and spliced RNAs for efficient virus replication. The 5'ss region of different HIV-1 isolates and the related SIVmac239 can fold a similar structure. This evolutionary conservation supports the importance of this structure in viral replication.

Comparative analyses between retained introns and constitutively spliced introns in Arabidopsis thaliana using random forest and support vector machine

One of the important modes of pre-mRNA post-transcriptional modification is alternative splicing. Alternative splicing allows creation of many distinct mature mRNA transcripts from a single gene by utilizing different splice sites. In plants like Arabidopsis thaliana, the most common type of alternative splicing is intron retention. Many studies in the past focus on positional distribution of retained introns (RIs) among different genic regions and their expression regulations, while little systematic classification of RIs from constitutively spliced introns (CSIs) has been conducted using machine learning approaches. We used random forest and support vector machine (SVM) with radial basis kernel function (RBF) to differentiate these two types of introns in Arabidopsis. By comparing coordinates of introns of all annotated mRNAs from TAIR10, we obtained our high-quality experimental data. To distinguish RIs from CSIs, We investigated the unique characteristics of RIs in comparison with CSIs and finally extracted 37 quantitative features: local and global nucleotide sequence features of introns, frequent motifs, the signal strength of splice sites, and the similarity between sequences of introns and their flanking regions. We demonstrated that our proposed feature extraction approach was more accurate in effectively classifying RIs from CSIs in comparison with other four approaches. The optimal penalty parameter C and the RBF kernel parameter in SVM were set based on particle swarm optimization algorithm (PSOSVM). Our classification performance showed F-Measure of 80.8% (random forest) and 77.4% (PSOSVM). Not only the basic sequence features and positional distribution characteristics of RIs were obtained, but also putative regulatory motifs in intron splicing were predicted based on our feature extraction approach. Clearly, our study will facilitate a better understanding of underlying mechanisms involved in intron retention.

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.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

Eukaryotic TPP riboswitch regulation of alternative splicing involving long-distance base pairing

Thiamin pyrophosphate (TPP) riboswitches are found in organisms from all three domains of life. Examples in bacteria commonly repress gene expression by terminating transcription or by blocking ribosome binding, whereas most eukaryotic TPP riboswitches are predicted to regulate gene expression by modulating RNA splicing. Given the widespread distribution of eukaryotic TPP riboswitches and the diversity of their locations in precursor messenger RNAs (pre-mRNAs), we sought to examine the mechanism of alternative splicing regulation by a fungal TPP riboswitch from Neurospora crassa, which is mostly located in a large intron separating protein-coding exons. Our data reveal that this riboswitch uses a long-distance (∼530-nt separation) base-pairing interaction to regulate alternative splicing. Specifically, a portion of the TPP-binding aptamer can form a base-paired structure with a conserved sequence element (α) located near a 5′ splice site, which greatly increases use of this 5′ splice site and promotes gene expression. Comparative sequence analyses indicate that many fungal species carry a TPP riboswitch with similar intron architecture, and therefore the homologous genes in these fungi are likely to use the same mechanism. Our findings expand the scope of genetic control mechanisms relying on long-range RNA interactions to include riboswitches.

Genetic variations and alternative splicing: the Glioma associated oncogene 1, GLI1

Alternative splicing is a post-transcriptional regulatory process that is attaining stronger recognition as a modulator of gene expression. Alternative splicing occurs when the primary RNA transcript is differentially processed into more than one mature RNAs. This is the result of a variable definition/inclusion of the exons, the sequences that are excised from the primary RNA to form the mature RNAs. Consequently, RNA expression can generate a collection of differentially spliced RNAs, which may distinctly influence subsequent biological events, such as protein synthesis or other biomolecular interactions. Still the mechanisms that control exon definition and exon inclusion are not fully clarified. This mini-review highlights advances in this field as well as the impact of single nucleotide polymorphisms in affecting splicing decisions. The Glioma-associated oncogene 1, GLI1, is taken as an example in addressing the role of nucleotide substitutions for splicing regulation.

Evidence for widespread association of mammalian splicing and conserved long-range RNA structures

Pre-mRNA structure impacts many cellular processes, including splicing in genes associated with disease. The contemporary paradigm of RNA structure prediction is biased toward secondary structures that occur within short ranges of pre-mRNA, although long-range base-pairings are known to be at least as important. Recently, we developed an efficient method for detecting conserved RNA structures on the genome-wide scale, one that does not require multiple sequence alignments and works equally well for the detection of local and long-range base-pairings. Using an enhanced method that detects base-pairings at all possible combinations of splice sites within each gene, we now report RNA structures that could be involved in the regulation of splicing in mammals. Statistically, we demonstrate strong association between the occurrence of conserved RNA structures and alternative splicing, where local RNA structures are generally more frequent at alternative donor splice sites, while long-range structures are more associated with weak alternative acceptor splice sites. As an example, we validated the RNA structure in the human SF1 gene using minigenes in the HEK293 cell line. Point mutations that disrupted the base-pairing of two complementary boxes between exons 9 and 10 of this gene altered the splicing pattern, while the compensatory mutations that reestablished the base-pairing reverted splicing to that of the wild-type. There is statistical evidence for a Dscam-like class of mammalian genes, in which mutually exclusive RNA structures control mutually exclusive alternative splicing. In sum, we propose that long-range base-pairings carry an important, yet unconsidered part of the splicing code, and that, even by modest estimates, there must be thousands of such potentially regulatory structures conserved throughout the evolutionary history of mammals.

Post-transcriptional exon shuffling events in humans can be evolutionarily conserved and abundant

In silico analyses have established that transcripts from some genes can be processed into RNAs with rearranged exon order relative to genomic structure (post-transcriptional exon shuffling, or PTES). Although known to contribute to transcriptome diversity in some species, to date the structure, distribution, abundance, and functional significance of human PTES transcripts remains largely unknown. Here, using high-throughput transcriptome sequencing, we identify 205 putative human PTES products from 176 genes. We validate 72 out of 112 products analyzed using RT-PCR, and identify additional PTES products structurally related to 61% of validated targets. Sequencing of these additional products reveals GT-AG dinucleotides at >95% of the splice junctions, confirming that they are processed by the spliceosome. We show that most PTES transcripts are expressed in a wide variety of human tissues, that they can be polyadenylated, and that some are conserved in mouse. We also show that they can extend into 5' and 3' UTRs, consistent with formation via trans-splicing of independent pre-mRNA molecules. Finally, we use real-time PCR to compare the abundance of PTES exon junctions relative to canonical exon junctions within the transcripts from seven genes. PTES exon junctions are present at <0.01% to >90% of the levels of canonical junctions, with transcripts from MAN1A2, PHC3, TLE4, and CDK13 exhibiting the highest levels. This is the first systematic experimental analysis of PTES in human, and it suggests both that the phenomenon is much more widespread than previously thought and that some PTES transcripts could be functional.

A conditional random fields method for RNA sequence-structure relationship modeling and conformation sampling

Accurate tertiary structures are very important for the functional study of non-coding RNA molecules. However, predicting RNA tertiary structures is extremely challenging, because of a large conformation space to be explored and lack of an accurate scoring function differentiating the native structure from decoys. The fragment-based conformation sampling method (e.g. FARNA) bears shortcomings that the limited size of a fragment library makes it infeasible to represent all possible conformations well. A recent dynamic Bayesian network method, BARNACLE, overcomes the issue of fragment assembly. In addition, neither of these methods makes use of sequence information in sampling conformations. Here, we present a new probabilistic graphical model, conditional random fields (CRFs), to model RNA sequence–structure relationship, which enables us to accurately estimate the probability of an RNA conformation from sequence. Coupled with a novel tree-guided sampling scheme, our CRF model is then applied to RNA conformation sampling. Experimental results show that our CRF method can model RNA sequence–structure relationship well and sequence information is important for conformation sampling. Our method, named as TreeFolder, generates a much higher percentage of native-like decoys than FARNA and BARNACLE, although we use the same simple energy function as BARNACLE.

Contact:zywang@ttic.edu; j3xu@ttic.edu

Supplementary Information:Supplementary data are available at Bioinformatics online.

Secondary structure is required for 3' splice site recognition in yeast

Higher order RNA structures can mask splicing signals, loop out exons, or constitute riboswitches all of which contributes to the complexity of splicing regulation. We identified a G to A substitution between branch point (BP) and 3′ splice site (3′ss) of Saccharomyces cerevisiae COF1 intron, which dramatically impaired its splicing. RNA structure prediction and in-line probing showed that this mutation disrupted a stem in the BP-3′ss region. Analyses of various COF1 intron modifications revealed that the secondary structure brought about the reduction of BP to 3′ss distance and masked potential 3′ss. We demonstrated the same structural requisite for the splicing of UBC13 intron. Moreover, RNAfold predicted stable structures for almost all distant BP introns in S. cerevisiae and for selected examples in several other Saccharomycotina species. The employment of intramolecular structure to localize 3′ss for the second splicing step suggests the existence of pre-mRNA structure-based mechanism of 3′ss recognition.

Stoichiometry of a regulatory splicing complex revealed by single-molecule analyses

Splicing is regulated by complex interactions of numerous RNA-binding proteins. The molecular mechanisms involved remain elusive, in large part because of ignorance regarding the numbers of proteins in regulatory complexes. Polypyrimidine tract-binding protein (PTB), which regulates tissue-specific splicing, represses exon 3 of alpha-tropomyosin through distant pyrimidine-rich tracts in the flanking introns. Current models for repression involve either PTB-mediated looping or the propagation of complexes between tracts. To test these models, we used single-molecule approaches to count the number of bound PTB molecules both by counting the number of bleaching steps of GFP molecules linked to PTB within complexes and by analysing their total emissions. Both approaches showed that five or six PTB molecules assemble. Given the domain structures, this suggests that the molecules occupy primarily multiple overlapping potential sites in the polypyrimidine tracts, excluding propagation models. As an alternative to direct looping, we propose that repression involves a multistep process in which PTB binding forms small local loops, creating a platform for recruitment of other proteins that bring these loops into close proximity.

Limited complementarity between U1 snRNA and a retroviral 5' splice site permits its attenuation via RNA secondary structure

Multiple types of regulation are used by cells and viruses to control alternative splicing. In murine leukemia virus, accessibility of the 5′ splice site (ss) is regulated by an upstream region, which can fold into a complex RNA stem–loop structure. The underlying sequence of the structure itself is negligible, since most of it could be functionally replaced by a simple heterologous RNA stem–loop preserving the wild-type splicing pattern. Increasing the RNA duplex formation between U1 snRNA and the 5′ss by a compensatory mutation in position +6 led to enhanced splicing. Interestingly, this mutation affects splicing only in the context of the secondary structure, arguing for a dynamic interplay between structure and primary 5′ss sequence. The reduced 5′ss accessibility could also be counteracted by recruiting a splicing enhancer domain via a modified MS2 phage coat protein to a single binding site at the tip of the simple RNA stem–loop. The mechanism of 5′ss attenuation was revealed using hyperstable U1 snRNA mutants, showing that restricted U1 snRNP access is the cause of retroviral alternative splicing.

Insulin Receptor and its Relationship with Different Forms of Insulin Resistance

Insulin plays an important role in maintaining the whole organism’s homeostasis. The presence of insulin receptors in all vertebrates and invertebrates cells reflects the diversity of regulatory processes in which this hormone is involved. Furthermore, many different factors may influence the level of insulin receptor expression. These factors include e.g. the sole insulin or stage of development. Mutations in the receptor may lead to the development of insulin resistance. These mutations differ in the level of severity and are frequently associated with diabetes mellitus, hypertension, cardiovascular disorders, heart failure, metabolic syndrome and infertility in women. More than 50 mutations in insulin receptor gene have already been characterized. These mutations are associated with rare forms of insulin resistance like leprechaunism, insulin resistance type A or Rabson-Mendenhall syndrome. Molecular analysis of insulin receptor gene may lead to a better understanding of molecular mechanisms underlying various types of insulin resistance and help to develop more efficient treatment.

Alu exonization events reveal features required for precise recognition of exons by the splicing machinery

Despite decades of research, the question of how the mRNA splicing machinery precisely identifies short exonic islands within the vast intronic oceans remains to a large extent obscure. In this study, we analyzed Alu exonization events, aiming to understand the requirements for correct selection of exons. Comparison of exonizing Alus to their non-exonizing counterparts is informative because Alus in these two groups have retained high sequence similarity but are perceived differently by the splicing machinery. We identified and characterized numerous features used by the splicing machinery to discriminate between Alu exons and their non-exonizing counterparts. Of these, the most novel is secondary structure: Alu exons in general and their 5′ splice sites (5′ss) in particular are characterized by decreased stability of local secondary structures with respect to their non-exonizing counterparts. We detected numerous further differences between Alu exons and their non-exonizing counterparts, among others in terms of exon–intron architecture and strength of splicing signals, enhancers, and silencers. Support vector machine analysis revealed that these features allow a high level of discrimination (AUC = 0.91) between exonizing and non-exonizing Alus. Moreover, the computationally derived probabilities of exonization significantly correlated with the biological inclusion level of the Alu exons, and the model could also be extended to general datasets of constitutive and alternative exons. This indicates that the features detected and explored in this study provide the basis not only for precise exon selection but also for the fine-tuned regulation thereof, manifested in cases of alternative splicing.

A typical human gene consists of 9 exons around 150 nucleotides in length, separated by introns that are ∼3,000 nucleotides long. The challenge of the splicing machinery is to precisely identify and ligate the exons, while removing the introns. We aimed to understand how the splicing machinery meets this momentous challenge, based on Alu exonization events. Alus are transposable elements, of which approximately one million copies exist in the human genome, a large portion of which within introns. Throughout evolution, some intronic Alus accumulated mutations and became recognized by the splicing machinery as exons, a process termed exonization. Such Alus remain highly similar to their non-exonizing counterparts but are perceived as different by the splicing machinery. By comparing exonizing Alus to their non-exonizing counterparts, we were able to identify numerous features in which they differ and which presumably lead to the recognition only of the former by the splicing machinery. Our findings reveal insights regarding the role of local RNA secondary structures, exon–intron architecture constraints, and splicing regulatory signals. We integrated these features in a computational model, which was able to successfully mimic the function of the splicing machinery and discriminate between true Alu exons and their intronic counterparts, highlighting the functional importance of these features.

Intronic Alus influence alternative splicing

Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called Alu. A large fraction of Alus are found in intronic sequences, implying extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these intronic Alus on splicing of the flanking exons is largely unknown. Here, we show that more Alus flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that Alu insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two Alu elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic Alus in influencing the splicing of flanking exons, further emphasizing the role of Alus in shaping of the human transcriptome.

Ligand-induced sequestering of branchpoint sequence allows conditional control of splicing

Background

Despite tremendous progress in understanding the mechanisms of constitutive and alternative splicing, an important and widespread step along the gene expression pathway, our ability to deliberately regulate gene expression at this step remains rudimentary. The present study was performed to investigate whether a theophylline-dependent "splice switch" that sequesters the branchpoint sequence (BPS) within RNA-theophylline complex can regulate alternative splicing.

Results

We constructed a series of pre-mRNAs in which the BPS was inserted within theophylline aptamer. We show that theophylline-induced sequestering of BPS inhibits pre-mRNA splicing both in vitro and in vivo in a dose-dependent manner. Several lines of evidence suggest that theophylline-dependent inhibition of splicing is highly specific, and thermodynamic stability of RNA-theophylline complex as well as the location of BPS within this complex affects the efficiency of splicing inhibition. Finally, we have constructed an alternative splicing model pre-mRNA substrate in which theophylline caused exon skipping both in vitro and in vivo, suggesting that a small molecule-RNA interaction can modulate alternative splicing.

Conclusion

These findings provide the ability to control splicing pattern at will and should have important implications for basic, biotechnological, and biomedical research.

Association of polymorphisms in odorant-binding protein genes with variation in olfactory response to benzaldehyde in Drosophila

Adaptive evolution of animals depends on behaviors that are essential for their survival and reproduction. The olfactory system of Drosophila melanogaster has emerged as one of the best characterized olfactory systems, which in addition to a family of odorant receptors, contains an approximately equal number of odorant-binding proteins (OBPs), encoded by a multigene family of 51 genes. Despite their abundant expression, little is known about their role in chemosensation, largely due to the lack of available mutations in these genes. We capitalized on naturally occurring mutations (polymorphisms) to gain insights into their functions. We analyzed the sequences of 13 Obp genes in two chromosomal clusters in a population of wild-derived inbred lines, and asked whether polymorphisms in these genes are associated with variation in olfactory responsiveness. Four polymorphisms in 3 Obp genes exceeded the statistical permutation threshold for association with responsiveness to benzaldehyde, suggesting redundancy and/or combinatorial recognition by these OBPs of this odorant. Model predictions of alternative pre-mRNA secondary structures associated with polymorphic sites suggest that alterations in Obp mRNA structure could contribute to phenotypic variation in olfactory behavior.

Transposable element Ds2 of Zea mays influences polyadenylation and splice site selection

In the allele adh1-2F11 of Zea mays the 1.3 kb transposable element Ds2 is inserted in the fourth exon of Adh1. Two major RNAs of 3.0 and 1.6 kb are transcribed from this allele. While the 3.0 kb transcript also contains the Ds2 sequences, the 1.6 kb transcript has lost the Ds2 sequences by an alternative splicing process. In this process, a normal 5' splice site is joined to a cryptic 3' splice site located within exon 4. This cryptic 3' splice site is not used in a wild-type Adh1F allele. Other minor transcripts of adh1-2F11 are prematurely terminated and polyadenylated.

An artificial riboswitch for controlling pre-mRNA splicing

Riboswitches, as previously reported, are natural RNA aptamers that regulate the expression of numerous bacterial metabolic genes in response to small molecule ligands. It has recently been shown that these RNA genetic elements are also present near the splice site junctions of plant and fungal introns, thus raising the possibility of their involvement in regulating mRNA splicing. Here it is shown for the first time that a riboswitch can be engineered to regulate pre-mRNA splicing in vitro. We show that insertion of a high-affinity theophylline binding aptamer into the 3' splice site (3' ss) region of a model pre-mRNA (AdML-Theo29AG) enables its splicing to be repressed by the addition theophylline. Our results indicate that the location of 3' ss AG within the aptamer plays a crucial role in conferring theophylline-dependent control of pre-mRNA splicing. We also show that theophylline-mediated control of pre-mRNA splicing is highly specific by first demonstrating that a small molecule ligand similar in shape and size to theophylline had no effect on the splicing of AdML-Theo29AG pre-mRNA. Second, theophylline failed to exert any influence on the splicing of a pre-mRNA that does not contain its binding site. Third, theophylline specifically blocks the step II of the splicing reaction. Finally, we provide evidence that theophylline-dependent control of pre-mRNA splicing is functionally relevant.

Alternative splicing of type II procollagen exon 2 is regulated by the combination of a weak 5' splice site and an adjacent intronic stem-loop cis element

Alternative splicing of the type II procollagen gene (COL2A1) is developmentally regulated during chondrogenesis. Chondroprogenitor cells produce the type IIA procollagen isoform by splicing (including) exon 2 during pre-mRNA processing, whereas differentiated chondrocytes synthesize the type IIB procollagen isoform by exon 2 skipping (exclusion). Using a COL2A1 mini-gene and chondrocytes at various stages of differentiation, we identified a non-classical consensus splicing sequence in intron 2 adjacent to the 5' splice site, which is essential in regulating exon 2 splicing. RNA mapping confirmed this region contains secondary structure in the form of a stem-loop. Mutational analysis identified three cis elements within the conserved double-stranded stem region that are functional only in the context of the natural weak 5' splice site of exon 2; they are 1) a uridine-rich enhancer element in all cell types tested except differentiated chondrocytes; 2) an adenine-rich silencer element, and 3) an enhancer cis element functional in the context of secondary structure. This is the first report identifying key cis elements in the COL2A1 gene that modulate the cell type-specific alternative splicing switch of exon 2 during cartilage development.

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

BACKGROUND

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

RESULTS

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

CONCLUSION

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.

Ultraconserved elements in insect genomes: a highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Construction of a herpes simplex virus type 1 mutant with only a three-nucleotide change in the branchpoint region of the latency-associated transcript (LAT) and the stability of its two-kilobase LAT intron

Previous studies using a eukaryotic expression system indicated that the unusual stability of the latency-associated transcript (LAT) intron was due to its nonconsensus branchpoint sequence (T.-T Wu, Y.-H. Su, T. M. Block, and J. M. Taylor, Virology, 243:140-149, 1998). The present study investigated the role of the branchpoint sequence in the stability of the intron expressed from the herpes simplex virus type 1 (HSV-1) genome and the role of LAT intron stability in the HSV-1 life cycle. A branchpoint mutant called Sy2 and the corresponding rescued viruses, SyRA and SyRB, were constructed. To preserve the coding sequence of the immediate early gene icp0, which overlaps with the branchpoint region of the 2-kb LAT, a 3-nucleotide mutation into the branchpoint region of the 2-kb LAT was introduced, resulting in a branchpoint that is 85% identical to the consensus intron branchpoint sequence of eukaryotic cells. As anticipated, there was a 90- to 96-fold reduction in 2-kb LAT accumulation following productive infection in tissue culture and latent infection in mice with Sy2, as determined by Northern blot analysis. These results clearly suggest that the accumulation of the 2-kb intron in tissue culture and in vivo is, at least in part, due to the nonconsensus branchpoint sequence of the LAT intron. Interestingly, a failure to accumulate LAT was associated with greater progeny production of Sy2 at a low multiplicity of infection (0.01) in tissue culture, but not in mice. However, the ability of mutant Sy2 to reactivate from trigeminal ganglia (TG) derived from latently infected mice was indistinguishable from that of wild-type virus, as assayed in the mouse TG explant reactivation system.

Analysis of the LAMB3 gene in a junctional epidermolysis bullosa patient reveals exonic splicing and allele-specific nonsense-mediated mRNA decay

How splicing, the process of intron removal in pre-messenger RNA (mRNA), is carried out with such fidelity in human cells is still not understood, although some general rules are being proposed mainly by in vitro experiments. These rules are currently being redefined by analysis of splicing mechanisms in patients presenting splicing defects. We analysed material of a patient suffering from junctional epidermolysis bullosa, a heritable blistering skin disease. Absence of laminin-5 protein together with hypoplastic hemidesmosomes at the dermo-epidermal junction in the patient's skin was shown by immunohistochemical analysis and immunoelectron microscopy. Subsequent DNA analysis revealed heterozygosity for the mutations R635X and 3009C-->T in the LAMB3 gene. The latter did not alter codon translation, but introduced an exonic splice site in exon 20. Interestingly, this exonic splice site, which presented a splice score of only 68.6, was preferentially used by the spliceosome over the wild-type splice site at the exon 20-intron 20 border, which showed a splice score of 92.2. LAMB3 mRNA was still detectable in RT-PCR analysis although the aberrantly spliced mRNA leads to a stop codon in exon 21, 5' of the commonly assumed 3' border for nonsense-mediated mRNA decay. These results describe an exception to the proposed rules of pre-mRNA splicing and RNA degradation.

Solution structure of the pseudo-5' splice site of a retroviral splicing suppressor

Control of Rous sarcoma virus RNA splicing depends in part on the interaction of U1 and U11 snRNPs with an intronic RNA element called the negative regulator of splicing (NRS). A 23mer RNA hairpin (NRS23) of the NRS directly binds U1 and U11 snRNPs. Mutations that disrupt base-pairing between the loop of NRS23 and U1 snRNA abolish its negative control of splicing. We have determined the solution structure of NRS23 using NOEs, torsion angles, and residual dipolar couplings that were extracted from multidimensional heteronuclear NMR spectra. Our structure showed that the 6-bp stem of NRS23 adopts a nearly A-form duplex conformation. The loop, which consists of 11 residues according to secondary structure probing, was in a closed conformation. U913, the first residue in the loop, was bulged out or dynamic, and loop residues G914-C923, G915-U922, and U916-A921 were base-paired. The remaining UUGU tetraloop sequence did not adopt a stable structure and appears flexible in solution. This tetraloop differs from the well-known classes of tetraloops (GNRA, CUYG, UNCG) in terms of its stability, structure, and function. Deletion of the bulged U913, which is not complementary to U1 snRNA, increased the melting temperature of the RNA hairpin. This hyperstable hairpin exhibited a significant decrease in binding to U1 snRNP. Thus, the structure of the NRS RNA, as well as its sequence, is important for interaction with U1 snRNP and for splicing suppression.

An intronic enhancer regulates splicing of the twintron of Drosophila melanogaster prospero pre-mRNA by two different spliceosomes

We have examined the alternative splicing of the Drosophila melanogaster prospero twintron, which contains splice sites for both the U2- and U12-type spliceosome and generates two forms of mRNA, pros-L (U2-type product) and pros-S (U12-type product). We find that twintron splicing is developmentally regulated: pros-L is abundant in early embryogenesis while pros-S displays the opposite pattern. We have established a Kc cell in vitro splicing system that accurately splices a minimal pros substrate containing the twintron and have examined the sequence requirements for pros twintron splicing. Systematic deletion and mutation analysis of intron sequences established that twintron splicing requires a 46-nucleotide purine-rich element located 32 nucleotides downstream of the U2-type 5' splice site. While this element regulates both splicing pathways, its alteration showed the severest effects on the U2-type splicing pathway. Addition of an RNA competitor containing the wild-type purine-rich element to the Kc extract abolished U2-type splicing and slightly repressed U12-type splicing, suggesting that a trans-acting factor(s) binds the enhancer element to stimulate twintron splicing. Thus, we have identified an intron region critical for prospero twintron splicing as a first step towards elucidating the molecular mechanism of splicing regulation involving competition between the two kinds of spliceosomes.

A Non-sequence-specific double-stranded RNA structural element regulates splicing of two mutually exclusive exons of fibroblast growth factor receptor 2 (FGFR2)

Alternative splicing of fibroblast growth factor receptor 2 (FGFR2) mutually exclusive exons IIIb and IIIc represents a tightly regulated and functionally relevant example of post-transcriptional gene regulation. Rat prostate cancer DT3 and AT3 cell lines demonstrate exclusive selection of either exon IIIb or exon IIIc, respectively, and have been used to characterize regulatory FGFR2 RNA cis-elements that are required for splicing regulation. Two sequences termed ISE-2 and ISAR are located in the intron between exons IIIb and IIIc and are required for cell-type specific exon IIIb. Previous studies suggest that the function of these elements involves formation of an RNA stem structure, even though they are separated by more than 700 nucleotides. Using transfected minigenes, we performed a systematic analysis of the sequence and structural components of ISE-2 and ISAR that are required for their ability to regulate FGFR2 splicing. We found that the primary sequence of these elements can be replaced by completely unrelated sequences, provided that they are also predicted to form an RNA stem structure. Thus, a nonsequence-specific double stranded RNA stem constitutes a functional element required for FGFR2 splicing; suggesting that a double-stranded RNA binding protein is a component of the splicing regulatory machinery.

Stimulatory effect of splicing factors on transcriptional elongation

Transcription and pre-mRNA splicing are tightly coupled gene expression events in eukaryotic cells. An interaction between the carboxy-terminal domain of the largest subunit of RNA polymerase (Pol) II and components of the splicing machinery is postulated to mediate this coupling. Here, we show that splicing factors function directly to promote transcriptional elongation, demonstrating that transcription is more intimately coupled to splicing than previously thought. The spliceosomal U small nuclear ribonucleoproteins (snRNPs) interact with human transcription elongation factor TAT-SF1 (refs 6,7,8,9) and strongly stimulate polymerase elongation when directed to an intron-free human immunodeficiency virus-1 (HIV-1) template. This effect is likely to be mediated through the binding of TAT-SF1 to elongation factor P-TEFb, a proposed component of the transcription elongation complex. Inclusion of splicing signals in the nascent transcript further stimulates transcription, supporting the notion that the recruitment of U snRNPs near the elongating polymerase is important for transcription. Because the TAT-SF1-U snRNP complex also stimulates splicing in vitro, it may serve as a dual-function factor to couple transcription and splicing and to facilitate their reciprocal activation.

Five novel alternatively spliced transcripts of DNA (cytosine-5) methyltransferase 2 in human peripheral blood leukocytes

Alternative splicing of RNA molecules transcribed from DNA (cytosine-5) methyltransferases has been proposed as a mechanism by which methylation is able to effect diverse biological processes in higher eukaryotes. This study has investigated transcriptional versatility of DNA (cytosine-5) methyltransferase 2, which may methylate cytosine residues within 5'-CCTGG-3' pentanucleotides in regions of the human genome devoid of 5'-CG-3' methylation. Five novel splice variants of DNA (cytosine-5) methyltransferase 2 were identified in the peripheral blood leukocytes of healthy subjects following cloning and sequencing of RT-PCR products amplified using gene specific oligodeoxyribonucleotide primers. The generation of some of these splice variants may be influenced by the formation of secondary structures within pre-mRNA due to the repetition of sequences flanking alternatively spliced exons in a reverse and complementary orientation on the same strand. These findings enable novel approaches to investigate the role of RNA secondary structures in alternative splicing. The DNA (cytosine-5) methyltransferase 2 splice variants are generated in all the major cell types of peripheral blood, as well as in neoplastic lymphoid cells indicating that they are unlikely to generate proteins involved in control of the cell cycle or cellular differentiation. Interestingly, the gene products generated by some splice variants completely or partially lack highly conserved amino acid motifs shown to be important for the catalysis of cytosine methylation. The possibility cannot be excluded, therefore, that alternative splicing of DNA (cytosine-5) methyltransferase 2 pre-mRNA may generate protein isoforms which have different methylating capabilities or which are involved in biological processes other than the catalysis of cytosine methylation.

Mutational analysis of 3' splice site selection during trans-splicing

trans-Splicing is essential for mRNA maturation in trypanosomatids. A conserved AG dinucleotide serves as the 3' splice acceptor site, and analysis of native processing sites suggests that selection of this site is determined according to a 5'-3' scanning model. A series of stable gene replacement lines were generated that carried point mutations at or near the 3' splice site within the intergenic region separating CUB2.65, the calmodulin-ubiquitin associated gene, and FUS1, the ubiquitin fusion gene of Trypanosoma cruzi. In one stable line, the elimination of the native 3' splice acceptor site led to the accumulation of Y-branched splicing intermediates, which served as templates for mapping the first trans-splicing branch points in T. cruzi. In other lines, point mutations shifted the position of the first consensus AG dinucleotide either upstream or downstream of the wild-type 3' splice acceptor site in this intergenic region. Consistent with the scanning model, the first AG dinucleotide downstream of the branch points was used as the predominant 3' splice acceptor site. In all of the stable lines, the point mutations affected splicing efficiency in this region.

Aberrant splicing of tau pre-mRNA caused by intronic mutations associated with the inherited dementia frontotemporal dementia with parkinsonism linked to chromosome 17

Frontotemporal dementia accounts for a significant fraction of dementia cases. Frontotemporal dementia with parkinsonism linked to chromosome 17 is associated with either exonic or intronic mutations in the tau gene. This highlights the involvement of aberrant pre-mRNA splicing in the pathogenesis of neurodegenerative disorders. Little is known about the molecular mechanisms of the splicing defects underlying these diseases. To establish a model system for studying the role of pre-mRNA splicing in neurodegenerative diseases, we have constructed a tau minigene that reproduces tau alternative splicing in both cultured cells and in vitro biochemical assays. We demonstrate that mutations in a nonconserved intronic region of the human tau gene lead to increased splicing between exon 10 and exon 11. Systematic biochemical analyses indicate the importance of U1 snRNP and, to a lesser extent, U6 snRNP in differentially recognizing wild-type versus intron mutant tau pre-mRNAs. Gel mobility shift assays with purified U1 snRNP and oligonucleotide-directed RNase H cleavage experiments support the idea that the intronic mutations destabilize a stem-loop structure that sequesters the 5' splice site downstream of exon 10 in tau pre-mRNA, leading to increases in U1 snRNP binding and in splicing between exon 10 and exon 11. Thus, mutations in nonconserved intronic regions that increase rather than decrease alternative splicing can be an important pathogenic mechanism for the development of human diseases.