The A-rich RNA sequences of HIV-1 pol are important for the synthesis of viral cDNA

The sound of silence: mechanisms and implications of HUSH complex function

The vertebrate genome is under constant threat of invasion by genetic parasites. Whether the host can immediately recognize and respond to invading elements has been unclear. The discovery of the human silencing hub (HUSH) complex, and the finding that it provides immediate protection from genome invasion by silencing products of reverse transcription, have important implications for mammalian genome evolution. In this review, we summarize recent insights into HUSH function and describe how cellular introns provide a novel means of self-nonself discrimination, allowing HUSH to recognize and transcriptionally repress a broad range of intronless genetic elements. We discuss how HUSH contributes to genome evolution, and highlight studies reporting the critical role of HUSH in development and implicating HUSH in the control of immune signaling and cancer progression.

Genome surveillance by HUSH-mediated silencing of intronless mobile elements

All life forms defend their genome against DNA invasion. Eukaryotic cells recognize incoming DNA and limit its transcription through repressive chromatin modifications. The human silencing hub (HUSH) complex transcriptionally represses long interspersed element-1 retrotransposons (L1s) and retroviruses through histone H3 lysine 9 trimethylation (H3K9me3)1-3. How HUSH recognizes and initiates silencing of these invading genetic elements is unknown. Here we show that HUSH is able to recognize and transcriptionally repress a broad range of long, intronless transgenes. Intron insertion into HUSH-repressed transgenes counteracts repression, even in the absence of intron splicing. HUSH binds transcripts from the target locus, prior to and independent of H3K9me3 deposition, and target transcription is essential for both initiation and propagation of HUSH-mediated H3K9me3. Genomic data reveal how HUSH binds and represses a subset of endogenous intronless genes generated through retrotransposition of cellular mRNAs. Thus intronless cDNA-the hallmark of reverse transcription-provides a versatile way to distinguish invading retroelements from host genes and enables HUSH to protect the genome from 'non-self' DNA, despite there being no previous exposure to the invading element. Our findings reveal the existence of a transcription-dependent genome-surveillance system and explain how it provides immediate protection against newly acquired elements while avoiding inappropriate repression of host genes.

Interplay between Host tRNAs and HIV-1: A Structural Perspective

The cellular metabolism of host tRNAs and life cycle of HIV-1 cross paths at several key virus-host interfaces. Emerging data suggest a multi-faceted interplay between host tRNAs and HIV-1 that plays essential roles, both structural and regulatory, in viral genome replication, genome packaging, and virion biogenesis. HIV-1 not only hijacks host tRNAs and transforms them into obligatory reverse transcription primers but further commandeers tRNAs to regulate the localization of its major structural protein, Gag, via a specific interface. This review highlights recent advances in understanding tRNA-HIV-1 interactions, primarily from a structural perspective, which start to elucidate their underlying molecular mechanisms, intrinsic specificities, and biological significances. Such understanding may provide new avenues toward developing HIV/AIDS treatments and therapeutics including small molecules and RNA biologics that target these host-virus interfaces.

RNA Structures and Their Role in Selective Genome Packaging

To generate infectious viral particles, viruses must specifically select their genomic RNA from milieu that contains a complex mixture of cellular or non-genomic viral RNAs. In this review, we focus on the role of viral encoded RNA structures in genome packaging. We first discuss how packaging signals are constructed from local and long-range base pairings within viral genomes, as well as inter-molecular interactions between viral and host RNAs. Then, how genome packaging is regulated by the biophysical properties of RNA. Finally, we examine the impact of RNA packaging signals on viral evolution.

Biased Mutation and Selection in RNA Viruses

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based data set spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases toward U rather than A in genomes of all viruses. In +ssRNA viruses, we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses, the genomic mutational bias toward U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.

A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection

In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.

Role of host tRNAs and aminoacyl-tRNA synthetases in retroviral replication

The lifecycle of retroviruses and retrotransposons includes a reverse transcription step, wherein dsDNA is synthesized from genomic RNA for subsequent insertion into the host genome. Retroviruses and retrotransposons commonly appropriate major components of the host cell translational machinery, including cellular tRNAs, which are exploited as reverse transcription primers. Nonpriming functions of tRNAs have also been proposed, such as in HIV-1 virion assembly, and tRNA-derived fragments may also be involved in retrovirus and retrotransposon replication. Moreover, host cellular proteins regulate retroviral replication by binding to tRNAs and thereby affecting various steps in the viral lifecycle. For example, in some cases, tRNA primer selection is facilitated by cognate aminoacyl-tRNA synthetases (ARSs), which bind tRNAs and ligate them to their corresponding amino acids, but also have many known nontranslational functions. Multi-omic studies have revealed that ARSs interact with both viral proteins and RNAs and potentially regulate retroviral replication. Here, we review the currently known roles of tRNAs and their derivatives in retroviral and retrotransposon replication and shed light on the roles of tRNA-binding proteins such as ARSs in this process.

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.

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.

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.

Codon usage of HIV regulatory genes is not determined by nucleotide composition

Codon usage bias can be a result of either mutational bias or selection for translational efficiency and/or accuracy. Previous data has suggested that nucleotide composition constraint was the main determinant of HIV codon usage, and that nucleotide composition and codon usage were different between the regulatory genes, tat and rev, and other viral genes. It is not clear whether translational selection contributed to the codon usage difference and how nucleotide composition and translational selection interact to determine HIV codon usage. In this study, a model of codon bias due to GC composition with modification for the A-rich third codon position was used to calculate predicted HIV codon frequencies based on its nucleotide composition. The predicted codon usage of each gene was compared with the actual codon frequency. The predicted codon usage based on GC composition matched well with the actual codon frequencies for the structural genes (gag, pol and env). However, the codon usage of the regulatory genes (tat and rev) could not be predicted. Codon usage of the regulatory genes was also relatively unbiased showing the highest effective number of codons (ENC). Moreover, the codon adaptation index (CAI) of the regulatory genes showed better adaptation to human codons when compared to other HIV genes. Therefore, the early expressed genes responsible for regulation of the replication cycle, tat and rev, were more similar to humans in terms of codon usage and GC content than other HIV genes. This may help these genes to be expressed efficiently during the early stages of infection.

HIV-1 tolerates changes in A-count in a small segment of the pol gene

Background

The HIV-1 RNA genome has a biased nucleotide composition with a surplus of As. Several hypotheses have been put forward to explain this striking phenomenon, but the A-count of the HIV-1 genome has thus far not been systematically manipulated. The reason for this reservation is the likelihood that known and unknown sequence motifs will be affected by such a massive mutational approach, thus resulting in replication-impaired virus mutants. We present the first attempt to increase and decrease the A-count in a relatively small polymerase (pol) gene segment of HIV-1 RNA.

Results

To minimize the mutational impact, a new mutational approach was developed that is inspired by natural sequence variation as present in HIV-1 isolates. This phylogeny-instructed mutagenesis allowed us to create replication-competent HIV-1 mutants with a significantly increased or decreased local A-count. The local A-count of the wild-type (wt) virus (40.2%) was further increased to 46.9% or reduced to 31.7 and 26.3%. These HIV-1 variants replicate efficiently in vitro, despite the fact that the pol changes cause a quite profound move in HIV–SIV sequence space.

Conclusions

Extrapolating these results to the complete 9 kb RNA genome, we may cautiously suggest that the A-rich signature does not have to be maintained. This survey also provided clues that silent codon changes, in particular from G-to-A, determine the subtype-specific sequence signatures.

Regulation of human immunodeficiency virus type 1 (HIV-1) mRNA translation

Human immunodeficiency virus type 1 (HIV-1) mRNA translation is a complex process that uses the host translation machinery to synthesise viral proteins. Several mechanisms for HIV-1 mRNA translation initiation have been proposed including (1) cap-dependent, eIF4E-dependent, (2) cap-dependent, cap-binding complex-dependent, (3) internal ribosome entry sites, and (4) ribosome shunting. While these mechanisms promote HIV-1 mRNA translation in the context of in vitro systems and subgenomic constructs, there are substantial knowledge gaps in understanding how they regulate viral protein production in the context of full-length virus infection. In this review, we will summarise the different translation mechanisms used by HIV-1 mRNAs and the challenges in understanding how they regulate protein synthesis during viral infection.

Biased Nucleotide Composition and Differential Codon Usage Pattern in HIV-1 and HIV-2

HIV-1 and HIV-2 are closely related retroviruses with differences in pathogenicity and geographic distribution. HIV-2 infection is associated with slower disease progression and transmission, longer latency period, low or undetectable plasmatic viral loads, and reduced likelihood of progression to AIDS, compared to HIV-1. In this investigation, we analyzed HIV-2 genes and genomes and compared them with that of HIV-1 belonging to various subtypes. Comparative analysis of the effective number of codons (ENC) for each of the nine genes of the two viruses revealed that the tat gene of HIV-2 had a higher ENC value compared to HIV-1 tat, reflecting lower levels of expression of HIV-2 tat. Lower levels of tat protein particularly during the early stages of infection could result in a lower viral load, lower viral set point, and delayed progression of disease in HIV-2-infected individuals compared to HIV-1-infected subjects. Furthermore, the GC3 composition of the regulatory genes of HIV-2 was ≥50%, suggesting a firm effort by these viruses to adapt themselves to evolutionary survival. We hypothesize that differential codon usage could be one of the possible factors that could contribute to the diminished pathogenicity of HIV-2 in the host as compared to HIV-1.

The thermodynamics of Pr55Gag-RNA interaction regulate the assembly of HIV

The interactions that occur during HIV Pr55^Gag oligomerization and genomic RNA packaging are essential elements that facilitate HIV assembly. However, mechanistic details of these interactions are not clearly defined. Here, we overcome previous limitations in producing large quantities of full-length recombinant Pr55^Gag that is required for isothermal titration calorimetry (ITC) studies, and we have revealed the thermodynamic properties of HIV assembly for the first time. Thermodynamic analysis showed that the binding between RNA and HIV Pr55^Gag is an energetically favourable reaction (ΔG<0) that is further enhanced by the oligomerization of Pr55^Gag. The change in enthalpy (ΔH) widens sequentially from: (1) Pr55^Gag-Psi RNA binding during HIV genome selection; to (2) Pr55^Gag-Guanosine Uridine (GU)-containing RNA binding in cytoplasm/plasma membrane; and then to (3) Pr55^Gag-Adenosine(A)-containing RNA binding in immature HIV. These data imply the stepwise increments of heat being released during HIV biogenesis may help to facilitate the process of viral assembly. By mimicking the interactions between A-containing RNA and oligomeric Pr55^Gag in immature HIV, it was noted that a p6 domain truncated Pr50^{Gag Δp6} is less efficient than full-length Pr55^Gag in this thermodynamic process. These data suggest a potential unknown role of p6 in Pr55^Gag-Pr55^Gag oligomerization and/or Pr55^Gag-RNA interaction during HIV assembly. Our data provide direct evidence on how nucleic acid sequences and the oligomeric state of Pr55^Gag regulate HIV assembly.

Formation of any virus particle will require energy, yet the precise biophysical properties that drive the formation of HIV particles remain undefined. Isothermal titration calorimetry (ITC) is a biophysical technique that is the gold standard to reveal parameters governing biochemical and biophysical reaction. However, ITC requires large amount of proteins for analysis. As large quantities of full-length recombinant HIV Pr55^Gag proteins have not been available in the past 30 years due to technical limitation, a comprehensive thermodynamic analysis of full-length HIV Pr55^Gag has not been possible. Here, we have generated sufficient amount of full-length recombinant HIV Pr55^Gag protein for isothermal titration calorimetry analysis. Our analyses have shown that the major interactions amongst HIV proteins and RNA sequences during viral assembly are energetically favourable reactions. In other words, HIV Pr55^Gag proteins and viral RNA have evolved to overcome the energy barrier for virus formation by utilising energy obtained from protein-RNA interactions in order to facilitate the viral assembly process. Furthermore, HIV also use the oligomeric states of HIV Pr55^Gag proteins and the RNA sequences as means to regulate the viral assembly process.

Sequence-Specific Sensing of Nucleic Acids

Innate immune cells are endowed with many nucleic acid receptors, but the role of sequence in the detection of foreign organisms remains unclear. Can sequence patterns influence recognition? In addition, how can we infer those patterns from sequence data? Here, we detail recent computational and experimental evidence associated with sequence-specific sensing. We review the mechanisms underlying the detection and discrimination of foreign sequences from self. We also describe quantitative approaches used to infer the stimulatory capacity of a given pathogen nucleic acid species, and the influence of sequence-specific sensing on host-pathogen coevolution, including endogenous sequences of foreign origin. Finally, we speculate how further studies of sequence-specific sensing will be useful to improve vaccine design, gene therapy and cancer treatment.

Structure-Based Alignment and Consensus Secondary Structures for Three HIV-Related RNA Genomes

HIV and related primate lentiviruses possess single-stranded RNA genomes. Multiple regions of these genomes participate in critical steps in the viral replication cycle, and the functions of many RNA elements are dependent on the formation of defined structures. The structures of these elements are still not fully understood, and additional functional elements likely exist that have not been identified. In this work, we compared three full-length HIV-related viral genomes: HIV-1NL4-3, SIVcpz, and SIVmac (the latter two strains are progenitors for all HIV-1 and HIV-2 strains, respectively). Model-free RNA structure comparisons were performed using whole-genome structure information experimentally derived from nucleotide-resolution SHAPE reactivities. Consensus secondary structures were constructed for strongly correlated regions by taking into account both SHAPE probing structural data and nucleotide covariation information from structure-based alignments. In these consensus models, all known functional RNA elements were recapitulated with high accuracy. In addition, we identified multiple previously unannotated structural elements in the HIV-1 genome likely to function in translation, splicing and other replication cycle processes; these are compelling targets for future functional analyses. The structure-informed alignment strategy developed here will be broadly useful for efficient RNA motif discovery.

On the nucleotide composition and structure of retroviral RNA genomes

Retroviral RNA genomes display a rich variety in their nucleotide composition. For instance, the single-stranded RNA genome of human T cell leukemia virus (HTLV-1) is C-rich and G-poor and that of the human immunodeficiency virus (HIV-1) is A-rich and C-poor. Animal retroviruses add further variation to this unexplained, but many times remarkable virus-specific property. We previously described that the nucleotide bias is even more extreme in the unpaired regions of the structured HIV-1 RNA genome, which has been probed by SHAPE technology. We now document that the same trend is apparent for the MFold-predicted RNA structure of HIV-1 RNA and subsequently investigated the predicted structures of the RNA genomes of other retroviruses. We conclude that all virus-specific signatures are enhanced for the unpaired nucleotides in the RNA genome. Consequently, the differences in nucleotide count between the diverse human and animal retroviruses are further exposed in the single stranded genome regions. We used a skew analysis to visualize these striking differences in nucleotide usage. Evolutionary events responsible for these nucleotide signatures will be discussed.

Large-scale nucleotide optimization of simian immunodeficiency virus reduces its capacity to stimulate type I interferon in vitro

Lentiviral RNA genomes present a strong bias in their nucleotide composition with extremely high frequencies of A nucleotide in human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV). Based on the observation that human optimization of RNA virus gene fragments may abolish their ability to stimulate the type I interferon (IFN-I) response, we identified the most biased sequences along the SIV genome and showed that they are the most potent IFN-I stimulators. With the aim of designing an attenuated SIV genome based on a reduced capacity to activate the IFN-I response, we synthesized artificial SIV genomes whose biased sequences were optimized toward macaque average nucleotide composition without altering their regulatory elements or amino acid sequences. A synthetic SIV optimized with 169 synonymous mutations in gag and pol genes showed a 100-fold decrease in replicative capacity. Interestingly, a synthetic SIV optimized with 70 synonymous mutations in pol had a normal replicative capacity. Its ability to stimulate IFN-I was reduced when infected cells were cocultured with reporter cells. IFN regulatory factor 3 (IRF3) transcription factor was required for IFN-I stimulation, implicating cytosolic sensors in the detection of SIV-biased RNA in infected cells. No reversion of introduced mutations was observed for either of the optimized viruses after 10 serial passages. In conclusion, we have designed large-scale nucleotide-modified SIVs that may display attenuated pathogenic potential.

IMPORTANCE

In this study, we synthesized artificial SIV genomes in which the most hyperbiased sequences were optimized to bring them closer to the nucleotide composition of the macaque SIV host. Interestingly, we generated a stable synthetic SIV optimized with 70 synonymous mutations in pol gene, which had a normal replicative capacity but a reduced ability to stimulate type I IFN. This demonstrates the possibility to rationally change viral nucleotide composition to design replicative and genetically stable lentiviruses with attenuated pathogenic potentials.

Quantifying selection against synonymous mutations in HIV-1 env evolution

Intrapatient evolution of human immunodeficiency virus type 1 (HIV-1) is driven by the adaptive immune system resulting in rapid change of HIV-1 proteins. When cytotoxic CD8(+) T cells or neutralizing antibodies target a new epitope, the virus often escapes via nonsynonymous mutations that impair recognition. Synonymous mutations do not affect this interplay and are often assumed to be neutral. We test this assumption by tracking synonymous mutations in longitudinal intrapatient data from the C2-V5 part of the env gene. We find that most synonymous variants are lost even though they often reach high frequencies in the viral population, suggesting a cost to the virus. Using published data from SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) assays, we find that synonymous mutations that disrupt base pairs in RNA stems flanking the variable loops of gp120 are more likely to be lost than other synonymous changes: these RNA hairpins might be important for HIV-1. Computational modeling indicates that, to be consistent with the data, a large fraction of synonymous mutations in this genomic region need to be deleterious with a cost on the order of 0.002 per day. This weak selection against synonymous substitutions does not result in a strong pattern of conservation in cross-sectional data but slows down the rate of evolution considerably. Our findings are consistent with the notion that large-scale patterns of RNA structure are functionally relevant, whereas the precise base pairing pattern is not.

Changes in codon-pair bias of human immunodeficiency virus type 1 have profound effects on virus replication in cell culture

Background

Human immunodeficiency virus type 1 (HIV-1) has a biased nucleotide composition different from human genes. This raises the question of how evolution has chosen the nucleotide sequence of HIV-1 that is observed today, or to what extent the actual encoding contributes to virus replication capacity, evolvability and pathogenesis. Here, we applied the previously described synthetic attenuated virus engineering (SAVE) approach to HIV-1.

Results

Using synonymous codon pairs, we rationally recoded and codon pair–optimized and deoptimized different moieties of the HIV-1 gag and pol genes. Deoptimized viruses had significantly lower viral replication capacity in MT-4 and peripheral blood mononuclear cells (PBMCs). Varying degrees of ex vivo attenuation were obtained, depending upon both the specific deoptimized region and the number of deoptimized codons. A protease optimized virus carrying 38 synonymous mutations was not attenuated and displayed a replication capacity similar to that of the wild-type virus in MT-4 cells and PBMCs. Although attenuation is based on several tens of nucleotide changes, deoptimized HIV-1 reverted to wild-type virulence after serial passages in MT-4 cells. Remarkably, no reversion was observed in the optimized virus.

Conclusion

These data demonstrate that SAVE is a useful strategy to phenotypically affect the replicative properties of HIV-1.

The A-nucleotide preference of HIV-1 in the context of its structured RNA genome

A bipartition of HIV-1 RNA genome sequences into single- and double-stranded nucleotides is possible based on the secondary structure model of a complete 9 kb genome. Subsequent analysis revealed that the well-known lentiviral property of A-accumulation is profoundly present in single-stranded domains, yet absent in double-stranded domains. Mutational rate analysis by means of an unrestricted model of nucleotide substitution suggests the presence of an evolutionary equilibrium to preserve this biased nucleotide distribution.

HIV-1 suppressive sequences are modulated by Rev transport of unspliced RNA and are required for efficient HIV-1 production

The unspliced human immunodeficiency virus type 1 (HIV-1) RNAs are translated as Gag and Gag-Pol polyproteins or packaged as genomes into viral particles. Efficient translation is necessary before the transition to produce infective virions. The viral protein Rev exports all intron-containing viral RNAs; however, it also appears to enhance translation. Cellular microRNAs target cellular and viral mRNAs to silence their translation and enrich them at discrete cytoplasmic loci that overlap with the putative interim site of Gag and the genome. Here, we analyzed how Rev-mediated transport and the splicing status of the mRNA influenced the silencing status imposed by microRNA. Through identification and mutational analysis of the silencing sites in the HIV-1 genome, we elucidated the effect of silencing on virus production. Renilla luciferase mRNA, which contains a let-7 targeting site in its 3' untranslated region, was mediated when it was transported by Rev and not spliced, but it was either not mediated when it was spliced even in a partial way or it was Rev-independent. The silencing sites in the pol and env-nef regions of the HIV-1 genome, which were repressed in T cells and other cell lines, were Drosha-dependent and could also be modulated by Rev in an unspliced state. Mutant viruses that contained genomic mutations that reflect alterations to show more derepressive effects in the 3' untranslated region of the Renilla luciferase gene replicated more slowly than wild-type virus. These findings yield insights into the HIV-1 silencing sites that might allow the genome to avoid translational machinery and that might be utilized in coordinating virus production during initial virus replication. However, the function of Rev to modulate the silencing sites of unspliced RNAs would be advantageous for the efficient translation that is required to support protein production prior to viral packaging and particle production.

The biased nucleotide composition of the HIV genome: a constant factor in a highly variable virus

Viruses often deviate from their hosts in the nucleotide composition of their genomes. The RNA genome of the lentivirus family of retroviruses, including human immunodeficiency virus (HIV), contains e.g. an above average percentage of adenine (A) nucleotides, while being extremely poor in cytosine (C). Such a deviant base composition has implications for the amino acids that are encoded by the open reading frames (ORFs), both in the requirement of specific tRNA species and in the preference for amino acids encoded by e.g. A-rich codons. Nucleotide composition does obviously affect the secondary and tertiary structure of the RNA genome and its biological functions, but it does also influence phylogenetic analysis of viral genome sequences, and possibly the activity of the integrated DNA provirus. Over time, the nucleotide composition of the HIV-1 genome is exceptionally conserved, varying by less than 1% per base position per isolate within either group M, N, or O during 1983–2009. This extreme stability of the nucleotide composition may possibly be achieved by negative selection, perhaps conserving semi-stable RNA secondary structure as reverse transcription would be significantly affected for a less A-rich genome where secondary structures are expected to be more stable and thus more difficult to unfold.

This review will discuss all aspects of the lentiviral genome composition, both of the RNA and of its derived double-stranded DNA genome, with a focus on HIV-1, the nucleotide composition over time, the effects of artificially humanized codons as well as contributions of immune system pressure on HIV nucleotide bias.

HIV-1 modulates the tRNA pool to improve translation efficiency

Despite its poorly adapted codon usage, HIV-1 replicates and is expressed extremely well in human host cells. HIV-1 has recently been shown to package non-lysyl transfer RNAs (tRNAs) in addition to the tRNA(Lys) needed for priming reverse transcription and integration of the HIV-1 genome. By comparing the codon usage of HIV-1 genes with that of its human host, we found that tRNAs decoding codons that are highly used by HIV-1 but avoided by its host are overrepresented in HIV-1 virions. In particular, tRNAs decoding A-ending codons, required for the expression of HIV's A-rich genome, are highly enriched. Because the affinity of Gag-Pol for all tRNAs is nonspecific, HIV packaging is most likely passive and reflects the tRNA pool at the time of viral particle formation. Codon usage of HIV-1 early genes is similar to that of highly expressed host genes, but codon usage of HIV-1 late genes was better adapted to the selectively enriched tRNA pool, suggesting that alterations in the tRNA pool are induced late in viral infection. If HIV-1 genes are adapting to an altered tRNA pool, codon adaptation of HIV-1 may be better than previously thought.

Patterns of Human Immunodeficiency Virus type 1 recombination ex vivo provide evidence for coadaptation of distant sites, resulting in purifying selection for intersubtype recombinants during replication

High-frequency recombination is a hallmark of HIV-1 replication. Recombination can occur between two members of the same subtype or between viruses from two different subtypes, generating intra- or intersubtype recombinants, respectively. Many intersubtype recombinants have been shown to circulate in human populations. We hypothesize that sequence diversity affects the emergence of viable recombinants by decreasing recombination events and reducing the ability of the recombinants to replicate. To test our hypothesis, we compared recombination between two viruses containing subtype B pol genes (B/B) and between viruses with pol genes from subtype B or F (B/F). Recombination events generated during a single cycle of infection without selection pressure on pol gene function were analyzed by single-genome sequencing. We found that recombination occurred slightly ( approximately 30%) less frequently in B/F than in B/B viruses, and the overall distribution of crossover junctions in pol was similar for the two classes of recombinants. We then examined the emergence of recombinants in a multiple cycle assay, so that functional pol gene products were selected. We found that the emerging B/B recombinants had complex patterns, and the crossover junctions were distributed throughout the pol gene. In contrast, selected B/F recombinants had limited recombination patterns and restricted crossover junction distribution. These results provide evidence for the evolved coadapted sites in variants from different subtypes; these sites may be segregated by recombination events, causing the newly generated intersubtype recombinants to undergo purifying selection. Therefore, the ability of the recombinants to replicate is the major barrier for many of these viruses.

Accurately measuring recombination between closely related HIV-1 genomes

Retroviral recombination is thought to play an important role in the generation of immune escape and multiple drug resistance by shuffling pre-existing mutations in the viral population. Current estimates of HIV-1 recombination rates are derived from measurements within reporter gene sequences or genetically divergent HIV sequences. These measurements do not mimic the recombination occurring in vivo, between closely related genomes. Additionally, the methods used to measure recombination make a variety of assumptions about the underlying process, and often fail to account adequately for issues such as co-infection of cells or the possibility of multiple template switches between recombination sites. We have developed a HIV-1 marker system by making a small number of codon modifications in gag which allow recombination to be measured over various lengths between closely related viral genomes. We have developed statistical tools to measure recombination rates that can compensate for the possibility of multiple template switches. Our results show that when multiple template switches are ignored the error is substantial, particularly when recombination rates are high, or the genomic distance is large. We demonstrate that this system is applicable to other studies to accurately measure the recombination rate and show that recombination does not occur randomly within the HIV genome.