Influence of the small leader exons 2 and 3 on human immunodeficiency virus type 1 gene expression

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.

Genome-Wide Analysis of Heterogeneous Nuclear Ribonucleoprotein (hnRNP) Binding to HIV-1 RNA Reveals a Key Role for hnRNP H1 in Alternative Viral mRNA Splicing

Alternative splicing of HIV-1 mRNAs increases viral coding potential and controls the levels and timing of gene expression. HIV-1 splicing is regulated in part by heterogeneous nuclear ribonucleoproteins (hnRNPs) and their viral target sequences, which typically repress splicing when studied outside their native viral context. Here, we determined the location and extent of hnRNP binding to HIV-1 mRNAs and their impact on splicing in a native viral context. Notably, hnRNP A1, hnRNP A2, and hnRNP B1 bound to many dispersed sites across viral mRNAs. Conversely, hnRNP H1 bound to a few discrete purine-rich sequences, a finding that was mirrored in vitro hnRNP H1 depletion and mutation of a prominent viral RNA hnRNP H1 binding site decreased the use of splice acceptor A1, causing a deficit in Vif expression and replicative fitness. This quantitative framework for determining the regulatory inputs governing alternative HIV-1 splicing revealed an unexpected splicing enhancer role for hnRNP H1 through binding to its target element.IMPORTANCE Alternative splicing of HIV-1 mRNAs is an essential yet quite poorly understood step of virus replication that enhances the coding potential of the viral genome and allows the temporal regulation of viral gene expression. Although HIV-1 constitutes an important model system for general studies of the regulation of alternative splicing, the inputs that determine the efficiency with which splice sites are utilized remain poorly defined. Our studies provide an experimental framework to study an essential step of HIV-1 replication more comprehensively and in much greater detail than was previously possible and reveal novel cis-acting elements regulating HIV-1 splicing.

Behind the scenes of HIV-1 replication: Alternative splicing as the dependency factor on the quiet

Alternative splicing plays a key role in the HIV-1 life cycle and is essential to maintain an equilibrium of mRNAs that encode viral proteins and polyprotein-isoforms. In particular, since all early HIV-1 proteins are expressed from spliced intronless and late enzymatic and structural proteins from intron containing, i.e. splicing repressed viral mRNAs, cellular splicing factors and splicing regulatory proteins are crucial for the replication capacity. In this review, we will describe the complex network of cis-acting splicing regulatory elements (SREs), which are mainly localized in the neighbourhoods of all HIV-1 splice sites and warrant the proper ratio of individual transcript isoforms. Since SREs represent binding sites for trans-acting cellular splicing factors interacting with the cellular spliceosomal apparatus we will review the current knowledge of interactions between viral RNA and cellular proteins as well as their impact on viral replication. Finally, we will discuss potential therapeutic approaches targeting HIV-1 alternative splicing.

Characterizing HIV-1 Splicing by Using Next-Generation Sequencing

Full-length human immunodeficiency virus type 1 (HIV-1) RNA serves as the genome or as an mRNA, or this RNA undergoes splicing using four donors and 10 acceptors to create over 50 physiologically relevant transcripts in two size classes (1.8 kb and 4 kb). We developed an assay using Primer ID-tagged deep sequencing to quantify HIV-1 splicing. Using the lab strain NL4-3, we found that A5 (env/nef) is the most commonly used acceptor (about 50%) and A3 (tat) the least used (about 3%). Two small exons are made when a splice to acceptor A1 or A2 is followed by activation of donor D2 or D3, and the high-level use of D2 and D3 dramatically reduces the amount of vif and vpr transcripts. We observed distinct patterns of temperature sensitivity of splicing to acceptors A1 and A2. In addition, disruption of a conserved structure proximal to A1 caused a 10-fold reduction in all transcripts that utilized A1. Analysis of a panel of subtype B transmitted/founder viruses showed that splicing patterns are conserved, but with surprising variability of usage. A subtype C isolate was similar, while a simian immunodeficiency virus (SIV) isolate showed significant differences. We also observed transsplicing from a downstream donor on one transcript to an upstream acceptor on a different transcript, which we detected in 0.3% of 1.8-kb RNA reads. There were several examples of splicing suppression when the env intron was retained in the 4-kb size class. These results demonstrate the utility of this assay and identify new examples of HIV-1 splicing regulation. IMPORTANCE During HIV-1 replication, over 50 conserved spliced RNA variants are generated. The splicing assay described here uses new developments in deep-sequencing technology combined with Primer ID-tagged cDNA primers to efficiently quantify HIV-1 splicing at a depth that allows even low-frequency splice variants to be monitored. We have used this assay to examine several features of HIV-1 splicing and to identify new examples of different mechanisms of regulation of these splicing patterns. This splicing assay can be used to explore in detail how HIV-1 splicing is regulated and, with moderate throughput, could be used to screen for structural elements, small molecules, and host factors that alter these relatively conserved splicing patterns.

An intronic G run within HIV-1 intron 2 is critical for splicing regulation of vif mRNA

Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3' splice site (3'ss) A1 but lack splicing at 5'ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3'ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3'ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5'ss D2. Here we show that an intronic G run (G(I2)-1) represses the use of a second 5'ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of G(I2)-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

Chapter 1. Regulation of HIV-1 alternative RNA splicing and its role in virus replication

Over 40 different human immunodeficiency virus type 1 (HIV-1) mRNA species, both completely and incompletely spliced, are produced by alternative splicing of the primary viral RNA transcript. In addition, about half of the viral RNA remains unspliced and is transported to the cytoplasm where it is used both as mRNA and as genomic RNA. In general, the identities of the completely and incompletely spliced HIV-1 mRNA species are determined by the proximity of the open reading frames to the 5'-end of the mRNAs. The relative abundance of the mRNAs encoding the HIV-1 gene products is determined by the frequency of splicing at the different alternative 3'-splice sites. This chapter will highlight studies showing how HIV-1 uses exon definition to control the level of splicing at each of its 3'-splice sites through a combination of positively acting exonic splicing enhancer (ESE) elements, negatively acting exonic and intronic splicing silencer elements (ESS and ISS elements, respectively), and the 5'-splice sites of the regulated exons. Each of these splicing elements represent binding sites for cellular factors whose levels in the infected cell can determine the dominance of the positive or negative elements on HIV-1 alternative splicing. Both mutations of HIV-1 splicing elements and overexpression or inhibition of cellular splicing factors that bind to these elements have been used to show that disruption of regulated splicing inhibits HIV-1 replication. These studies have provided strong rationale for the investigation and development of antiviral drugs that specifically inhibit HIV-1 RNA splicing.

Analysis of Jembrana disease virus replication dynamics in vivo reveals strain variation and atypical responses to infection

Jembrana disease virus (JDV) is an acute lentiviral infection of Bali cattle in Indonesia. Data generated during a series of cattle infection experiments was examined and significant differences were identified in the mean plasma viral load on the first and second days of the febrile response in cattle infected with JDV(TAB/87) compared to those infected with JDV(PUL/01). The peak and total viral loads >or=10(6) genome copies/ml during the acute stage of the disease were significantly higher in JDV(TAB/87) infected cattle. JDV(PUL/01) infected cattle developed peak rectal temperatures earlier than the JDV(TAB/87) cattle but there were no differences in the duration of the febrile responses observed for the 2 groups of animals. The plasma viremia was above 10(6) genome copies/ml for almost 3 days longer in JDV(TAB/87) compared to JDV(PUL/01) infected cattle. Atypical responses to infection occurred in approximately 15% of experimentally infected animals, characterized by reduced viral loads, lower or absent febrile responses and absence of p26-specific antibody responses. Most of these cattle developed normal Tm-specific antibody responses between 4-12 weeks post-infection.

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.

A minimal uORF within the HIV-1 vpu leader allows efficient translation initiation at the downstream env AUG

The HIV-1 Vpu and Env proteins are translated from 16 alternatively spliced bicistronic mRNA isoforms. Translation of HIV-1 mRNAs generally follows the ribosome scanning mechanism. However, by using subgenomic env expression vectors, we found that translation of glycoprotein from polycistronic mRNAs was inconsistent with leaky scanning. Instead a conserved minimal upstream open reading frame (uORF) consisting only of a start and stop codon that overlaps with the vpu start site, appears to augment access to the env start codon downstream. Mutating the translational start and stop codons of this uORF resulted in up to fivefold reduction in Env expression. Removing the vpu uORF and increasing the strength of the authentic vpu initiation sequence abolished Env expression from subgenomic constructs and replication of HIV-1, whereas an identical increase in the strength of the minimal uORF initiation site did not alter Env expression.

Both linear and discontinuous ribosome scanning are used for translation initiation from bicistronic human immunodeficiency virus type 1 env mRNAs

Human immunodeficiency virus type 1 (HIV-1) generates 16 alternatively spliced isoforms of env mRNA that contain the same overlapping open reading frames for Vpu and Env proteins but differ in their 5' untranslated regions (UTR). A subset of env mRNAs carry the extra upstream Rev initiation codon in the 5' UTR. We explored the effect of the alternative UTR on the translation of Vpu and Env proteins from authentic env mRNAs expressed from cDNA constructs. Vpu expression from the subset of env mRNA isoforms with exons containing an upstream Rev AUG codon was minimal. However, every env mRNA isoform expressed similar levels of Env protein. Mutations that removed, altered the strength of, or introduced upstream AUG codons dramatically altered Vpu expression but had little impact on the consistent expression of Env. These data show that the different isoforms of env mRNA are not redundant but instead regulate Vpu production in HIV-1-infected cells. Furthermore, while the initiation of Vpu translation conforms to the leaky ribosome-scanning model, the consistent Env synthesis infers a novel, discontinuous ribosome-scanning mechanism to translate Env.

TCRzeta mRNA splice variant forms observed in the peripheral blood T cells from systemic lupus erythematosus patients

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease of unknown etiology. Tyrosine phosphorylation and protein expression of the T-cell receptor zeta chain (zeta) have been reported to be significantly decreased in SLE T cells. In addition, zeta mRNA with alternatively spliced 3' untranslated region (zetamRNA/as-3'UTR) is detected predominantly in SLE T cells, and aberrant zeta mRNA accompanied by the mutations in the open reading frame including zeta mRNA lacking exon7 (zetamRNA/exon7-) is observed in SLE T cells. These zeta mRNA splice variant forms exhibit a reduction in the expression of TCR/CD3 complex and zeta protein on their cell surface due to the instability of zeta mRNA splice variant forms as well as the reduction in interleukin (IL)-2 production after stimulating with anti-CD3 antibody. Data from cDNA microarray showed that 36 genes encoding cytokines and chemokines, including IL-2, IL-15, IL-18, and TGF-beta2, were down-regulated in the MA5.8 cells transfected with the zeta mRNA splice variant forms. Another 16 genes were up-regulated and included genes associated with membranous proteins and cell damage granules, including the genes encoding poliovirus-receptor-related 2, syndecan-1, and granzyme A.

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 hepatitis B virus PRE contains a splicing regulatory element

The posttranscriptional regulatory element (PRE) is considered to enhance hepatitis B virus (HBV) gene expression by facilitating the nuclear export of intronless viral subgenomic RNAs. Its role in the RNA metabolism of the viral pregenomic RNA (pgRNA) is currently unknown. We identified a positively cis-acting splicing regulatory element (SRE-1) and present two lines of evidence for its functionality. Firstly, in a heterologous context SRE-1 functionally substitutes for a retroviral bidirectional exonic splicing enhancer (ESE). As expected, SRE-1 is a splicing enhancer also in its natural viral sequence context, since deletion of SRE-1 reduces splicing of pgRNA in cell culture experiments. Secondly, we show that stimulation of HBV RNA splicing by the splicing factor PSF was repressed by the PRE. Analysis of a variety of PSF mutants indicated that RNA-binding and protein-protein interaction were required to enhance splicing. In addition, we show that the PRE contributed to pgRNA stability, but has little influence on its nuclear export. Herein, we report for the first time that the PRE harbors splicing stimulating and inhibiting regulatory elements controlling processing of the viral pregenome. We discuss a model in which the regulation of pgRNA splicing depends on cellular factors interacting with the PRE.

Disruption of the putative splice acceptor site for SIV(mac239)Vif reveals tight control of SIV splicing and impaired replication in Vif non-permissive cells

Vif is dispensable for simian immunodeficiency virus (SIV) replication in some cells, termed permissive (i.e., CEM-SS), but not in others, termed non-permissive (i.e., H9, CEMx174, and peripheral blood lymphocytes). Non-permissive cells express the RNA editing enzyme, APOBEC3G. To determine whether vif mRNA could be alternatively spliced, a mutation altering the putative vif splice acceptor site (SA1) was introduced into SIV(mac239) (SIV(Deltavif-SA)). Despite three consensus splice acceptor sites nearby SA1, SIV(Deltavif-SA) did not efficiently generate alternatively spliced vif mRNA. SIV(Deltavif-SA) was growth attenuated in CEMx174 and H9 cells but not in CEM-SS cells. Following SIV(Deltavif-SA), but not SIV(mac239), infection in either H9 or CEMx174 cells viral cDNA contained numerous G to A mutations; no such differences were observed in CEM-SS cells. This pattern is consistent with mutations generated by APOBEC3G in the absence of Vif. Therefore, efficient splicing of SIV vif mRNA is tightly controlled and requires the SA1 site.

Extended base pair complementarity between U1 snRNA and the 5' splice site does not inhibit splicing in higher eukaryotes, but rather increases 5' splice site recognition

Spliceosome formation is initiated by the recognition of the 5′ splice site through formation of an RNA duplex between the 5′ splice site and U1 snRNA. We have previously shown that RNA duplex formation between U1 snRNA and the 5′ splice site can protect pre-mRNAs from degradation prior to splicing. This initial RNA duplex must be disrupted to expose the 5′ splice site sequence for base pairing with U6 snRNA and to form the active spliceosome. Here, we investigated whether hyperstabilization of the U1 snRNA/5′ splice site duplex interferes with splicing efficiency in human cell lines or nuclear extracts. Unlike observations in Saccharomyces cerevisiae, we demonstrate that an extended U1 snRNA/5′ splice site interaction does not decrease splicing efficiency, but rather increases 5′ splice site recognition and exon inclusion. However, low complementarity of the 5′ splice site to U1 snRNA significantly increases exon skipping and RNA degradation. Although the splicing mechanisms are conserved between human and S.cerevisiae, these results demonstrate that distinct differences exist in the activation of the spliceosome.

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.

A splice variant of the TCR zeta mRNA lacking exon 7 leads to the down-regulation of TCR zeta, the TCR/CD3 complex, and IL-2 production in systemic lupus erythematosus T cells

The reduction or absence of TCR zeta-chain (zeta) expression in patients with systemic lupus erythematosus (SLE) is thought to be a factor in the pathogenesis of SLE. We previously reported a splice variant of zeta mRNA that lacks the 36-bp exon 7 (zeta mRNA/exon 7(-)) and is accompanied by the down-regulation of zeta protein in T cells from SLE patients. In this study, we show that EX7- mutants (MA5.8 cells deficient in zeta protein that have been transfected with zeta mRNA/exon 7(-)) exhibit a reduction in the expression of TCR/CD3 complex and zeta protein on their cell surface as well as a reduction in the production of IL-2 after stimulation with anti-CD3 Ab, compared with that in wild-type (WT) mutants (MA5.8 cells transfected with the WT zeta mRNA). Furthermore, real-time PCR analyses demonstrated that zeta mRNA/exon 7(-) in EX7- mutants was easily degraded compared with zeta mRNA by the WT mutants. Pulse-chase experiment showed zeta protein produced by this EX7- mutants was more rapidly decreased compared with the WT mutants. Thus, the lower stability of zeta mRNA/exon 7(-) might also be responsible for the reduced expression of the TCR/CD3 complex, including zeta protein, in SLE T cells.

Function of alternative splicing

Alternative splicing is one of the most important mechanisms to generate a large number of mRNA and protein isoforms from the surprisingly low number of human genes. Unlike promoter activity, which primarily regulates the amount of transcripts, alternative splicing changes the structure of transcripts and their encoded proteins. Together with nonsense-mediated decay (NMD), at least 25% of all alternative exons are predicted to regulate transcript abundance. Molecular analyses during the last decade demonstrate that alternative splicing determines the binding properties, intracellular localization, enzymatic activity, protein stability and posttranslational modifications of a large number of proteins. The magnitude of the effects range from a complete loss of function or acquisition of a new function to very subtle modulations, which are observed in the majority of cases reported. Alternative splicing factors regulate multiple pre-mRNAs and recent identification of physiological targets shows that a specific splicing factor regulates pre-mRNAs with coherent biological functions. Therefore, evidence is now accumulating that alternative splicing coordinates physiologically meaningful changes in protein isoform expression and is a key mechanism to generate the complex proteome of multicellular organisms.