More than one way to splice an RNA: branching without a bulge and splicing without branching in group II introns.

Domain 6 (D6) of group II introns contains a bulged adenosine that serves as the branch-site during self-splicing. In addition to this adenosine, other structural features in D6 are likely to contribute to the efficiency of branching. To understand their role in promoting self-splicing, the branch-site and surrounding nucleotides were mutagenized. Detailed kinetic analysis on the self-splicing efficiency of the mutants revealed several interesting features. First, elimination of the branch-site does not preclude efficient splicing, which takes place instead through a hydrolytic first step. Second, pairing of the branch-site does not eliminate branching, particularly if the adenosine is involved in a mispair. Third, the G-U pairs that often surround group II intron branch-points contribute to the efficiency of branching. These results suggest that there is a strong driving force for promoting self-splicing by group II introns, which employ a versatile set of different mechanisms for ensuring that splicing is successful. In addition, the behavior of these mutants indicates that a bulged adenosine per se is not the important determinant for branch-site recognition in group II introns. Rather, the data suggest that the branch-site adenosine is recognized as a flipped base, a conformation that can be promoted by a variety of different substructures in RNA and DNA.     

The receptor for branch-site docking within a group II intron active site

The distinguishing feature of group II introns, and the property that links them with spliceosomal catalysis, is their ability to undergo splicing through branching. In this reaction, the 2'-hydroxyl group of a specific adenosine within intron domain 6 serves as the nucleophile for attack on the 5' splice site. We know less about branching than any other feature of group II intron catalysis, largely because the receptor structure for activating the branch site is unknown. Here, we identify the intronic region that binds the branch site of a group IIB intron. Located in domain 1, close to receptors for intron domain 5 and both splice sites, we demonstrate that the branch-site receptor is a functional element required for transesterification. Furthermore, we show that crosslinked branch sites can carry out both steps of splicing, suggesting that the conformational state of the intron core is set early and that it persists throughout the entire splicing process.     

Adenosine is inherently favored as the branch-site RNA nucleotide in a structural context that resembles natural RNA splicing

We previously used in vitro selection to identify the 7S11 deoxyribozyme, which catalyzes formation of 2',5'-branched RNA using a branch-site adenosine nucleophile and a 5'-triphosphate electrophile. An unanswered question is whether the use of branch-site adenosine is inherently preferred or a chance event during the particular selection experiment. Here we have found that deoxyribozymes newly selected to use uridine as the branch-site RNA nucleotide in a structural context that resembles natural RNA splicing instead prefer a branch-site adenosine, although adenosine was never available during the selection itself. Our results support a chemical basis for nature's choice of the branch-site nucleotide, which is almost always adenosine in group II introns and the spliceosome.     

Lariat formation and a hydrolytic pathway in plant chloroplast group II intron splicing

Lariat formation has been studied intensively only with a few self-splicing group II introns, and little is known about how the numerous diverse introns in plant organelles are excised. Several of these introns have branch-points that are not a single bulge but are adjoined by A:A, A:C, A:G and G:G pairs. Using a highly sensitive in vivo approach, we demonstrate that all but one of the barley chloroplast introns splice via the common pathway that produces a branched product. RNA editing does not improve domain 5 and 6 structures of these introns. The conserved branch-point in tobacco rpl16 is chosen even if an adjacent unpaired adenosine is available, suggesting that spatial arrangements in domain 6 determine correct branch-point selection. Lariats were not detected for the chloroplast trnV intron, which lacks an unpaired adenosine in domain 6. Instead, this intron is released as linear molecules that undergo further polyadenylation. trnV, which is conserved throughout plant evolution, constitutes the first example of naturally occurring hydrolytic group II intron splicing in vivo.     

A conserved pseudouridine modification in eukaryotic U2 snRNA induces a change in branch-site architecture.

The removal of noncoding sequences (introns) from eukaryotic precursor mRNA is catalyzed by the spliceosome, a dynamic assembly involving specific and sequential RNA-RNA and RNA-protein interactions. An essential RNA-RNA pairing between the U2 small nuclear (sn)RNA and a complementary consensus sequence of the intron, called the branch site, results in positioning of the 2'OH of an unpaired intron adenosine residue to initiate nucleophilic attack in the first step of splicing. To understand the structural features that facilitate recognition and chemical activity of the branch site, duplexes representing the paired U2 snRNA and intron sequences from Saccharomyces cerevisiae were examined by solution NMR spectroscopy. Oligomers were synthesized with pseudouridine (psi) at a conserved site on the U2 snRNA strand (opposite an A-A dinucleotide on the intron strand, one of which forms the branch site) and with uridine, the unmodified analog. Data from NMR spectra of nonexchangeable protons demonstrated A-form helical backbone geometry and continuous base stacking throughout the unmodified molecule. Incorporation of psi at the conserved position, however, was accompanied by marked deviation from helical parameters and an extrahelical orientation for the unpaired adenosine. Incorporation of psi also stabilized the branch-site interaction, contributing -0.7 kcal/mol to duplex deltaG degrees 37. These findings suggest that the presence of this conserved U2 snRNA pseudouridine induces a change in the structure and stability of the branch-site sequence, and imply that the extrahelical orientation of the branch-site adenosine may facilitate recognition of this base during splicing.     

Mimicking the first step of RNA splicing: an artificial DNA enzyme can synthesize branched RNA using an oligonucleotide leaving group as a 5'-exon analogue

The 7S11 deoxyribozyme synthesizes 2',5'-branched RNA by mediating the nucleophilic attack of an internal 2'-hydroxyl group of one RNA substrate into the 5'-triphosphate of a second RNA substrate, with pyrophosphate as the leaving group. Here we comprehensively examined the role of the leaving group in the 7S11-catalyzed reaction by altering the 5'-phosphorylation state and the length of the second RNA substrate. When the leaving group is the less stabilized phosphate or hydroxide anion as provided by a 5'-diphosphate or 5'-monophosphate, the same 2',5'-branched product is formed as when pyrophosphate is the leaving group, but with an approximately 50- or approximately 1000-fold lower rate (Br?nsted beta(LG) = -0.40). When the 5'-end of the RNA substrate that bears the leaving group is longer by one or more nucleotides, either the new 5'-terminal alpha-phosphate or the original alpha-phosphate can be attacked by the branch-site 2'-hydroxyl group; in the latter case, the leaving group is an oligonucleotide. The choice between these alpha-phosphate reaction sites is determined by the subtle balance between the length of the single-stranded 5'-extension and the stability of the leaving group. Because the branch-site adenosine is a bulged nucleotide flanked by Watson-Crick duplex regions, we earlier concluded that 7S11 structurally mimics the first step of natural RNA splicing. The observation of 7S11-catalyzed branch formation with an oligonucleotide leaving group strengthens this resemblance to natural RNA splicing, with the oligonucleotide playing the role of the 5'-exon in the first step. These findings reinforce the notion that splicing-related catalysis can be achieved by artificial nucleic acid enzymes that are much smaller than the spliceosome and group II introns.     

Group II introns in wheat mitochondria have degenerate structural features and varied splicing pathways

Mitochondrial introns in flowering plant genes are virtually all classified as members of the group II ribozyme family although certain structural features have degenerated to varying degrees over evolutionary time. We are interested in the impact that unconventional intron architecture might have on splicing biochemistry in vivo and we have focused in particular on intronic domains V and VI, which for self-splicing introns provide a key component of the catalytic core and the bulged branchpoint adenosine, respectively. Notably, the two transesterification steps in classical group II splicing are the same as for nuclear spliceosomal introns and release the intron as a lariat. Using RT-PCR and circularized RT-PCR, we had previously demonstrated that several wheat mitochondrial introns which lack a branchpoint adenosine have atypical splicing pathways, and we have now extended this analysis to the full set of wheat introns, namely six trans-splicing and sixteen cis-splicing ones. A number of introns are excised using non-lariat pathways and interestingly, we find that several introns which do have a conventional domain VI also use pathways that appear to exploit other internal or external nucleophiles, with the lariat form being relatively minor. Somewhat surprisingly, several introns with weakly-structured domain V/VI helices still exhibit classical lariat splicing, suggesting that accessory factors aid in restoring a splicing-competent conformation. Our observations illustrate that the loss of conventional group II features during evolution is correlated with altered splicing biochemistry in an intron-distinctive manner.     

Group II intron in Bacillus cereus has an unusual 3' extension and splices 56 nucleotides downstream of the predicted site

All group II introns known to date fold into six functional domains. However, we recently identified an intron in Bacillus cereus ATCC 10987, B.c.I4, that splices 56nt downstream of the expected 3′ splice site in vivo (Tourasse et al. 2005, J. Bacteriol., 187, 5437–5451). In this study, we confirmed by ribonuclease protection assay that the 56-bp segment is part of the intron RNA molecule, and computational prediction suggests that it might form a stable stem-loop structure downstream of domain VI. The splicing of B.c.I4 was further investigated both in vivo and in vitro. Lariat formation proceeded primarily by branching at the ordinary bulged adenosine in domain VI without affecting the fidelity of splicing. In addition, the splicing efficiency of the wild-type intron was better than that of a mutant construct deleted of the 56-bp 3′ extension. These results indicate that the intron has apparently adapted to the extra segment, possibly through conformational adjustments. The extraordinary group II intron B.c.I4 harboring an unprecedented extra 3′ segment constitutes a dramatic example of the flexibility and adaptability of group II introns.     

Group II introns deleted for multiple substructures retain self-splicing activity.

Group II introns can be folded into highly conserved secondary structures with six major substructures or domains. Domains 1 and 5 are known to play key roles in self-splicing, while the roles of domains 2, 3, 4, and 6 are less clear. A trans assay for domain 5 function has been developed which indicates that domain 5 has a binding site on the precursor RNA that is not predicted from any secondary structure element. In this study, the self-splicing group II intron 5 gamma of the coxI gene of yeast mitochondrial DNA was deleted for various intron domains, singly and in combinations. Those mutant introns were characterized for self-splicing reactions in vitro as a means of locating the domain 5 binding site. A single deletion of domain 2, 3, 4, or 6 does not block in vitro reactions at either splice junction, though the deletion of domain 6 reduces the fidelity of 3' splice site selection somewhat. Even the triple deletion lacking domains 2, 4, and 6 retains some self-splicing activity. The deletion of domains 2, 3, 4, and 6 blocks the reaction at the 3' splice junction but not at the 5' junction. From these results, we conclude that the binding site for domain 5 is within domain 1 and that the complex of 5' exon, domain 1, and domain 5 (plus short connecting sequences) constitutes the essential catalytic core of this intron.     

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5′ splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5′ and 3′ splice sites, and the activation of cryptic U12-dependent branch/3′ splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3′ splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.     

A random effects branch-site model for detecting episodic diversifying selection

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.     

Three recognition events at the branch-site adenine.

An adenosine at the branch site, the nucleophile for the first transesterification step of splicing, is nearly invariant in mammalian pre-mRNA introns. The chemical groups on the adenine base were varied systematically and assayed for formation of early spliceosome complexes and execution of the first and second steps of splicing. Recognition of constituents of the adenine is critical in formation of a U2 snRNP-containing complex on a minimal branch-site oligonucleotide. Furthermore, the efficiencies of the first and second chemical steps have different dependencies on the functional groups of the adenine. In total, the chemical groups on the adenine base at the branch site are differentially recognized during at least three different processes in the splicing of pre-mRNA. Moreover, a protein, p14, interacts with the adenine in a base-specific fashion and may mediate early recognition of this base.     

X-ray structures of U2 snRNA-branchpoint duplexes containing conserved pseudouridines

A pseudouridine-modified region of the U2 small nuclear (sn)RNA anneals with the intronic branchpoint sequence and positions a bulged adenosine to serve as the nucleophile in the first chemical step of pre-mRNA splicing. We have determined three X-ray structures of RNA oligonucleotides containing the pseudouridylated U2 snRNA and the branchpoint consensus sequences. The expected adenosine branchpoint is extrahelical in a 1.65 A resolution structure containing the mammalian consensus sequence variant and in a 2.10 A resolution structure containing a shortened Saccharomyces cerevisiae consensus sequence. The adenosine adjacent to the expected branchpoint is extrahelical in a third structure, which contains the intact yeast consensus sequence at 1.57 A resolution. The hydration and base stacking interactions mediated by the U2 snRNA pseudouridines correlate with the identity of the unpaired adenosine. The expected adenosine bulge is associated with a well-stacked pseudouridine, which is linked via an ordered water molecule to a neighboring nucleotide. In contrast, the bulge of the adjacent adenosine shifts the base stacking and disrupts the water-mediated interactions of the pseudouridine. These structural differences may contribute to the ability of the pseudouridine modification to promote the bulged conformation of the branch site adenosine and to enhance catalysis by snRNAs. Furthermore, iodide binding sites are identified adjacent to the unconventional bulged adenosine, and the structure of the mammalian consensus sequence variant provides a high-resolution view of a hydrated magnesium ion bound in a similar manner to a divalent cation binding site of the group II intron.     

False-positive results obtained from the branch-site test of positive selection

Natural selection operating at the amino acid sequence level can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitutions, where r(N)/r(S) (omega) > 1 and omega < 1 suggest positive and negative selection, respectively. The branch-site test has been developed for detecting positive selection operating at a group of amino acid sites for a pre-specified (foreground) branch of a phylogenetic tree by taking into account the heterogeneity of omega among sites and branches. Here the performance of the branch-site test was examined by computer simulation, with special reference to the false-positive rate when the divergence of the sequences analyzed was small. The false-positive rate was found to inflate when the assumptions made on the omega values for the foreground and other (background) branches in the branch-site test were violated. In addition, under a similar condition, false-positive results were often obtained even when Bonferroni correction was conducted and the false-discovery rate was controlled in a large-scale analysis. False-positive results were also obtained even when the number of nonsynonymous substitutions for the foreground branch was smaller than the minimum value required for detecting positive selection. The existence of a codon site with a possibility of occurrence of multiple nonsynonymous substitutions for the foreground branch often caused the branch-site test to falsely identify positive selection. In the re-analysis of orthologous trios of protein-coding genes from humans, chimpanzees, and macaques, most of the genes previously identified to be positively selected for the human or chimpanzee branch by the branch-site test contained such a codon site, suggesting a possibility that a significant fraction of these genes are false-positives.     

Relevance of the branch point adenosine, coordination loop, and 3' exon binding site for in vivo excision of the Sinorhizobium meliloti group II intron RmInt1

Excision of the bacterial group II intron RmInt1 has been demonstrated in vivo, resulting in the formation of both intron lariat and putative intron RNA circles. We show here that the bulged adenosine in domain VI of RmInt1 is required for splicing via the branching pathway, but branch site mutants produce small numbers of RNA molecules in which the first G residue of the intron is linked to the last C residue. Mutations in the coordination loop in domain I reduced splicing efficiency, but branched templates clearly predominated among splicing products. We also found that a single substitution at the EBS3 position (G329C), preventing EBS3-IBS3 pairing, resulted in the production of 50 to 100 times more RNA molecules in which the 5' and 3' extremities were joined. We provide evidence that these intron molecules may correspond to both, intron circles linked by a 2'-5' phosphodiester bond, and tandem, head-to-tail intron copies.     

Linking the branchpoint helix to a newly found receptor allows lariat formation by a group II intron

Like spliceosomal introns, the ribozyme-containing group II introns are excised as branched, lariat structures: a 2'-5' bond is created between the first nucleotide of the intron and an adenosine in domain VI, a component which is missing from available crystal structures of the ribozyme. Comparative sequence analysis, modelling and nucleotide substitutions point to the existence, and probable location, of a specific RNA receptor for the section of domain VI that lies just distal to the branchpoint adenosine. By designing oligonucleotides that tether domain VI to this novel binding site, we have been able to specifically activate lariat formation in an engineered, defective group II ribozyme. The location of the newly identified receptor implies that prior to exon ligation, the distal part of domain VI undergoes a major translocation, which can now be brought under control by the system of anchoring oligonucleotides we have developed. Interestingly, these oligonucleotides, which link the branchpoint helix and the binding site for intron nucleotides 3-4, may be viewed as counterparts of U2-U6 helix III in the spliceosome.     

Self-splicing of group II introns in vitro: mapping of the branch point and mutational inhibition of lariat formation

Group II intron bl1 from yeast mitochondria can undergo self-splicing in vitro. Exons become correctly ligated, and the excised intron has a lariat structure similar to that of introns from nuclear mRNA. The branch point of the bl1 lariat is located eight or nine nucleotides upstream of the 3' end of the intron and is part of a hairpin structure that is well conserved among group II introns. Several mutations next to the branch point and in other parts of the core structure of group II introns are shown to affect lariat formation. One of them, carried by strain M4873, abolishes splicing in vivo and in vitro, apparently by changing the architecture of the hairpin structure containing the branch point. Similarities between group II introns and nuclear pre-mRNA introns are discussed in terms of evolutionary relatedness.     

An RNA conformational change between the two chemical steps of group II self-splicing.

As for nuclear pre-mRNA introns, the splicing pathway of group II self-splicing introns proceeds by two successive transesterifications involving substrates with different chemical configurations. These two reactions have been proposed to be catalysed by two active sites, or alternatively by a single active site rearranging its components to accommodate the successive substrates. Here we show that the structural elements specific for the second splicing step are clustered in peripheral structures of domains II and VI. We show that these structures are not required for catalysis of the second chemical step but, instead, take part in a conformational change that occurs between the two catalytic steps. This rearrangement involves the formation of a tertiary contact between part of domain II and a GNRA tetraloop at the tip of domain VI. The fact that domain VI, which carries the branched structure, is involved in this structural rearrangement and the fact that modifications affecting the structures involved have almost no effect when splicing proceeds without branch formation, suggest that the conformational change results in the displacement of the first-step product out of the active site. These observations give further support to the existence of a single active site in group II introns.     

Identification of a U2/U6 helix la mutant that influences 3' splice site selection during nuclear pre-mRNA splicing

Base substitutions in U2/U6 helix I, a conserved base-pairing interaction between the U6 and U2 snRNAs, have previously been found to specifically block the second catalytic step of nuclear pre-mRNA splicing. To further assess the role of U2/U6 helix I in the second catalytic step, we have screened mutations in U2/U6 helix I to identify those that influence 3' splice site selection using a derivative of the yeast actin pre-mRNA. In these derivatives, the spacing between the branch site adenosine and 3' splice site has been reduced from 43 to 12 nt and this results in enhanced splicing of mutants in the conserved 3' terminal intron residue. In this context, mutation of the conserved 3' intron terminal G to a C also results in the partial activation of a nearby cryptic 3' splice site with U as the 3' terminal intron nucleotide. Using this highly sensitive mutant substrate, we have identified a mutation in the U6 snRNA (U57A) that significantly increases the selection of the cryptic 3' splice site over the normal 3' splice site and augments its utilization relative to that observed with the wild-type U2 or U6 snRNAs. In a previous study, we found that the same U6 mutation suppressed the effects of an A-to-G branch site mutation in an allele-specific fashion. The ability of U6-U57 mutants to influence the fidelity of both branch site and 3' splice site recognition suggests that this nucleotide may participate in the formation of the active site(s) of the spliceosome.