Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Splicing of vertebrate introns involves recognition of three consensus elements at the 3′ end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3′ splice site (3′ ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3′ ss is the absence of intervening AG dinucleotides, leading to its designation as the “AG exclusion zone” (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.     

Restriction of an intron size en route to endothermy

Ca²⁺-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5′ splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5′ splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca²⁺ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca²⁺-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.     

Branch point selection in alternative splicing of tropomyosin pre-mRNAs.

The rat tropomyosin 1 gene gives rise to two mRNAs encoding rat fibroblast TM-1 and skeletal muscle beta-tropomyosin via an alternative splicing mechanism. The gene is comprised of 11 exons. Exons 1 through 5 and exons 8 and 9 are common to all mRNAs expressed from this gene. Exons 6 and 11 are used in fibroblasts as well as smooth muscle whereas exons 7 and 10 are used exclusively in skeletal muscle. In the present studies we have focused on the mutually exclusive internal alternative splice choice involving exon 6 (fibroblast-type splice) and exon 7 (skeletal muscle-type splice). To study the mechanism and regulation of alternative splice site selection we have characterized the branch points used in processing of the tropomyosin pre-mRNAs in vitro using nuclear extracts obtained from HeLa cells. Splicing of exon 5 to exon 6 (fibroblast-type splice) involves the use of three branch points located 25, 29, and 36 nucleotides upstream of the 3' splice site of exon 6. Splicing of exon 6 (fibroblast-type splice) or exon 7 (skeletal muscle type-splice) to exon 8 involves the use of the same branch point located 24 nucleotides upstream of this shared 3' splice site. In contrast, the splicing of exon 5 to exon 7 (skeletal muscle-type splice) involves the use of three branch sites located 144, 147 and 153 nucleotides, upstream of the 3' splice site of exon 7. In addition, the pyrimidine content of the region between these unusual branch points and the 3' splice site of exon 7 was found to be greater than 80%. These studies raise the possibility that the use of branch points located a long distance from a 3' splice site may be an essential feature of some alternatively spliced exons. The possible significance of these unusual branch points as well as a role for the polypyrimidine stretch in intron 6 in splice site selection are discussed.     

A sequence compilation and comparison of exons that are alternatively spliced in neurons.

Alternative splicing is an important regulatory mechanism to create protein diversity. In order to elucidate possible regulatory elements common to neuron specific exons, we created and statistically analysed a database of exons that are alternatively spliced in neurons. The splice site comparison of alternatively and constitutively spliced exons reveals that some, but not all alternatively spliced exons have splice sites deviating from the consensus sequence, implying diverse patterns of regulation. The deviation from the consensus is most evident at the -3 position of the 3' splice site and the +4 and -3 position of the 5' splice site. The nucleotide composition of alternatively and constitutively spliced exons is different, with alternatively spliced exons being more AU rich. We performed overlapping k-tuple analysis to identify common motifs. We found that alternatively and constitutively spliced exons differ in the frequency of several trinucleotides that cannot be explained by the amino acid composition and may be important for splicing regulation.     

Accumulation of GC donor splice signals in mammals

Abstract

The GT dinucleotide in the first two intron positions is the most conserved element of the U2 donor splice signals. However, in a small fraction of donor sites, GT is replaced by GC. A substantial enrichment of GC in donor sites of alternatively spliced genes has been observed previously in human, nematode and Arabidopsis, suggesting that GC signals are important for regulation of alternative splicing. We used parsimony analysis to reconstruct evolution of donor splice sites and inferred 298 GT > GC conversion events compared to 40 GC > GT conversion events in primate and rodent genomes. Thus, there was substantive accumulation of GC donor splice sites during the evolution of mammals. Accumulation of GC sites might have been driven by selection for alternative splicing.

Reviewers

This article was reviewed by Jerzy Jurka and Anton Nekrutenko. For the full reviews, please go to the Reviewers' Reports section.     

Both the polypyrimidine tract and the 3' splice site function prior to the first step of splicing in fission yeast.

While it is known that several trans -acting splicing factors are highly conserved between Schizosaccharomyces pombe and mammals, the roles of cis -acting signals have received comparatively little attention. In Saccharomyces cerevisiae, sequences downstream from the branch point are not required prior to the first transesterification reaction, whereas in mammals the polypyrimidine tract and, in some introns, the 3' AG dinucleotide are critical for initial recognition of an intron. We have investigated the contribution of these two sequence elements to splicing in S.pombe. To determine the stage at which the polypyrimidine tract functions, we analyzed the second intron of the cdc2 gene (cdc 2-Int2), in which pyrimidines span the entire interval between the branch point and 3' splice site. Our data indicate that substitution of a polypurine tract results in accumulation of linear pre-mRNA, while expanding the polypyrimidine tract enhances splicing efficiency, as in mammals. To examine the role of the AG dinucleotide in cdc 2-Int2 splicing, we mutated the 3' splice junction in both the wild-type and pyrimidine tract variant RNAs. These changes block the first transesterification reaction, as in a subset of mammalian introns. However, in contrast to the situation in mammals, we were unable to rescue the first step of splicing in a 3' splice site mutant by expanding the polypyrimidine tract. Mutating the terminal G in the third intron of the nda 3 gene (nda 3-Int3) also blocks the first transesterification reaction, suggesting that early recognition of the 3' splice site is a general property of fission yeast introns. Counter to earlier work with an artificial intron, it is not possible to restore the first step of splicing in cdc 2-Int2 and nda 3-Int3 3' splice site mutants by introducing compensatory changes in U1 snRNA. These results highlight the diversity and probable redundancy of mechanisms for identifying the 3' ends of introns.     

The conserved dinucleotide AG of the 3′ splice site may be recognized twice during in vitro splicing of mammalian mRNA precursors

We have previously used site-directed mutagenesis to introduce an additional branch site into the first intron of the human β-globin gene at nt −24 between the natural branch site (nt−37) and the normal 3′ splice site (nt−1). We found that either the upstream or downstream branch site could be used during in vitro splicing, depending on which site best matched the mammalian branch site consensus YURAC (R = purine; Y = pyrimidine). Here we show that introduction of an additional AG dinucleotide at nt −20 between the downstream branch site and the normal 3′ splice site results in alternative 3′ splicing. Splicing to the new AG uses the upstream branch site exclusively, presumably because the downstream branch site is only 4 nt from this 3′ splice site. We were surprised, however, to find that the presence of the new AG also prevents the use of the upstream branch site for splicing to the normal 3′ splice site. Analysis of additional mutants confirmed earlier work [Krainer et al.: Mechanisms of human β-globin pre-mRNA splicing. In Berg, P. (Ed.), The Robert A. Welch Foundation Conferences on Chemical Research XXIX. Genetic Chemistry: The Molecular Basis of Heredity. Welch Foundation, Houston, TX, 1985, pp. 353–382] that the new AG cannot function by itself as a complete 3′ splice site; rather, it appears that alternative 3′ splicing initiates at the normal 3′ splices site but then searches, once the reaction is underway, for the first AG downstream from the chosen branch site. Taken together, our data suggest that the conserved AG dinucleotide at the 3′ splice site may be recognized twice during mammalian mRNA splicing in vitro.     

Intron sequences involved in lariat formation during pre-mRNA splicing

We have shown that lariat formation during in vitro splicing of several RNA precursors, from Drosophila to man, occurs at a unique and identifiable but weakly conserved site, 18 to 37 nucleotides proximal to the 3' splice site. Lariat formation within an artificial intron lacking a normal branch-point sequence occurs at a cryptic site a conserved distance (approximately 23 nucleotides) from the 3' splice site. Analysis of beta-thalassemia splicing mutations revealed that lariat formation in the first intron of the human beta-globin gene occurs at the same site in normal and mutant precursors, even though alternate 5' and 3' splice sites are utilized in the mutants. Remarkably, cleavage at the 5' splice site and lariat formation do not occur when the precursor contains a beta-thalassemia deletion removing the polypyrimidine stretch and AG dinucleotide at the 3' splice site. In contrast, a single base substitution in the AG dinucleotide blocks cleavage at the 3' splice site but not at the 5' site.     

In vitro splicing of mutually exclusive exons from the chicken beta-tropomyosin gene: role of the branch point location and very long pyrimidine stretch.

The chicken beta-tropomyosin gene contains 11 exons, two of which are spliced into mRNA only in skeletal muscle. One pair of alternative exons, 6A and 6B, is found in the middle of the gene; they are spliced in a mutually exclusive manner. The non-muscle splice 6A-7 is by far the predominant in vitro reaction in a HeLa cell nuclear extract. A minor product is the 6A-6B splice, which is excluded in all tissues. This minor product results from the use of a branch point located 105 nt upstream of the 3' end of the intron separating exons 6A and 6B. The region between the branch point sequence and the final AG contains a stretch of approximately 80 pyrimidines. We have examined the role of the distance of the branchpoint to the 3' splice site and of the sequences between these two elements. Our results suggest that at least two cis-acting elements contribute to the mutual exclusivity of exons 6A and 6B. The intron between exons 6A and 6B is intrinsically poorly 'spliceable' both because the branch point is too far upstream of the 3' end of the intron to give efficient splicing and because of the particular sequence lying between this branch point and the 3' splice site.     

Functional analysis of the polypyrimidine tract in pre-mRNA splicing.

The polypyrimidine tract is one of the important cis-acting sequence elements directing intron removal in pre-mRNA splicing. Progressive deletions of the polypyrimidine tract have been found to abolish correct lariat formation, spliceosome assembly and splicing. In addition, the polypyrimidine tract can alter 3'-splice site selection by promoting alternative branch site selection. However, there appears to be great flexibility in the specific sequence of a given tract. Not only the optimal composition of the polypyrimidine tract, but also the role of the tract in introns with no apparent polypyrimidine tracts or where changes in the tract are apparently harmless are uncertain. Accordingly, we have designed a series of cis-competition splicing constructs to test the functional competitive efficiency of a variety of systematically mutated polypyrimidine tracts. An RT/PCR assay was used to detect spliced product formation as a result of differential branch point selection dependent on direct competition between two opposing polypyrimidine tracts. We found that pyrimidine tracts containing 11 continuous uridines are the strongest pyrimidine tracts. In such cases, the position of the uridine stretch between the branch point and 3'-splice site AG is unimportant. In contrast, decreasing the continuous uridine stretch to five or six residues requires that the tract be located immediately adjacent to the AG for optimal competitive efficiency. The block to splicing with decreasing polypyrimidine tract strength is primarily prior to the first step of splicing. While lengthy continuous uridine tracts are the most competitive, tracts with decreased numbers of consecutive uridines and even tracts with alternating purine/pyrimidine residues can still function to promote branch point selection, but are far less effective competitors in 3'-splice site selection assays.     

Splicing of internal large exons is defined by novel cis-acting sequence elements

Human internal exons have an average size of 147 nt, and most are <300 nt. This small size is thought to facilitate exon definition. A small number of large internal exons have been identified and shown to be alternatively spliced. We identified 1115 internal exons >1000 nt in the human genome; these were found in 5% of all protein-coding genes, and most were expressed and translated. Surprisingly, 40% of these were expressed at levels similar to the flanking exons, suggesting they were constitutively spliced. While all of the large exons had strong splice sites, the constitutively spliced large exons had a higher ratio of splicing enhancers/silencers and were more conserved across mammals than the alternatively spliced large exons. We asked if large exons contain specific sequences that promote splicing and identified 38 sequences enriched in the large exons relative to small exons. The consensus sequence is C-rich with a central invariant CA dinucleotide. Mutation of these sequences in a candidate large exon indicated that these are important for recognition of large exons by the splicing machinery. We propose that these sequences are large exon splicing enhancers (LESEs).     

The exon sequence TAGG can inhibit splicing.

The fibroblast growth factor receptor-2 gene contains a pair of alternative exons, K-SAM and BEK, which are spliced in a cell type specific manner. We have shown previously that a 10 nucleotide sequence within the K-SAM exon exerts a negative effect on K-SAM exon splicing independent of cell type. We demonstrate here that this sequence works autonomously, as it can repress splicing of a heterologous exon, the EIIIb alternative exon of the rat fibronectin gene. By introducing point mutations into the 10 nucleotide sequence, we have shown that the functional portion is limited to 4 nucleotides, TAGG, the dinucleotide AG of which is particularly important. This short sequence may participate in the control of splicing of exons carrying it, provided that they carry weak splice sites.     

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.     

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.