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

HIV-1 RNA genome packaging: it's G-rated

A member of the Retroviridae, human immunodeficiency virus type 1 (HIV-1), uses the RNA genome packaged into nascent virions to transfer genetic information to its progeny. The genome packaging step is a highly regulated and extremely efficient process as a vast majority of virus particles contain two copies of full-length unspliced HIV-1 RNA that form a dimer. Thus, during virus assembly HIV-1 can identify and selectively encapsidate HIV-1 unspliced RNA from an abundant pool of cellular RNAs and various spliced HIV-1 RNAs. Several "G" features facilitate the packaging of a dimeric RNA genome. The viral polyprotein Gag orchestrates virus assembly and mediates RNA genome packaging. During this process, Gag preferentially binds unpaired guanosines within the highly structured 5' untranslated region (UTR) of HIV-1 RNA. In addition, the HIV-1 unspliced RNA provides a scaffold that promotes Gag:Gag interactions and virus assembly, thereby ensuring its packaging. Intriguingly, recent studies have shown that the use of different guanosines at the junction of U3 and R as transcription start sites results in HIV-1 unspliced RNA species with 99.9% identical sequences but dramatically distinct 5' UTR conformations. Consequently, one species of unspliced RNA is preferentially packaged over other nearly identical RNAs. These studies reveal how conformations affect the functions of HIV-1 RNA elements and the complex regulation of HIV-1 replication. In this review, we summarize cis- and trans-acting elements critical for HIV-1 RNA packaging, locations of Gag:RNA interactions that mediate genome encapsidation, and the effects of transcription start sites on the structure and packaging of HIV-1 RNA.

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.

Synonymous Variants of Uncertain Silence

Synonymous variants, traditionally regarded as silent mutations due to their lack of impact on protein sequence, structure and function, have been the subject of increasing scrutiny. This commentary explores the emerging evidence challenging the notion of synonymous variants as functionally inert. Analysis of the activity of 70 synonymous variants in the HIV Tat transcription factor revealed that 50% of the variants exhibited significant deviations from wild-type activity. Our analysis supports previous work and raises important questions about the broader impact of non-silent synonymous variants in human genes. Considering the potential functional implications, the authors propose classifying such variants as "synonymous variants of uncertain silence" (sVUS), highlighting the need for cautious interpretation and further investigations in clinical and genetic testing settings.

A Synthetic Biology Approach for Vaccine Candidate Design against Delta Strain of SARS-CoV-2 Revealed Disruption of Favored Codon Pair as a Better Strategy over Using Rare Codons

The SARS-CoV-2 delta variant (B.1.617.2) appeared for the first time in December 2020 and later spread worldwide. Currently available vaccines are not so efficacious in curbing the viral pathogenesis of the delta strain of COVID; therefore, the development of a safe and effective vaccine is required. In the present study, we envisaged molecular patterns in the structural genes' spike, nucleoprotein, membrane, and envelope of the SARS-CoV-2 delta variant. The study was based on determining compositional features, dinucleotide odds ratio, synonymous codon usage, positive and negative codon contexts, rare codons, and insight into relatedness between the human host isoacceptor tRNA and preferred codons from the structural genes. We found specific patterns, including a significant abundance of T nucleotide over all other three nucleotides. The underrepresentation of GpA, GpG, CpC, and CpG dinucleotides and the overrepresentation of TpT, ApA, CpT, and TpG were observed. A preference towards ACT- (Thr), AAT- (Asn), TTT- (Phe), and TTG- (Leu) initiated codons and aversion towards CGG (Arg), CCG (Pro), and CAC (His) was present in the structural genes of the delta strain. The interaction between the host tRNA pool and preferred codons of the envisaged structural genes revealed that the virus preferred the codons for those suboptimal numbers of isoacceptor tRNA were present. We see this as a strategy adapted by the virus to keep the translation rate low to facilitate the correct folding of viral proteins. The information generated in the study helps design the attenuated vaccine candidate against the SARS-CoV-2 delta variant using a synthetic biology approach. Three strategies were tested: changing TpT to TpA, introducing rare codons, and disrupting favored codons. It found that disrupting favored codons is a better approach to reducing virus fitness and attenuating SARS-CoV-2 delta strain using structural genes.

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

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

Protein-protein interaction and non-interaction predictions using gene sequence natural vector

Predicting protein–protein interaction and non-interaction are two important different aspects of multi-body structure predictions, which provide vital information about protein function. Some computational methods have recently been developed to complement experimental methods, but still cannot effectively detect real non-interacting protein pairs. We proposed a gene sequence-based method, named NVDT (Natural Vector combine with Dinucleotide and Triplet nucleotide), for the prediction of interaction and non-interaction. For protein–protein non-interactions (PPNIs), the proposed method obtained accuracies of 86.23% for Homo sapiens and 85.34% for Mus musculus, and it performed well on three types of non-interaction networks. For protein-protein interactions (PPIs), we obtained accuracies of 99.20, 94.94, 98.56, 95.41, and 94.83% for Saccharomyces cerevisiae, Drosophila melanogaster, Helicobacter pylori, Homo sapiens, and Mus musculus, respectively. Furthermore, NVDT outperformed established sequence-based methods and demonstrated high prediction results for cross-species interactions. NVDT is expected to be an effective approach for predicting PPIs and PPNIs.

Protein-protein non-interactions and interactions are distinguished and predicted by gene sequence using single nucleotide and contiguous nucleotides combined with machine learning models.

Altered HIV-1 mRNA Splicing Due to Drug-Resistance-Associated Mutations in Exon 2/2b

The underlying molecular mechanism and their general effect on the replication capacity of HIV 1 drug-resistance-associated mutations is often poorly understood. To elucidate the effect of two such mutations located in a region with a high density of spicing regulatory elements on the HIV-1-splicing outcome, bioinformatic predictions were combined with transfection and infection experiments. Results show that the previously described R263K drug-resistance-associated integrase mutation has additionally a severe effect on the ESE2b splicing regulatory element (SRE) in exon 2b, which causes loss of SD2b recognition. This was confirmed by an R263R silent mutation with a similar predicted effect on the exon 2b SRE. In contrast, a V260I mutation and its silent counterpart with a lower effect on ESS2b did not exhibit any differences in the splicing pattern. Since HIV-1 highly relies on a balanced splicing reaction, changes in the splicing outcome can contribute to changes in viral replication and might add to the effect of escape mutations toward antiviral drugs. Thus, a classification of mutations purely addressing proteins is insufficient.

An ancient retroviral RNA element hidden in mammalian genomes and its involvement in co-opted retroviral gene regulation

Background

Retroviruses utilize multiple unique RNA elements to control RNA processing and translation. However, it is unclear what functional RNA elements are present in endogenous retroviruses (ERVs). Gene co-option from ERVs sometimes entails the conservation of viral cis-elements required for gene expression, which might reveal the RNA regulation in ERVs.

Results

Here, we characterized an RNA element found in ERVs consisting of three specific sequence motifs, called SPRE. The SPRE-like elements were found in different ERV families but not in any exogenous viral sequences examined. We observed more than a thousand of copies of the SPRE-like elements in several mammalian genomes; in human and marmoset genomes, they overlapped with lineage-specific ERVs. SPRE was originally found in human syncytin-1 and syncytin-2. Indeed, several mammalian syncytin genes: mac-syncytin-3 of macaque, syncytin-Ten1 of tenrec, and syncytin-Car1 of Carnivora, contained the SPRE-like elements. A reporter assay revealed that the enhancement of gene expression by SPRE depended on the reporter genes. Mutation of SPRE impaired the wild-type syncytin-2 expression while the same mutation did not affect codon-optimized syncytin-2, suggesting that SPRE activity depends on the coding sequence.

Conclusions

These results indicate multiple independent invasions of various mammalian genomes by retroviruses harboring SPRE-like elements. Functional SPRE-like elements are found in several syncytin genes derived from these retroviruses. This element may facilitate the expression of viral genes, which were suppressed due to inefficient codon frequency or repressive elements within the coding sequences. These findings provide new insights into the long-term evolution of RNA elements and molecular mechanisms of gene expression in retroviruses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-021-00580-2.

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.

Epistatic Net allows the sparse spectral regularization of deep neural networks for inferring fitness functions

Despite recent advances in high-throughput combinatorial mutagenesis assays, the number of labeled sequences available to predict molecular functions has remained small for the vastness of the sequence space combined with the ruggedness of many fitness functions. While deep neural networks (DNNs) can capture high-order epistatic interactions among the mutational sites, they tend to overfit to the small number of labeled sequences available for training. Here, we developed Epistatic Net (EN), a method for spectral regularization of DNNs that exploits evidence that epistatic interactions in many fitness functions are sparse. We built a scalable extension of EN, usable for larger sequences, which enables spectral regularization using fast sparse recovery algorithms informed by coding theory. Results on several biological landscapes show that EN consistently improves the prediction accuracy of DNNs and enables them to outperform competing models which assume other priors. EN estimates the higher-order epistatic interactions of DNNs trained on massive sequence spaces-a computational problem that otherwise takes years to solve.

MrHAMER yields highly accurate single molecule viral sequences enabling analysis of intra-host evolution

Technical challenges remain in the sequencing of RNA viruses due to their high intra-host diversity. This bottleneck is particularly pronounced when interrogating long-range co-evolved genetic interactions given the read-length limitations of next-generation sequencing platforms. This has hampered the direct observation of these genetic interactions that code for protein-protein interfaces with relevance in both drug and vaccine development. Here we overcome these technical limitations by developing a nanopore-based long-range viral sequencing pipeline that yields accurate single molecule sequences of circulating virions from clinical samples. We demonstrate its utility in observing the evolution of individual HIV Gag-Pol genomes in response to antiviral pressure. Our pipeline, called Multi-read Hairpin Mediated Error-correction Reaction (MrHAMER), yields >1000s of viral genomes per sample at 99.9% accuracy, maintains the original proportion of sequenced virions present in a complex mixture, and allows the detection of rare viral genomes with their associated mutations present at <1% frequency. This method facilitates scalable investigation of genetic correlates of resistance to both antiviral therapy and immune pressure and enables the identification of novel host-viral and viral-viral interfaces that can be modulated for therapeutic benefit.

Synonymous Codon Pair Recoding of the HIV-1 env Gene Affects Virus Replication Capacity

Synonymous codon pair deoptimization is an efficient strategy for virus attenuation; however, the underlying mechanism remains controversial. Here, we optimized and deoptimized the codon pair bias (CPB) of the human immunodeficiency virus type 1 (HIV-1) envelope (env) gene to investigate the influence of env synonymous CPB recoding on virus replication capacity, as well as the potential mechanism. We found that env CPB deoptimization did not always generate attenuation, whereas CPB optimization attenuated virus replication in MT-4 cells. Furthermore, virus attenuation correlated with reduced Env protein production but not with decreased viral RNA synthesis. Remarkably, in our model, increasing the number of CpG dinucleotides in the 5′ end of env did not reduce the replication capacity of HIV-1. These results indicate that factors other than CPB or CpG content may have impacted the viral fitness of the synonymously recoded study variants. Our findings provide evidence that CPB recoding-associated attenuation can affect translation efficiency. Moreover, we demonstrated that an increased number of CpGs in the 5′ end of HIV-1 env is not always associated with reduced virus replication capacity.

Targeted Restriction of Viral Gene Expression and Replication by the ZAP Antiviral System

The zinc finger antiviral protein (ZAP) restricts the replication of a broad range of RNA and DNA viruses. ZAP directly binds viral RNA, targeting it for degradation and inhibiting its translation. While the full scope of RNA determinants involved in mediating selective ZAP activity are unclear, ZAP binds CpG dinucleotides, dictating at least part of its target specificity. ZAP interacts with many cellular proteins, although only a few have been demonstrated to be essential for its antiviral activity, including the 3'-5' exoribonuclease exosome complex, TRIM25, and KHNYN. In addition to inhibiting viral gene expression, ZAP also directly and indirectly targets a subset of cellular messenger RNAs to regulate the innate immune response. Overall, ZAP protects a cell from viral infection by restricting viral replication and regulating cellular gene expression. Further understanding of the ZAP antiviral system may allow for novel viral vaccine and anticancer therapy development. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The HIV 5' Gag Region Displays a Specific Nucleotide Bias Regulating Viral Splicing and Infectivity

Alternative splicing and the expression of intron-containing mRNAs is one hallmark of HIV gene expression. To facilitate the otherwise hampered nuclear export of non-fully processed mRNAs, HIV encodes the Rev protein, which recognizes its intronic response element and fuels the HIV RNAs into the CRM-1-dependent nuclear protein export pathway. Both alternative splicing and Rev-dependency are regulated by the primary HIV RNA sequence. Here, we show that these processes are extremely sensitive to sequence alterations in the 5’coding region of the HIV genomic RNA. Increasing the GC content by insertion of either GFP or silent mutations activates a cryptic splice donor site in gag, entirely deregulates the viral splicing pattern, and lowers infectivity. Interestingly, an adaptation of the inserted GFP sequence toward an HIV-like nucleotide bias reversed these phenotypes completely. Of note, the adaptation yielded completely different primary sequences although encoding the same amino acids. Thus, the phenotypes solely depend on the nucleotide composition of the two GFP versions. This is a strong indication of an HIV-specific mRNP code in the 5′ gag region wherein the primary RNA sequence bias creates motifs for RNA-binding proteins and controls the fate of the HIV-RNA in terms of viral gene expression and infectivity.

Impact of Synonymous Genome Recoding on the HIV Life Cycle

Synonymous mutations within protein coding regions introduce changes in DNA or messenger (m) RNA, without mutating the encoded proteins. Synonymous recoding of virus genomes has facilitated the identification of previously unknown virus biological features. Moreover, large-scale synonymous recoding of the genome of human immunodeficiency virus type 1 (HIV-1) has elucidated new antiviral mechanisms within the innate immune response, and has improved our knowledge of new functional virus genome structures, the relevance of codon usage for the temporal regulation of viral gene expression, and HIV-1 mutational robustness and adaptability. Continuous improvements in our understanding of the impacts of synonymous substitutions on virus phenotype - coupled with the decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments - have enhanced our ability to identify potential HIV-1 and host factors and other aspects involved in the infection process. In this review, we address how silent mutagenesis impacts HIV-1 phenotype and replication capacity. We also discuss the general potential of synonymous recoding of the HIV-1 genome to elucidate unknown aspects of the virus life cycle, and to identify new therapeutic targets.

Plasmacytoid dendritic cells sense HIV replication before detectable viremia following treatment interruption

Plasmacytoid dendritic cells (pDCs) are robust producers of IFNα and one of the first immune cells to respond to SIV infection. To elucidate responses to early HIV-1 replication, we studied blood pDCs in 29 HIV-infected participants who initiated antiretroviral therapy during acute infection and underwent analytic treatment interruption (ATI). We observed an increased frequency of partially activated pDCs in the blood before detection of HIV RNA. Concurrent with peak pDC frequency, we detected a transient decline in the ability of pDCs to produce IFNα in vitro, which correlated with decreased phosphorylation of IFN regulatory factory 7 (IRF7) and NF-κB. The levels of phosphorylated IRF7 and NF-κB inversely correlated with plasma IFNα2 levels, implying that pDCs were refractory to in vitro stimulation after IFNα production in vivo. After ATI, decreased expression of IFN genes in pDCs inversely correlated with the time to viral detection, suggesting that pDC IFN loss is part of an effective early immune response. These data from a limited cohort provide a critical first step in understanding the earliest immune response to HIV-1 and suggest that changes in blood pDC frequency and function can be used as an indicator of viral replication before detectable plasma viremia.

HIV-1: To Splice or Not to Splice, That Is the Question

The transcription of the HIV-1 provirus results in only one type of transcript-full length genomic RNA. To make the mRNA transcripts for the accessory proteins Tat and Rev, the genomic RNA must completely splice. The mRNA transcripts for Vif, Vpr, and Env must undergo splicing but not completely. Genomic RNA (which also functions as mRNA for the Gag and Gag/Pro/Pol precursor polyproteins) must not splice at all. HIV-1 can tolerate a surprising range in the relative abundance of individual transcript types, and a surprising amount of aberrant and even odd splicing; however, it must not over-splice, which results in the loss of full-length genomic RNA and has a dramatic fitness cost. Cells typically do not tolerate unspliced/incompletely spliced transcripts, so HIV-1 must circumvent this cell policing mechanism to allow some splicing while suppressing most. Splicing is controlled by RNA secondary structure, cis-acting regulatory sequences which bind splicing factors, and the viral protein Rev. There is still much work to be done to clarify the combinatorial effects of these splicing regulators. These control mechanisms represent attractive targets to induce over-splicing as an antiviral strategy. Finally, splicing has been implicated in latency, but to date there is little supporting evidence for such a mechanism. In this review we apply what is known of cellular splicing to understand splicing in HIV-1, and present data from our newer and more sensitive deep sequencing assays quantifying the different HIV-1 transcript types.

RanDeL-Seq: a High-Throughput Method to Map Viral cis- and trans-Acting Elements

Recent studies have renewed interest in developing novel antiviral therapeutics and vaccines based on defective interfering particles (DIPs)—a subset of viral deletion mutants that conditionally replicate. Identifying and engineering DIPs require that viral cis- and trans-acting elements be accurately mapped.

It has long been known that noncoding genomic regions can be obligate cis elements acted upon in trans by gene products. In viruses, cis elements regulate gene expression, encapsidation, and other maturation processes, but mapping these elements relies on targeted iterative deletion or laborious prospecting for rare spontaneously occurring mutants. Here, we introduce a method to comprehensively map viral cis and trans elements at single-nucleotide resolution by high-throughput random deletion. Variable-size deletions are randomly generated by transposon integration, excision, and exonuclease chewback and then barcoded for tracking via sequencing (i.e., random deletion library sequencing [RanDeL-seq]). Using RanDeL-seq, we generated and screened >23,000 HIV-1 variants to generate a single-base resolution map of HIV-1’s cis and trans elements. The resulting landscape recapitulated HIV-1’s known cis-acting elements (i.e., long terminal repeat [LTR], Ψ, and Rev response element [RRE]) and, surprisingly, indicated that HIV-1’s central DNA flap (i.e., central polypurine tract [cPPT] to central termination sequence [CTS]) is as critical as the LTR, Ψ, and RRE for long-term passage. Strikingly, RanDeL-seq identified a previously unreported ∼300-bp region downstream of RRE extending to splice acceptor 7 that is equally critical for sustained viral passage. RanDeL-seq was also used to construct and screen a library of >90,000 variants of Zika virus (ZIKV). Unexpectedly, RanDeL-seq indicated that ZIKV’s cis-acting regions are larger than the untranscribed (UTR) termini, encompassing a large fraction of the nonstructural genes. Collectively, RanDeL-seq provides a versatile framework for generating viral deletion mutants, enabling discovery of replication mechanisms and development of novel antiviral therapeutics, particularly for emerging viral infections.

HIV-1 Lethality and Loss of Env Protein Expression Induced by Single Synonymous Substitutions in the Virus Genome Intronic-Splicing Silencer

Synonymous genome recoding has been widely used to study different aspects of virus biology. Codon usage affects the temporal regulation of viral gene expression. In this study, we performed synonymous codon mutagenesis to investigate whether codon usage affected HIV-1 Env protein expression and virus viability. We replaced the codons AGG, GAG, CCU, ACU, CUC, and GGG of the HIV-1 env gene with the synonymous codons CGU, GAA, CCG, ACG, UUA, and GGA, respectively. We found that recoding the Env protein gp120 coding region (excluding the Rev response element [RRE]) did not significantly affect virus replication capacity, even though we introduced 15 new CpG dinucleotides. In contrast, changing a single codon (AGG to CGU) located in the gp41 coding region (HXB2 env position 2125 to 2127), which was included in the intronic splicing silencer (ISS), completely abolished virus replication and Env expression. Computational analyses of this mutant revealed a severe disruption in the ISS RNA secondary structure. A variant that restored ISS secondary RNA structure also reestablished Env production and virus viability. Interestingly, this codon variant prevented both virus replication and Env translation in a eukaryotic expression system. These findings suggested that disrupting mRNA splicing was not the only means of inhibiting translation. Our findings indicated that synonymous gp120 recoding was not always deleterious to HIV-1 replication. Importantly¸ we found that disrupting an external ISS loop strongly affected HIV-1 replication and Env translation.IMPORTANCE Synonymous substitutions can influence virus phenotype, replication capacity, and virulence. In this study, we explored how synonymous codon mutations impacted HIV-1 Env protein expression and virus replication capacity. We changed a single codon, AGG to CGU, which was located in the gp41 coding region (env nucleotide residues 2125 to 2127) and was included in the HIV-1 intronic splicing silencer. This change completely abolished virus replication and Env expression. We also found that changing codon usage in the gp120 region by including an increased number of CpG dinucleotides did not significantly affect Env expression or virus viability. Our findings showed that synonymous recoding was useful for altering viral phenotype and exploring virus biology.

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.

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.

Determination of RNA structural diversity and its role in HIV-1 RNA splicing

Human immunodeficiency virus-1 (HIV-1) is a retrovirus with a 10-kb single-stranded RNA genome. HIV-1 must express all of its gene products from the same primary transcript, which undergoes alternative splicing to produce diverse protein products, including structural proteins and regulatory factors^1,2. Despite the critical role of alternative splicing, the mechanisms driving splice-site choice are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements³. We use DMS-MaPseq to probe the structure of HIV-1 RNA in cells and develop an algorithm called Detection of RNA folding Ensembles using Expectation-Maximization (DREEM), which reveals alternative conformations assumed by the same RNA sequence. Contrary to previous models, which analyzed population averages⁴, our results reveal the widespread heterogeneous nature of HIV-1 RNA structure. In addition to confirming that in vitro characterized alternative structures for the HIV-1 Rev Responsive Element (RRE) exist in cells, we discover alternative conformations at critical splice sites that influence the ratio of transcript isoforms. Our simultaneous measurement of splicing and intracellular RNA structure provides evidence for the long-standing hypothesis^5–7 that RNA conformation heterogeneity regulates splice site usage and viral gene expression.

CpG Dinucleotides Inhibit HIV-1 Replication through Zinc Finger Antiviral Protein (ZAP)-Dependent and -Independent Mechanisms

CpG dinucleotides are suppressed in the genomes of many vertebrate RNA viruses, including HIV-1. The cellular antiviral protein ZAP (zinc finger antiviral protein) binds CpGs and inhibits HIV-1 replication when CpGs are introduced into the viral genome. However, it is not known if ZAP-mediated restriction is the only mechanism driving CpG suppression. To determine how CpG dinucleotides affect HIV-1 replication, we increased their abundance in multiple regions of the viral genome and analyzed the effect on RNA expression, protein abundance, and infectious-virus production. We found that the antiviral effect of CpGs was not correlated with their abundance. Interestingly, CpGs inserted into some regions of the genome sensitize the virus to ZAP antiviral activity more efficiently than insertions into other regions, and this sensitivity can be modulated by interferon treatment or ZAP overexpression. Furthermore, the sensitivity of the virus to endogenous ZAP was correlated with its sensitivity to the ZAP cofactor KHNYN. Finally, we show that CpGs in some contexts can also inhibit HIV-1 replication by ZAP-independent mechanisms, and one of these is the activation of a cryptic splice site at the expense of a canonical splice site. Overall, we show that the location and sequence context of the CpG in the viral genome determines its antiviral activity.IMPORTANCE Some RNA virus genomes are suppressed in the nucleotide combination of a cytosine followed by a guanosine (CpG), indicating that they are detrimental to the virus. The antiviral protein ZAP binds viral RNA containing CpGs and prevents the virus from multiplying. However, it remains unknown how the number and position of CpGs in viral genomes affect restriction by ZAP and whether CpGs have other antiviral mechanisms. Importantly, manipulating the CpG content in viral genomes could help create new vaccines. HIV-1 shows marked CpG suppression, and by introducing CpGs into its genome, we show that ZAP efficiently targets a specific region of the viral genome, that the number of CpGs does not predict the magnitude of antiviral activity, and that CpGs can inhibit HIV-1 gene expression through a ZAP-independent mechanism. Overall, the position of CpGs in the HIV-1 genome determines the magnitude and mechanism through which they inhibit the virus.

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.

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.

HIV-1 is a very rapidly mutating organism, however some parts of its genetic material change more than others. We looked for coding regions of HIV-1 that change relatively little, by turning the problem of finding such regions into a problem in signal processing, and solving this using a novel analytical approach that we recently described. We investigated why the regions we identified change less, including using the genetic code in the regions we found to prime an algorithm to predict their structures. In some cases there are already known functions for the features we found, in others they provide new insights into the properties of known regions, and in some cases we identify new regions that vary less for as yet unknown functional reasons.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

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.

Integrative Analysis of Zika Virus Genome RNA Structure Reveals Critical Determinants of Viral Infectivity

Zika virus (ZIKV) strains can be classified into the ancestral African and contemporary Asian lineages, with the latter responsible for recent epidemics associated with neurological conditions. To understand how Asian strains lead to exacerbated disease, a crucial step is identifying genomic variations that affect infectivity and pathogenicity. Here we use two high-throughput sequencing approaches to assess RNA secondary structures and intramolecular RNA-RNA interactions in vivo for the RNA genomes of Asian and African ZIKV lineages. Our analysis identified functional RNA structural elements and a functional long-range intramolecular interaction specific for the Asian epidemic strains. Mutants that disrupt this extended RNA interaction between the 5' UTR and the E protein coding region reduce virus infectivity, which is partially rescued with compensatory mutants, restoring this RNA-RNA interaction. These findings illuminate the structural basis of ZIKV regulation and provide a resource for the discovery of RNA structural elements important for ZIKV infection.

RNA Structure-A Neglected Puppet Master for the Evolution of Virus and Host Immunity

The central dogma of molecular biology describes the flow of genetic information from DNA to protein via an RNA intermediate. For many years, RNA has been considered simply as a messenger relaying information between DNA and proteins. Recent advances in next generation sequencing technology, bioinformatics, and non-coding RNA biology have highlighted the many important roles of RNA in virtually every biological process. Our understanding of RNA biology has been further enriched by a number of significant advances in probing RNA structures. It is now appreciated that many cellular and viral biological processes are highly dependent on specific RNA structures and/or sequences, and such reliance will undoubtedly impact on the evolution of both hosts and viruses. As a contribution to this special issue on host immunity and virus evolution, it is timely to consider how RNA sequences and structures could directly influence the co-evolution between hosts and viruses. In this manuscript, we begin by stating some of the basic principles of RNA structures, followed by describing some of the critical RNA structures in both viruses and hosts. More importantly, we highlight a number of available new tools to predict and to evaluate novel RNA structures, pointing out some of the limitations readers should be aware of in their own analyses.

HIV-1 Protease Evolvability Is Affected by Synonymous Nucleotide Recoding

One unexplored aspect of HIV-1 genetic architecture is how codon choice influences population diversity and evolvability. Here we compared the levels of development of HIV-1 resistance to protease inhibitors (PIs) between wild-type (WT) virus and a synthetic virus (MAX) carrying a codon-pair-reengineered protease sequence including 38 (13%) synonymous mutations. The WT and MAX viruses showed indistinguishable replication in MT-4 cells or peripheral blood mononuclear cells (PBMCs). Both viruses were subjected to serial passages in MT-4 cells, with selective pressure from the PIs atazanavir (ATV) and darunavir (DRV). After 32 successive passages, both the WT and MAX viruses developed phenotypic resistance to PIs (50% inhibitory concentrations [IC₅₀s] of 14.6 ± 5.3 and 21.2 ± 9 nM, respectively, for ATV and 5.9 ± 1.0 and 9.3 ± 1.9, respectively, for DRV). Ultradeep sequence clonal analysis revealed that both viruses harbored previously described mutations conferring resistance to ATV and DRV. However, the WT and MAX virus proteases showed different resistance variant repertoires, with the G16E and V77I substitutions observed only in the WT and the L33F, S37P, G48L, Q58E/K, and L89I substitutions detected only in the MAX virus. Remarkably, the G48L and L89I substitutions are rarely found in vivo in PI-treated patients. The MAX virus showed significantly higher nucleotide and amino acid diversity of the propagated viruses with and without PIs (P < 0.0001), suggesting a higher selective pressure for change in this recoded virus. Our results indicate that the HIV-1 protease position in sequence space delineates the evolution of its mutant spectrum. Nevertheless, the investigated synonymously recoded variant showed mutational robustness and evolvability similar to those of the WT virus.IMPORTANCE Large-scale synonymous recoding of virus genomes is a new tool for exploring various aspects of virus biology. Synonymous virus genome recoding can be used to investigate how a virus's position in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenesis. In this study, we evaluated how synonymous recoding of the human immunodeficiency virus type 1 (HIV-1) protease affects the development of protease inhibitor (PI) resistance. HIV-1 protease is a main target of current antiretroviral therapies. Our present results demonstrate that the wild-type (WT) virus and a virus with recoded protease exhibited different patterns of resistance mutations after PI treatment. Nevertheless, the developed PI resistance phenotypes were indistinguishable between the recoded virus and the WT virus, suggesting that the HIV-1 strain with synonymously recoded protease and the WT virus are equally robust and evolvable.

In cell mutational interference mapping experiment (in cell MIME) identifies the 5' polyadenylation signal as a dual regulator of HIV-1 genomic RNA production and packaging

Non-coding RNA regulatory elements are important for viral replication, making them promising targets for therapeutic intervention. However, regulatory RNA is challenging to detect and characterise using classical structure-function assays. Here, we present in cell Mutational Interference Mapping Experiment (in cell MIME) as a way to define RNA regulatory landscapes at single nucleotide resolution under native conditions. In cell MIME is based on (i) random mutation of an RNA target, (ii) expression of mutated RNA in cells, (iii) physical separation of RNA into functional and non-functional populations, and (iv) high-throughput sequencing to identify mutations affecting function. We used in cell MIME to define RNA elements within the 5' region of the HIV-1 genomic RNA (gRNA) that are important for viral replication in cells. We identified three distinct RNA motifs controlling intracellular gRNA production, and two distinct motifs required for gRNA packaging into virions. Our analysis reveals the 73AAUAAA78 polyadenylation motif within the 5' PolyA domain as a dual regulator of gRNA production and gRNA packaging, and demonstrates that a functional polyadenylation signal is required for viral packaging even though it negatively affects gRNA production.