The external domains of the HIV-1 envelope are a mutational cold spot

Evolutionary Dynamics of Accelerated Antiviral Resistance Development in Hypermutator Herpesvirus

Antiviral therapy is constantly challenged by the emergence of resistant pathogens. At the same time, experimental approaches to understand and predict resistance are limited by long periods required for evolutionary processes. Here, we present a herpes simplex virus 1 mutant with impaired proofreading capacity and consequently elevated mutation rates. Comparing this hypermutator to parental wild type virus, we study the evolution of antiviral drug resistance in vitro. We model resistance development and elucidate underlying genetic changes against three antiviral substances. Our analyzes reveal no principle difference in the evolutionary behavior of both viruses, adaptive processes are overall similar, however significantly accelerated for the hypermutator. We conclude that hypermutator viruses are useful for modeling adaptation to antiviral therapy. They offer the benefit of expedited adaptation without introducing apparent bias and can therefore serve as an accelerator to predict natural evolution.

Mapping the Evolutionary Space of SARS-CoV-2 Variants to Anticipate Emergence of Subvariants Resistant to COVID-19 Therapeutics

New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with mutations in the spike glycoprotein. In most cases, the sublineage-defining mutations vary between the VOCs. It is unclear whether these differences reflect lineage-specific likelihoods for mutations at each spike position or the stochastic nature of their appearance. Here we show that SARS-CoV-2 lineages have distinct evolutionary spaces (a probabilistic definition of the sequence states that can be occupied by expanding virus subpopulations). This space can be accurately inferred from the patterns of amino acid variability at the whole-protein level. Robust networks of co-variable sites identify the highest-likelihood mutations in new VOC sublineages and predict remarkably well the emergence of subvariants with resistance mutations to COVID-19 therapeutics. Our studies reveal the contribution of low frequency variant patterns at heterologous sites across the protein to accurate prediction of the changes at each position of interest.

Influence of viral genome properties on polymerase fidelity

The first step of viral evolution takes place during genome replication via the error-prone viral polymerase. Among the mutants that arise through this process, only a few well-adapted variants will be selected by natural selection, renewing the viral genome population. Viral polymerase-mediated errors are thought to occur stochastically. However, accumulating evidence suggests that viral polymerase-mediated mutations are heterogeneously distributed throughout the viral genome. Here, we review work that supports this concept and provides mechanistic insights into how specific features of the viral genome could modulate viral polymerase-mediated errors. A predisposition to accumulate viral polymerase-mediated errors at specific loci in the viral genome may guide evolution to specific pathways, thus opening new directions of research to better understand viral evolutionary dynamics.

Host Molecule Incorporation into HIV Virions, Potential Influences in HIV Pathogenesis

During the last phase of HIV viral production, nascent HIV virions acquire a fraction of the cellular lipid membrane to create the external lipid envelope, a process by which cellular proteins present on the surface of the infected cell can be incorporated along with Env trimers. Interestingly, several studies indicated that these incorporated host molecules could conserve their biological activity and consequently contribute to HIV pathogenesis either by enhancing the infectivity of HIV virions, their tissue tropism or by affecting immune cell functions. The following review will describe the main approaches used to characterize membrane bound host molecule incorporation into HIV virions, the proposed mechanisms involved, and the role of a non-exhaustive list of incorporated molecules.

Mechanistic Interplay between HIV-1 Reverse Transcriptase Enzyme Kinetics and Host SAMHD1 Protein: Viral Myeloid-Cell Tropism and Genomic Mutagenesis

Human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) has been the primary interest among studies on antiviral discovery, viral replication kinetics, drug resistance, and viral evolution. Following infection and entry into target cells, the HIV-1 core disassembles, and the viral RT concomitantly converts the viral RNA into double-stranded proviral DNA, which is integrated into the host genome. The successful completion of the viral life cycle highly depends on the enzymatic DNA polymerase activity of RT. Furthermore, HIV-1 RT has long been known as an error-prone DNA polymerase due to its lack of proofreading exonuclease properties. Indeed, the low fidelity of HIV-1 RT has been considered as one of the key factors in the uniquely high rate of mutagenesis of HIV-1, which leads to efficient viral escape from immune and therapeutic antiviral selective pressures. Interestingly, a series of studies on the replication kinetics of HIV-1 in non-dividing myeloid cells and myeloid specific host restriction factor, SAM domain, and HD domain-containing protein, SAMHD1, suggest that the myeloid cell tropism and high rate of mutagenesis of HIV-1 are mechanistically connected. Here, we review not only HIV-1 RT as a key antiviral target, but also potential evolutionary and mechanistic crosstalk among the unique enzymatic features of HIV-1 RT, the replication kinetics of HIV-1, cell tropism, viral genetic mutation, and host SAMHD1 protein.

Mutability Patterns Across the Spike Glycoprotein Reveal the Diverging and Lineage-specific Evolutionary Space of SARS-CoV-2

Mutations in the spike glycoprotein of SARS-CoV-2 allow the virus to probe the sequence space in search of higher-fitness states. New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with such mutations. Interestingly, the sites of mutation in these sublineages vary between the VOCs. Whether such differences reflect the random nature of mutation appearance or distinct evolutionary spaces of spike in the VOCs is unclear. Here we show that each position of spike has a lineage-specific likelihood for mutations to appear and dominate descendent sublineages. This likelihood can be accurately estimated from the lineage-specific mutational profile of spike at a protein-wide level. The mutability environment of each position, including adjacent sites on the protein structure and neighboring sites on the network of comutability, accurately forecast changes in descendent sublineages. Mapping of imminent changes within the VOCs can contribute to the design of immunogens and therapeutics that address future forms of SARS-CoV-2.

Mutational spectrum of hepatitis C virus in patients with chronic hepatitis C determined by single molecule real-time sequencing

The emergence of hepatitis C virus (HCV) with resistance-associated substitution (RAS), produced by mutations in the HCV genome, is a major problem in direct acting antivirals (DAA) treatment. This study aimed to clarify the mutational spectrum in HCV-RNA and the substitution pattern for the emergence of RASs in patients with chronic HCV infection. HCV-RNA from two HCV replicon cell lines and the serum HCV-RNA of four non-liver transplant and four post-liver transplant patients with unsuccessful DAA treatment were analyzed using high-accuracy single-molecule real-time long-read sequencing. Transition substitutions, especially A>G and U>C, occurred prominently under DAAs in both non-transplant and post-transplant patients, with a mutational bias identical to that occurring in HCV replicon cell lines during 10-year culturing. These mutational biases were reproduced in natural courses after DAA treatment. RASs emerged via both transition and transversion substitutions. NS3-D168 and NS5A-L31 RASs resulted from transversion mutations, while NS5A-Y93 RASs was caused by transition substitutions. The fidelity of the RNA-dependent RNA polymerase, HCV-NS5B, produces mutational bias in the HCV genome, characterized by dominant transition mutations, notably A>G and U>C substitutions. However, RASs are acquired by both transition and transversion substitutions, and the RASs-positive HCV clones are selected and proliferated under DAA treatment pressure.

Collaborative Mining of Whole Genome Sequences for Intelligent HIV-1 Sub-Strain(s) Discovery

BACKGROUND

Effective global antiretroviral vaccines and therapeutic strategies depend on the diversity, evolution, and epidemiology of their various strains as well as their transmission and pathogenesis. Most viral disease-causing particles are clustered into a taxonomy of subtypes to suggest pointers toward nucleotide-specific vaccines or therapeutic applications of clinical significance sufficient for sequence-specific diagnosis and homologous viral studies. These are very useful to formulate predictors to induce cross-resistance to some retroviral control drugs being used across study areas.

OBJECTIVE

This research proposed a collaborative framework of hybridized (Machine Learning and Natural Language Processing) techniques to discover hidden genome patterns and feature predictors, for HIV-1 genome sequences mining.

METHOD

630 human HIV-1 genome sequences above 8500 bps were excavated from the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov) for 21 countries across different continents, Antarctica exempt. These sequences were transformed and learned using a self-organizing map (SOM). To discriminate emerging/new sub-strain(s), the HIV-1 reference genome was included as part of the input isolates/samples during the training. After training the SOM, component planes defining pattern clusters of the input datasets were generated, for cognitive knowledge mining and subsequent labelling of the datasets. Additional genome features including dinucleotide transmission recurrences, codon recurrences, and mutation recurrences, were finally extracted from the raw genomes to construct output classification targets for supervised learning.

RESULTS

SOM training explains the inherent pattern diversity of HIV-1 genomes as well as inter- and intra-country transmissions in which mobility might play an active role, as corroborated by literature. Nine sub-strains were discovered after disassembling the SOM correlation hunting matrix space attributed to disparate clusters. Cognitive knowledge mining separated similar pattern clusters bounded by a certain degree of correlation range, discovered by the SOM. A Kruskal-Wallis rank-sum test and Wilcoxon rank-sum test showed statistically significant variations in dinucleotide, codon, and mutation patterns.

CONCLUSION

Results of the discovered sub-strains and response clusters visualizations corroborate existing literature, with significant haplotype variations. The proposed framework would assist in the development of decision support systems for easy contact tracing, infectious disease surveillance, and studying the progressive evolution of the reference HIV-1 genome.

Phylodynamic Study of the Conserved RNA Structure Encompassing the Hemagglutinin Cleavage Site Encoding Region of H5 and H7 Low Pathogenic Avian Influenza Viruses

Highly Pathogenic Avian Influenza Viruses (HPAIV) evolve from Low Pathogenic Avian Influenza Viruses (LPAIV) of the H5 and H7 subtypes. This evolution is characterized by the acquisition of a multi-basic cleavage site (MBCS) motif in the hemagglutinin (HA) that leads to an extended viral tropism and severe disease in poultry. One key unanswered question is whether the risk of transition to HPAIV is similar for all LPAIV H5 or H7 strains, or whether specific determinants in the HA sequence of some H5 or H7 LPAIV strains correlate with a higher risk of transition to HPAIV. Here we determined if specific features of the conserved RNA stem loop located at the hemagglutinin cleavage site-encoding region could be detected along the LPAIV to HPAIV evolutionary pathway. Analysis of the thermodynamic stability of the predicted RNA structures showed no specific patterns common to HA sequences leading to HPAIV and distinct from those remaining LPAIV. However, RNA structure clustering analysis revealed that most of the American lineage ancestors leading to H7 emergences via recombination shared the same vRNA structure topology at the HA1/HA2 boundary region. Our study thus identified predicted secondary RNA structures present in the HA of H7 viruses, which could promote genetic recombination and acquisition of a MBCS.

HIV-1 Capsid Core: A Bullet to the Heart of the Target Cell

The first step of the intracellular phase of retroviral infection is the release of the viral capsid core in the cytoplasm. This structure contains the viral genetic material that will be reverse transcribed and integrated into the genome of infected cells. Up to recent times, the role of the capsid core was considered essentially to protect this genetic material during the earlier phases of this process. However, increasing evidence demonstrates that the permanence inside the cell of the capsid as an intact, or almost intact, structure is longer than thought. This suggests its involvement in more aspects of the infectious cycle than previously foreseen, particularly in the steps of viral genomic material translocation into the nucleus and in the phases preceding integration. During the trip across the infected cell, many host factors are brought to interact with the capsid, some possessing antiviral properties, others, serving as viral cofactors. All these interactions rely on the properties of the unique component of the capsid core, the capsid protein CA. Likely, the drawback of ensuring these multiple functions is the extreme genetic fragility that has been shown to characterize this protein. Here, we recapitulate the busy agenda of an HIV-1 capsid in the infectious process, in particular in the light of the most recent findings.

The Case for S2: The Potential Benefits of the S2 Subunit of the SARS-CoV-2 Spike Protein as an Immunogen in Fighting the COVID-19 Pandemic

As COVID-19 cases continue to rise, it is imperative to learn more about antibodies and T-cells produced against the causative virus, SARS-CoV-2, in order to guide the rapid development of therapies and vaccines. While much of the current antibody and vaccine research focuses on the receptor-binding domain of S1, a less-recognized opportunity is to harness the potential benefits of the more conserved S2 subunit. Similarities between the spike proteins of both SARS-CoV-2 and HIV-1 warrant exploring S2. Possible benefits of employing S2 in therapies and vaccines include the structural conservation of S2, extant cross-reactive neutralizing antibodies in populations (due to prior exposure to common cold coronaviruses), the steric neutralization potential of antibodies against S2, and the stronger memory B-cell and T-cell responses. More research is necessary on the effect of glycans on the accessibility and stability of S2, SARS-CoV-2 mutants that may affect infectivity, the neutralization potential of antibodies produced by memory B-cells, cross-reactive T-cell responses, antibody-dependent enhancement, and antigen competition. This perspective aims to highlight the evidence for the potential advantages of using S2 as a target of therapy or vaccine design.

Spontaneous Mutations in HIV-1 Gag, Protease, RT p66 in the First Replication Cycle and How They Appear: Insights from an In Vitro Assay on Mutation Rates and Types

While drug resistant mutations in HIV-1 are largely credited to its error prone HIV-1 RT, the time point in the infection cycle that these mutations can arise and if they appear spontaneously without selection pressures both remained enigmatic. Many HIV-1 RT mutational in vitro studies utilized reporter genes (LacZ) as a template to investigate these questions, thereby not accounting for the possible contribution of viral codon usage. To address this gap, we investigated HIV-1 RT mutation rates and biases on its own Gag, protease, and RT p66 genes in an in vitro selection pressure free system. We found rare clinical mutations with a general avoidance of crucial functional sites in the background mutations rates for Gag, protease, and RT p66 at 4.71 × 10-5, 6.03 × 10-5, and 7.09 × 10-5 mutations/bp, respectively. Gag and p66 genes showed a large number of 'A to G' mutations. Comparisons with silently mutated p66 sequences showed an increase in mutation rates (1.88 × 10-4 mutations/bp) and that 'A to G' mutations occurred in regions reminiscent of ADAR neighbor sequence preferences. Mutational free energies of the 'A to G' mutations revealed an avoidance of destabilizing effects, with the natural p66 gene codon usage providing barriers to disruptive amino acid changes. Our study demonstrates the importance of studying mutation emergence in HIV genes in a RT-PCR in vitro selection pressure free system to understand how fast drug resistance can emerge, providing transferable applications to how new viral diseases and drug resistances can emerge.

Evolutive Temporal Footprint of an HIV-1 Envelope Protein in an Epidemiologically Linked Cluster

Background:

The C2V3C3 region of gp 120, encoded by the HIV-1 envelope gene (env), is an important antigenic target, a key determinant for viral evolution and essential for determining epitopes for vaccines.

Methods:

The relationships among genetic sequence diversity, selective pressure, constraints on HIV-1 envelope protein were explored and also correlated this analysis with information entropy; hypermutation; HIV tropism; CD4+ T cell counts or HIV viral load. A total of 179 HIV-1 C2V3C3 sequences derived from cell-free plasma were used, determined from serial samples, in four epidemiologically linked individuals (one infected blood donor, two transfusion recipients and a sexual partner infected by one of the recipients) over a maximum period of 8 years. This study is important because it considers the analysis of patterns in genomic sequences, without drugs and over time.

Results:

A temporal relationship among information entropy, hypermutation, tropism switch, viral load, and CD4+ T cell count was determined. Changes in information entropy were time-dependent, and an increase in entropy was observed in the C2V3C3 region at amino acids G313 and F317-I320 (related to the GPGR-motif and coreceptor tropism), and at amino acids A281 in C2 and A346 in C3, related to immune escape.

Conclusion:

The increase of information entropy over time was correlated with hypermutation and the emergence of nonR5- strains, which are both associated with more variable genomes.

The frequency of defective genes in vif and vpr genes in 20 hemophiliacs is associated with Korean Red Ginseng and highly active antiretroviral therapy: the impact of lethal mutations in vif and vpr genes on HIV-1 evolution

Background

We have reported that internal deletions in the nef, gag, and pol genes in HIV-1–infected patients are induced in those treated with Korean Red Ginseng (KRG). KRG delays the development of resistance mutations to antiretroviral drugs.

Methods

The vif-vpr genes over 26 years in 20 hemophiliacs infected with HIV-1 from a single source were sequenced to investigate whether vif-vpr genes were affected by KRG and KRG plus highly active antiretroviral therapy (ART) (hereafter called GCT) and compared the results with our previous data.

Results

A significantly higher number of in-frame small deletions were found in the vif-vpr genes of KRG-treated patients than at the baseline, in control patients, and in ART-alone patients (p < 0.001). These were significantly reduced in GCT patients (p < 0.05). In contrast, sequences harboring a premature stop codon (SC) were more significant in GCT patients (10.1%) than in KRG-alone patients, control (p < 0.01), and ART-alone patients (p = 0.078 for peripheral blood mononuclear cells). The proportion of SC in Vpr was similar to that in Vif, whereas the proportion of sequences revealing SC in the env-nef genes was significantly lower than that in the pol-vif-vpr genes (p < 0.01). The genetic distance was 1.8 times higher in the sequences harboring SC than in the sequences without SC (p < 0.001). Q135P in the vif gene is significantly associated with rapid progression to AIDS (p < 0.01).

Conclusion

Our data show that KRG might induce sΔ in the vif-vpr genes and that vif-vpr genes are similarly affected by lethal mutations.

The Determination of HIV-1 RT Mutation Rate, Its Possible Allosteric Effects, and Its Implications on Drug Resistance

The high mutation rate of the human immunodeficiency virus type 1 (HIV-1) plays a major role in treatment resistance, from the development of vaccines to therapeutic drugs. In addressing the crux of the issue, various attempts to estimate the mutation rate of HIV-1 resulted in a large range of 10⁻⁵–10⁻³ errors/bp/cycle due to the use of different types of investigation methods. In this review, we discuss the different assay methods, their findings on the mutation rates of HIV-1 and how the locations of mutations can be further analyzed for their allosteric effects to allow for new inhibitor designs. Given that HIV is one of the fastest mutating viruses, it serves as a good model for the comprehensive study of viral mutations that can give rise to a more horizontal understanding towards overall viral drug resistance as well as emerging viral diseases.

Chromatin maturation of the HIV-1 provirus in primary resting CD4+ T cells

Human immunodeficiency virus type 1 (HIV-1) infection is a chronic condition, where viral DNA integrates into the genome. Latently infected cells form a persistent, heterogeneous reservoir that at any time can reactivate the integrated HIV-1. Here we confirmed that latently infected cells from HIV-1 positive study participants exhibited active HIV-1 transcription but without production of mature spliced mRNAs. To elucidate the mechanisms behind this we employed primary HIV-1 latency models to study latency establishment and maintenance. We characterized proviral transcription and chromatin development in cultures of resting primary CD4⁺ T-cells for four months after ex vivo HIV-1 infection. As heterochromatin (marked with H3K9me3 or H3K27me3) gradually stabilized, the provirus became less accessible with reduced activation potential. In a subset of infected cells, active marks (e.g. H3K27ac) and elongating RNAPII remained detectable at the latent provirus, despite prolonged proviral silencing. In many aspects, latent HIV-1 resembled an active enhancer in a subset of resting cells. The enhancer chromatin actively promoted latency and the enhancer-specific CBP/P300-inhibitor GNE049 was identified as a new latency reversal agent. The division of the latent reservoir according to distinct chromatin compositions with different reactivation potential enforces the notion that even though a relatively large set of cells contains the HIV-1 provirus, only a discrete subset is readily able to reactivate the provirus and spread the infection.

HIV infection is a devastating disease affecting 35 million people worldwide. Current anti-retroviral treatment is highly effective and has made the HIV infection chronic. However, despite more effective treatments, the prospects of a cure are distant. The problem for an HIV cure is that, even though the virus particles are eradicated, the infected cells maintain the information of remake the virus. This information is integrated in the host cell as a provirus. The provirus switches between active and inactive states. Thereby, the infected cells evade both the immune system and death associated with massive viral production. We have characterized the composition of proviral chromatin and how it connects with transcription and viral production. In resting primary CD4⁺ T-cells, we follow the fate of the provirus starting at infection until latency is firmly established. Only in a fraction of intact proviruses were we able to reverse latency and that this was highly regulated by the chromatin composition. Whereas the proviruses encompassed in heterochromatin were refractory to activation, latent proviruses with “enhancer” characteristics were readily activated. Our study provides key insights as to detect the remaining HIV-1 infected cells capable of reseeding the infection, and the mechanisms whereby they are maintained.

Reviewing HIV-1 Gag Mutations in Protease Inhibitors Resistance: Insights for Possible Novel Gag Inhibitor Designs

HIV protease inhibitors against the viral protease are often hampered by drug resistance mutations in protease and in the viral substrate Gag. To overcome this drug resistance and inhibit viral maturation, targeting Gag alongside protease rather than targeting protease alone may be more efficient. In order to successfully inhibit Gag, understanding of its drug resistance mutations and the elicited structural changes on protease binding needs to be investigated. While mutations on Gag have already been mapped to protease inhibitor resistance, there remain many mutations, particularly the non-cleavage mutations, that are not characterized. Through structural studies to unravel how Gag mutations contributes to protease drug resistance synergistically, it is thus possible to glean insights to design novel Gag inhibitors. In this review, we discuss the structural role of both novel and previously reported Gag mutations in PI resistance, and how new Gag inhibitors can be designed.

APOBEC3G Regulation of the Evolutionary Race Between Adaptive Immunity and Viral Immune Escape Is Deeply Imprinted in the HIV Genome

APOBEC3G (A3G) is a host enzyme that mutates the genomes of retroviruses like HIV. Since A3G is expressed pre-infection, it has classically been considered an agent of innate immunity. We and others previously showed that the impact of A3G-induced mutations on the HIV genome extends to adaptive immunity also, by generating cytotoxic T cell (CTL) escape mutations. Accordingly, HIV genomic sequences encoding CTL epitopes often contain A3G-mutable “hotspot” sequence motifs, presumably to channel A3G action toward CTL escape. Here, we studied the depths and consequences of this apparent viral genome co-evolution with A3G. We identified all potential CTL epitopes in Gag, Pol, Env, and Nef restricted to several HLA class I alleles. We simulated A3G-induced mutations within CTL epitope-encoding sequences, and flanking regions. From the immune recognition perspective, we analyzed how A3G-driven mutations are predicted to impact CTL-epitope generation through modulating proteasomal processing and HLA class I binding. We found that A3G mutations were most often predicted to result in diminishing/abolishing HLA-binding affinity of peptide epitopes. From the viral genome evolution perspective, we evaluated enrichment of A3G hotspots at sequences encoding CTL epitopes and included control sequences in which the HIV genome was randomly shuffled. We found that sequences encoding immunogenic epitopes exhibited a selective enrichment of A3G hotspots, which were strongly biased to translate to non-synonymous amino acid substitutions. When superimposed on the known mutational gradient across the entire length of the HIV genome, we observed a gradient of A3G hotspot enrichment, and an HLA-specific pattern of the potential of A3G hotspots to lead to CTL escape mutations. These data illuminate the depths and extent of the co-evolution of the viral genome to subvert the host mutator A3G.

Evolutionary analysis of the Chikungunya virus epidemic in Mexico reveals intra-host mutational hotspots in the E1 protein

The epidemic potential of Chikungunya virus (CHIKV) was recently made evident by its introduction and rapid expansion in the Caribbean and the Americas. We sought to gain a detailed understanding of the dynamics of the epidemic in Mexico, the country with the highest number of confirmed CHIKV cases in the Americas, and to characterise viral evolution at the population and intra-host level. Analysis of the spatiotemporal distribution of 2,739 diagnosed cases in Mexico from December 2014 to December 2015 showed a rapid nationwide expansion of the epidemic with focalisation in the South West of the country. We sequenced the envelope glycoprotein 1 gene (E1) from 25 patients using the Illumina MiSeq platform and report synonymous and non-synonymous consensus mutations. Bayesian phylogenetic analysis using 249 Asian lineage E1 sequences gave updated estimates of nucleotide substitution rates for E1 and time to most recent common ancestor of major lineages. The analysis indicates phylogenetically-related emergent Latin American clusters in South Western Mexico, Nicaragua and Honduras and transmission of American strains in the Pacific islands. Detailed analysis showed that intra-host changes in E1 mainly occurred in two variable regions (E1:189-220 and E1:349-358) in domains II and III, respectively, in residues involved in inter and intra-envelope spike interactions. At the population level, this study sheds light on the introduction and evolutionary dynamics of CHIKV in the Americas. At the intra-host level, this study identifies mutational hotspots of the E1 protein with implications for understanding the relationship between the CHIKV quasispecies, viral fitness and pathogenesis.

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.

Complexities of Viral Mutation Rates

Many viruses evolve rapidly. This is due, in part, to their high mutation rates. Mutation rate estimates for over 25 viruses are currently available. Here, we review the population genetics of virus mutation rates. We specifically cover the topics of mutation rate estimation, the forces that drive the evolution of mutation rates, and how the optimal mutation rate can be context-dependent.

Viral subversion of APOBEC3s: Lessons for anti-tumor immunity and tumor immunotherapy

APOBEC3s (A3) are endogenous DNA-editing enzymes that are expressed in immune cells including T lymphocytes. A3s target and mutate the genomes of retroviruses that infect immune tissues such as the human immunodeficiency virus (HIV). Therefore, A3s were classically defined as host anti-viral innate immune factors. In contrast, we and others showed that A3s can also benefit the virus by mediating escape from adaptive immune recognition and drugs. Crucially, whether A3-mediated mutations help or hinder HIV, is not up to chance. Rather, the virus has evolved multiple mechanisms to actively and maximally subvert A3 activity. More recently, extensive A3 mutational footprints in tumor genomes have been observed in many different cancers. This suggests a role for A3s in cancer initiation and progression. On the other hand, multiple anti-tumor activities of A3s have also come to light, including impact on immune checkpoint molecules and possible generation of tumor neo-antigens. Here, we review the studies that reshaped the view of A3s from anti-viral innate immune agents to host factors exploited by HIV to escape from immune recognition. Viruses and tumors share many attributes, including rapid evolution and adeptness at exploiting mutations. Given this parallel, we then discuss the pro- and anti-tumor roles of A3s, and suggest that lessons learned from studying A3s in the context of anti-viral immunity can be applied to tumor immunotherapy.

Time Intervals in Sequence Sampling, Not Data Modifications, Have a Major Impact on Estimates of HIV Escape Rates

The ability of human immunodeficiency virus (HIV) to avoid recognition by humoral and cellular immunity (viral escape) is well-documented, but the strength of the immune response needed to cause such a viral escape remains poorly quantified. Several previous studies observed a more rapid escape of HIV from CD8 T cell responses in the acute phase of infection compared to chronic infection. The rate of HIV escape was estimated with the help of simple mathematical models, and results were interpreted to suggest that CD8 T cell responses causing escape in acute HIV infection may be more efficient at killing virus-infected cells than responses that cause escape in chronic infection, or alternatively, that early escapes occur in epitopes mutations in which there is minimal fitness cost to the virus. However, these conclusions were challenged on several grounds, including linkage and interference of multiple escape mutations due to a low population size and because of potential issues associated with modifying the data to estimate escape rates. Here we use a sampling method which does not require data modification to show that previous results on the decline of the viral escape rate with time since infection remain unchanged. However, using this method we also show that estimates of the escape rate are highly sensitive to the time interval between measurements, with longer intervals biasing estimates of the escape rate downwards. Our results thus suggest that data modifications for early and late escapes were not the primary reason for the observed decline in the escape rate with time since infection. However, longer sampling periods for escapes in chronic infection strongly influence estimates of the escape rate. More frequent sampling of viral sequences in chronic infection may improve our understanding of factors influencing the rate of HIV escape from CD8 T cell responses.

Population Diversity and Collective Interactions during Influenza Virus Infection

Influenza A virus (IAV) continues to pose an enormous and unpredictable global public health threat, largely due to the continual evolution of escape from preexisting immunity and the potential for zoonotic emergence. Understanding how the unique genetic makeup and structure of IAV populations influences their transmission and evolution is essential for developing more-effective vaccines, therapeutics, and surveillance capabilities. Owing to their mutation-prone replicase and unique genome organization, IAV populations exhibit enormous amounts of diversity both in terms of sequence and functional gene content. Here, I review what is currently known about the genetic and genomic diversity present within IAV populations and how this diversity may shape the replicative and evolutionary dynamics of these viruses.

The substitution rate of HIV-1 subtypes: a genomic approach

HIV-1M causes most infections in the AIDS pandemic. Its genetic diversity is defined by nine pure subtypes and more than sixty recombinant forms. We have performed a comparative analysis of the evolutionary rate of five pure subtypes (A1, B, C, D, and G) and two circulating recombinant forms (CRF01_AE and CRF02 AG) using data obtained from nearly complete genome coding sequences. Times to the most recent common ancestor (tMRCA) and substitution rates of these HIV genomes, and their genomic partitions, were estimated by Bayesian coalescent analyses. Genomic substitution rate estimates were compared between the HIV-1 datasets analyzed by means of randomization tests. Significant differences in the rate of evolution were found between subtypes, with subtypes C and A1 and CRF01_AE displaying the highest rates. On the other hand, CRF02_AG and subtype D were the slowest evolving types. Using a different molecular clock model for each genomic partition led to more precise tMRCA estimates than when linking the same clock along the HIV genome. Overall, the earliest tMRCA corresponded to subtype A1 (median = 1941, 95% HPD = 1943-55), whereas the most recent tMRCA corresponded to subtype G and CRF01_AE subset 3 (median = 1971, 95% HPD = 1967-75 and median = 1972, 95% HPD = 1970-75, respectively). These results suggest that both biological and epidemiological differences among HIV-1M subtypes are reflected in their evolutionary dynamics. The estimates obtained for tMRCAs and substitution rates provide information that can be used as prior distributions in future Bayesian coalescent analyses of specific HIV-1 subtypes/CRFs and genes.

Detection of Replication Competent Lentivirus Using a qPCR Assay for VSV-G

Lentiviral vectors are a common tool used to introduce new and corrected genes into cell therapy products for treatment of human diseases. Although lentiviral vectors are ideal for delivery and stable integration of genes of interest into the host cell genome, they potentially pose risks to human health, such as integration-mediated transformation and generation of a replication competent lentivirus (RCL) capable of infecting non-target cells. In consideration of the latter risk, all cell-based products modified by lentiviral vectors and intended for patient use must be tested for RCL prior to treatment of the patient. Current Food and Drug Administration (FDA) guidelines recommend use of cell-based assays to this end, which can take up to 6 weeks for results. However, qPCR-based assays are a quick alternative for rapid assessment of RCL in products intended for fresh infusion. We describe here the development and qualification of a qPCR assay based on detection of envelope gene sequences (vesicular stomatitis virus G glycoprotein [VSV-G]) for RCL in accordance with Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines. Our results demonstrate the sensitivity, linearity, specificity, and reproducibility of detection of VSV-G sequences, with a low false-positive rate. These procedures are currently being used in our phase 1 clinical investigations.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

Concave gold bipyramids bound with multiple high-index facets: improved Raman and catalytic activities

Concave nanocrystals usually exhibit a large electromagnetic-field enhancement and superior catalytic performance due to their sharp corners, negative curvature and high-index facets. Conventional gold bipyramids (AuBPs) possess intriguing plasmonic properties which are attractive for various applications while the surface curvature of the reported bipyramids has not been fine-tuned to concave or convex structures to date. Additionally, the longitudinal surface plasmon resonance (LSPR) wavelengths of conventional AuBPs are mostly located in the range of 650-1350 nm and the sizes of these nanoparticles are usually not beyond 350 nm, which are not facilitated to some nano-focusing and nanophotonic applications. Herein, we reported a facile and robust approach for fabricating concave AuBPs (CAuBPs) with multiple high-index facets which are distinct from the conventional AuBPs and nanojavelin structures. The length of the as-prepared CAuBPs can even extend up to 800 nm. The CAuBP nanoparticles exhibit a strikingly pronounced broader plasmonic tuning range (even exceeding 1800 nm) and provide much higher electromagnetic-field enhancements in comparison to the conventional AuBPs, which broaden the promising applications of CAuBPs for many single-particle analyses. More importantly, the surface-enhanced Raman scattering (SERS) signals of CAuBPs on the single-particle or aqueous solution both displayed an enhanced intensity compared to conventional AuBPs. The CAuBP nanoparticles also exhibited improved catalytic activity due to the incredible abundance of uncoordinated atoms as active sites.

Survivors Remorse: antibody-mediated protection against HIV-1

It is clear that antibodies can play a pivotal role in preventing the transmission of HIV-1 and large efforts to identify an effective antibody-based vaccine to quell the epidemic. Shortly after HIV-1 was discovered as the cause of AIDS, the search for epitopes recognized by neutralizing antibodies became the driving strategy for an antibody-based vaccine. Neutralization escape variants were discovered shortly thereafter, and, after almost three decades of investigation, it is now known that autologous neutralizing antibody responses and their selection of neutralization resistant HIV-1 variants can lead to broadly neutralizing antibodies in some infected individuals. This observation drives an intensive effort to identify a vaccine to elicit broadly neutralizing antibodies. In contrast, there has been less systematic study of antibody specificities that must rely mainly or exclusively on other protective mechanisms, although non-human primate (NHP) studies as well as the RV144 vaccine trial indicate that non-neutralizing antibodies can contribute to protection. Here we propose a novel strategy to identify new epitope targets recognized by these antibodies for which viral escape is unlikely or impossible.

Mechanisms of viral mutation

The remarkable capacity of some viruses to adapt to new hosts and environments is highly dependent on their ability to generate de novo diversity in a short period of time. Rates of spontaneous mutation vary amply among viruses. RNA viruses mutate faster than DNA viruses, single-stranded viruses mutate faster than double-strand virus, and genome size appears to correlate negatively with mutation rate. Viral mutation rates are modulated at different levels, including polymerase fidelity, sequence context, template secondary structure, cellular microenvironment, replication mechanisms, proofreading, and access to post-replicative repair. Additionally, massive numbers of mutations can be introduced by some virus-encoded diversity-generating elements, as well as by host-encoded cytidine/adenine deaminases. Our current knowledge of viral mutation rates indicates that viral genetic diversity is determined by multiple virus- and host-dependent processes, and that viral mutation rates can evolve in response to specific selective pressures.

Buffering deleterious polymorphisms in highly constrained parts of HIV-1 envelope by flexible regions

Background

Covariation is an essential process that leads to coevolution of parts of proteins and genomes. In organisms subject to strong selective pressure, coevolution is central to keep the balance between the opposite requirements of antigenic variation and retention of functionality. Being the viral component most exposed to the external environment, the HIV-1 glycoprotein gp120 constitutes the main target of the immune response. Accordingly its more external portions are characterised by extensive sequence heterogeneity fostering constant antigenic variation.

Results

We report that a single polymorphism, present at the level of the viral population in the conserved internal region C2, was sufficient to totally abolish Env functionality when introduced in an exogenous genetic context. The prominent defect of the non-functional protein is a block occurring after recognition of the co-receptor CCR5, likely due to an interference with the subsequent conformational changes that lead to membrane fusion. We also report that the presence of compensatory polymorphisms at the level of the external and hypervariable region V3 fully restored the functionality of the protein. The functional revertant presents different antigenic profiles and sensitivity to the entry inhibitor TAK 779.

Conclusions

Our data suggest that variable regions, besides harbouring intrinsic extensive antigenic diversity, can also contribute to sequence diversification in more structurally constrained parts of the gp120 by buffering the deleterious effect of polymorphisms, further increasing the genetic flexibility of the protein and the antigenic repertoire of the viral population.

A step forward understanding HIV-1 diversity

Human immunodeficiency virus (HIV) populations are characterized by extensive genetic diversity. Antigenic diversification is essential for escape from immune selection and therapy, and remains one of the major obstacles for the development of an efficient vaccine strategy. Even if intensive efforts have been made for understanding the molecular mechanisms responsible for genetic diversity in HIV, conclusive data in vivo is still lacking. Recent works have addressed this issue, focusing on the identification of the sources of genetic diversity during in vivo infections and on the estimate of the pervasiveness of genetic recombination during replication in vivo. Surprisingly, it appears that despite the error-prone nature of the viral polymerase, the bulk of mutations found in patients are indeed due to the effect of a cellular restriction factor. This factor tends to hypermutate the viral genome abolishing viral infectivity. When hypermutation is incomplete, the virus retains infectivity and converts the effect of the cellular factor to its advantage by exploiting it to generate genetic diversity that is beneficial for viral propagation. This view contrasts the long-standing dogma that viral diversity is due to the intrinsic error-prone nature of the viral replication cycle. Besides hypermutations and mutations, recombination is also a pervasive source of genetic diversity. The estimate of the frequency at which this process takes place in vivo has remained elusive, despite extensive efforts in this sense. Now, using single genome amplification, and starting from publically available datasets, it has been obtained a confirmation of the estimates previously made using tissue culture studies. These recent findings are presented here and their implications for the development of future researches are discussed.

Highly heterogeneous mutation rates in the hepatitis C virus genome

Spontaneous mutations are the ultimate source of genetic variation and have a prominent role in evolution. RNA viruses such as hepatitis C virus (HCV) have extremely high mutation rates, but these rates have been inferred from a minute fraction of genome sites, limiting our view of how RNA viruses create diversity. Here, by applying high-fidelity ultradeep sequencing to a modified replicon system, we scored >15,000 spontaneous mutations, encompassing more than 90% of the HCV genome. This revealed >1,000-fold differences in mutability across genome sites, with extreme variations even between adjacent nucleotides. We identify base composition, the presence of high- and low-mutation clusters and transition/transversion biases as the main factors driving this heterogeneity. Furthermore, we find that mutability correlates with the ability of HCV to diversify in patients. These data provide a site-wise baseline for interrogating natural selection, genetic load and evolvability in HCV, as well as for evaluating drug resistance and immune evasion risks.