Genetic selection for balanced retroviral splicing: novel regulation involving the second step can be mediated by transitions in the polypyrimidine tract

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

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

Finding, Conducting, and Nurturing Science: A Virologist's Memoir

My laboratory investigations have been driven by an abiding interest in understanding the consequences of genetic rearrangement in evolution and disease, and in using viruses to elucidate fundamental mechanisms in biology. Starting with bacteriophages and moving to the retroviruses, my use of the tools of genetics, molecular biology, biochemistry, and biophysics has spanned more than half a century-from the time when DNA structure was just discovered to the present day of big data and epigenetics. Both riding and contributing to the successive waves of technology, my laboratory has elucidated fundamental mechanisms in DNA replication, repair, and recombination. We have made substantial contributions in the area of retroviral oncogenesis, delineated mechanisms that control retroviral gene expression, and elucidated critical details of the structure and function of the retroviral enzymes-reverse transcriptase, protease, and integrase-and have had the satisfaction of knowing that the fundamental knowledge gained from these studies contributed important groundwork for the eventual development of antiviral drugs to treat AIDS. While pursuing laboratory research as a principal investigator, I have also been a science administrator-moving from laboratory head to department chair and, finally, to institute director. In addition, I have undertaken a number of community service, science-related "extracurricular" activities during this time. Filling all of these roles, while being a wife and mother, has required family love and support, creative management, and, above all, personal flexibility-with not too much long-term planning. I hope that this description of my journey, with various roles, obstacles, and successes, will be both interesting and informative, especially to young female scientists.

Sequencing of lariat termini in S. cerevisiae reveals 5' splice sites, branch points, and novel splicing events

Pre-mRNA splicing is a central step in the shaping of the eukaryotic transcriptome and in the regulation of gene expression. Yet, due to a focus on fully processed mRNA, common approaches for defining pre-mRNA splicing genome-wide are suboptimal-especially with respect to defining the branch point sequence, a key cis-element that initiates the chemistry of splicing. Here, we report a complementary intron-centered approach designed to more efficiently, simply, and directly define splicing events genome-wide. Specifically, we developed a method distinguished by deep sequencing of lariat intron termini (LIT-seq). In a test of LIT-seq using the budding yeast Saccharomyces cerevisiae, we not only successfully captured the majority of annotated, expressed splicing events but also uncovered 45 novel splicing events, establishing the sensitivity of LIT-seq. Moreover, our libraries were highly enriched with reads that reported on splice sites; by a simple and direct inspection of sequencing reads, we empirically defined both 5' splice sites and branch sites, as well as their consensus sequences, with nucleotide resolution. Additionally, our study revealed that the 3' termini of lariat introns are subject to nontemplated addition of adenosines, characteristic of signals sensed by 3' to 5' RNA turnover machinery. Collectively, this work defines a novel, genome-wide approach for analyzing splicing with unprecedented depth, specificity, and resolution.

Neurofibromatosis type 1 alternative splicing is a key regulator of Ras signaling in neurons

Neurofibromatosis type I (Nf1) is a GTPase-activating protein (GAP) that inactivates the oncoprotein Ras and plays important roles in nervous system development and learning. Alternative exon 23a falls within the Nf1 GAP domain coding sequence and is tightly regulated in favor of skipping in neurons; however, its biological function is not fully understood. Here we generated mouse embryonic stem (ES) cells with a constitutive endogenous Nf1 exon 23a inclusion, termed Nf1 23aIN/23aIN cells, by mutating the splicing signals surrounding the exon to better match consensus sequences. We also made Nf1 23aΔ/23aΔ cells lacking the exon. Active Ras levels are high in wild-type (WT) and Nf1 23aIN/23aIN ES cells, where the Nf1 exon 23a inclusion level is high, and low in Nf1 23aΔ/23aΔ cells. Upon neuronal differentiation, active Ras levels are high in Nf1 23aIN/23aIN cells, where the exon inclusion level remains high, but Ras activation is low in the other two genotypes, where the exon is skipped. Signaling downstream of Ras is significantly elevated in Nf1 23aIN/23aIN neurons. These results suggest that exon 23a suppresses the Ras-GAP activity of Nf1. Therefore, regulation of Nf1 exon 23a inclusion serves as a mechanism for providing appropriate levels of Ras signaling and may be important in modulating Ras-related neuronal functions.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

RNA processing control in avian retroviruses

Upon integration into the host chromosome, retroviral gene expression requires transcription by the host RNA polymerase II, and viral messages are subject RNA processing events including 5'-end capping, pre-mRNA splicing, and polyadenylation. At a minimum, RNA splicing is required to generate the env mRNA, but viral replication requires substantial amounts of unspliced RNA to serve as mRNA and for incorporation into progeny virions as genomic RNA. Therefore, splicing has to be controlled to preserve the large unspliced RNA pool. Considering the current view that splicing and polyadenylation are coupled, the question arises as to how genome-length viral RNA is efficiently polyadenylated in the absence of splicing. Polyadenylation of many retroviral mRNAs is inefficient; in avian retroviruses, approximately 15 percent of viral transcripts extend into and are polyadenylated at downstream host genes, which often has profound biological consequences. Retroviruses have served as important models to study RNA processing and this review summarizes a body of work using avian retroviruses that has led to the discovery of novel RNA splicing and polyadenylation control mechanisms.

Large-scale comparative analysis of splicing signals and their corresponding splicing factors in eukaryotes

Introns are among the hallmarks of eukaryotic genes. Splicing of introns is directed by three main splicing signals: the 5' splice site (5'ss), the branch site (BS), and the polypyrimdine tract/3'splice site (PPT-3'ss). To study the evolution of these splicing signals, we have conducted a systematic comparative analysis of these signals in over 1.2 million introns from 22 eukaryotes. Our analyses suggest that all these signals have dramatically evolved: The PPT is weak among most fungi, intermediate in plants and protozoans, and strongest in metazoans. Within metazoans it shows a gradual strengthening from Caenorhabditis elegans to human. The 5'ss and the BS were found to be degenerate among most organisms, but highly conserved among some fungi. A maximum parsimony-based algorithm for reconstructing ancestral position-specific scoring matrices suggested that the ancestral 5'ss and BS were degenerate, as in metazoans. To shed light on the evolutionary variation in splicing signals, we have analyzed the evolutionary changes in the factors that bind these signals. Our analysis reveals coevolution of splicing signals and their corresponding splicing factors: The strength of the PPT is correlated to changes in key residues in its corresponding splicing factor U2AF2; limited correlation was found between changes in the 5'ss and U1 snRNA that binds it; but not between the BS and U2 snRNA. Thus, although the basic ability of eukaryotes to splice introns has remained conserved throughout evolution, the splicing signals and their corresponding splicing factors have considerably evolved, uniquely shaping the splicing mechanisms of different organisms.

Adaptive evolution of a tagged chimeric gammaretrovirus: identification of novel cis-acting elements that modulate splicing

Retroviruses are well known for their ability to incorporate envelope (Env) proteins from other retroviral strains and genera, and even from other virus families. This characteristic has been widely exploited for the generation of replication-defective retroviral vectors, including those derived from murine leukemia virus (MLV), bearing heterologous Env proteins. We investigated the possibility of "genetically pseudotyping" replication-competent MLV by replacing the native env gene in a full-length viral genome with that of another gammaretrovirus. Earlier, we developed replication-competent versions of MLV that stably transmit and express transgenes inserted into the 3' untranslated region of the viral genome. In one such tagged MLV expressing green fluorescent protein, we replaced the native env sequence with that of gibbon ape leukemia virus (GALV). Although the GALV Env protein is commonly used to make high-titer pseudotypes of MLV vectors, we found that the env replacement greatly attenuated viral replication. However, extended cultivation of cells exposed to the chimeric virus resulted in selection of mutants exhibiting rapid replication kinetics and different variants arose in different infections. Two of these variants had acquired mutations at or adjacent to the splice acceptor site, and three others had acquired dual mutations within the long terminal repeat. Analysis of the levels of unspliced and spliced viral RNA produced by the parental and adapted viruses showed that the mutations gained by each of these variants functioned to reverse an imbalance in splicing caused by the env gene substitution. Our results reveal the presence of previously unknown cis-acting sequences in MLV that modulate splicing of the viral transcript and demonstrate that tagging of the retroviral genome with an easily assayed transgene can be combined with in vitro evolution as an approach to efficiently generating and screening for replicating mutants of replication-impaired recombinant viruses.

Murine leukemia virus regulates alternative splicing through sequences upstream of the 5' splice site

Alternative splicing of the primary transcript plays a key role in retroviral gene expression. In contrast to all known mechanisms that mediate alternative splicing in retroviruses, we found that in murine leukemia virus, distinct elements located upstream of the 5' splice site either inhibited or activated splicing of the genomic RNA. Detailed analysis of the first untranslated exon showed that the primer binding site (PBS) activates splicing, whereas flanking sequences either downstream or upstream of the PBS are inhibitory. This new function of the PBS was independent of its orientation and primer binding but associated with a particular destabilizing role in a proposed secondary structure. On the contrary, all sequences surrounding the PBS that are involved in stem formation of the first exon were found to suppress splicing. Targeted mutations that destabilized the central stem and compensatory mutations of the counter strand clearly validated the concept that murine leukemia virus attenuates its 5' splice site by forming an inhibitory stem-loop in its first exon. Importantly, this mode of splice regulation was conserved in a complete proviral clone. Some of the mutants that increase splicing revealed an opposite effect on translation, implying that the first exon also regulates this process. Together, these findings suggest that sequences upstream of the 5' splice site play an important role in splice regulation of simple retroviruses, directly or indirectly attenuating the efficiency of splicing.

Characterization and prediction of alternative splice sites

Human alternative isoform, cryptic, skipped, and constitutive splice sites from the ALTEXTRON database were analysed regarding splice site strength, composition, GC content, position and binding site strength of polypyrimidine tract and branch site. Several features were identified which distinguish alternative isoform and cryptic splice sites, but not skipped splice sites from constitutive ones. These include splice site strength, introns GC content, U2AF35 binding site score, and oligonucleotide frequencies. For the predictive classification of splice sites, pattern recognition models for different splicing factor binding sites and oligonucleotide frequency models (OFMs) were combined using backpropagation networks. 67.45% of acceptor sites and 71.23% of donor sites are correctly classified by networks trained for classification of constitutive and alternative isoform/cryptic splice sites. A web-application for the prediction of alternative splice sites is available at http://es.embnet.org/~mwang/assp.html .

Biased exon/intron distribution of cryptic and de novo 3' splice sites

We compiled sequences of previously published aberrant 3′ splice sites (3′ss) that were generated by mutations in human disease genes. Cryptic 3′ss, defined here as those resulting from a mutation of the 3′YAG consensus, were more frequent in exons than in introns. They clustered in ∼20 nt region adjacent to authentic 3′ss, suggesting that their under-representation in introns is due to a depletion of AG dinucleotides in the polypyrimidine tract (PPT). In contrast, most aberrant 3′ss that were induced by mutations outside the 3′YAG consensus (designated ‘de novo’) were in introns. The activation of intronic de novo 3′ss was largely due to AG-creating mutations in the PPT. In contrast, exonic de novo 3′ss were more often induced by mutations improving the PPT, branchpoint sequence (BPS) or distant auxiliary signals, rather than by direct AG creation. The Shapiro–Senapathy matrix scores had a good prognostic value for cryptic, but not de novo 3′ss. Finally, AG-creating mutations in the PPT that produced aberrant 3′ss upstream of the predicted BPS in vivo shared a similar ‘BPS-new AG’ distance. Reduction of this distance and/or the strength of the new AG PPT in splicing reporter pre-mRNAs improved utilization of authentic 3′ss, suggesting that AG-creating mutations that are located closer to the BPS and are preceded by weaker PPT may result in less severe splicing defects.

The exon splicing silencer in human immunodeficiency virus type 1 Tat exon 3 is bipartite and acts early in spliceosome assembly

Inefficient splicing of human immunodeficiency virus type 1 (HIV-1) RNA is necessary to preserve unspliced and singly spliced viral RNAs for transport to the cytoplasm by the Rev-dependent pathway. Signals within the HIV-1 genome that control the rate of splicing include weak 3' splice sites, exon splicing enhancers (ESE), and exon splicing silencers (ESS). We have previously shown that an ESS present within tat exon 2 (ESS2) and a suboptimal 3' splice site together act to inhibit splicing at the 3' splice site flanking tat exon 2. This occurs at an early step in spliceosome assembly. Splicing at the 3' splice site flanking tat exon 3 is regulated by a bipartite element composed of an ESE and an ESS (ESS3). Here we show that ESS3 is composed of two smaller elements (AGAUCC and UUAG) that can inhibit splicing independently. We also show that ESS3 is more active in the context of a heterologous suboptimal splice site than of an optimal 3' splice site. ESS3 inhibits splicing by blocking the formation of a functional spliceosome at an early step, since A complexes are not detected in the presence of ESS3. Competitor RNAs containing either ESS2 or ESS3 relieve inhibition of splicing of substrates containing ESS3 or ESS2. This suggests that a common cellular factor(s) may be required for the inhibition of tat mRNA splicing mediated by ESS2 and ESS3.

Role of the constitutive splicing factors U2AF65 and SAP49 in suboptimal RNA splicing of novel retroviral mutants

Retroviruses display a unique form of alternative splicing in which both spliced and unspliced RNAs accumulate in the cytoplasm. Simple retroviruses, such as avian sarcoma virus, do not encode regulatory proteins that affect splicing; this process is controlled solely through interactions between the viral RNA and the host cell splicing machinery. Previously, we described the selection and characterization of novel avian sarcoma virus mutants. These viruses were separated into two classes based upon analysis of splicing intermediates produced in infected cells and in a cell-free system. One class, which included mutants with altered polypyrimidine tract or branch point sequences, showed significant accumulation of intermediates, suggesting that splicing was regulated in step 2. The other class, which included mutants with deletions of exonic enhancer sequences, did not accumulate splicing intermediates, suggesting that splicing was regulated before step 1 of the splicing reaction. In this report, we show that a mutant blocked at step 1 fails to form a stable spliceosomal complex, whereas one blocked at step 2 shows a defect in its ability to transit through the last spliceosomal complex. Using UV cross-linking methods, we show that regulation at each step is associated with specific changes in the binding of cellular splicing factors. Regulation at step 1 is correlated with decreased cross-linking of the factor U2AF65, whereas regulation at step 2 is correlated with enhanced cross-linking of the factor SAP49. Because these mutations were isolated by selection for replication-competent viruses, we conclude that retroviral splicing may be regulated in vivo through altered binding of constitutive splicing factors.

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.

The Moloney murine sarcoma virus ts110 5' splice site signal contributes to the regulation of splicing efficiency and thermosensitivity

The 5' splice site signal (5'ss) in Moloney murine sarcoma virus ts110 (MuSVts110) RNA was found to participate in the regulation of its splicing phenotype. This 5'ss (CAG/GUAGGA) departs from the mammalian consensus (CAG/GURAGU) at positions +4 and +6, both of which base pair with U1 and U6 small nuclear RNAs during splicing. A doubling in splicing efficiency and near elimination of the splicing thermosensitivity characteristic of MuSVts110 were observed in 5'ss mutants containing a U at position +6 (termed 5' A6U), even in those in which U1-5'ss complementarity had been reduced. At the permissive temperature (28 degrees C), the 5' A6U mutation increased the efficiency of the second splicing reaction, while at the nonpermissive temperature (39 degrees C), both splicing reactions were positively affected.

Selection of the bovine papillomavirus type 1 nucleotide 3225 3' splice site is regulated through an exonic splicing enhancer and its juxtaposed exonic splicing suppressor

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. However, one 3' splice site, located at nucleotide (nt) 3225, is used for the processing of most BPV-1 pre-mRNAs in BPV-1-transformed C127 cells and at early to intermediate times in productively infected warts. At late stages of the viral life cycle, an alternative 3' splice site at nt 3605 is used for the processing of the late pre-mRNA. In this study, we used in vitro splicing in HeLa cell nuclear extracts to identify cis elements which regulate BPV-1 3' splice site selection. Two purine-rich exonic splicing enhancers were identified downstream of nt 3225. These sequences, designated SE1 (nt 3256 to 3305) and SE2 (nt 3477 to 3526), were shown to strongly stimulate the splicing of a chimeric Drosophila doublesex pre-mRNA, which contains a weak 3' splice site. A BPV-1 late pre-mRNA containing the nt 3225 3' splice site but lacking both SE1 and SE2 was spliced poorly, indicating that this 3' splice site is inherently weak. Analysis of the 3' splice site suggested that this feature is due to both a nonconsensus branch point sequence and a suboptimal polypyrimidine tract. Addition of SE1 to the late pre-mRNA dramatically stimulated splicing, indicating that SE1 also functions as an exonic splicing enhancer in its normal context. However, a late pre-mRNA containing both SE1 and SE2 as well as the sequence in between was spliced inefficiently. Further mapping studies demonstrated that a 48-nt pyrimidine-rich region immediately downstream of SE1 was responsible for this suppression of splicing. Thus, these data suggest that selection of the BPV-1 nt 3225 3' splice site is regulated by both positive and negative exonic sequences.

Sequence of the polypyrimidine tract of the 3'-terminal 3' splicing signal can affect intron-dependent pre-mRNA processing in vivo

Most pre-mRNAs require an intron for efficient processing in higher eukaryotes. However, not all introns can provide this function. For example, transcripts synthesized from a variant of the human beta-globin gene lacking its second intervening sequence (IVS2), yet retaining its first intervening sequence (IVS1), exhibit multiple defects in mRNA biogenesis. To investigate why, we transfected into monkey cells plasmids containing the human beta-globin gene and variants of it altered in (i) IVS1, (ii) the 3'-terminal exon, and (iii) the polyadenylation signal. The beta-globin RNAs accumulated in these cells were analyzed by quantitative S1 nuclease mapping for nuclear accumulation, intron excision, polyadenylation and cytoplasmic accumulation. We found that the 3' splicing signal of IVS1, with multiple purines interrupting its polypyrimidine tract, could efficiently function as an internal 3' splicing signal; however, it could not efficiently function as the 3'-terminal 3' splicing signal for any of these steps in intron-dependent mRNA biogenesis unless (i) its polypyrimidine tract was made uninterrupted in pyrimidines, or (ii) specific sequences were deleted from the 3'-terminal exon. We conclude that whether an intron can provide the function necessary for efficient processing of intron-dependent pre-mRNA is dependent upon the ability of its 3' splicing signal to define the 3'-terminal exon. On the practical side, this finding means one needs to consider both the sequence and location of the intron to be included in an intron-dependent gene to obtain efficient expression in vivo.

Branchpoint and polypyrimidine tract mutations mediating the loss and partial recovery of the Moloney murine sarcoma virus MuSVts110 thermosensitive splicing phenotype

Balanced splicing of retroviral RNAs is mediated by weak signals at the 3' splice site (ss) acting in concert with other cis elements. Moloney murine sarcoma virus MuSVts110 shows a similar balance between unspliced and spliced RNAs, differing only in that the splicing of its RNA is, in addition, growth temperature sensitive. We have generated N-nitroso-N-methylurea (NMU)-treated MuSVts110 revertants in which splicing was virtually complete at all temperatures and have investigated the molecular basis of this reversion on the assumption that the findings would reveal cis-acting elements controlling MuSVts110 splicing thermosensitivity. In a representative revertant (NMU-20), we found that complete splicing was conferred by a G-to-A substitution generating a consensus branchpoint (BP) signal (-CCCUGGC- to -CCCUGAC- [termed G(-25)A]) at -25 relative to the 3' ss. Weakening this BP to -CCCGAC- [G(-25)A,U(-27)C] moderately reduced splicing at the permissive temperature and sharply inhibited splicing at the originally nonpermissive temperature, arguing that MuSVts110 splicing thermosensitivity depends on a suboptimal BP-U2 small nuclear RNA interaction. This conclusion was supported by results indicating that lengthening the short MuSVts110 polypyrimidine tract and altering its uridine content doubled splicing efficiency at permissive temperatures and nearly abrogated splicing thermosensitivity. In vitro splicing experiments showed that MuSVts110 G(-25)A RNA intermediates were far more efficiently ligated than RNAs carrying the wild-type BP, the G(-25)A,U (-27)C BP, or the extended polypyrimidine tract. The efficiency of ligation in vitro roughly paralleled splicing efficiency in vivo [G(-25)A BP > extended polypyrimidine tract > G(-25)A,U(-27)C BP > wild-type BP]. These results suggest that MuSVts110 RNA splicing is balanced by cis elements similar to those operating in other retroviruses and, in addition, that its splicing thermosensitivity is a response to the presence of multiple suboptimal splicing signals.