Comparison of SIV and HIV-1 genomic RNA structures reveals impact of sequence evolution on conserved and non-conserved structural motifs

Minimally Modified HIV-1 Infection of Macaques: Development, Utility, and Limitations of Current Models

Nonhuman primate (NHP) studies that utilize simian immunodeficiency virus (SIV) to model human immunodeficiency virus (HIV-1) infection have proven to be powerful, highly informative research tools. However, there are substantial differences between SIV and HIV-1. Accordingly, there are numerous research questions for which SIV-based models are not well suited, including studies of certain aspects of basic HIV-1 biology, and pre-clinical evaluations of many proposed HIV-1 treatment, prevention, and vaccination strategies. To overcome these limitations of NHP models of HIV-1 infection, several groups have pursued the derivation of a minimally modified HIV-1 (mmHIV-1) capable of establishing pathogenic infection in macaques that authentically recapitulates key features of HIV-1 in humans. These efforts have focused on three complementary objectives: (1) engineering HIV-1 to circumvent species-specific cellular restriction factors that otherwise potently inhibit HIV-1 in macaques, (2) introduction of a C chemokine receptor type 5 (CCR5)-tropic envelope, ideally that can efficiently engage macaque CD4, and (3) correction of gene expression defects inadvertently introduced during viral genome manipulations. While some progress has been made toward development of mmHIV-1 variants for use in each of the three macaque species (pigtail, cynomolgus, and rhesus), model development progress has been most promising in pigtail macaques (PTMs), which do not express an HIV-1-restricting tripartite motif-containing protein 5 α (TRIM5α). In our work, we have derived a CCR5-tropic mmHIV-1 clone designated stHIV-A19 that comprises 94% HIV-1 genome sequence and replicates to high acute-phase titers in PTMs. In animals treated with a cell-depleting CD8α antibody at the time of infection, stHIV-A19 maintains chronically elevated plasma viral loads with progressive CD4+ T-cell loss and the development of acquired immune-deficiency syndrome (AIDS)-defining clinical endpoints. However, in the absence of CD8α+ cell depletion, no mmHIV-1 model has yet displayed high levels of chronic viremia or AIDS-like pathogenesis. Here, we review mmHIV-1 development approaches, the phenotypes, features, limitations, and potential utility of currently available mmHIV-1s, and propose future directions to further advance these models.

Nef defect attenuates HIV viremia and immune dysregulation in the bone marrow-liver-thymus-spleen (BLTS) humanized mouse model

In vitro studies have shown that deletion of nef and deleterious mutation in the Nef dimerization interface attenuates HIV replication and associated pathogenesis. Humanized rodents with human immune cells and lymphoid tissues are robust in vivo models for investigating the interactions between HIV and the human immune system. Here, we demonstrate that nef deletion impairs HIV replication and HIV-induced immune dysregulation in the blood and human secondary lymphoid tissue (human spleen) in bone marrow-liver-thymus-spleen (BLTS) humanized mice. Furthermore, we also show that nef defects (via deleterious mutations in the dimerization interface) impair HIV replication and HIV-induced immune dysregulation in the blood and human spleen in BLTS-humanized mice. We demonstrate that the reduced replication of nef-deleted and nef-defective HIV is associated with robust antiviral innate immune response, and T helper 1 response. Our results support the proposition that Nef may be a therapeutic target for adjuvants in HIV cure strategies.

Mapping the structural landscape of the yeast Ty3 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons are significant contributors to the evolution and diversity of eukaryotic genomes. Their RNA genomes (gRNA) serve as a template for protein synthesis and reverse transcription to a DNA copy, which can integrate into the host genome. Here, we used the SHAPE-MaP strategy to explore Ty3 retrotransposon gRNA structure in yeast and under cell-free conditions. Our study reveals the structural dynamics of Ty3 gRNA and the well-folded core, formed independently of the cellular environment. Based on the detailed map of Ty3 gRNA structure, we characterized the structural context of cis-acting sequences involved in reverse transcription and frameshifting. We also identified a novel functional sequence as a potential initiator for Ty3 gRNA dimerization. Our data indicate that the dimer is maintained by direct interaction between short palindromic sequences at the 5' ends of the two Ty3 gRNAs, resembling the model characteristic for other retroelements like HIV-1 and Ty1. This work points out a range of cell-dependent and -independent Ty3 gRNA structural changes that provide a solid background for studies on RNA structure-function relationships important for retroelement biology.

Toward a Categorization of Virus-ncRNA Interactions in the World of RNA to Disentangle the Tiny Secrets of Dengue Virus

In recent years, the function of noncoding RNAs (ncRNAs) as regulatory molecules of cell physiology has begun to be better understood. Advances in viral molecular biology have shown that host ncRNAs, cellular factors, and virus-derived ncRNAs and their interplay are strongly disturbed during viral infections. Nevertheless, the folding of RNA virus genomes has also been identified as a critical factor in regulating canonical and non-canonical functions. Due to the influence of host ncRNAs and the structure of RNA viral genomes, complex molecular and cellular processes in infections are modulated. We propose three main categories to organize the current information about RNA-RNA interactions in some well-known human viruses. The first category shows examples of host ncRNAs associated with the immune response triggered in viral infections. Even though miRNAs introduce a standpoint, they are briefly presented to keep researchers moving forward in uncovering other RNAs. The second category outlines interactions between virus-host ncRNAs, while the third describes how the structure of the RNA viral genome serves as a scaffold for processing virus-derived RNAs. Our grouping may provide a comprehensive framework to classify ncRNA-host-cell interactions for emerging viruses and diseases. In this sense, we introduced them to organize DENV-host-cell interactions.

Involvement of a Rarely Used Splicing SD2b Site in the Regulation of HIV-1 vif mRNA Production as Revealed by a Growth-Adaptive Mutation

We have previously reported an HIV-1 mutant designated NL-Y226tac that expresses Vif at an ultra-low level, being replication-defective in high-APOBEC3G cells, such as H9. It carries a synonymous mutation within the splicing SA1 site relative to its parental clone. In order to determine whether a certain mutant(s) emerges during multi-infection cycles, we maintained H9 cells infected with a relatively low or high input of NL-Y226tac for extended time periods. Unexpectedly, we reproducibly identified a g5061a mutation in the SD2b site in the two independent long-term culture experiments that partially increases Vif expression and replication ability. Importantly, the adaptive mutation g5061a was demonstrated to enhance vif mRNA production by activation of the SA1 site mediated through increasing usage of a rarely used SD2b site. In the long-term culture initiated by a high virus input, we additionally found a Y226Fttc mutation at the original Y226tac site in SA1 that fully restores Vif expression and replication ability. As expected, the adaptive mutation Y226Fttc enhances vif mRNA production through increasing the splicing site usage of SA1. Our results here revealed the importance of the SD2b nucleotide sequence in producing vif mRNA involved in the HIV-1 adaptation and of mutual antagonism between Vif and APOBEC3 proteins in HIV-1 adaptation/evolution and survival.

Structural diversity of the region encompassing DIS, SD and Psi hairpins in HIV and SIV genomes

We investigated in silico the secondary structure of the region encompassing DIS, SD and Psi hairpins in HIV-1 genomes of rare groups N, O and P, HIV-2 genomes and SIV genomes from chimpanzees, gorillas and monkeys. We found that the structure of this region in SIVcpzptt genomes of the 1st and the 2nd clusters is similar to that in HIV-1 genomes of groups M and N, respectively. Further, the structure of the region encompassing DIS, SD and Psi hairpins is similar in HIV-1 genomes of groups O and P and SIVgor genomes. Here we report that the DIS hairpin and truncated Psi hairpin are conserved in all HIV-1 and SIVcpz/gor genomes studied, while only the sequence of the splice donor site, but not the architecture of the SD hairpin involving this signal is conserved in HIV-1N/O/P and SIVcpz/gor genomes. A study on the 5' leader structure in genomes of 28 different SIV lineages infecting monkeys showed that the domain closed by U5-AUG duplex can form in all these genomes. This domain mainly consists of 2 subdomains, one of which includes the signal PBS (PBS subdomain) and another contains a putative DIS hairpin (DIS subdomain). DIS subdomains contain 1-8 hairpins. None of them is similar to those in HIV-1 or SIVcpz/gor genomes. The palindrome GUGCAC was found only in SIVdrl/mnd-2, the GACGC-GCGUC duplex (Sakuragi et al., 2012) - only in SIVrcm/drl/mnd-2 and a putative 5' G-quadruplex - in SIVdeb/drl/rcm/stm genomes. In genomes of eight SIV lineages, DIS hairpin has palindrome UGCGCA. Studies on the 5' leader in 64 HIV-2 genomes of different subtypes showed, in particular, that this region has sequences of a putative 5' G-quadruplex and a putative duplex similar to the GACGC-GCGUC duplex. The secondary structures of the region encompassing DIS, SD and Psi hairpins in HIV-2 genomes of subtype B and recombinant 01_AB are similar and differ from that in genomes of subtype A.

Structure-first identification of RNA elements that regulate dengue virus genome architecture and replication

The genomes of RNA viruses encode the information required for replication in host cells both in their linear sequence and in complex higher-order structures. A subset of these RNA genome structures show clear sequence conservation, and have been extensively described for well-characterized viruses. However, the extent to which viral RNA genomes contain functional structural elements-unable to be detected by sequence alone-that nonetheless are critical to viral fitness is largely unknown. Here, we devise a structure-first experimental strategy and use it to identify 22 structure-similar motifs across the coding sequences of the RNA genomes for the four dengue virus serotypes. At least 10 of these motifs modulate viral fitness, revealing a significant unnoticed extent of RNA structure-mediated regulation within viral coding sequences. These viral RNA structures promote a compact global genome architecture, interact with proteins, and regulate the viral replication cycle. These motifs are also thus constrained at the levels of both RNA structure and protein sequence and are potential resistance-refractory targets for antivirals and live-attenuated vaccines. Structure-first identification of conserved RNA structure enables efficient discovery of pervasive RNA-mediated regulation in viral genomes and, likely, other cellular RNAs.

A potential long-range RNA-RNA interaction in the HIV-1 RNA

It is well-established that viral and cellular mRNAs alike harbour functional long-range intra-molecular RNA-RNA interactions. Despite the biological importance of such interactions, their identification and characterization remain challenging. Here we present a computational method for the identification of certain kinds of long-range intra-molecular RNA-RNA interactions involving the loop nucleotides of a hairpin loop. Using the computational method, we analysed 4272 HIV-1 genomic mRNAs. A potential long-range intra-molecular RNA-RNA interaction within the HIV-1 genomic RNA was identified. The long-range interaction is mediated by a kissing loop structure between two stem-loops of the previously reported SHAPE-based secondary structure of the entire HIV-1 genome. Structural modelling studies were carried out to show that the kissing loop structure not only is sterically feasible, but also contains a conserved RNA structural motif often found in compact RNA pseudoknots. The computational method should be generally applicable to the identification of potential long-range intra-molecular RNA-RNA interactions in any viral or cellular mRNA sequence.Communicated by Ramaswamy H. Sarma.

Genetic stability of SIV Gag/Tat gene inserted into Del-II in modified vaccinia virus ankara after serial passage of recombinant vector in pCEFs cells

Modified vaccinia virus Ankara (MVA) is an attenuated vaccinia virus with restricted replication in human cells. The virus serves as an ideal vaccine vector suitable for safe use even in immune-compromised individuals. With its inherently large packaging capacity, expression cassettes encoding bulky genes can be inserted into deletion regions within the MVA genome. These deletion sites develop during the process of the attenuation of the virus by passage in Chicken Embryo Fibroblasts (pCEFs). Transgene stability in MVA is important to assure immunogenicity and efficacy. In the present study, we assessed the effect of substantial passage of recombinant MVA vectors on the stability of expression cassette encoding SIV Gag/Tat genes inserted at the Del-II site, as part of generating a vaccine to protect from HIV. Our data indicated that after 15 passages there was a significant loss or mutation of the inserted genes.

Genome sequence analysis suggests coevolution of the DIS, SD, and Psi hairpins in HIV-1 genomes

HIV-1 RNA dimerization is a critical step in viral life cycle. It is a prerequisite for genome packaging and plays an important role in reverse transcription and recombination. Dimerization is promoted by the DIS (dimerization initiation site) hairpin located in the 5' leader of HIV-1 genome. Despite the high genetic diversity in HIV-1 group M, only five apical loops (AAGCGCGCA, AAGUGCGCA, AAGUGCACA, AGGUGCACA and AGUGCAC) are commonly found in DIS hairpins. We refer to the parent DISes with these apical loops as DIS_Lai, DIS_Trans, DIS_F, DIS_Mal, and DIS_C, respectively. Based on identity or similarity of DIS hairpins to parent DISes, we distributed HIV-1 M genomes into five dimerization groups. Comparison of the primary and secondary structures of DIS, SD and Psi hairpins in about 3000 HIV-1 M genomes showed that the mutation frequencies at particular nucleotide positions of these hairpins differ among the dimerization groups, and DIS_F may be an origin of other parent DISes. We found that DIS, SD and Psi hairpins have hundreds of variants, only some of them occurring rather frequently. The lower part of DIS hairpin with G x AGG internal loop is highly conserved in both HIV-1 and SIV genomes. We supposed that the G-quadruplex, located 56 nts downstream of the Gag start codon, may participate in switching of HIV-1 leader RNA from BMH (branched multiple hairpins) to LDI (long distance interaction) conformation.

Animal models for studies of HIV-1 brain reservoirs

The HIV-1 often evades a robust antiretroviral-mediated immune response, leading to persistent infection within anatomically privileged sites including the CNS. Continuous low-level infection occurs in the presence of effective antiretroviral therapy (ART) in CD4+ T cells and mononuclear phagocytes (MP; monocytes, macrophages, microglia, and dendritic cells). Within the CNS, productive viral infection is found exclusively in microglia and meningeal, perivascular, and choroidal macrophages. MPs serve as the principal viral CNS reservoir. Animal models have been developed to recapitulate natural human HIV-1 infection. These include nonhuman primates, humanized mice, EcoHIV, and transgenic rodent models. These models have been used to study disease pathobiology, antiretroviral and immune modulatory agents, viral reservoirs, and eradication strategies. However, each of these models are limited to specific component(s) of human disease. Indeed, HIV-1 species specificity must drive therapeutic and cure studies. These have been studied in several model systems reflective of latent infections, specifically in MP (myeloid, monocyte, macrophages, microglia, and histiocyte cell) populations. Therefore, additional small animal models that allow productive viral replication to enable viral carriage into the brain and the virus-susceptible MPs are needed. To this end, this review serves to outline animal models currently available to study myeloid brain reservoirs and highlight areas that are lacking and require future research to more effectively study disease-specific events that could be useful for viral eradication studies both in and outside the CNS.

Computational and Experimental Approaches to Study the RNA Secondary Structures of RNA Viruses

Most pandemics of recent decades can be traced to RNA viruses, including HIV, SARS, influenza, dengue, Zika, and SARS-CoV-2. These RNA viruses impose considerable social and economic burdens on our society, resulting in a high number of deaths and high treatment costs. As these RNA viruses utilize an RNA genome, which is important for different stages of the viral life cycle, including replication, translation, and packaging, studying how the genome folds is important to understand virus function. In this review, we summarize recent advances in computational and high-throughput RNA structure-mapping approaches and their use in understanding structures within RNA virus genomes. In particular, we focus on the genome structures of the dengue, Zika, and SARS-CoV-2 viruses due to recent significant outbreaks of these viruses around the world.

A Stretch of Unpaired Purines in the Leader Region of Simian Immunodeficiency Virus (SIV) Genomic RNA is Critical for its Packaging into Virions

Simian immunodeficiency virus (SIV) is an important lentivirus used as a non-human primate model to study HIV replication, and pathogenesis of human AIDS, as well as a potential vector for human gene therapy. This study investigated the role of single-stranded purines (ssPurines) as potential genomic RNA (gRNA) packaging determinants in SIV replication. Similar ssPurines have been implicated as important motifs for gRNA packaging in many retroviruses like, HIV-1, MPMV, and MMTV by serving as Gag binding sites during virion assembly. In examining the secondary structure of the SIV 5' leader region, as recently deduced using SHAPE methodology, we identified four specific stretches of ssPurines (I-IV) in the region that harbors major packaging determinants of SIV. The significance of these ssPurine motifs were investigated by mutational analysis coupled with a biologically relevant single round of replication assay. These analyses revealed that while ssPurine II was essential, the others (ssPurines I, III, & IV) did not significantly contribute to SIV gRNA packaging. Any mutation in the ssPurine II, such as its deletion or substitution, or other mutations that caused base pairing of ssPurine II loop resulted in near abrogation of RNA packaging, further substantiating the crucial role of ssPurine II and its looped conformation in SIV gRNA packaging. Structure prediction analysis of these mutants further corroborated the biological results and further revealed that the unpaired nature of ssPurine II is critical for its function during SIV RNA packaging perhaps by enabling it to function as a specific binding site for SIV Gag.

The Expression Level of HIV-1 Vif Is Optimized by Nucleotide Changes in the Genomic SA1D2prox Region during the Viral Adaptation Process

HIV-1 Vif plays an essential role in viral replication by antagonizing anti-viral cellular restriction factors, a family of APOBEC3 proteins. We have previously shown that naturally-occurring single-nucleotide mutations in the SA1D2prox region, which surrounds the splicing acceptor 1 and splicing donor 2 sites of the HIV-1 genome, dramatically alter the Vif expression level, resulting in variants with low or excessive Vif expression. In this study, we investigated how these HIV-1 variants with poor replication ability adapt and evolve under the pressure of APOBEC3 proteins. Adapted clones obtained through adaptation experiments exhibited an altered replication ability and Vif expression level compared to each parental clone. While various mutations were present throughout the viral genome, all replication-competent adapted clones with altered Vif expression levels were found to bear them within SA1D2prox, without exception. Indeed, the mutations identified within SA1D2prox were responsible for changes in the Vif expression levels and altered the splicing pattern. Moreover, for samples collected from HIV-1-infected patients, we showed that the nucleotide sequences of SA1D2prox can be chronologically changed and concomitantly affect the Vif expression levels. Taken together, these results demonstrated the importance of the SA1D2prox nucleotide sequence for modulating the Vif expression level during HIV-1 replication and adaptation.

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.

High-throughput interrogation of programmed ribosomal frameshifting in human cells

Programmed ribosomal frameshifting (PRF) is the controlled slippage of the translating ribosome to an alternative frame. This process is widely employed by human viruses such as HIV and SARS coronavirus and is critical for their replication. Here, we developed a high-throughput approach to assess the frameshifting potential of a sequence. We designed and tested >12,000 sequences based on 15 viral and human PRF events, allowing us to systematically dissect the rules governing ribosomal frameshifting and discover novel regulatory inputs based on amino acid properties and tRNA availability. We assessed the natural variation in HIV gag-pol frameshifting rates by testing >500 clinical isolates and identified subtype-specific differences and associations between viral load in patients and the optimality of PRF rates. We devised computational models that accurately predict frameshifting potential and frameshifting rates, including subtle differences between HIV isolates. This approach can contribute to the development of antiviral agents targeting PRF.

Evolutionary Analyses of Base-Pairing Interactions in DNA and RNA Secondary Structures

Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.

Expression Level of HIV-1 Vif Can Be Fluctuated by Natural Nucleotide Variations in the vif-Coding and Regulatory SA1D2prox Sequences of the Proviral Genome

Vif is required for HIV-1 replication in natural target cells by counteracting host restriction factors, APOBEC3 (A3) proteins. We recently demonstrated that Vif expression level can be changed by naturally occurring single-nucleotide variations within SA1D2prox of the HIV-1 genome. We also found that levels for vif/vpr mRNAs are inversely correlated. While amino acid sequence per se is critical for functionality, Vif expression level modulated by signal sequences in its coding region is likely to be important as well. There are two splicing sites in the region involved in vpr expression. To reveal possible fluctuations of Vif-expression level, we examined SA1D2prox and vif gene by chimeric approaches using HIV-1 subtypes B and C with distinct anti-A3 activity. In this report, recombinant clones in subtype B backbone carrying chimeric sequences with respect to SA1D2prox/vif and those within the vif-coding region were generated. Of these, clones containing vif-coding sequence of subtype C, especially its 3′ region, expressed vif/Vif at a decreased level but did at an increased level for vpr/Vpr. Clones with reduced vif/Vif level grew similarly or slightly better than a parental clone in weakly A3G-positive cells but more poorly in highly A3G-expressing cells. Three clones with this property were also tested for their A3-degrading activity. One of the clones appeared to have some defect in addition to the poor ability to express vif/Vif. Taken all together, our results show that natural variations in the SA1D2prox and vif-coding region can change the Vif-expression level and affect the HIV-1 replication potential.

Physical and Functional Analysis of Viral RNA Genomes by SHAPE

RNA viruses encode the information required to usurp cellular metabolism and gene regulation and to enable their own replication in two ways: in the linear sequence of their RNA genomes and in higher-order structures that form when the genomic RNA strand folds back on itself. Application of high-resolution SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) structure probing to viral RNA genomes has identified numerous new regulatory elements, defined new principles by which viral RNAs interact with the cellular host and evade host immune responses, and revealed relationships between virus evolution and RNA structure. This review summarizes our current understanding of genome structure-function interrelationships for RNA viruses, as informed by SHAPE structure probing, and outlines opportunities for future studies.

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.

On the generation of the MSD-Ѱ class of defective HIV proviruses

Antiretroviral therapy (ART) can effectively suppress ongoing HIV replication and block disease progression, but the infection is never cured due to the persistence of a small pool of latently infected cells hosting integrated replication-competent HIV proviruses. However, the vast majority of HIV proviruses in ART-treated patients are replication-incompetent due to a variety of genetic defects. Most defective proviruses (around 90%) contain large internal deletions or are G-to-A hypermutated, resulting in destruction of most if not all viral open reading frames, which is consistent with the idea that cytotoxic T cells (CTLs) effectively remove cells that produce viral antigens. An intriguing subclass of defective proviruses (around 10%) that are consistently detected in such patients carry a small deletion or a point mutation in a relatively precise and well conserved region near the 5ʹ end of the HIV genome, in the area that encodes the major splice donor (MSD) site and the packaging signal Ѱ in the viral RNA genome. Why this subclass of proviruses is defective has never been properly understood. We now propose a mechanistic scenario for how these MSD-Ѱ mutations can prevent viral protein expression. Based on ample results in literature, we argue that MSD inactivation triggers the activity of the 5ʹ-polyadenylation site, resulting in the production of ultra-short non-protein-coding HIV transcripts.

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts

The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5′ untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5′ UTR. As expected, mutations in double-stranded (DS) regions of the 5′ UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5′ cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.

IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Guidelines for SHAPE Reagent Choice and Detection Strategy for RNA Structure Probing Studies

Chemical probing is an important tool for characterizing the complex folded structures of RNA molecules, many of which play key cellular roles. Electrophilic SHAPE reagents create adducts at the 2'-hydroxyl position on the RNA backbone of flexible ribonucleotides with relatively little dependence on nucleotide identity. Strategies for adduct detection such as mutational profiling (MaP) allow accurate, automated calculation of relative adduct frequencies for each nucleotide in a given RNA or group of RNAs. A number of alternative reagents and adduct detection strategies have been proposed, especially for use in living cells. Here we evaluate five SHAPE reagents: three previously well-validated reagents 1M7 (1-methyl-7-nitroisatoic anhydride), 1M6 (1-methyl-6-nitroisatoic anhydride), and NMIA ( N-methylisatoic anhydride), one more recently proposed NAI (2-methylnicotinic acid imidazolide), and one novel reagent 5NIA (5-nitroisatoic anhydride). We clarify the importance of carefully designed software in reading out SHAPE experiments using massively parallel sequencing approaches. We examine SHAPE modification in living cells in diverse cell lines, compare MaP and reverse transcription-truncation as SHAPE adduct detection strategies, make recommendations for SHAPE reagent choice, and outline areas for future development.

Relationship between genome packaging and Gag translation/AUG of primate lentiviruses

About the relationship between retroviral genome packaging and translation, three possible modes (random-, trans-, and cis-) of packaging process could be assumed. In this report, we developed an assay system based on the RT-qPCR to measure the packaging efficiency of primate lentiviruses. With this system, we analyzed the genome packaging modes of primate lentiviruses such as HIV-1, 2, SIVmac and SIVagm. The data suggested that the modes of all viruses analyzed were very similar. In addition, we observed that the Gag-AUG sequences of them played important roles for maintaining efficient packaging, other than the initiation of translation.

Adeno-associated virus gene delivery of broadly neutralizing antibodies as prevention and therapy against HIV-1

Vectored gene delivery of HIV-1 broadly neutralizing antibodies (bNAbs) using recombinant adeno-associated virus (rAAV) is a promising alternative to conventional vaccines for preventing new HIV-1 infections and for therapeutically suppressing established HIV-1 infections. Passive infusion of single bNAbs has already shown promise in initial clinical trials to temporarily decrease HIV-1 load in viremic patients, and to delay viral rebound from latent reservoirs in suppressed patients during analytical treatment interruptions of antiretroviral therapy. Long-term, continuous, systemic expression of such bNAbs could be achieved with a single injection of rAAV encoding antibody genes into muscle tissue, which would bypass the challenges of eliciting such bNAbs through traditional vaccination in naïve patients, and of life-long repeated passive transfers of such biologics for therapy. rAAV delivery of single bNAbs has already demonstrated protection from repeated HIV-1 vaginal challenge in humanized mouse models, and phase I clinical trials of this approach are underway. Selection of which individual, or combination of, bNAbs to deliver to counter pre-existing resistance and the rise of escape mutations in the virus remains a challenge, and such choices may differ depending on use of this technology for prevention versus therapy.

N6-Methyladenosine-binding proteins suppress HIV-1 infectivity and viral production

The internal N6-methyladenosine (m6A) modification of cellular mRNA regulates post-transcriptional gene expression. The YTH domain family proteins (YTHDF1-3 or Y1-3) bind to m6A-modified cellular mRNAs and modulate their metabolism and processing, thereby affecting cellular protein translation. We previously reported that HIV-1 RNA contains the m6A modification and that Y1-3 proteins inhibit HIV-1 infection by decreasing HIV-1 reverse transcription activity. Here, we investigated the mechanisms of Y1-3-mediated inhibition of HIV-1 infection in target cells and the effect of Y1-3 on viral production levels in virus-producing cells. We found that Y1-3 protein overexpression in HIV-1 target cells decreases viral genomic RNA (gRNA) levels and inhibits both early and late reverse transcription. Purified recombinant Y1-3 proteins preferentially bound to the m6A-modified 5' leader sequence of gRNA compared with its unmodified RNA counterpart, consistent with the strong binding of Y1-3 proteins to HIV-1 gRNA in infected cells. HIV-1 mutants with two altered m6A modification sites in the 5' leader sequence of gRNA exhibited significantly lower infectivity than WT, replication-competent HIV-1, confirming that these sites alter viral infection. HIV-1 produced from cells in which endogenous Y1, Y3, or Y1-3 proteins were knocked down singly or together had increased viral infectivity compared with HIV-1 produced in control cells. Interestingly, we found that Y1-3 proteins and HIV-1 Gag protein formed a complex with RNA in HIV-1-producing cells. Overall, these results indicate that Y1-3 proteins inhibit HIV-1 infection and provide new insights into the mechanisms by which the m6A modification of HIV-1 RNA affects viral replication.

Neutral Theory and Rapidly Evolving Viral Pathogens

The evolution of viral pathogens is shaped by strong selective forces that are exerted during jumps to new hosts, confrontations with host immune responses and antiviral drugs, and numerous other processes. However, while undeniably strong and frequent, adaptive evolution is largely confined to small parts of information-packed viral genomes, and the majority of observed variation is effectively neutral. The predictions and implications of the neutral theory have proven immensely useful in this context, with applications spanning understanding within-host population structure, tracing the origins and spread of viral pathogens, predicting evolutionary dynamics, and modeling the emergence of drug resistance. We highlight the multiple ways in which the neutral theory has had an impact, which has been accelerated in the age of high-throughput, high-resolution genomics.

Structural divergence creates new functional features in alphavirus genomes

Alphaviruses are mosquito-borne pathogens that cause human diseases ranging from debilitating arthritis to lethal encephalitis. Studies with Sindbis virus (SINV), which causes fever, rash, and arthralgia in humans, and Venezuelan equine encephalitis virus (VEEV), which causes encephalitis, have identified RNA structural elements that play key roles in replication and pathogenesis. However, a complete genomic structural profile has not been established for these viruses. We used the structural probing technique SHAPE-MaP to identify structured elements within the SINV and VEEV genomes. Our SHAPE-directed structural models recapitulate known RNA structures, while also identifying novel structural elements, including a new functional element in the nsP1 region of SINV whose disruption causes a defect in infectivity. Although RNA structural elements are important for multiple aspects of alphavirus biology, we found the majority of RNA structures were not conserved between SINV and VEEV. Our data suggest that alphavirus RNA genomes are highly divergent structurally despite similar genomic architecture and sequence conservation; still, RNA structural elements are critical to the viral life cycle. These findings reframe traditional assumptions about RNA structure and evolution: rather than structures being conserved, alphaviruses frequently evolve new structures that may shape interactions with host immune systems or co-evolve with viral proteins.

Structural analyses of NEAT1 lncRNAs suggest long-range RNA interactions that may contribute to paraspeckle architecture

Paraspeckles are nuclear bodies that regulate multiple aspects of gene expression. The long non-coding RNA (lncRNA) NEAT1 is essential for paraspeckle formation. NEAT1 has a highly ordered spatial organization within the paraspeckle, such that its 5' and 3' ends localize on the periphery of paraspeckle, while central sequences of NEAT1 are found within the paraspeckle core. As such, the structure of NEAT1 RNA may be important as a scaffold for the paraspeckle. In this study, we used SHAPE probing and computational analyses to investigate the secondary structure of human and mouse NEAT1. We propose a secondary structural model of the shorter (3,735 nt) isoform hNEAT1_S, in which the RNA folds into four separate domains. The secondary structures of mouse and human NEAT1 are largely different, with the exception of several short regions that have high structural similarity. Long-range base-pairing interactions between the 5' and 3' ends of the long isoform NEAT1 (NEAT1_L) were predicted computationally and verified using an in vitro RNA-RNA interaction assay. These results suggest that the conserved role of NEAT1 as a paraspeckle scaffold does not require extensively conserved RNA secondary structure and that long-range interactions among NEAT1 transcripts may have an important architectural function in paraspeckle formation.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Global synonymous mutagenesis identifies cis-acting RNA elements that regulate HIV-1 splicing and replication

The ~9.5 kilobase HIV-1 genome contains RNA sequences and structures that control many aspects of viral replication, including transcription, splicing, nuclear export, translation, packaging and reverse transcription. Nonetheless, chemical probing and other approaches suggest that the HIV-1 genome may contain many more RNA secondary structures of unknown importance and function. To determine whether there are additional, undiscovered cis-acting RNA elements in the HIV-1 genome that are important for viral replication, we undertook a global silent mutagenesis experiment. Sixteen mutant proviruses containing clusters of ~50 to ~200 synonymous mutations covering nearly the entire HIV-1 protein coding sequence were designed and synthesized. Analyses of these mutant viruses resulted in their division into three phenotypic groups. Group 1 mutants exhibited near wild-type replication, Group 2 mutants exhibited replication defects accompanied by perturbed RNA splicing, and Group 3 mutants had replication defects in the absence of obvious splicing perturbation. The three phenotypes were caused by mutations that exhibited a clear regional bias in their distribution along the viral genome, and those that caused replication defects all caused reductions in the level of unspliced RNA. We characterized in detail the underlying defects for Group 2 mutants. Second-site revertants that enabled viral replication could be derived for Group 2 mutants, and generally contained point mutations that reduced the utilization of proximal splice sites. Mapping of the changes responsible for splicing perturbations in Group 2 viruses revealed the presence of several RNA sequences that apparently suppressed the use of cryptic or canonical splice sites. Some sequences that affected splicing were diffusely distributed, while others could be mapped to discrete elements, proximal or distal to the affected splice site(s). Overall, our data indicate complex negative regulation of HIV-1 splicing by RNA elements in various regions of the HIV-1 genome that enable balanced splicing and viral replication.

In addition to encoding viral proteins, the HIV-1 genome contains sequence elements that act at the level of RNA to enable replication. We undertook an experiment to discover new RNA elements that act in this way by altering nearly the entire coding sequence of the viral genome so as to change the RNA sequence without changing protein sequences. This experiment uncovered two classes of defective mutants. One class had profound defects in RNA splicing, the other had no obvious defects in splicing. Through an analysis of the splicing-defective mutants, we found several previously RNA sequences in the viral genome that affected splicing, enabling a nearly complete catalogue of signals that regulate HIV-1 alternative splicing in infected cells to be derived. Because these newly described sequences lack sequence motifs that are known to bind to canonical splicing-regulatory proteins, they may function through novel mechanisms.

New windows into retroviral RNA structures

Background

The multiple roles of both viral and cellular RNAs have become increasingly apparent in recent years, and techniques to model them have become significantly more powerful, enabling faster and more accurate visualization of RNA structures.

Main body

Techniques such as SHAPE (selective 2’OH acylation analysed by primer extension) have revolutionized the field, and have been used to examine RNAs belonging to many and diverse retroviruses. Secondary structure probing reagents such as these have been aided by the development of faster methods of analysis either via capillary or next-generation sequencing, allowing the analysis of entire genomes, and of retroviral RNA structures within virions. Techniques to model the three-dimensional structures of these large RNAs have also recently developed.

Conclusions

The flexibility of retroviral RNAs, both structural and functional, is clear from the results of these new experimental techniques. Retroviral RNA structures and structural changes control many stages of the lifecycle, and both the RNA structures themselves and their interactions with ligands are potential new drug targets. In addition, our growing understanding of retroviral RNA structures is aiding our knowledge of cellular RNA form and function.

Probing the Structures of Viral RNA Regulatory Elements with SHAPE and Related Methodologies

Viral RNAs were selected by evolution to possess maximum functionality in a minimal sequence. Depending on the classification of the virus and the type of RNA in question, viral RNAs must alternately be replicated, spliced, transcribed, transported from the nucleus into the cytoplasm, translated and/or packaged into nascent virions, and in most cases, provide the sequence and structural determinants to facilitate these processes. One consequence of this compact multifunctionality is that viral RNA structures can be exquisitely complex, often involving intermolecular interactions with RNA or protein, intramolecular interactions between sequence segments separated by several thousands of nucleotides, or specialized motifs such as pseudoknots or kissing loops. The fluidity of viral RNA structure can also present a challenge when attempting to characterize it, as genomic RNAs especially are likely to sample numerous conformations at various stages of the virus life cycle. Here we review advances in chemoenzymatic structure probing that have made it possible to address such challenges with respect to cis-acting elements, full-length viral genomes and long non-coding RNAs that play a major role in regulating viral gene expression.

Evaluation of global HIV/SIV envelope gp120 RNA structure and evolution within and among infected hosts

Lentiviral RNA genomes contain structural elements that play critical roles in viral replication. Although structural features of 5'-untranslated regions have been well characterized, attempts to identify important structures in other genomic regions by Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) have led to conflicting structural and mechanistic conclusions. Previous approaches accounted neither for sequence heterogeneity that is ubiquitous in viral populations, nor for selective constraints operating at the protein level. We developed an approach that augments SHAPE with phylogenetic analyses and applied it to investigate structure in coding regions (cRNA) within the HIV and SIV envelope genes. Analysis of single-genome SHAPE data with phylogenetic information from diverse lentiviral sequences argues against the conservation of a putative global gp120 RNA structure but points to the existence of core RNA sub-structures. Our findings establish a framework for considering sequence heterogeneity and protein function in de novo RNA structure inference approaches.

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.

Increasing the CpG dinucleotide abundance in the HIV-1 genomic RNA inhibits viral replication

BACKGROUND

The human immunodeficiency virus type 1 (HIV-1) structural protein Gag is necessary and sufficient to form viral particles. In addition to encoding the amino acid sequence for Gag, the underlying RNA sequence could encode cis-acting elements or nucleotide biases that are necessary for viral replication. Furthermore, RNA sequences that inhibit viral replication could be suppressed in gag. However, the functional relevance of RNA elements and nucleotide biases that promote or repress HIV-1 replication remain poorly understood.

RESULTS

To characterize if the RNA sequence in gag controls HIV-1 replication, the matrix (MA) region was codon modified, allowing the RNA sequence to be altered without affecting the protein sequence. Codon modification of nucleotides (nt) 22-261 or 22-378 in gag inhibited viral replication by decreasing genomic RNA (gRNA) abundance, gRNA stability, Gag expression, virion production and infectivity. Comparing the effect of these point mutations to deletions of the same region revealed that the mutations inhibited infectious virus production while the deletions did not. This demonstrated that codon modification introduced inhibitory sequences. There is a much lower than expected frequency of CpG dinucleotides in HIV-1 and codon modification introduced a substantial increase in CpG abundance. To determine if they are necessary for inhibition of HIV-1 replication, codons introducing CpG dinucleotides were mutated back to the wild type codon, which restored efficient Gag expression and infectious virion production. To determine if they are sufficient to inhibit viral replication, CpG dinucleotides were inserted into gag in the absence of other changes. The increased CpG dinucleotide content decreased HIV-1 infectivity and viral replication.

CONCLUSIONS

The HIV-1 RNA sequence contains low abundance of CpG dinucleotides. Increasing the abundance of CpG dinucleotides inhibits multiple steps of the viral life cycle, providing a functional explanation for why CpG dinucleotides are suppressed in HIV-1.

RNA structure inference through chemical mapping after accidental or intentional mutations

Despite the critical roles RNA structures play in regulating gene expression, sequencing-based methods for experimentally determining RNA base pairs have remained inaccurate. Here, we describe a multidimensional chemical-mapping method called "mutate-and-map read out through next-generation sequencing" (M2-seq) that takes advantage of sparsely mutated nucleotides to induce structural perturbations at partner nucleotides and then detects these events through dimethyl sulfate (DMS) probing and mutational profiling. In special cases, fortuitous errors introduced during DNA template preparation and RNA transcription are sufficient to give M2-seq helix signatures; these signals were previously overlooked or mistaken for correlated double-DMS events. When mutations are enhanced through error-prone PCR, in vitro M2-seq experimentally resolves 33 of 68 helices in diverse structured RNAs including ribozyme domains, riboswitch aptamers, and viral RNA domains with a single false positive. These inferences do not require energy minimization algorithms and can be made by either direct visual inspection or by a neural-network-inspired algorithm called M2-net. Measurements on the P4-P6 domain of the Tetrahymena group I ribozyme embedded in Xenopus egg extract demonstrate the ability of M2-seq to detect RNA helices in a complex biological environment.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

Comparative and integrative analysis of RNA structural profiling data: current practices and emerging questions

BACKGROUND

Structure profiling experiments provide single-nucleotide information on RNA structure. Recent advances in chemistry combined with application of high-throughput sequencing have enabled structure profiling at transcriptome scale and in living cells, creating unprecedented opportunities for RNA biology. Propelled by these experimental advances, massive data with ever-increasing diversity and complexity have been generated, which give rise to new challenges in interpreting and analyzing these data.

RESULTS

We review current practices in analysis of structure profiling data with emphasis on comparative and integrative analysis as well as highlight emerging questions. Comparative analysis has revealed structural patterns across transcriptomes and has become an integral component of recent profiling studies. Additionally, profiling data can be integrated into traditional structure prediction algorithms to improve prediction accuracy.

CONCLUSIONS

To keep pace with experimental developments, methods to facilitate, enhance and refine such analyses are needed. Parallel advances in analysis methodology will complement profiling technologies and help them reach their full potential.

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.

In vivo mutation rates and the landscape of fitness costs of HIV-1

Mutation rates and fitness costs of deleterious mutations are difficult to measure in vivo but essential for a quantitative understanding of evolution. Using whole genome deep sequencing data from longitudinal samples during untreated HIV-1 infection, we estimated mutation rates and fitness costs in HIV-1 from the dynamics of genetic variation. At approximately neutral sites, mutations accumulate with a rate of 1.2 × 10-5 per site per day, in agreement with the rate measured in cell cultures. We estimated the rate from G to A to be the largest, followed by the other transitions C to T, T to C, and A to G, while transversions are less frequent. At other sites, mutations tend to reduce virus replication. We estimated the fitness cost of mutations at every site in the HIV-1 genome using a model of mutation selection balance. About half of all non-synonymous mutations have large fitness costs (>10 percent), while most synonymous mutations have costs <1 percent. The cost of synonymous mutations is especially low in most of pol where we could not detect measurable costs for the majority of synonymous mutations. In contrast, we find high costs for synonymous mutations in important RNA structures and regulatory regions. The intra-patient fitness cost estimates are consistent across multiple patients, indicating that the deleterious part of the fitness landscape is universal and explains a large fraction of global HIV-1 group M diversity.

Mechanistic differences between HIV-1 and SIV nucleocapsid proteins and cross-species HIV-1 genomic RNA recognition

BACKGROUND

The nucleocapsid (NC) domain of HIV-1 Gag is responsible for specific recognition and packaging of genomic RNA (gRNA) into new viral particles. This occurs through specific interactions between the Gag NC domain and the Psi packaging signal in gRNA. In addition to this critical function, NC proteins are also nucleic acid (NA) chaperone proteins that facilitate NA rearrangements during reverse transcription. Although the interaction with Psi and chaperone activity of HIV-1 NC have been well characterized in vitro, little is known about simian immunodeficiency virus (SIV) NC. Non-human primates are frequently used as a platform to study retroviral infection in vivo; thus, it is important to understand underlying mechanistic differences between HIV-1 and SIV NC.

RESULTS

Here, we characterize SIV NC chaperone activity for the first time. Only modest differences are observed in the ability of SIV NC to facilitate reactions that mimic the minus-strand annealing and transfer steps of reverse transcription relative to HIV-1 NC, with the latter displaying slightly higher strand transfer and annealing rates. Quantitative single molecule DNA stretching studies and dynamic light scattering experiments reveal that these differences are due to significantly increased DNA compaction energy and higher aggregation capability of HIV-1 NC relative to the SIV protein. Using salt-titration binding assays, we find that both proteins are strikingly similar in their ability to specifically interact with HIV-1 Psi RNA. In contrast, they do not demonstrate specific binding to an RNA derived from the putative SIV packaging signal.

CONCLUSIONS

Based on these studies, we conclude that (1) HIV-1 NC is a slightly more efficient NA chaperone protein than SIV NC, (2) mechanistic differences between the NA interactions of highly similar retroviral NC proteins are revealed by quantitative single molecule DNA stretching, and (3) SIV NC demonstrates cross-species recognition of the HIV-1 Psi RNA packaging signal.

Changes in the Plasticity of HIV-1 Nef RNA during the Evolution of the North American Epidemic

Because of a high mutation rate, HIV exists as a viral swarm of many sequence variants evolving under various selective pressures from the human immune system. Although the Nef gene codes for the most immunogenic of HIV accessory proteins, which alone makes it of great interest to HIV research, it also encodes an RNA structure, whose contribution to HIV virulence has been largely unexplored. Nef RNA helps HIV escape RNA interference (RNAi) through nucleotide changes and alternative folding. This study examines Historic and Modern Datasets of patient HIV-1 Nef sequences during the evolution of the North American epidemic for local changes in RNA plasticity. By definition, RNA plasticity refers to an RNA molecule’s ability to take alternative folds (i.e., alternative conformations). Our most important finding is that an evolutionarily conserved region of the HIV-1 Nef gene, which we denote by R2, recently underwent a statistically significant increase in its RNA plasticity. Thus, our results indicate that Modern Nef R2 typically accommodates an alternative fold more readily than Historic Nef R2. Moreover, the increase in RNA plasticity resides mostly in synonymous nucleotide changes, which cannot be a response to selective pressures on the Nef protein. R2 may therefore be of interest in the development of antiviral RNAi therapies.

Analysis of HIV-1 intersubtype recombination breakpoints suggests region with high pairing probability may be a more fundamental factor than sequence similarity affecting HIV-1 recombination

Background

With increasing data on HIV-1, a more relevant molecular model describing mechanism details of HIV-1 genetic recombination usually requires upgrades. Currently an incomplete structural understanding of the copy choice mechanism along with several other issues in the field that lack elucidation led us to perform an analysis of the correlation between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarity to further explore structural mechanisms.

Methods

Near full length sequences of URFs from Asia, Europe, and Africa (one sequence/patient), and representative sequences of worldwide CRFs were retrieved from the Los Alamos HIV database. Their recombination patterns were analyzed by jpHMM in detail. Then the relationships between breakpoint distributions and (1) the probability of base pairing, and (2) intersubtype genetic similarities were investigated.

Results

Pearson correlation test showed that all URF groups and the CRF group exhibit the same breakpoint distribution pattern. Additionally, the Wilcoxon two-sample test indicated a significant and inexplicable limitation of recombination in regions with high pairing probability. These regions have been found to be strongly conserved across distinct biological states (i.e., strong intersubtype similarity), and genetic similarity has been determined to be a very important factor promoting recombination. Thus, the results revealed an unexpected disagreement between intersubtype similarity and breakpoint distribution, which were further confirmed by genetic similarity analysis. Our analysis reveals a critical conflict between results from natural HIV-1 isolates and those from HIV-1-based assay vectors in which genetic similarity has been shown to be a very critical factor promoting recombination.

Conclusions

These results indicate the region with high-pairing probabilities may be a more fundamental factor affecting HIV-1 recombination than sequence similarity in natural HIV-1 infections. Our findings will be relevant in furthering the understanding of HIV-1 recombination mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0616-1) contains supplementary material, which is available to authorized users.

A Phylogenetic Survey on the Structure of the HIV-1 Leader RNA Domain That Encodes the Splice Donor Signal

RNA splicing is a critical step in the human immunodeficiency virus type 1 (HIV-1) replication cycle because it controls the expression of the complex viral proteome. The major 5′ splice site (5′ss) that is positioned in the untranslated leader of the HIV-1 RNA transcript is of particular interest because it is used for the production of the more than 40 differentially spliced subgenomic mRNAs. HIV-1 splicing needs to be balanced tightly to ensure the proper levels of all viral proteins, including the Gag-Pol proteins that are translated from the unspliced RNA. We previously presented evidence that the major 5′ss is regulated by a repressive local RNA structure, the splice donor (SD) hairpin, that masks the 11 nucleotides (nts) of the 5′ss signal for recognition by U1 small nuclear RNA (snRNA) of the spliceosome machinery. A strikingly different multiple-hairpin RNA conformation was recently proposed for this part of the HIV-1 leader RNA. We therefore inspected the sequence of natural HIV-1 isolates in search for support, in the form of base pair (bp) co-variations, for the different RNA conformations.

Transcending the prediction paradigm: novel applications of SHAPE to RNA function and evolution

Selective 2′‐hydroxyl acylation analyzed by primer extension (SHAPE) provides information on RNA structure at single‐nucleotide resolution. It is most often used in conjunction with RNA secondary structure prediction algorithms as a probabilistic or thermodynamic restraint. With the recent advent of ultra‐high‐throughput approaches for collecting SHAPE data, the applications of this technology are extending beyond structure prediction. In this review, we discuss recent applications of SHAPE data in the transcriptomic context and how this new experimental paradigm is changing our understanding of these experiments and RNA folding in general. SHAPE experiments probe both the secondary and tertiary structure of an RNA, suggesting that model‐free approaches for within and comparative RNA structure analysis can provide significant structural insight without the need for a full structural model. New methods incorporating SHAPE at different nucleotide resolutions are required to parse these transcriptomic data sets to transcend secondary structure modeling with global structural metrics. These ‘multiscale’ approaches provide deeper insights into RNA global structure, evolution, and function in the cell. WIREs RNA 2017, 8:e1374. doi: 10.1002/wrna.1374

For further resources related to this article, please visit the WIREs website.

3D RNA and Functional Interactions from Evolutionary Couplings

Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences far outpaces the research on the structure and functional interactions of these RNA gene sequences. We mine the evolutionary sequence record to derive precise information about the function and structure of RNAs and RNA-protein complexes. As in protein structure prediction, we use maximum entropy global probability models of sequence co-variation to infer evolutionarily constrained nucleotide-nucleotide interactions within RNA molecules and nucleotide-amino acid interactions in RNA-protein complexes. The predicted contacts allow all-atom blinded 3D structure prediction at good accuracy for several known RNA structures and RNA-protein complexes. For unknown structures, we predict contacts in 160 non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings help identify important functional interactions-e.g., at switch points in riboswitches and at a complex nucleation site in HIV. Aided by increasing sequence accumulation, evolutionary coupling analysis can accelerate the discovery of functional interactions and 3D structures involving RNA.

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.

Innate immune system activation by viral RNA: How to predict it?

The immune system is able to identify foreign pathogens via different pathways. In the case of viral infection, recognition of the viral RNA is a crucial step, and many efforts have been made to understand which features of viral RNA are detected by the immune system. The biased viral RNA composition, measured as host-virus nucleotidic divergence, or CpG enrichment, has been proposed as salient signal. Peculiar structural features of these RNA could also be related to the immune system activation. Here, we gather multiple datasets and proceed to a meta-analysis to uncover the best predictors of immune system activation by viral RNA. "A" nucleotide content and Minimum Folding Energy are good predictors, and are more easily generalized than more complex indicators suggested previously. As RNA composition and structure are highly correlated, we suggest further experiments on synthetic sequences to identify the viral RNA sensing mechanisms by immune system receptors.

Full-length RNA structure prediction of the HIV-1 genome reveals a conserved core domain

A distance constrained secondary structural model of the ≈10 kb RNA genome of the HIV-1 has been predicted but higher-order structures, involving long distance interactions, are currently unknown. We present the first global RNA secondary structure model for the HIV-1 genome, which integrates both comparative structure analysis and information from experimental data in a full-length prediction without distance constraints. Besides recovering known structural elements, we predict several novel structural elements that are conserved in HIV-1 evolution. Our results also indicate that the structure of the HIV-1 genome is highly variable in most regions, with a limited number of stable and conserved RNA secondary structures. Most interesting, a set of long distance interactions form a core organizing structure (COS) that organize the genome into three major structural domains. Despite overlapping protein-coding regions the COS is supported by a particular high frequency of compensatory base changes, suggesting functional importance for this element. This new structural element potentially organizes the whole genome into three major domains protruding from a conserved core structure with potential roles in replication and evolution for the virus.