Conserved stem-loop structures in the HIV-1 RNA region containing the A3 3' splice site and its cis-regulatory element: possible involvement in RNA splicing

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

Cooperative Binding of SRSF3 to Structured 3'ss-α Exon RNA during α Exon Inclusion in the ZO-1 mRNA

ZO-1α+ and ZO-1α- proteins are expressed in hermetic and leaky tight junctions, respectively. Two cis-acting distant exonic elements partly activate the 240 nucleotide-long α exon producing the ZO-1α+ isoform. However, the elements within and around the α exon and their respective factors involved in its splicing are unknown. To study the dynamic interaction between SRSF3 and its bioinformatically predicted target sites around the 3'ss upstream of the α exon during its activation, we performed EMSA, crosslinking, and in vivo splicing assays by ZO-1 minigene expression and siRNA-mediated silencing in transfected cells. Using V1 RNase, we probed the possible formation of a hairpin RNA structure between the intronic and proximal exonic SRSF3 binding sites. The hairpin sufficed for complex formations in the EMSA. The interaction of SRSF3 with the intronic site promoted the cooperative binding of SRSF3 to the exonic site. Finally, SRSF3 restored α exon activation in SRSF3 knockdown transfectants. Altogether, our results show that SRSF3-hairpin RNA interaction is crucial in the early recognition of 3'ss for α exon activation. It remains to be explored whether SRSF3 recruits or stabilizes the binding of other factors or brings separate splice sites into proximity.

Structural and Functional RNA Motifs of SARS-CoV-2 and Influenza A Virus as a Target of Viral Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic, whereas the influenza A virus (IAV) causes seasonal epidemics and occasional pandemics. Both viruses lead to widespread infection and death. SARS-CoV-2 and the influenza virus are RNA viruses. The SARS-CoV-2 genome is an approximately 30 kb, positive sense, 5' capped single-stranded RNA molecule. The influenza A virus genome possesses eight single-stranded negative-sense segments. The RNA secondary structure in the untranslated and coding regions is crucial in the viral replication cycle. The secondary structure within the RNA of SARS-CoV-2 and the influenza virus has been intensively studied. Because the whole of the SARS-CoV-2 and influenza virus replication cycles are dependent on RNA with no DNA intermediate, the RNA is a natural and promising target for the development of inhibitors. There are a lot of RNA-targeting strategies for regulating pathogenic RNA, such as small interfering RNA for RNA interference, antisense oligonucleotides, catalytic nucleic acids, and small molecules. In this review, we summarized the knowledge about the inhibition of SARS-CoV-2 and influenza A virus propagation by targeting their RNA secondary structure.

Alternative RNA Conformations: Companion or Combatant

RNA molecules, in one form or another, are involved in almost all aspects of cell physiology, as well as in disease development. The diversity of the functional roles of RNA comes from its intrinsic ability to adopt complex secondary and tertiary structures, rivaling the diversity of proteins. The RNA molecules form dynamic ensembles of many interconverting conformations at a timescale of seconds, which is a key for understanding how they execute their cellular functions. Given the crucial role of RNAs in various cellular processes, we need to understand the RNA molecules from a structural perspective. Central to this review are studies aimed at revealing the regulatory role of conformational equilibria in RNA in humans to understand genetic diseases such as cancer and neurodegenerative diseases, as well as in pathogens such as bacteria and viruses so as to understand the progression of infectious diseases. Furthermore, we also summarize the prior studies on the use of RNA structures as platforms for the rational design of small molecules for therapeutic applications.

Encoded Conformational Dynamics of the HIV Splice Site A3 Regulatory Locus: Implications for Differential Binding of hnRNP Splicing Auxiliary Factors

Alternative splicing of the HIV transcriptome is controlled through cis regulatory elements functioning as enhancers or silencers depending on their context and the type of host RNA binding proteins they recruit. Splice site acceptor A3 (ssA3) is one of the least used acceptor sites in the HIV transcriptome and its activity determines the levels of tat mRNA. Splice acceptor 3 is regulated by a combination of cis regulatory sequences, auxiliary splicing factors, and presumably RNA structure. The mechanisms by which these multiple regulatory components coordinate to determine the frequency in which ssA3 is utilized is poorly understood. By NMR spectroscopy and phylogenetic analysis, we show that the ssA3 regulatory locus is conformationally heterogeneous and that the sequences that encompass the locus are conserved across most HIV isolates. Despite the conformational heterogeneity, the major stem loop (A3SL1) observed in vitro folds to base pair the Polypyrimdine Tract (PPyT) to the Exon Splicing Silencer 2p (ESS2p) element and to a conserved downstream linker. The 3D structure as determined by NMR spectroscopy further reveals that the A3 consensus cleavage site is embedded within a unique stereochemical environment within the apical loop, where it is surrounded by alternating base-base interactions. Despite being described as a receptor for hnRNP H, the ESS2p element is sequestered by base pairing to the 3' end of the PPyT and within this context it cannot form a stable complex with hnRNP H. By comparison, hnRNP A1 directly binds to the A3 consensus cleavage site located within the apical loop, suggesting that it can directly modulate U2AF assembly. Sequence mutations designed to destabilize the PPyT:ESS2p helix results in an increase usage of ssA3 within HIV-infected cells, consistent with the PPyT becoming more accessible for U2AF recognition. Additional mutations introduced into the downstream ESS2 element synergize with ESS2p to cause further increases in ssA3 usage. When taken together, our work provides a unifying picture by which cis regulatory sequences, splicing auxiliary factors and RNA structure cooperate to provide stringent control over ssA3. We describe this as the pair-and-lock mechanism to restrict access of the PPyT, and posit that it operates to regulate a subset of the heterogenous structures encompassing the ssA3 regulatory locus.

Plasma-Derived HIV-1 Virions Contain Considerable Levels of Defective Genomes

Genetically-characterizing full-length HIV-1 RNA is critical for identifying genetically-intact genomes and for comparing these RNA genomes to proviral DNA. We have developed a method for sequencing plasma-derived RNA using long-range sequencing (PRLS assay; ∼8.3 kb from gag to the 3′ end or ∼5 kb from integrase to the 3′ end). We employed the gag-3′ PRLS assay to sequence HIV-1 RNA genomes from ART-naive participants during acute/early infection (n = 6) or chronic infection (n = 2). On average, only 65% of plasma-derived genomes were genetically-intact. Defects were found in all genomic regions but were concentrated in env and pol. We compared these genomes to near-full-length proviral sequences from paired peripheral blood mononuclear cell (PBMC) samples for the acute/early group and found that near-identical (>99.98% identical) sequences were identified only during acute infection. For three participants who initiated therapy during acute infection, we used the int-3′ PRLS assay to sequence plasma-derived genomes from an analytical treatment interruption and identified 100% identical genomes between pretherapy and rebound time points. The PRLS assay provides a new level of sensitivity for understanding the genetic composition of plasma-derived HIV-1 RNA from viremic individuals either pretherapy or after treatment interruption, which will be invaluable in assessing possible HIV-1 curative strategies.

IMPORTANCE We developed novel plasma-derived RNA using long-range sequencing assays (PRLS assay; 8.3 kb, gag-3′, and 5.0 kb, int-3′). Employing the gag-3′ PRLS assay, we found that 26% to 51% of plasma-derived genomes are genetically-defective, largely as a result of frameshift mutations and deletions. These genetic defects were concentrated in the env region compared to gag and pol, likely a reflection of viral immune escape in env during untreated HIV-1 infection. Employing the int-3′ PRLS assay, we found that analytical treatment interruption (ATI) plasma-derived sequences were identical and genetically-intact. Several sequences from the ATI plasma samples were identical to viral sequences from pretherapy plasma and PBMC samples, indicating that HIV-1 reservoirs established prior to therapy contribute to viral rebound during an ATI. Therefore, near-full-length sequencing of HIV-1 particles is required to gain an accurate picture of the genetic landscape of plasma HIV-1 virions in studies of HIV-1 replication and persistence.

Conserved long-range base pairings are associated with pre-mRNA processing of human genes

The ability of nucleic acids to form double-stranded structures is essential for all living systems on Earth. Current knowledge on functional RNA structures is focused on locally-occurring base pairs. However, crosslinking and proximity ligation experiments demonstrated that long-range RNA structures are highly abundant. Here, we present the most complete to-date catalog of conserved complementary regions (PCCRs) in human protein-coding genes. PCCRs tend to occur within introns, suppress intervening exons, and obstruct cryptic and inactive splice sites. Double-stranded structure of PCCRs is supported by decreased icSHAPE nucleotide accessibility, high abundance of RNA editing sites, and frequent occurrence of forked eCLIP peaks. Introns with PCCRs show a distinct splicing pattern in response to RNAPII slowdown suggesting that splicing is widely affected by co-transcriptional RNA folding. The enrichment of 3'-ends within PCCRs raises the intriguing hypothesis that coupling between RNA folding and splicing could mediate co-transcriptional suppression of premature pre-mRNA cleavage and polyadenylation.

Dynamic nanopore long-read sequencing analysis of HIV-1 splicing events during the early steps of infection

Background

Alternative splicing is a key step in Human Immunodeficiency Virus type 1 (HIV-1) replication that is tightly regulated both temporally and spatially. More than 50 different transcripts can be generated from a single HIV-1 unspliced pre-messenger RNA (pre-mRNA) and a balanced proportion of unspliced and spliced transcripts is critical for the production of infectious virions. Understanding the mechanisms involved in the regulation of viral RNA is therefore of potential therapeutic interest. However, monitoring the regulation of alternative splicing events at a transcriptome-wide level during cell infection is challenging. Here we used the long-read cDNA sequencing developed by Oxford Nanopore Technologies (ONT) to explore in a quantitative manner the complexity of the HIV-1 transcriptome regulation in infected primary CD4+ T cells.

Results

ONT reads mapping to the viral genome proved sufficiently long to span all possible splice junctions, even distant ones, and to be assigned to a total of 150 exon combinations. Fifty-three viral RNA isoforms, including 14 new ones were further considered for quantification. Relative levels of viral RNAs determined by ONT sequencing showed a high degree of reproducibility, compared favourably to those produced in previous reports and highly correlated with quantitative PCR (qPCR) data. To get further insights into alternative splicing regulation, we then compiled quantifications of splice site (SS) usage and transcript levels to build “splice trees”, a quantitative representation of the cascade of events leading to the different viral isoforms. This approach allowed visualizing the complete rewiring of SS usages upon perturbation of SS D2 and its impact on viral isoform levels. Furthermore, we produced the first dynamic picture of the cascade of events occurring between 12 and 24 h of viral infection. In particular, our data highlighted the importance of non-coding exons in viral RNA transcriptome regulation.

Conclusion

ONT sequencing is a convenient and reliable strategy that enabled us to grasp the dynamic of the early splicing events modulating the viral RNA landscape in HIV-1 infected cells.

Tat IRES modulator of tat mRNA (TIM-TAM): a conserved RNA structure that controls Tat expression and acts as a switch for HIV productive and latent infection

Tat protein is essential to fully activate HIV transcription and processing of viral mRNA, and therefore determines virus expression in productive replication and the establishment and maintenance of latent infection. Here, we used thermodynamic and structure analyses to define a highly conserved sequence-structure in tat mRNA that functions as Tat IRES modulator of tat mRNA (TIM-TAM). By impeding cap-dependent ribosome progression during authentic spliced tat mRNA translation, TIM-TAM stable structure impacts on timing and level of Tat protein hence controlling HIV production and infectivity along with promoting latency. TIM-TAM also adopts a conformation that mediates Tat internal ribosome entry site (IRES)-dependent translation during the early phases of infection before provirus integration. Our results document the critical role of TIM-TAM in Tat expression to facilitate virus reactivation from latency, with implications for HIV treatment and drug development.

A scale-free analysis of the HIV-1 genome demonstrates multiple conserved regions of structural and functional importance

HIV-1 replicates via a low-fidelity polymerase with a high mutation rate; strong conservation of individual nucleotides is highly indicative of the presence of critical structural or functional properties. Identifying such conservation can reveal novel insights into viral behaviour. We analysed 3651 publicly available sequences for the presence of nucleic acid conservation beyond that required by amino acid constraints, using a novel scale-free method that identifies regions of outlying score together with a codon scoring algorithm. Sequences with outlying score were further analysed using an algorithm for producing local RNA folds whilst accounting for alignment properties. 11 different conserved regions were identified, some corresponding to well-known cis-acting functions of the HIV-1 genome but also others whose conservation has not previously been noted. We identify rational causes for many of these, including cis functions, possible additional reading frame usage, a plausible mechanism by which the central polypurine tract primes second-strand DNA synthesis and a conformational stabilising function of a region at the 5' end of env.

Production of HIV-1 vif mRNA Is Modulated by Natural Nucleotide Variations and SLSA1 RNA Structure in SA1D2prox Genomic Region

Genomic RNA of HIV-1 contains localized structures critical for viral replication. Its structural analysis has demonstrated a stem-loop structure, SLSA1, in a nearby region of HIV-1 genomic splicing acceptor 1 (SA1). We have previously shown that the expression level of vif mRNA is considerably altered by some natural single-nucleotide variations (nSNVs) clustering in SLSA1 structure. In this study, besides eleven nSNVs previously identified by us, we totally found nine new nSNVs in the SLSA1-containing sequence from SA1, splicing donor 2, and through to the start codon of Vif that significantly affect the vif mRNA level, and designated the sequence SA1D2prox (142 nucleotides for HIV-1 NL4-3). We then examined by extensive variant and mutagenesis analyses how SA1D2prox sequence and SLSA1 secondary structure are related to vif mRNA level. While the secondary structure and stability of SLSA1 was largely changed by nSNVs and artificial mutations introduced to restore the original NL4-3 form from altered ones by nSNVs, no clear association of the two SLSA1 properties with vif mRNA level was observed. In contrast, when naturally occurring SA1D2prox sequences that contain multiple nSNVs were examined, we attained significant inverse correlation between the vif level and SLSA1 stability. These results may suggest that SA1D2prox sequence adapts over time, and also that the altered SA1D2prox sequence, SLSA1 stability, and vif level are mutually related. In total, we show here that the entire SA1D2prox sequence and SLSA1 stability critically contribute to the modulation of vif mRNA level.

HIV-1 Vpr N-terminal tagging affects alternative splicing of the viral genome

To facilitate studies on Vpr function in replicating HIV-1, we aimed to tag the protein in an infectious virus. First we showed that N-, but not C-terminal HA/FLAG tagging of Vpr protein preserves Vpr cytopathicity. Cloning the tags into proviral DNA however ablated viral production and replication. By construction of additional viral variants we could show this defect was not protein- but RNA-dependent and sequence specific, and characterized by oversplicing of the genomic RNA. Simulation of genomic RNA folding suggested that introduction of the tag sequence induced an alternative folding structure in a region enriched in splice sites and splicing regulatory sequences. In silico predictions identified the HA/His₆-Vpr tagging in HIV-1 to affect mRNA folding less than HA/FLAG-Vpr tagging. In vitro infectivity and mRNA splice pattern improved but did not reach wild-type values. Thus, sequence-specific insertions may interfere with mRNA splicing, possibly due to altered RNA folding. Our results point to the complexity of viral RNA genome sequence interactions. This should be taken into consideration when designing viral manipulation strategies, for both research as for biological interventions.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

Sequence Analysis of In Vivo-Expressed HIV-1 Spliced RNAs Reveals the Usage of New and Unusual Splice Sites by Viruses of Different Subtypes

HIV-1 RNAs are generated through a complex splicing mechanism, resulting in a great diversity of transcripts, which are classified in three major categories: unspliced, singly spliced (SS), and doubly spliced (DS). Knowledge on HIV-1 RNA splicing in vivo and by non-subtype B viruses is scarce. Here we analyze HIV-1 RNA splice site usage in CD4+CD25+ lymphocytes from HIV-1-infected individuals through pyrosequencing. HIV-1 DS and SS RNAs were amplified by RT-PCR in 19 and 12 samples, respectively. 13,108 sequences from HIV-1 spliced RNAs, derived from viruses of five subtypes (A, B, C, F, G), were identified. In four samples, three of non-B subtypes, five 3' splice sites (3'ss) mapping to unreported positions in the HIV-1 genome were identified. Two, designated A4i and A4j, were used in 22% and 25% of rev RNAs in two viruses of subtypes B and A, respectively. Given their close proximity (one or two nucleotides) to A4c and A4d, respectively, they could be viewed as variants of these sites. Three 3'ss, designated A7g, A7h, and A7i, located 20, 32, and 18 nucleotides downstream of A7, respectively, were identified in a subtype C (A7g, A7h) and a subtype G (A7i) viruses, each in around 2% of nef RNAs. The new splice sites or variants of splice sites were associated with the usual sequence features of 3'ss. Usage of unusual 3'ss A4d, A4e, A5a, A7a, and A7b was also detected. A4f, previously identified in two subtype C viruses, was preferentially used by rev RNAs of a subtype C virus. These results highlight the great diversity of in vivo splice site usage by HIV-1 RNAs. The fact that four of five newly identified splice sites or variants of splice sites were detected in non-subtype B viruses allows anticipating an even greater diversity of HIV-1 splice site usage than currently known.

Natural Single-Nucleotide Variations in the HIV-1 Genomic SA1prox Region Can Alter Viral Replication Ability by Regulating Vif Expression Levels

We previously found that natural single-nucleotide variations located within a proximal region of splicing acceptor 1 (SA1prox) in the HIV-1 genome could alter the viral replication potential and mRNA expression pattern, especially the vif mRNA level. Here, we studied the virological and molecular basis of nucleotide sequence variations in SA1prox for alterations of viral replication ability. Consistent with our previous findings, variant clones indeed expressed Vif at different levels and grew distinctively in cells with various APOBEC3G expression levels. Similar effects were observed for natural variations found in HIV-2 SA1prox, suggesting the importance of the SA1prox sequence. To define nucleotides critical for the regulation of HIV-1 Vif expression, effects of natural SA1prox variations newly found in the HIV Sequence Compendium database on vif mRNA/Vif protein levels were examined. Seven out of nine variations were found to produce Vif at lower, higher, or more excessive levels than wild-type NL4-3. Combination experiments of variations giving distinct Vif levels suggested that the variations mutually affected vif transcript production. While low and high producers of Vif grew in an APOBEC3G-dependent manner, excessive expressers always showed an impeded growth phenotype due to defects in single-cycle infectivity and/or virion production levels. The phenotype of excessive expressers was not due primarily to inadequate expression of Tat or Rev, although SA1prox variations altered the overall HIV-1 mRNA expression pattern. Collectively, our results demonstrate that HIV SA1prox regulates Vif expression levels and suggest a relationship between SA1prox and viral adaptation/evolution given that variations occurred naturally.

IMPORTANCE

While human cells possess restriction factors to inhibit HIV-1 replication, HIV-1 encodes antagonists to overcome these barriers. Conflicts between host restriction factors and viral counterparts are critical driving forces behind mutual evolution. The interplay of cellular APOBEC3G and viral Vif proteins is a typical example. Here, we demonstrate that naturally occurring single-nucleotide variations in the proximal region of splicing acceptor 1 (SA1prox) of the HIV-1 genome frequently alter Vif expression levels, thereby modulating viral replication potential in cells with various ABOBEC3G levels. The results of the present study reveal a previously unidentified and important way for HIV-1 to compete with APOBEC3G restriction by regulating its Vif expression levels. We propose that SA1prox plays a regulatory role in Vif counteraction against APOBEC3G in order to contribute to HIV-1 replication and evolution, and this may be applicable to other primate lentiviruses.

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.

Identification of U2AF(35)-dependent exons by RNA-Seq reveals a link between 3' splice-site organization and activity of U2AF-related proteins

The auxiliary factor of U2 small nuclear RNA (U2AF) is a heterodimer consisting of 65- and 35-kD proteins that bind the polypyrimidine tract (PPT) and AG dinucleotides at the 3′ splice site (3′ss). The gene encoding U2AF35 (U2AF1) is alternatively spliced, giving rise to two isoforms U2AF35a and U2AF35b. Here, we knocked down U2AF35 and each isoform and characterized transcriptomes of HEK293 cells with varying U2AF35/U2AF65 and U2AF35a/b ratios. Depletion of both isoforms preferentially modified alternative RNA processing events without widespread failure to recognize 3′ss or constitutive exons. Over a third of differentially used exons were terminal, resulting largely from the use of known alternative polyadenylation (APA) sites. Intronic APA sites activated in depleted cultures were mostly proximal whereas tandem 3′UTR APA was biased toward distal sites. Exons upregulated in depleted cells were preceded by longer AG exclusion zones and PPTs than downregulated or control exons and were largely activated by PUF60 and repressed by CAPERα. The U2AF(35) repression and activation was associated with a significant interchange in the average probabilities to form single-stranded RNA in the optimal PPT and branch site locations and sequences further upstream. Although most differentially used exons were responsive to both U2AF subunits and their inclusion correlated with U2AF levels, a small number of transcripts exhibited distinct responses to U2AF35a and U2AF35b, supporting the existence of isoform-specific interactions. These results provide new insights into function of U2AF and U2AF35 in alternative RNA processing.

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.

Strong epistatic selection on the RNA secondary structure of HIV

A key question in evolutionary genomics is how populations navigate the adaptive landscape in the presence of epistasis, or interactions among loci. This problem can be directly addressed by studying the evolution of RNA secondary structures, for which there is constraint to maintain pairing between Watson-Crick (WC) sites. Replacement of a nucleotide at one site of a WC pair reduces fitness by disrupting binding, which can be restored via a compensatory replacement at the interacting site. Here, I present the first genome-scale analysis of epistasis on the RNA secondary structure of human immunodeficiency virus type 1 (HIV-1). Comparison of polymorphism frequencies at ancestrally conserved sites reveals that selection against replacements is ∼2.7 times stronger at WC than at non-WC sites, such that nearly 50% of constraint can be attributed to epistasis. However, almost all epistatic constraint is due to selection against conversions of WC pairs to unpaired (UP) nucleotides, whereas conversions to GU wobbles are only slightly deleterious. This disparity is also evident in pairs with second-site compensatory replacements; conversions from UP nucleotides to WC pairs increase median fitness by ∼4.2%, whereas conversions from GU wobbles to WC pairs only increase median fitness by ∼0.3%. Moreover, second-site replacements that convert UP nucleotides to GU wobbles also increase median fitness by ∼4%, indicating that such replacements are nearly as compensatory as those that restore WC pairing. Thus, WC peaks of the HIV-1 epistatic adaptive landscape are connected by high GU ridges, enabling the viral population to rapidly explore distant peaks without traversing deep UP valleys.

Epistasis is an evolutionary process in which the effect of a nucleotide at one site in the genome is dependent on the presence or absence of particular nucleotides at other sites in the genome. One of the simplest types of epistasis occurs between Watson-Crick (WC) nucleotides in RNA secondary structures, which are under constraint to maintain base-pairing. In this study, I examine the effects of mutations at WC sites in the RNA secondary structure of HIV-1. I show that while epistasis plays a major role in the evolution of the HIV-1 secondary structure, different types of mutations have variable effects on fitness. Therefore, by favoring certain mutational trajectories, HIV-1 can evolve rapidly despite strong epistatic constraint on its RNA secondary structure.

Impact of human pathogenic micro-insertions and micro-deletions on post-transcriptional regulation

Small insertions/deletions (INDELs) of ≤21 bp comprise 18% of all recorded mutations causing human inherited disease and are evident in 24% of documented Mendelian diseases. INDELs affect gene function in multiple ways: for example, by introducing premature stop codons that either lead to the production of truncated proteins or affect transcriptional efficiency. However, the means by which they impact post-transcriptional regulation, including alternative splicing, have not been fully evaluated. In this study, we collate disease-causing INDELs from the Human Gene Mutation Database (HGMD) and neutral INDELs from the 1000 Genomes Project. The potential of these two types of INDELs to affect binding-site affinity of RNA-binding proteins (RBPs) was then evaluated. We identified several sequence features that can distinguish disease-causing INDELs from neutral INDELs. Moreover, we built a machine-learning predictor called PinPor (predicting pathogenic small insertions and deletions affecting post-transcriptional regulation, http://watson.compbio.iupui.edu/pinpor/) to ascertain which newly observed INDELs are likely to be pathogenic. Our results show that disease-causing INDELs are more likely to ablate RBP-binding sites and tend to affect more RBP-binding sites than neutral INDELs. Additionally, disease-causing INDELs give rise to greater deviations in binding affinity than neutral INDELs. We also demonstrated that disease-causing INDELs may be distinguished from neutral INDELs by several sequence features, such as their proximity to splice sites and their potential effects on RNA secondary structure. This predictor showed satisfactory performance in identifying numerous pathogenic INDELs, with a Matthews correlation coefficient (MCC) value of 0.51 and an accuracy of 0.75.

The influenza A segment 7 mRNA 3' splice site pseudoknot/hairpin family

The 3′ splice site of the influenza A segment 7 transcript is utilized to produce mRNA for the critical M2 ion-channel protein. In solution a 63 nt fragment that includes this region can adopt two conformations: a pseudoknot and a hairpin. In each conformation, the splice site, a binding site for the SF2/ASF exonic splicing enhancer and a polypyrimidine tract, each exists in a different structural context. The most dramatic difference occurs for the splice site. In the hairpin the splice site is between two residues that are involved in a 2 by 2 nucleotide internal loop. In the pseudoknot, however, these bases are canonically paired within one of the pseudoknotted helices. The conformational switching observed in this region has implications for the regulation of splicing of the segment 7 mRNA. A measure of stability of the structures also shows interesting trends with respect to host specificity: avian strains tend to be the most stable, followed by swine and then human.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Solution structure of the HIV-1 exon splicing silencer 3

Alternative splicing of the human immunodeficiency virus type 1 (HIV-1) genomic RNA is necessary to produce the complete viral protein complement, and aberrations in the splicing pattern impair HIV-1 replication. Genome splicing in HIV-1 is tightly regulated by the dynamic assembly/disassembly of trans host factors with cis RNA control elements. The host protein, heterogeneous nuclear ribonucleoprotein (hnRNP) A1, regulates splicing at several highly conserved HIV-1 3' splice sites by binding 5'-UAG-3' elements embedded within regions containing RNA structure. The physical determinants of hnRNP A1 splice site recognition remain poorly defined in HIV-1, thus precluding a detailed understanding of the molecular basis of the splicing pattern. Here, the three-dimensional structure of the exon splicing silencer 3 (ESS3) from HIV-1 has been determined using NMR spectroscopy. ESS3 adopts a 27-nucleotide hairpin with a 10-bp A-form stem that contains a pH-sensitive A(+)C wobble pair. The seven-nucleotide hairpin loop contains the high-affinity hnRNP-A1-responsive 5'-UAGU-3' element and a proximal 5'-GAU-3' motif. The NMR structure shows that the heptaloop adopts a well-organized conformation stabilized primarily by base stacking interactions reminiscent of a U-turn. The apex of the loop is quasi-symmetric with UA dinucleotide steps from the 5'-GAU-3' and 5'-UAGU-3' motifs stacking on opposite sides of the hairpin. As a step towards understanding the binding mechanism, we performed calorimetric and NMR titrations of several hnRNP A1 subdomains into ESS3. The data show that the UP1 domain forms a high-affinity (K(d)=37.8±1.1 nM) complex with ESS3 via site-specific interactions with the loop.

HMGA1a is involved in specific splice site regulation of human immunodeficiency virus type 1

Human immunodeficiency virus type 1 (HIV-1) utilizes a highly complex splice site regulation system, taking advantage of host proteins, to express its own viral protein in an orderly way. We show here that one of the host proteins, high mobility group A protein 1a (HMGA1a), is involved in splice site regulation of 3' splice site 2 (A2) and 5'splice site 3 (D3) of HIV-1 genomic RNA. shRNA knockdown of HMGA1 in HeLa cells resulting in a decrease of HMGA1 showed a significant decrease of Vpr mRNA. RNA electrophoretic mobility shift assays showed HMGA1a specifically binds to a sequence adjacently upstream D3. In vitro splicing using heterologous pre-mRNA with A2 and D3, showed HMGA1a induced a splicing intermediate which decreased when an RNA decoy of the HMGA1a binding site was added. RT-PCR of in vitro splicing products revealed that HMGA1a induced an incomplete splicing product resulting from usage of A2 but inhibition of D3, which is reminiscent of the splicing pattern necessary for Vpr mRNA formation. HMGA1a interacted with hnRNPA1 shown by coimmunoprecipitation and supershifted U1 snRNP in an RNA electrophoretic mobility shift assay. We conclude that HMGA1a anchors U1 snRNP to inhibit D3 function, and that HMGA1a inhibits hnRNPA1 function on exon splicing silencer of Vpr (ESSV) to activate A2 function. We show here for the first time that HMGA1a is involved in specific splice site regulation of HIV-1.

GC content around splice sites affects splicing through pre-mRNA secondary structures

Background

Alternative splicing increases protein diversity by generating multiple transcript isoforms from a single gene through different combinations of exons or through different selections of splice sites. It has been reported that RNA secondary structures are involved in alternative splicing. Here we perform a genomic study of RNA secondary structures around splice sites in humans (Homo sapiens), mice (Mus musculus), fruit flies (Drosophila melanogaster), and nematodes (Caenorhabditis elegans) to further investigate this phenomenon.

Results

We observe that GC content around splice sites is closely associated with the splice site usage in multiple species. RNA secondary structure is the possible explanation, because the structural stability difference among alternative splice sites, constitutive splice sites, and skipped splice sites can be explained by the GC content difference. Alternative splice sites tend to be GC-enriched and exhibit more stable RNA secondary structures in all of the considered species. In humans and mice, splice sites of first exons and long exons tend to be GC-enriched and hence form more stable structures, indicating the special role of RNA secondary structures in promoter proximal splicing events and the splicing of long exons. In addition, GC-enriched exon-intron junctions tend to be overrepresented in tissue-specific alternative splice sites, indicating the functional consequence of the GC effect. Compared with regions far from splice sites and decoy splice sites, real splice sites are GC-enriched. We also found that the GC-content effect is much stronger than the nucleotide-order effect to form stable secondary structures.

Conclusion

All of these results indicate that GC content is related to splice site usage and it may mediate the splicing process through RNA secondary structures.

Structural and functional analysis of the Rous Sarcoma virus negative regulator of splicing and demonstration of its activation by the 9G8 SR protein

Retroviruses require both spliced and unspliced RNAs for replication. Accumulation of Rous Sarcoma virus (RSV) unspliced RNA depends upon the negative regulator of splicing (NRS). Its 5′-part is considered as an ESE binding SR proteins. Its 3′-part contains a decoy 5′-splice site (ss), which inhibits splicing at the bona fide 5′-ss. Only the 3D structure of a small NRS fragment had been experimentally studied. Here, by chemical and enzymatic probing, we determine the 2D structure of the entire RSV NRS. Structural analysis of other avian NRSs and comparison with all sequenced avian NRSs is in favour of a phylogenetic conservation of the NRS 2D structure. By combination of approaches: (i) in vitro and in cellulo splicing assays, (ii) footprinting assays and (iii) purification and analysis of reconstituted RNP complex, we define a small NRS element retaining splicing inhibitory property. We also demonstrate the capability of the SR protein 9G8 to increase NRS activity in vitro and in cellulo. Altogether these data bring new insights on how NRS fine tune splicing activity.

Excessive RNA splicing and inhibition of HIV-1 replication induced by modified U1 small nuclear RNAs

HIV-1 RNA undergoes a complex splicing process whereby over 40 different mRNA species are produced by alternative splicing. In addition, approximately half of the RNA transcripts remain unspliced and either are used to encode Gag and Gag-Pol proteins or are packaged into virions as genomic RNA. It has previously been shown that HIV-1 splicing is regulated by cis elements that bind to cellular factors. These factors either enhance or repress definition of exons that are flanked by the HIV-1 3' splice sites. Here we report that expression of modified U1 snRNPs with increased affinity to HIV-1 downstream 5' splice sites and to sequences within the first tat coding exon act to selectively increase splicing at the upstream 3' splice sites in cotransfected 293T cells. This results in a decrease of unspliced viral RNA levels and an approximately 10-fold decrease in virus production. In addition, excessive splicing of viral RNA is concomitant with a striking reduction in the relative amounts of Gag processing intermediates and products. We also show that T cell lines expressing modified U1 snRNAs exhibit reduced HIV-1 replication. Our results suggest that induction of excessive HIV-1 RNA splicing may be a novel strategy to inhibit virus replication in human patients.

Influenza virus RNA structure: unique and common features

The influenza A virus genome consists of eight negative-sense RNA segments. Here we review the currently available data on structure-function relationships in influenza virus RNAs. Various ideas and hypotheses about the roles of influenza virus RNA folding in the virus replication are also discussed in relation to other viruses.

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.

Antagonistic factors control the unproductive splicing of SC35 terminal intron

Alternative splicing is regulated in part by variations in the relative concentrations of a variety of factors, including serine/arginine-rich (SR) proteins. The SR protein SC35 self-regulates its expression by stimulating unproductive splicing events in the 3' untranslated region of its own pre-mRNA. Using various minigene constructs containing the terminal retained intron and flanking exons, we identified in the highly conserved last exon a number of exonic splicing enhancer elements responding specifically to SC35, and showed an inverse correlation between affinity of SC35 and enhancer strength. The enhancer region, which is included in a long stem loop, also contains repressor elements, and is recognized by other RNA-binding proteins, notably hnRNP H protein and TAR DNA binding protein (TDP-43). Finally, in vitro and in cellulo experiments indicated that hnRNP H and TDP-43 antagonize the binding of SC35 to the terminal exon and specifically repress the use of SC35 terminal 3' splice site. Our study provides new information about the molecular mechanisms of SC35-mediated splicing activation. It also highlights the existence of a complex network of self- and cross-regulatory mechanisms between splicing regulators, which controls their homeostasis and offers many ways of modulating their concentration in response to the cellular environment.

Role of RNA structure in regulating pre-mRNA splicing

Pre-mRNA splicing involves removing non-coding introns from RNA transcripts. It is carried out by the spliceosome, along with other auxiliary factors. In general, research in splicing has focused on the sequences within the pre-mRNA, without considering the structures that these sequences might form. We propose that the role of RNA structure deserves more consideration when thinking about splicing mechanisms. RNA structures can inhibit or aid binding of spliceosomal components to the pre-mRNA, or can increase splicing efficiency by bringing important sequences into close proximity. Recent reports have identified proteins and small molecules that can regulate splicing by modulating RNA structures, thereby expanding our knowledge of the mechanisms used to regulate splicing.

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.

Inhibition of human immunodeficiency virus (HIV-1) infection in human peripheral blood leucocytes-SCID reconstituted mice by rapamycin

The capacity of the immunomodulatory drug rapamycin (RAPA) to inhibit replication of the CCR5 strain of human immunodeficiency virus (HIV) in vitro prompted us to test its effects in a murine preclinical model of HIV infection. RAPA (0.6 or 6 mg/kg body weight) or its vehicle were administered daily, per os, to SCID mice reconstituted with human peripheral blood leucocytes (hu-PBL) starting 2 days before the intraperitoneal challenge with the R5 tropic SF162 strain of HIV-1 (1000 50% tissue culture infective dose/ml). Relative to hu-PBL-SCID mice that received no treatment, HIV-infected hu-PBL-SCID mice treated with the vehicle control for 3 weeks exhibited a severe depletion of CD4(+) cells (90%), an increase in CD8(+) cells and an inversion of the CD4(+)/CD8(+) cell ratio. In contrast, treatment of HIV-infected mice with RAPA prevented a decrease in CD4(+) cells and the increase of CD8(+) cells, thereby preserving the original CD4(+):CD8(+) cell ratio. Viral infection also resulted in the detection of HIV-DNA within peritoneal cells and spleen, and lymph node tissues of the vehicle-treated mice within 3 weeks of the viral challenge. In contrast, treatment with RAPA decreased cellular provirus integration and reduced HIV-RNA levels in the blood. Furthermore, in co-cultivation assays, spleens from RAPA-treated mice exhibited a reduced capacity for infecting allogeneic T cells which was dose-dependent. These data show that RAPA possesses powerful anti-viral activity against R5 strains of HIV in vivo and support the use of additional studies to evaluate the potential application of this drug in the management of HIV patients.

Extensive purifying selection acting on synonymous sites in HIV-1 Group M sequences

Background

Positive selection pressure acting on protein-coding sequences is usually inferred when the rate of nonsynonymous substitution is greater than the synonymous rate. However, purifying selection acting directly on the nucleotide sequence can lower the synonymous substitution rate. This could result in false inference of positive selection because when synonymous changes at some sites are under purifying selection, the average synonymous rate is an underestimate of the neutral rate of evolution. Even though HIV-1 coding sequences contain a number of regions that function at the nucleotide level, and are thus likely to be affected by purifying selection, studies of positive selection assume that synonymous substitutions can be used to estimate the neutral rate of evolution.

Results

We modelled site-to-site variation in the synonymous substitution rate across coding regions of the HIV-1 genome. Synonymous substitution rates were found to vary significantly within and between genes. Surprisingly, regions of the genome that encode proteins in more than one frame had significantly higher synonymous substitution rates than regions coding in a single frame. We found evidence of strong purifying selection pressure affecting synonymous mutations in fourteen regions with known functions. These included an exonic splicing enhancer, the rev-responsive element, the poly-purine tract and a transcription factor binding site. A further five highly conserved regions were located within known functional domains. We also found four conserved regions located in env and vpu which have not been characterized previously.

Conclusion

We provide the coordinates of genomic regions with markedly lower synonymous substitution rates, which are putatively under the influence of strong purifying selection pressure at the nucleotide level as well as regions encoding proteins in more than one frame. These regions should be excluded from studies of positive selection acting on HIV-1 coding regions.

The secondary structure of the human immunodeficiency virus type 1 transcript modulates viral splicing and infectivity

Splicing of human immunodeficiency virus type 1 (HIV-1) exon 6D is regulated by the presence of a complex splicing regulatory element (SRE) sequence that interacts with the splicing factors hnRNP H and SC35. In this work, we show that, in the context of the wild-type viral sequence, hnRNP H acts as a repressor of exon 6D inclusion independent of its binding to the SRE. However, hnRNP H binding to the SRE acts as an enhancer of exon 6D inclusion in the presence of a critical T-to-C mutation. These seemingly contrasting functional properties of hnRNP H appear to be caused by a change in the RNA secondary structure induced by the T-to-C mutation that affects the spatial location of bound hnRNP H with respect to the exon 6D splicing determinants. We propose a new regulatory mechanism mediated by RNA folding that may also explain the dual properties of hnRNP H in splicing regulation.

RNA structure modulates splicing efficiency at the human immunodeficiency virus type 1 major splice donor

The untranslated leader of the human immunodeficiency virus type 1 (HIV-1) RNA genome encodes essential sequence and structural motifs that control various replication steps. The 5' splice site or splice donor (SD) is embedded in a semistable hairpin, but the function of this structure is unknown. We stabilized this SD hairpin by creating an additional base pair and demonstrated a severe HIV-1 replication defect. A splicing defect was apparent in RNA analyses of virus-infected cells and cells transfected with appropriate reporter constructs. We selected multiple virus revertants in search for interesting second-site escape pathways. Most revertants acquired an additional mutation that modulated the stability of the mutant SD hairpin. One revertant acquired a single nucleotide change in the upstream DIS hairpin. We demonstrate that a novel SD site is created by this upstream mutation, which obviously reduces the number of leader nucleotides that are included in spliced HIV-1 transcripts. These results suggest a novel role of RNA structure in the regulation of HIV-1 splicing.

(1)H, (13)C and (15)N resonance assignment of the first N-terminal RNA recognition motif (RRM) of the human heterogeneous nuclear ribonucleoprotein H (hnRNP H)

Human heterogeneous nuclear ribonucleoprotein H (hnRNP H) regulates alternative splicing of HIV-1 Tat pre-mRNA. The structure of the first N-terminal domain (residues 1-104) of hnRNP H was solved and its binding to an exonic splicing silencer (pESS2) studied. For this, all backbone and 85% of side-chain resonance frequencies were assigned.

Evaluation of the OPTC gene in primary open angle glaucoma: functional significance of a silent change

Background

We investigated the molecular basis of primary open-angle glaucoma (POAG) using Opticin (OPTC) as a candidate gene on the basis of its expression in the trabecular meshwork cells involved in the disease pathogenesis. Two hundred POAG patients and 100 controls were enrolled in this study. The coding sequence of OPTC was amplified by PCR from genomic DNA of POAG patients, followed by SSCP, DHPLC and DNA sequencing. Subsequent bioinformatic analysis, site-directed mutagenesis, quantitative RT-PCR and western blot experiments were performed to address the functional significance of a 'silent' change in the OPTC coding region while screening for mutations in POAG patients.

Results

We detected two missense (p.Glu66Gly & p.Ile89Thr) and one silent change (p.Phe162Phe; c.602 C>T) that was present in 3 different patients but in none of the 100 controls screened. The mutant (c.602T) mRNA was predicted to have remarkably different secondary structure compared to the wild-type transcript by in silico approaches. Subsequent wet-lab experiments showed lower expression of the gene both at the mRNA and protein levels.

Conclusion

Our study suggests OPTC as a candidate gene for POAG. Further, it highlights the importance of investigating the 'silent' variations for functional implication that might not be apparent from only in silico analysis.

Modulating role of RNA structure in alternative splicing of a critical exon in the spinal muscular atrophy genes

Humans have two nearly identical copies of the survival motor neuron (SMN ) gene, SMN1 and SMN2. Homozygous loss of SMN1 causes spinal muscular atrophy (SMA). SMN2 is unable to prevent the disease due to skipping of exon 7. Using a systematic approach of in vivo selection, we have previously demonstrated that a weak 5' splice site (ss) serves as the major cause of skipping of SMN2 exon 7. Here we show the inhibitory impact of RNA structure on the weak 5' ss of exon 7. We call this structure terminal stem-loop 2 (TSL2). Confirming the inhibitory nature of TSL2, point mutations that destabilize TSL2 promote exon 7 inclusion in SMN2, whereas strengthening of TSL2 promotes exon 7 skipping even in SMN1. We also demonstrate that TSL2 negatively affects the recruitment of U1snRNP at the 5' ss of exon 7. Using enzymatic structure probing, we confirm that the sequence at the junction of exon 7/intron 7 folds into TSL2 and show that mutations in TSL2 cause predicted structural changes in this region. Our findings reveal for the first time the critical role of RNA structure in regulation of alternative splicing of human SMN.

The strength of the HIV-1 3' splice sites affects Rev function

Background

The HIV-1 Rev protein is a key component in the early to late switch in HIV-1 splicing from early intronless (e.g. tat, rev) to late intron-containing Rev-dependent (e.g. gag, vif, env) transcripts. Previous results suggested that cis-acting sequences and inefficient 5' and 3' splice sites are a prerequisite for Rev function. However, we and other groups have shown that two of the HIV-1 5' splice sites, D1 and D4, are efficiently used in vitro and in vivo. Here, we focus on the efficiency of the HIV-1 3' splice sites taking into consideration to what extent their intrinsic efficiencies are modulated by their downstream cis-acting exonic sequences. Furthermore, we delineate their role in RNA stabilization and Rev function.

Results

In the presence of an efficient upstream 5' splice site the integrity of the 3' splice site is not essential for Rev function whereas an efficient 3' splice site impairs Rev function. The detrimental effect of a strong 3' splice site on the amount of Rev-dependent intron-containing HIV-1 glycoprotein coding (env) mRNA is not compensatable by weakening the strength of the upstream 5' splice site. Swapping the HIV-1 3' splice sites in an RRE-containing minigene, we found a 3' splice site usage which was variably dependent on the presence of the usual downstream exonic sequence. The most evident activation of 3' splice site usage by its usual downstream exonic sequence was observed for 3' splice site A1 which was turned from an intrinsic very weak 3' splice site into the most active 3' splice site, even abolishing Rev activity. Performing pull-down experiments with nuclear extracts of HeLa cells we identified a novel ASF/SF2-dependent exonic splicing enhancer (ESE) within HIV-1 exon 2 consisting of a heptameric sequence motif occurring twice (M1 and M2) within this short non-coding leader exon. Single point mutation of M1 within an infectious molecular clone is detrimental for HIV-1 exon 2 recognition without affecting Rev-dependent vif expression.

Conclusion

Under the conditions of our assay, the rate limiting step of retroviral splicing, competing with Rev function, seems to be exclusively determined by the functional strength of the 3' splice site. The bipartite ASF/SF2-dependent ESE within HIV-1 exon 2 supports cross-talk between splice site pairs across exon 2 (exon definition) which is incompatible with processing of the intron-containing vif mRNA. We propose that Rev mediates a switch from exon to intron definition necessary for the expression of all intron-containing mRNAs.

Biochemical and NMR Study on the Competition between Proteins SC35, SRp40, and Heterogeneous Nuclear Ribonucleoprotein A1 at the HIV-1 Tat Exon 2 Splicing Site

The human immunodeficiency virus, type 1, Tat protein plays a key role in virus multiplication. Because of its apoptotic property, its production is highly controlled. It depends upon the A3 splicing site utilization. A key control of site A3 activity is the ESS2 splicing silencer, which is located within the long stem-loop structure 3 (SLS3), far downstream from site A3. Here, by enzymatic footprints, we demonstrate the presence of several heterogeneous nuclear ribonucleoprotein (hnRNP) A1-binding sites on SLS3 and show the importance of the C-terminal Gly domain of hnRNP A1 in the formation of stable complexes containing several hnRNP A1 molecules bound on SLS3. Mutations in each of the UAG triplets in ESS2 strongly reduce the overall hnRNP A1 binding, showing the central role of ESS2 in hnRNP A1 assembly on SLS2-SLS3. Using NMR spectroscopy, we demonstrate the direct interaction of ESS2 with the RNA recognition motifs domains of hnRNP A1. This interaction has limited effect on the RNA two-dimensional structure. The SR proteins SC35 and SRp40 were found previously to be strong activators of site A3 utilization. By enzymatic and chemical footprints, we delineate their respective binding sites on SLS2 and SLS3 and find a strong similarity between the hnRNP A1-, SC35-, and SRp40-binding sites. The strongest SC35-binding site only has a modest contribution to site A3 activation. Hence, the main role of SR proteins at site A3 is to counteract hnRNP A1 binding on ESS2 and ESE2. Indeed, we found that ESE2 has inhibitory properties because of its ability to bind hnRNP A1.

A suboptimal 5' splice site downstream of HIV-1 splice site A1 is required for unspliced viral mRNA accumulation and efficient virus replication

Background

Inefficient alternative splicing of the human immunodeficiency virus type 1(HIV-1) primary RNA transcript results in greater than half of all viral mRNA remaining unspliced. Regulation of HIV-1 alternative splicing occurs through the presence of suboptimal viral 5' and 3' splice sites (5' and 3'ss), which are positively regulated by exonic splicing enhancers (ESE) and negatively regulated by exonic splicing silencers (ESS) and intronic splicing silencers (ISS). We previously showed that splicing at HIV-1 3'ss A2 is repressed by ESSV and enhanced by the downstream 5'ss D3 signal. Disruption of ESSV results in increased vpr mRNA accumulation and exon 3 inclusion, decreased accumulation of unspliced viral mRNA, and decreased virus production.

Results

Here we show that optimization of the 5'ss D2 signal results in increased splicing at the upstream 3'ss A1, increased inclusion of exon 2 into viral mRNA, decreased accumulation of unspliced viral mRNA, and decreased virus production. Virus production from the 5'ss D2 and ESSV mutants was rescued by transient expression of HIV-1 Gag and Pol. We further show that the increased inclusion of either exon 2 or 3 does not significantly affect the stability of viral mRNA but does result in an increase and decrease, respectively, in HIV-1 mRNA levels. The changes in viral mRNA levels directly correlate with changes in tat mRNA levels observed upon increased inclusion of exon 2 or 3.

Conclusion

These results demonstrate that splicing at HIV-1 3'ss A1 is regulated by the strength of the downstream 5'ss signal and that suboptimal splicing at 3'ss A1 is necessary for virus replication. Furthermore, the replication defective phenotype resulting from increased splicing at 3'ss A1 is similar to the phenotype observed upon increased splicing at 3'ss A2. Further examination of the role of 5'ss D2 and D3 in the alternative splicing of 3'ss A1 and A2, respectively, is necessary to delineate a role for non-coding exon inclusion in HIV-1 replication.

The Effect of a Single Nucleotide Substitution in the Splicing Silencer in the tat/rev Intron on HIV Type 1 Envelope Expression

A complex mRNA splicing pattern, which remains to be fully characterized, influences HIV-1 gene expression. In this study, poor envelope expression of a primary HIV-1 isolate was observed and linked to increased splicing of the two coding exons of tat/rev. The substitution of a nucleotide G, located 28 nucleotides upstream of the splice acceptor site SA7 in the recently identified intron splicing silencer sequence, was found to be responsible for the poor envelope expression. A single nucleotide substitution of G with A at this position results in a poor envelope expression phenotype. Moreover, substitution of the nucleotide G with any other nucleotide in an infectious HIV-1 proviral clone, HXB2RU3, results in poor envelope expression. The substitution of this nucleotide reduces the hnRNP A1 binding affinity but increases the splicing of env mRNA. The nucleotide G at this position is highly conserved among HIV-1 isolates and appears to play a critical role in HIV-1 splicing.

Characterization of the molecular mechanisms involved in the differential production of erythrose-4-phosphate dehydrogenase, 3-phosphoglycerate kinase and class II fructose-1,6-bisphosphate aldolase in Escherichia coli

A gapA-pgk gene tandem coding the glyceraldehyde 3-phosphate dehydrogenase and 3-phosphoglycerate kinase, is most frequently found in bacteria. However, in Enterobacteriaceae, gapA is replaced by an epd open reading frame (ORF) coding an erythrose-4-phosphate dehydrogenase and an fbaA ORF coding the class II fructose-1,6-bisphosphate aldolase follows pgk. Although epd expression is very low in Escherichia coli, we show that, in the presence of glucose, the 3 epd, pgk and fbaA ORFs are efficiently cotranscribed from promoter epd P0. Conservation of promoter epd P0 is likely due to its important role in modulation of the metabolic flux during glycolysis and gluconeogenesis. As a consequence, we found that the epd translation initiation region and ORF have been adapted in order to limit epd translation and to create an efficient RNase E entry site. We also show that fbaA is cotranscribed with pgk, from promoter epd P0 or an internal pgk P1 promoter of the extended -10 class. The differential expression of pgk and fbaA also depends upon an RNase E segmentation process, leading to individual mRNAs with different stabilities. The secondary structures of the RNA regions containing the RNase E sites were experimentally determined which brings important information on the structural features of RNase E ectopic sites.

An exonic splicing silencer downstream of the 3' splice site A2 is required for efficient human immunodeficiency virus type 1 replication

Alternative splicing of the human immunodeficiency virus type 1 (HIV-1) genomic mRNA produces more than 40 unique viral mRNA species, of which more than half remain incompletely spliced within an HIV-1-infected cell. Regulation of splicing at HIV-1 3' splice sites (3'ss) requires suboptimal polypyrimidine tracts, and positive or negative regulation of splicing occurs through binding of cellular factors to cis-acting splicing regulatory elements. We have previously shown that splicing at HIV-1 3'ss A2, which produces vpr mRNA and promotes inclusion of HIV-1 exon 3, is repressed by the hnRNP A/B-dependent exonic splicing silencer ESSV. Here we show that ESSV activity downstream of 3'ss A2 is localized to a 16-nucleotide element within HIV-1 exon 3. HIV-1 replication was reduced by 95% when ESSV was inactivated by mutagenesis. Reduced replication was concomitant with increased inclusion of exon 3 within spliced viral mRNA and decreased accumulation of unspliced viral mRNA, resulting in decreased cell-associated p55 Gag. Prolonged culture of ESSV mutant viruses resulted in two independent second-site reversions disrupting the splice sites that define exon 3, 3'ss A2 and 5' splice site D3. Either of these changes restored both HIV-1 replication and regulated viral splicing. Therefore, inhibition of HIV-1 3'ss A2 splicing is necessary for HIV-1 replication.

Dual effect of the SR proteins ASF/SF2, SC35 and 9G8 on HIV-1 RNA splicing and virion production

In HIV-1 infected cells transcription of the integrated provirus generates the single full length 9 kb viral RNA, a major fraction of which is spliced to produce the single-spliced 4 kb RNAs and the multiple-spliced 2 kb RNAs. These spliced RNAs are the messengers for the Env glycoproteins and the viral regulatory factors. The cellular SR and hnRNP proteins were shown in vitro to control alternative splicing by binding cis-regulatory elements on the viral RNA. To better understand in vivo the role of the SR proteins on HIV-1 genomic RNA splicing and virion production, we used a human cell line expressing high levels of complete HIV-1 and either one of the ASF/SF2, SC35, and 9G8 SR proteins. Results show that over-expressing SR proteins caused a large reduction of genomic RNA and that each SR protein modified the viral 9 kb RNA splicing pattern in a specific mode. In fact, ASF/SF2 increased the level of Vpr RNA while SC35 and 9G8 caused a large increase in Tat RNA. As expected, overexpressing SR proteins caused a strong reduction of total Gag made. However, we observed by immuno-confocal microscopy an accumulation of Gag at the plasma membrane and in intracellular compartments while there is a dramatic reduction of Env protein made in most cells. Due to the negative impact of the SR proteins on the levels of genomic RNA and HIV-1 structural proteins much less virions were produced which retained part of their infectivity. In conclusion, SR proteins can down-regulate the late steps of HIV-1 replication.

Evidence for the regulation of alternative splicing via complementary DNA sequence repeats

MOTIVATION

While the mechanism for regulating alternative splicing is poorly understood, secondary structure has been shown to be integral to this process. Due to their propensity for forming complementary hairpin loops and their elevated mutation rates, tandem repeated sequences have the potential to influence splicing regulation.

RESULTS

An analysis of human intronic sequences reveals a strong correlation between alternative splicing and the prevalence of mono- through hexanucleotide tandem repeats that may engage in complementary pairing in introns that flank alternatively spliced exons. While only 44% of the 18 173 genes in the Human Alternative Splicing Database are known to be alternatively spliced, they contain 84% of the 694 237 intronic complementary repeat pairs. Significantly, the normalized frequency and distribution of repeat sequences, independent of their potential for pairing, are indistinguishable between alternatively spliced and non-alternatively spliced genes. Thus, the increased prevalence of repeats with pairing potential in alternatively spliced genes is not merely a consequence of more repeats or repeat composition bias. These results suggest that complementary repeats may play a role in the regulation of alternative splicing.

CONTACT

harold.garner@utsouthwestern.edu.

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.

A novel approach to describe a U1 snRNA binding site

RNA duplex formation between U1 snRNA and a splice donor (SD) site can protect pre-mRNA from degradation prior to splicing and initiates formation of the spliceosome. This process was monitored, using sub-genomic HIV-1 expression vectors, by expression analysis of the glycoprotein env, whose formation critically depends on functional SD4. We systematically derived a hydrogen bond model for the complementarity between the free 5' end of U1 snRNA and 5' splice sites and numerous mutations following transient transfection of HeLa-T4+ cells with 5' splice site mutated vectors. The resulting model takes into account number, interdependence and neighborhood relationships of predicted hydrogen bond formation in a region spanning the three most 3' base pairs of the exon (-3 to -1) and the eight most 5' base pairs of the intron (+1 to +8). The model is represented by an algorithm classifying U1 snRNA binding sites which can or cannot functionally substitute SD4 with respect to Rev-mediated env expression. In a data set of 5' splice site mutations of the human ATM gene we found a significant correlation between the algorithmic classification and exon skipping (P = 0.018, chi2-test), showing that the applicability of the proposed model reaches far beyond HIV-1 splicing. However, the algorithmic classification must not be taken as an absolute measure of SD usage as it may be modified by upstream sequence elements. Upstream to SD4 we identified a fragment supporting ASF/SF2 binding. Mutating GAR nucleotide repeats within this site decreased the SD4-dependent Rev-mediated env expression, which could be balanced simply by artificially increasing the complementarity of SD4.