HIV evolution in early infection: selection pressures, patterns of insertion and deletion, and the impact of APOBEC

HIV-1 Vif global diversity and possible APOBEC-mediated response since 1980

HIV-1 Vif's principal function is to counter the antiretroviral activities of DNA-editing APOBEC3 cytidine deaminases. Unconstrained APOBEC3 activity introduces premature stop codons in HIV-1 genes and can lead to viral inactivation. To investigate the evolution and diversification of Vif over the HIV-1 pandemic and document evidence of APOBEC3-mediated pressure, we analyzed 4612 publicly available sequences derived from 10 dominant subtypes and circulating recombinant forms (CRFs) using the Hervé platform. We found widespread evidence of diversifying selection that was convergent across subtypes and CRFs, but remarkable stability in consensus sequences over time. Divergence and selection did not favor APOBEC3-interacting sites. We furthermore found that APOBEC3-induced substitutions in env and gag-pol genes increased over time and were positively associated with vif diversity. These results suggest that APOBEC3-driven adaptation in Vif is relatively rare and that permissiveness to human APOBEC3-induced substitution as a mechanism for generating diversity may be advantageous to HIV-1 evolution.

Archived HIV-1 Drug Resistance Mutations: Role of Proviral HIV-1 DNA Genotype for the Management of Virological Responder People Living with HIV

Despite its effectiveness in controlling plasma viremia, antiretroviral therapy (ART) cannot target proviral DNA, which remains an obstacle to HIV-1 eradication. When treatment is interrupted, the reservoirs can act as a source of viral rebound, highlighting the value of proviral DNA as an additional source of information on an individual's overall resistance burden. In cases where the viral load is too low for successful HIV-1 RNA genotyping, HIV-1 DNA can help identify resistance mutations in treated individuals. The absence of treatment history, the need to adjust ART despite undetectable viremia, or the presence of LLV further support the use of genotypic resistance tests (GRTs) on HIV-1 DNA. Conventionally, GRTs have been achieved through Sanger sequencing, but the advances in NGS are leading to an increase in its use, allowing the detection of minority variants present in less than 20% of the viral population. The clinical significance of these mutations remains under debate, with interpretations varying based on context. Additionally, proviral DNA is subject to APOBEC3-induced hypermutation, which can lead to defective, nonviable viral genomes, a factor that must be considered when performing GRTs on HIV-1 DNA.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels

Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.

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

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

Evolution of Antibody Responses in HIV-1 CRF01_AE Acute Infection: Founder Envelope V1V2 Impacts the Timing and Magnitude of Autologous Neutralizing Antibodies

Understanding the dynamics of early immune responses to HIV-1 infection, including the evolution of initial neutralizing and antibody-dependent cellular cytotoxicity (ADCC)-mediating antibodies, will inform HIV vaccine design. In this study, we assess the development of autologous neutralizing antibodies (ANAbs) against founder envelopes (Envs) from 18 participants with HIV-1 CRF01_AE acute infection. The timing of ANAb development directly associated with the magnitude of the longitudinal ANAb response. Participants that developed ANAbs within 6 months of infection had significantly higher ANAb responses at 1 year (50% inhibitory concentration [IC50] geometric mean titer [GMT] = 2,010 versus 184; P = 0.001) and 2 years (GMT = 3,479 versus 340; P = 0.015), compared to participants that developed ANAb responses after 6 months. Participants with later development of ANAb tended to develop an earlier, potent heterologous tier 1 (92TH023) neutralizing antibody (NAb) response (P = 0.049). CRF01_AE founder Env V1V2 loop lengths correlated indirectly with the timing (P = 0.002, r = -0.675) and directly with magnitude (P = 0.005, r = 0.635) of ANAb responses; Envs with longer V1V2 loop lengths elicited earlier and more potent ANAb responses. While ANAb responses did not associate with viral load, the viral load set point correlated directly with neutralization of the heterologous 92TH023 strain (P = 0.007, r = 0.638). In contrast, a striking inverse correlation was observed between viral load set point and peak ADCC against heterologous 92TH023 Env strain (P = 0.0005, r = -0.738). These data indicate that specific antibody functions can be differentially related to viral load set point and may affect HIV-1 pathogenesis. Exploiting Env properties, such as V1V2 length, could facilitate development of subtype-specific vaccines that elicit more effective immune responses and improved protection. IMPORTANCE Development of an effective HIV-1 vaccine will be facilitated by better understanding the dynamics between the founder virus and the early humoral responses. Variations between subtypes may influence the evolution of immune responses and should be considered as we strive to understand these dynamics. In this study, autologous founder envelope neutralization and heterologous functional humoral responses were evaluated after acute infection by HIV-1 CRF01_AE, a subtype that has not been thoroughly characterized. The evolution of these humoral responses was assessed in relation to envelope characteristics, magnitude of elicited immune responses, and viral load. Understanding immune parameters in natural infection will improve our understanding of protective responses and aid in the development of immunogens that elicit protective functional antibodies. Advancing our knowledge of correlates of positive clinical outcomes should lead to the design of more efficacious vaccines.

Cohort-Specific Peptide Reagents Broaden Depth and Breadth Estimates of the CD8 T Cell Response to HIV-1 Gag Potential T Cell Epitopes

An effective HIV vaccine will need to stimulate immune responses against the sequence diversity presented in circulating virus strains. In this study, we evaluate breadth and depth estimates of potential T-cell epitopes (PTEs) in transmitted founder virus sequence-derived cohort-specific peptide reagents against reagents representative of consensus and global sequences. CD8 T-cells from twenty-six HIV-1+ PBMC donor samples, obtained at 1-year post estimated date of infection, were evaluated. ELISpot assays compared responses to 15mer consensus (n = 121), multivalent-global (n = 320), and 10mer multivalent cohort-specific (n = 300) PTE peptides, all mapping to the Gag antigen. Responses to 38 consensus, 71 global, and 62 cohort-specific PTEs were confirmed, with sixty percent of common global and cohort-specific PTEs corresponding to consensus sequences. Both global and cohort-specific peptides exhibited broader epitope coverage compared to commonly used consensus reagents, with mean breadth estimates of 3.2 (global), 3.4 (cohort) and 2.2 (consensus) epitopes. Global or cohort peptides each identified unique epitope responses that would not be detected if these peptide pools were used alone. A peptide set designed around specific virologic and immunogenetic characteristics of a target cohort can expand the detection of CD8 T-cell responses to epitopes in circulating viruses, providing a novel way to better define the host response to HIV-1 with implications for vaccine development.

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Increased Frequency of Indels in Hypervariable Regions of SARS-CoV-2 Proteins—A Possible Signature of Adaptive Selection

Most attention in the surveillance of evolving SARS-CoV-2 genome has been centered on nucleotide substitutions in the spike glycoprotein. We show that, as the pandemic extends into its second year, the numbers and ratio of genomes with in-frame insertions and deletions (indels) increases significantly, especially among the variants of concern (VOCs). Monitoring of the SARS-CoV-2 genome evolution shows that co-occurrence (i.e., highly correlated presence) of indels, especially deletions on spike N-terminal domain and non-structural protein 6 (NSP6) is a shared feature in several VOCs such as Alpha, Beta, Delta, and Omicron. Indels distribution is correlated with spike mutations associated with immune escape and growth in the number of genomes with indels coincides with the increasing population resistance due to vaccination and previous infections. Indels occur most frequently in the spike, but also in other proteins, especially those involved in interactions with the host immune system. We also showed that indels concentrate in regions of individual SARS-CoV-2 proteins known as hypervariable regions (HVRs) that are mostly located in specific loop regions. Structural analysis suggests that indels remodel viral proteins’ surfaces at common epitopes and interaction interfaces, affecting the virus’ interactions with host proteins. We hypothesize that the increased frequency of indels, the non-random distribution of them and their independent co-occurrence in several VOCs is another mechanism of response to elevated global population immunity.

Genomic signatures of protease and reverse transcriptase genes from HIV-1 subtype C isolated from first-line ART patients in India

Genomic signatures of the protease and reverse transcriptase gene of HIV-1 from HIV infected North Indian patients who were under ART from 1 to ≤ 7 years were analyzed. The DNA from plasma samples of 9 patients and RNA from 57 patients were isolated and subjected to amplification for the protease and reverse transcriptase gene of HIV-1 subtype C. Then sequencing was carried out following the WHO dried blood spot protocol. The drug resistance mutation patterns were analyzed using the HIV Drug Resistance Database, Stanford University, USA. Lamivudine-associated drug-resistance mutations such as M184V/M184I, nevirapine-associated drug resistance mutations Y181C and H221Y, and efavirenz-associated drug resistance mutations M230I were observed in reverse transcriptase gene of archived DNA of two HIV-1 infected patients. No mutation was observed in the remaining 7 patients. Various computational tools and websites like viral epidemiological signature pattern analysis (VESPA), hyper mutation, SNAP version 2.1.1, and entropy were utilized for the analysis of the signature pattern of amino acids, hyper mutation, selection pressure, and Shannon entropy in the protease and reverse transcriptase gene sequences of the 9 archived DNA, 56 protease gene and 51 reverse transcriptase gene from the HIV-1 DNA amplified sequences of RNA. The HIV-1 Subtype-C (Gene bank accession number: AB023804) and first isolate HXB2 (Gene bank accession number: K03455.1) was taken as reference sequence. The signature amino acid sequences were identified in the protease and reverse transcriptase gene, no hyper mutation, highest entropy was marked in the amino acid positions and synonymous to non-synonymous nucleotide ratio was calculated in the protease and reverse transcriptase gene of 9 archived DNA sequences, 56 protease and 51 reverse transcriptase gene sequences of HIV-1 Subtype C isolates.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

Plasma-Derived HIV-1 Virions Contain Considerable Levels of Defective Genomes

Genetically-characterizing full-length HIV-1 RNA is critical for identifying genetically-intact genomes and for comparing these RNA genomes to proviral DNA. We have developed a method for sequencing plasma-derived RNA using long-range sequencing (PRLS assay; ∼8.3 kb from gag to the 3′ end or ∼5 kb from integrase to the 3′ end). We employed the gag-3′ PRLS assay to sequence HIV-1 RNA genomes from ART-naive participants during acute/early infection (n = 6) or chronic infection (n = 2). On average, only 65% of plasma-derived genomes were genetically-intact. Defects were found in all genomic regions but were concentrated in env and pol. We compared these genomes to near-full-length proviral sequences from paired peripheral blood mononuclear cell (PBMC) samples for the acute/early group and found that near-identical (>99.98% identical) sequences were identified only during acute infection. For three participants who initiated therapy during acute infection, we used the int-3′ PRLS assay to sequence plasma-derived genomes from an analytical treatment interruption and identified 100% identical genomes between pretherapy and rebound time points. The PRLS assay provides a new level of sensitivity for understanding the genetic composition of plasma-derived HIV-1 RNA from viremic individuals either pretherapy or after treatment interruption, which will be invaluable in assessing possible HIV-1 curative strategies.

IMPORTANCE We developed novel plasma-derived RNA using long-range sequencing assays (PRLS assay; 8.3 kb, gag-3′, and 5.0 kb, int-3′). Employing the gag-3′ PRLS assay, we found that 26% to 51% of plasma-derived genomes are genetically-defective, largely as a result of frameshift mutations and deletions. These genetic defects were concentrated in the env region compared to gag and pol, likely a reflection of viral immune escape in env during untreated HIV-1 infection. Employing the int-3′ PRLS assay, we found that analytical treatment interruption (ATI) plasma-derived sequences were identical and genetically-intact. Several sequences from the ATI plasma samples were identical to viral sequences from pretherapy plasma and PBMC samples, indicating that HIV-1 reservoirs established prior to therapy contribute to viral rebound during an ATI. Therefore, near-full-length sequencing of HIV-1 particles is required to gain an accurate picture of the genetic landscape of plasma HIV-1 virions in studies of HIV-1 replication and persistence.

Phylogenetic inference of changes in amino acid propensities with single-position resolution

Fitness conferred by the same allele may differ between genotypes and environments, and these differences shape variation and evolution. Changes in amino acid propensities at protein sites over the course of evolution have been inferred from sequence alignments statistically, but the existing methods are data-intensive and aggregate multiple sites. Here, we develop an approach to detect individual amino acids that confer different fitness in different groups of species from combined sequence and phylogenetic data. Using the fact that the probability of a substitution to an amino acid depends on its fitness, our method looks for amino acids such that substitutions to them occur more frequently in one group of lineages than in another. We validate our method using simulated evolution of a protein site under different scenarios and show that it has high specificity for a wide range of assumptions regarding the underlying changes in selection, while its sensitivity differs between scenarios. We apply our method to the env gene of two HIV-1 subtypes, A and B, and to the HA gene of two influenza A subtypes, H1 and H3, and show that the inferred fitness changes are consistent with the fitness differences observed in deep mutational scanning experiments. We find that changes in relative fitness of different amino acid variants within a site do not always trigger episodes of positive selection and therefore may not result in an overall increase in the frequency of substitutions, but can still be detected from changes in relative frequencies of different substitutions.

The mutational landscape of SARS-CoV-2 variants diversifies T cell targets in an HLA-supertype-dependent manner

The rapid, global dispersion of SARS-CoV-2 has led to the emergence of a diverse range of variants. Here, we describe how the mutational landscape of SARS-CoV-2 has shaped HLA-restricted T cell immunity at the population level during the first year of the pandemic. We analyzed a total of 330,246 high-quality SARS-CoV-2 genome assemblies, sampled across 143 countries and all major continents from December 2019 to December 2020 before mass vaccination or the rise of the Delta variant. We observed that proline residues are preferentially removed from the proteome of prevalent mutants, leading to a predicted global loss of SARS-CoV-2 T cell epitopes in individuals expressing HLA-B alleles of the B7 supertype family; this is largely driven by a dominant C-to-U mutation type at the RNA level. These results indicate that B7-supertype-associated epitopes, including the most immunodominant ones, were more likely to escape CD8⁺ T cell immunosurveillance during the first year of the pandemic.

RV144 vaccine imprinting constrained HIV-1 evolution following breakthrough infection

The scale of the HIV-1 epidemic underscores the need for a vaccine. The multitude of circulating HIV-1 strains together with HIV-1’s high evolvability hints that HIV-1 could adapt to a future vaccine. Here, we wanted to investigate the effect of vaccination on the evolution of the virus post-breakthrough infection. We analyzed 2,635 HIV-1 env sequences sampled up to a year post-diagnosis from 110 vaccine and placebo participants who became infected in the RV144 vaccine efficacy trial. We showed that the Env signature sites that were previously identified to distinguish vaccine and placebo participants were maintained over time. In addition, fewer sites were under diversifying selection in the vaccine group than in the placebo group. These results indicate that HIV-1 would possibly adapt to a vaccine upon its roll-out.

Structure-Based Design of First-Generation Small Molecule Inhibitors Targeting the Catalytic Pockets of AID, APOBEC3A, and APOBEC3B

Activation-induced cytidine deaminase (AID) initiates antibody diversification by mutating immunoglobulin loci in B lymphocytes. AID and related APOBEC3 (A3) enzymes also induce genome-wide mutations and lesions implicated in tumorigenesis and tumor progression. The most prevalent mutation signatures across diverse tumor genomes are attributable to the mistargeted mutagenic activities of AID/A3s. Thus, inhibiting AID/A3s has been suggested to be of therapeutic benefit. We previously used a computational-biochemical approach to gain insight into the structure of AID’s catalytic pocket, which resulted in the discovery of a novel type of regulatory catalytic pocket closure that regulates AID/A3s that we termed the “Schrodinger’s CATalytic pocket”. Our findings were subsequently confirmed by direct structural studies. Here, we describe our search for small molecules that target the catalytic pocket of AID. We identified small molecules that inhibit purified AID, AID in cell extracts, and endogenous AID of lymphoma cells. Analogue expansion yielded derivatives with improved potencies. These were found to also inhibit A3A and A3B, the two most tumorigenic siblings of AID. Two compounds exhibit low micromolar IC₅₀ inhibition of AID and A3A, exhibiting the strongest potency for A3A. Docking suggests key interactions between their warheads and residues lining the catalytic pockets of AID, A3A, and A3B and between the tails and DNA-interacting residues on the surface proximal to the catalytic pocket opening. Accordingly, mutants of these residues decreased inhibition potency. The chemistry and abundance of key stabilizing interactions between the small molecules and residues within and immediately outside the catalytic pockets are promising for therapeutic development.

Human APOBEC3 Variations and Viral Infection

Human APOBEC3 (apolipoprotein B mRNA-editing catalytic polypeptide-like 3) enzymes are capable of inhibiting a wide range of endogenous and exogenous viruses using deaminase and deaminase-independent mechanisms. These enzymes are essential components of our innate immune system, as evidenced by (a) their strong positive selection and expansion in primates, (b) the evolution of viral counter-defense mechanisms, such as proteasomal degradation mediated by HIV Vif, and (c) hypermutation and inactivation of a large number of integrated HIV-1 proviruses. Numerous APOBEC3 single nucleotide polymorphisms, haplotypes, and splice variants have been identified in humans. Several of these variants have been reported to be associated with differential antiviral immunity. This review focuses on the current knowledge in the field about these natural variations and their roles in infectious diseases.

HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens

Humanity is currently facing the challenge of two devastating pandemics caused by two very different RNA viruses: HIV-1, which has been with us for decades, and SARS-CoV-2, which has swept the world in the course of a single year. The same evolutionary strategies that drive HIV-1 evolution are at play in SARS-CoV-2. Single nucleotide mutations, multi-base insertions and deletions, recombination, and variation in surface glycans all generate the variability that, guided by natural selection, enables both HIV-1's extraordinary diversity and SARS-CoV-2's slower pace of mutation accumulation. Even though SARS-CoV-2 diversity is more limited, recently emergent SARS-CoV-2 variants carry Spike mutations that have important phenotypic consequences in terms of both antibody resistance and enhanced infectivity. We review and compare how these mutational patterns manifest in these two distinct viruses to provide the variability that fuels their evolution by natural selection.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Recapitulation of HIV-1 Env-antibody coevolution in macaques leading to neutralization breadth

Neutralizing antibodies elicited by HIV-1 coevolve with viral envelope proteins (Env) in distinctive patterns, in some cases acquiring substantial breadth. We report that primary HIV-1 envelope proteins-when expressed by simian-human immunodeficiency viruses in rhesus macaques-elicited patterns of Env-antibody coevolution very similar to those in humans, including conserved immunogenetic, structural, and chemical solutions to epitope recognition and precise Env-amino acid substitutions, insertions, and deletions leading to virus persistence. The structure of one rhesus antibody, capable of neutralizing 49% of a 208-strain panel, revealed a V2 apex mode of recognition like that of human broadly neutralizing antibodies (bNAbs) PGT145 and PCT64-35S. Another rhesus antibody bound the CD4 binding site by CD4 mimicry, mirroring human bNAbs 8ANC131, CH235, and VRC01. Virus-antibody coevolution in macaques can thus recapitulate developmental features of human bNAbs, thereby guiding HIV-1 immunogen design.

Simulating within host human immunodeficiency virus 1 genome evolution in the persistent reservoir

The complexities of viral evolution can be difficult to elucidate. Software simulating viral evolution provides powerful tools for exploring hypotheses of viral systems, especially in situations where thorough empirical data are difficult to obtain or parameters of interest are difficult to measure. Human immunodeficiency virus 1 (HIV-1) infection has no durable cure; this is primarily due to the virus’ ability to integrate into the genome of host cells, where it can remain in a transcriptionally latent state. An effective cure strategy must eliminate every copy of HIV-1 in this ‘persistent reservoir’ because proviruses can reactivate, even decades later, to resume an active infection. However, many features of the persistent reservoir remain unclear, including the temporal dynamics of HIV-1 integration frequency and the longevity of the resulting reservoir. Thus, sophisticated analyses are required to measure these features and determine their temporal dynamics. Here, we present software that is an extension of SANTA-SIM to include multiple compartments of viral populations. We used the resulting software to create a model of HIV-1 within host evolution that incorporates the persistent HIV-1 reservoir. This model is composed of two compartments, an active compartment and a latent compartment. With this model, we compared five different date estimation methods (Closest Sequence, Clade, Linear Regression, Least Squares, and Maximum Likelihood) to recover the integration dates of genomes in our model’s HIV-1 reservoir. We found that the Least Squares method performed the best with the highest concordance (0.80) between real and estimated dates and the lowest absolute error (all pairwise t tests: P < 0.01). Our software is a useful tool for validating bioinformatics software and understanding the dynamics of the persistent HIV-1 reservoir.

HIV infects astrocytes in vivo and egresses from the brain to the periphery

HIV invades the brain during acute infection. Yet, it is unknown whether long-lived infected brain cells release productive virus that can egress from the brain to re-seed peripheral organs. This understanding has significant implication for the brain as a reservoir for HIV and most importantly HIV interplay between the brain and peripheral organs. Given the sheer number of astrocytes in the human brain and their controversial role in HIV infection, we evaluated their infection in vivo and whether HIV infected astrocytes can support HIV egress to peripheral organs. We developed two novel models of chimeric human astrocyte/human peripheral blood mononuclear cells: NOD/scid-IL-2Rgc null (NSG) mice (huAstro/HuPBMCs) whereby we transplanted HIV (non-pseudotyped or VSVg-pseudotyped) infected or uninfected primary human fetal astrocytes (NHAs) or an astrocytoma cell line (U138MG) into the brain of neonate or adult NSG mice and reconstituted the animals with human peripheral blood mononuclear cells (PBMCs). We also transplanted uninfected astrocytes into the brain of NSG mice and reconstituted with infected PBMCs to mimic a biological infection course. As expected, the xenotransplanted astrocytes did not escape/migrate out of the brain and the blood brain barrier (BBB) was intact in this model. We demonstrate that astrocytes support HIV infection in vivo and egress to peripheral organs, at least in part, through trafficking of infected CD4+ T cells out of the brain. Astrocyte-derived HIV egress persists, albeit at low levels, under combination antiretroviral therapy (cART). Egressed HIV evolved with a pattern and rate typical of acute peripheral infection. Lastly, analysis of human cortical or hippocampal brain regions of donors under cART revealed that astrocytes harbor between 0.4-5.2% integrated HIV gag DNA and 2-7% are HIV gag mRNA positive. These studies establish a paradigm shift in the dynamic interaction between the brain and peripheral organs which can inform eradication of HIV reservoirs.

Structural Insights into APOBEC3-Mediated Lentiviral Restriction

Mammals have developed clever adaptive and innate immune defense mechanisms to protect against invading bacterial and viral pathogens. Human innate immunity is continuously evolving to expand the repertoire of restriction factors and one such family of intrinsic restriction factors is the APOBEC3 (A3) family of cytidine deaminases. The coordinated expression of seven members of the A3 family of cytidine deaminases provides intrinsic immunity against numerous foreign infectious agents and protects the host from exogenous retroviruses and endogenous retroelements. Four members of the A3 proteins-A3G, A3F, A3H, and A3D-restrict HIV-1 in the absence of virion infectivity factor (Vif); their incorporation into progeny virions is a prerequisite for cytidine deaminase-dependent and -independent activities that inhibit viral replication in the host target cell. HIV-1 encodes Vif, an accessory protein that antagonizes A3 proteins by targeting them for polyubiquitination and subsequent proteasomal degradation in the virus producing cells. In this review, we summarize our current understanding of the role of human A3 proteins as barriers against HIV-1 infection, how Vif overcomes their antiviral activity, and highlight recent structural and functional insights into A3-mediated restriction of lentiviruses.

Impact of Suboptimal APOBEC3G Neutralization on the Emergence of HIV Drug Resistance in Humanized Mice

HIV diversification facilitates immune escape and complicates antiretroviral therapy. In this study, we take advantage of a humanized-mouse model to probe the contribution of APOBEC3 mutagenesis to viral evolution. Humanized mice were infected with isogenic HIV molecular clones (HIV-WT, HIV-45G, and HIV-ΔSLQ) that differ in their abilities to counteract APOBEC3G (A3G). Infected mice remained naive or were treated with the reverse transcriptase (RT) inhibitor lamivudine (3TC). Viremia, emergence of drug-resistant variants, and quasispecies diversification in the plasma compartment were determined throughout infection. While both HIV-WT and HIV-45G achieved robust infection, over time, HIV-45G replication was significantly reduced compared to that of HIV-WT in the absence of 3TC treatment. In contrast, treatment responses differed significantly between HIV-45G- and HIV-WT-infected mice. Antiretroviral treatment failed in 91% of HIV-45G-infected mice, while only 36% of HIV-WT-infected mice displayed a similar negative outcome. Emergence of 3TC-resistant variants and nucleotide diversity were determined by analyzing 155,462 single HIV reverse transcriptase gene (RT) and 6,985 vif sequences from 33 mice. Prior to treatment, variants with genotypic 3TC resistance (RT-M184I/V) were detected at low levels in over a third of all the animals. Upon treatment, the composition of the plasma quasispecies rapidly changed, leading to a majority of circulating viral variants encoding RT-184I. Interestingly, increased viral diversity prior to treatment initiation correlated with higher plasma viremia in HIV-45G-infected animals, but not in HIV-WT-infected animals. Taken together, HIV variants with suboptimal anti-A3G activity were attenuated in the absence of selection but displayed a fitness advantage in the presence of antiretroviral treatment.IMPORTANCE Both viral (e.g., RT) and host (e.g., A3G) factors can contribute to HIV sequence diversity. This study shows that suboptimal anti-A3G activity shapes viral fitness and drives viral evolution in the plasma compartment in humanized mice.

Phylogenetic measures of indel rate variation among the HIV-1 group M subtypes

The transmission fitness and pathogenesis of HIV-1 is disproportionately influenced by evolution in the five variable regions (V1–V5) of the surface envelope glycoprotein (gp120). Insertions and deletions (indels) are a significant source of evolutionary change in these regions. However, the rate and composition of indels has not yet been quantified through a large-scale comparative analysis of HIV-1 sequences. Here, we develop and report results from a phylogenetic method to estimate indel rates for the gp120 variable regions across five major subtypes and two circulating recombinant forms (CRFs) of HIV-1 group M. We processed over 26,000 published HIV-1 gp120 sequences, from which we extracted 6,605 sequences for phylogenetic analysis. We reconstructed time-scaled phylogenies by maximum likelihood and fit a binomial-Poisson model to the observed distribution of indels between closely related pairs of sequences in each tree (cherries). By focusing on cherries in each tree, we obtained phylogenetically independent indel reconstructions, and the shorter time scales in cherries reduced the bias due to purifying selection. Rate estimates ranged from 3.0×10−5 to 1.5×10−3 indels/nt/year and varied significantly among variable regions and subtypes. Indel rates were significantly lower in V3 relative to V1, and were also lower in HIV-1 subtype B relative to the 01_AE reference. We also found that V1, V2, and V4 tended to accumulate significantly longer indels. Furthermore, we observed that the nucleotide composition of indels was distinct from the flanking sequence, with higher frequencies of G and lower frequencies of T. Indels affected N-linked glycosylation sites more often in V1 and V2 than expected by chance, consistent with positive selection on glycosylation patterns within these regions. These results represent the first comprehensive measures of indel rates in HIV-1 gp120 across multiple subtypes and CRFs, and identifies novel and unexpected patterns for further research in the molecular evolution of HIV-1.

Parallel evolution of HIV-1 in a long-term experiment

One of the most intriguing puzzles in biology is the degree to which evolution is repeatable. The repeatability of evolution, or parallel evolution, has been studied in a variety of model systems, but has rarely been investigated with clinically relevant viruses. To investigate parallel evolution of HIV-1, we passaged two replicate HIV-1 populations for almost 1 year in each of two human T-cell lines. For each of the four evolution lines, we determined the genetic composition of the viral population at nine time points by deep sequencing the entire genome. Mutations that were carried by the majority of the viral population accumulated continuously over 1 year in each evolution line. Many majority mutations appeared in more than one evolution line, that is, our experiments showed an extreme degree of parallel evolution. In one of the evolution lines, 62% of the majority mutations also occur in another line. The parallelism impairs our ability to reconstruct the evolutionary history by phylogenetic methods. We show that one can infer the correct phylogenetic topology by including minority mutations in our analysis. We also find that mutation diversity at the beginning of the experiment is predictive of the frequency of majority mutations at the end of the experiment.

Virus dynamics and phyloanatomy: Merging population dynamic and phylogenetic approaches

In evolutionary biology and epidemiology, phylodynamic methods are widely used to infer population biological characteristics, such as the rates of replication, death, migration, or, in the epidemiological context, pathogen spread. More recently, these methods have been used to elucidate the dynamics of viruses within their hosts. Especially the application of phylogeographic approaches has the potential to shed light on anatomical colonization pathways and the exchange of viruses between distinct anatomical compartments. We and others have termed this phyloanatomy. Here, we review the promise and challenges of phyloanatomy, and compare them to more classical virus dynamics and population genetic approaches. We argue that the extremely strong selection pressures that exist within the host may represent the main obstacle to reliable phyloanatomic analysis.

An insight into the dependence of the deamination rate of human APOBEC3F on the length of single-stranded DNA, which is affected by the concentrations of APOBEC3F and single-stranded DNA

BACKGROUND

APOBEC3F (A3F), a member of the human APOBEC3 (A3) family of cytidine deaminases, acts as an anti-HIV-1 factor by deaminating deoxycytidine in the complementary DNA of the viral genome. A full understanding of the deamination behavior of A3F awaits further investigation.

METHODS

The real-time NMR method and uracil-DNA glycosylase assay were used to track the activities of the C-terminal domain (CTD) of A3F at different concentrations of A3F-CTD and ssDNA. The steady-state fluorescence anisotropy measurement was used to examine the binding between A3F-CTD and ssDNA with different lengths. The use of the A3F-CTD N214H mutant, having higher activity than the wild-type, facilitated the tracking of the reactions.

RESULTS

A3F-CTD was found to efficiently deaminate the target deoxycytidine in long ssDNA in lower ssDNA concentration conditions ([A3F-CTD] ≫ [ssDNA]), while the target deoxycytidine in short ssDNA is deaminated efficiently in higher ssDNA concentration conditions ([A3F-CTD] ≪ [ssDNA]). This property is quite different from that of the previously studied A3 family member, A3B; the concentrations of the proteins and ssDNA had no effect.

CONCLUSIONS

The concentrations of A3F-CTD and ssDNA substrates affect the ssDNA-length-dependence of deamination rate of the A3F-CTD. This unique property of A3F is rationally interpreted on the basis of its binding characteristics with ssDNA.

GENERAL SIGNIFICANCE

The discovery of the unique property of A3F regarding the deamination rate deepens the understanding of its counteraction against HIV-1. Our strategy is applicable to investigate the other aspects of the A3 activities, such as those involved in the cancer development.

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 effects of the AID/APOBEC family of mutagenic enzymes on human gamma-herpesviruses

The human gamma-herpesviruses, Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, establish lifelong latency in B cells and are associated with multiple malignancies. Virus-host coevolution often drive changes in both host immunity and in the viral genome. We consider one host immune mechanism, the activation-induced deaminase (AID)/APOBEC family of cytidine deaminases, that induces mutations in viral DNA. AID, the ancestral gene in the family has a conserved role in somatic hypermutation, a key step in antibody affinity maturation. The APOBEC3 subfamily, of which there are seven genes in human, have evolved antiviral functions and have diversified in terms of their expression pattern, subcellular localization, and DNA mutation motifs (hotspots). In this study, we investigated how the human gamma-herpesviruses have evolved to avoid the action of the AID/APOBEC enzymes and determine if these enzymes are contributing to the ongoing evolution of the viruses. We used computational methods to evaluate observed versus expected frequency of AID/APOBEC hotspots in viral genomes and found that the viruses have evolved to limit the representation of AID and certain APOBEC3 motifs. At the same time, the remaining hotspots were highly likely to cause amino acid changes, suggesting prolonged evolutionary pressure of the enzymes on the viruses. To study current hypermutation, as opposed to historical mutation processes, we also analyzed putative mutations derived from alignments of published viral genomes and found again that AID and APOBEC3 appear to target the genome most frequently. New protein variants resulting from AID/APOBEC activity may have important consequences in health, including vaccine development (epitope evolution) and host immune evasion.

Inhibiting APOBEC3 Activity with Single-Stranded DNA Containing 2'-Deoxyzebularine Analogues

APOBEC3 enzymes form part of the innate immune system by deaminating cytosine to uracil in single-stranded DNA (ssDNA) and thereby preventing the spread of pathogenic genetic information. However, APOBEC mutagenesis is also exploited by viruses and cancer cells to increase rates of evolution, escape adaptive immune responses, and resist drugs. This raises the possibility of APOBEC3 inhibition as a strategy for augmenting existing antiviral and anticancer therapies. Here we show that, upon incorporation into short ssDNAs, the cytidine nucleoside analogue 2'-deoxyzebularine (dZ) becomes capable of inhibiting the catalytic activity of selected APOBEC variants derived from APOBEC3A, APOBEC3B, and APOBEC3G, supporting a mechanism in which ssDNA delivers dZ to the active site. Multiple experimental approaches, including isothermal titration calorimetry, fluorescence polarization, protein thermal shift, and nuclear magnetic resonance spectroscopy assays, demonstrate nanomolar dissociation constants and low micromolar inhibition constants. These dZ-containing ssDNAs constitute the first substrate-like APOBEC3 inhibitors and, together, comprise a platform for developing nucleic acid-based inhibitors with cellular activity.

HIV-1 immunogens and strategies to drive antibody responses towards neutralization breadth

Despite enormous efforts no HIV-1 vaccine has been developed that elicits broadly neutralizing antibodies (bNAbs) to protect against infection to date. The high antigenic diversity and dense N-linked glycan armor, which covers nearly the entire HIV-1 envelope protein (Env), are major roadblocks for the development of bNAbs by vaccination. In addition, the naive human antibody repertoire features a low frequency of exceptionally long heavy chain complementary determining regions (CDRH3s), which is a typical characteristic that many HIV-1 bNAbs use to penetrate the glycan armor. Native-like Env trimer immunogens can induce potent but strain-specific neutralizing antibody responses in animal models but how to overcome the many obstacles towards the development of bNAbs remains a challenge. Here, we review recent HIV-1 Env immunization studies and discuss strategies to guide strain-specific antibody responses towards neutralization breadth.

Estimating the mutational fitness effects distribution during early HIV infection

The evolution of HIV during acute infection is often considered a neutral process. Recent analysis of sequencing data from this stage of infection, however, showed high levels of shared mutations between independent viral populations. This suggests that selection might play a role in the early stages of HIV infection. We adapted an existing model for random evolution during acute HIV-infection to include selection. Simulations of this model were used to fit a global mutational fitness effects distribution to previously published sequencing data of the env gene of individuals with acute HIV infection. Measures of sharing between viral populations were used as summary statistics to compare the data to the simulations. We confirm that evolution during acute infection is significantly different from neutral. The distribution of mutational fitness effects is best fit by a distribution with a low, but significant fraction of beneficial mutations and a high fraction of deleterious mutations. While most mutations are neutral or deleterious in this model, about 5% of mutations are beneficial. These beneficial mutations will, on average, result in a small but significant increase in fitness. When assuming no epistasis, this indicates that, at the moment of transmission, HIV is near, but not on the fitness peak for early infection.

Characterization of the Plasmacytoid Dendritic Cell Response to Transmitted/Founder and Nontransmitted Variants of HIV-1

Human immunodeficiency virus type 1 (HIV-1) infection often arises from a single transmitted/founder (TF) viral variant among a large pool of viruses in the quasispecies in the transmitting partner. TF variants are typically nondominant in blood and genital secretions, indicating that they have unique traits. The plasmacytoid dendritic cell (pDC) is the primary alpha interferon (IFN-α)-producing cell in response to viral infections and is rapidly recruited to the female genital tract upon exposure to HIV-1. The impact of pDCs on transmission is unknown. We investigated whether evasion of pDC responses is a trait of TF viruses. pDCs from healthy donors were stimulated in vitro with a panel of 20 HIV-1 variants, consisting of one TF variant and three nontransmitted (NT) variants each from five transmission-linked donor pairs, and secretion of IFN-α and tumor necrosis factor alpha (TNF-α) was measured by enzyme-linked immunosorbent assay (ELISA). No significant differences in cytokine secretion in response to TF and NT viruses were observed, despite a trend toward enhanced IFN-α and TNF-α production in response to TF viruses. NT viruses demonstrated polarization toward production of either IFN-α or TNF-α, indicating possible dysregulation. Also, for NT viruses, IFN-α secretion was associated with increased resistance of the virus to inactivation by IFN-α in vitro, suggesting in vivo evolution. Thus, TF viruses do not appear to preferentially subvert pDC activation compared to that with nontransmitted HIV-1 variants. pDCs may, however, contribute to the in vivo evolution of HIV-1.IMPORTANCE The plasmacytoid dendritic cell (pDC) is the first cell type recruited to the site of HIV-1 exposure; however, its contribution to the viral bottleneck in HIV-1 transmission has not been explored previously. We hypothesized that transmitted/founder viruses are able to avoid the pDC response. In this study, we used previously established donor pair-linked transmitted/founder and nontransmitted (or chronic) variants of HIV-1 to stimulate pDCs. Transmitted/founder HIV-1, instead of suppressing pDC responses, induced IFN-α and TNF-α secretion to levels comparable to those induced by viruses from the transmitting partner. We noted several unique traits of chronic viruses, including polarization between IFN-α and TNF-α production as well as a strong relationship between IFN-α secretion and the resistance of the virus to neutralization. These data rule out the possibility that TF viruses preferentially suppress pDCs in comparison to the pDC response to nontransmitted HIV variants. pDCs may, however, be important drivers of viral evolution in vivo.

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.

Quasispecies characters of hepatitis B virus in immunoprophylaxis failure infants

Hepatitis B vaccination prevents 80-95% of transmission and reduces the incidence of HBV in children. The variations in the a determinant of HBV surface antigen (HBsAg) have been reported to be the most prevalent cause for vaccine or antibody escape. There is a conflicting evidence on as to whether escape mutants arise de novo in infected infants or whether the mutants, that have preexisted maternally, subsequently undergo selective replication in the infant under immune pressure. Here, we report that nearly 65% (55 of 85) vaccination failure in child patients has no amino acid substitution in a determinant as seen by Sanger sequencing. We further employed an Illumina sequencing platform-based method to detect HBV quasispecies in four immunoprophylaxis failure infants and their mothers. In our data, the substitution rate of amino acid located at a determinant is relatively low (< 10%), I/T126A, C124S, F134Y, K141Q, Q129H, D144A, G145V, and N146K, which showed no statistical difference to their mothers, proving that these vaccine escape mutants preexist maternally as minor variants. Besides that, bioinformatical analysis showed that the binding affinity of high variation epitopes (amino acid divergence in mother and their infants > 20%) to related HLA molecules was generally decreased, these traces of immune escape suggesting that immune pressure was present and was effective in all samples.

Biochemical Basis of APOBEC3 Deoxycytidine Deaminase Activity on Diverse DNA Substrates

The Apolipoprotein B mRNA editing complex (APOBEC) family of enzymes contains single-stranded polynucleotide cytidine deaminases. These enzymes catalyze the deamination of cytidine in RNA or single-stranded DNA, which forms uracil. From this 11 member enzyme family in humans, the deamination of single-stranded DNA by the seven APOBEC3 family members is considered here. The APOBEC3 family has many roles, such as restricting endogenous and exogenous retrovirus replication and retrotransposon insertion events and reducing DNA-induced inflammation. Similar to other APOBEC family members, the APOBEC3 enzymes are a double-edged sword that can catalyze deamination of cytosine in genomic DNA, which results in potential genomic instability due to the many mutagenic fates of uracil in DNA. Here, we discuss how these enzymes find their single-stranded DNA substrate in different biological contexts such as during human immunodeficiency virus (HIV) proviral DNA synthesis, retrotransposition of the LINE-1 element, and the "off-target" genomic DNA substrate. The enzymes must be able to efficiently deaminate transiently available single-stranded DNA during reverse transcription, replication, or transcription. Specific biochemical characteristics promote deamination in each situation to increase enzyme efficiency through processivity, rapid enzyme cycling between substrates, or oligomerization state. The use of biochemical data to clarify biological functions and alignment with cellular data is discussed. Models to bridge knowledge from biochemical, structural, and single molecule experiments are presented.

Perspective: APOBEC mutagenesis in drug resistance and immune escape in HIV and cancer evolution

The apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC) mutational signature has only recently been detected in a multitude of cancers through next-generation sequencing. In contrast, APOBEC has been a focus of virology research for over a decade. Many lessons learnt regarding APOBEC within virology are likely to be applicable to cancer. In this review, we explore the parallels between the role of APOBEC enzymes in HIV and cancer evolution. We discuss data supporting the role of APOBEC mutagenesis in creating HIV genome heterogeneity, drug resistance, and immune escape variants. We hypothesize similar functions of APOBEC will also hold true in cancer.

Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver

Studying the evolution of viruses and their molecular epidemiology relies on accurate viral sequence data, so that small differences between similar viruses can be meaningfully interpreted. Despite its higher throughput and more detailed minority variant data, next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of large between- and within-host diversity, including frequent indels, may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to pre-process reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We used shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read whole-genome data produced with the Illumina platform, for sixty-five existing publicly available samples and fifty new samples. We show the systematic superiority of mapping to shiver's constructed reference compared with mapping the same reads to the closest of 3,249 real references: median values of 13 bases called differently and more accurately, 0 bases called differently and less accurately, and 205 bases of missing sequence recovered. We also successfully applied shiver to whole-genome samples of Hepatitis C Virus and Respiratory Syncytial Virus. shiver is publicly available from https://github.com/ChrisHIV/shiver.

Evasion of adaptive immunity by HIV through the action of host APOBEC3G/F enzymes

APOBEC3G (A3G) and APOBEC3F (A3F) are DNA-mutating enzymes expressed in T cells, dendritic cells and macrophages. A3G/F have been considered innate immune host factors, based on reports that they lethally mutate the HIV genome in vitro. In vivo, A3G/F effectiveness is limited by viral proteins, entrapment in inactive complexes and filtration of mutations during viral life cycle. We hypothesized that the impact of sub-lethal A3G/F action could extend beyond the realm of innate immunity confined to the cytoplasm of infected cells. We measured recognition of wild type and A3G/F-mutated epitopes by cytotoxic T lymphocytes (CTL) from HIV-infected individuals and found that A3G/F-induced mutations overwhelmingly diminished CTL recognition of HIV peptides, in a human histocompatibility-linked leukocyte antigen (HLA)-dependent manner. Furthermore, we found corresponding enrichment of A3G/F-favored motifs in CTL epitope-encoding sequences within the HIV genome. These findings illustrate that A3G/F‐mediated mutations mediate immune evasion by HIV in vivo. Therefore, we suggest that vaccine strategies target T cell or antibody epitopes that are not poised for mutation into escape variants by A3G/F action.

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.

The Latent Reservoir for HIV-1: How Immunologic Memory and Clonal Expansion Contribute to HIV-1 Persistence

Combination antiretroviral therapy (ART) for HIV-1 infection reduces plasma virus levels to below the limit of detection of clinical assays. However, even with prolonged suppression of viral replication with ART, viremia rebounds rapidly after treatment interruption. Thus, ART is not curative. The principal barrier to cure is a remarkably stable reservoir of latent HIV-1 in resting memory CD4(+) T cells. In this review, we consider explanations for the remarkable stability of the latent reservoir. Stability does not appear to reflect replenishment from new infection events but rather normal physiologic processes that provide for immunologic memory. Of particular importance are proliferative processes that drive clonal expansion of infected cells. Recent evidence suggests that in some infected cells, proliferation is a consequence of proviral integration into host genes associated with cell growth. Efforts to cure HIV-1 infection by targeting the latent reservoir may need to consider the potential of latently infected cells to proliferate.

Rare HIV-1 transmitted/founder lineages identified by deep viral sequencing contribute to rapid shifts in dominant quasispecies during acute and early infection

In order to inform the rational design of HIV-1 preventive and cure interventions it is critical to understand the events occurring during acute HIV-1 infection (AHI). Using viral deep sequencing on six participants from the early capture acute infection RV217 cohort, we have studied HIV-1 evolution in plasma collected twice weekly during the first weeks following the advent of viremia. The analysis of infections established by multiple transmitted/founder (T/F) viruses revealed novel viral profiles that included: a) the low-level persistence of minor T/F variants, b) the rapid replacement of the major T/F by a minor T/F, and c) an initial expansion of the minor T/F followed by a quick collapse of the same minor T/F to low frequency. In most participants, cytotoxic T-lymphocyte (CTL) escape was first detected at the end of peak viremia downslope, proceeded at higher rates than previously measured in HIV-1 infection, and usually occurred through the exploration of multiple mutational pathways within an epitope. The rapid emergence of CTL escape variants suggests a strong and early CTL response. Minor T/F viral strains can contribute to rapid and varied profiles of HIV-1 quasispecies evolution during AHI. Overall, our results demonstrate that early, deep, and frequent sampling is needed to investigate viral/host interaction during AHI, which could help identify prerequisites for prevention and cure of HIV-1 infection.

Calculating site-specific evolutionary rates at the amino-acid or codon level yields similar rate estimates

Site-specific evolutionary rates can be estimated from codon sequences or from amino-acid sequences. For codon sequences, the most popular methods use some variation of the dN∕dS ratio. For amino-acid sequences, one widely-used method is called Rate4Site, and it assigns a relative conservation score to each site in an alignment. How site-wise dN∕dS values relate to Rate4Site scores is not known. Here we elucidate the relationship between these two rate measurements. We simulate sequences with known dN∕dS, using either dN∕dS models or mutation–selection models for simulation. We then infer Rate4Site scores on the simulated alignments, and we compare those scores to either true or inferred dN∕dS values on the same alignments. We find that Rate4Site scores generally correlate well with true dN∕dS, and the correlation strengths increase in alignments with greater sequence divergence and more taxa. Moreover, Rate4Site scores correlate very well with inferred (as opposed to true) dN∕dS values, even for small alignments with little divergence. Finally, we verify this relationship between Rate4Site and dN∕dS in a variety of empirical datasets. We conclude that codon-level and amino-acid-level analysis frameworks are directly comparable and yield very similar inferences.

Single-Strand Consensus Sequencing Reveals that HIV Type but not Subtype Significantly Impacts Viral Mutation Frequencies and Spectra

A long-standing question of human immunodeficiency virus (HIV) genetic variation and evolution has been whether differences exist in mutation rate and/or mutation spectra among HIV types (i.e., HIV-1 versus HIV-2) and among HIV groups (i.e., HIV-1 groups M-P and HIV-2 groups A-H) and HIV-1 Group M subtypes (i.e., subtypes A-D, F-H, and J-K). To address this, we developed a new single-strand consensus sequencing assay for the determination of HIV mutation frequencies and spectra using the Illumina sequencing platform. This assay enables parallel and standardized comparison of HIV mutagenesis among various viral vectors with lower background error than traditional methods of Illumina library preparation. We found significant differences in viral mutagenesis between HIV types but intriguingly no significant differences among HIV-1 Group M subtypes. More specifically, HIV-1 exhibited higher transition frequencies than HIV-2, due mostly to single G-to-A mutations and (to a lesser extent) G-to-A hypermutation. These data suggest that HIV-2 RT exhibits higher fidelity during viral replication, and taken together, these findings demonstrate that HIV type but not subtype significantly affects viral mutation frequencies and spectra. These differences may inform antiviral and vaccine strategies.

Viral reverse transcriptases

Reverse transcriptases (RTs) play a major role in the replication of Retroviridae, Metaviridae, Pseudoviridae, Hepadnaviridae and Caulimoviridae. RTs are enzymes that are able to synthesize DNA using RNA or DNA as templates (DNA polymerase activity), and degrade RNA when forming RNA/DNA hybrids (ribonuclease H activity). In retroviruses and LTR retrotransposons (Metaviridae and Pseudoviridae), the coordinated action of both enzymatic activities converts single-stranded RNA into a double-stranded DNA that is flanked by identical sequences known as long terminal repeats (LTRs). RTs of retroviruses and LTR retrotransposons are active as monomers (e.g. murine leukemia virus RT), homodimers (e.g. Ty3 RT) or heterodimers (e.g. human immunodeficiency virus type 1 (HIV-1) RT). RTs lack proofreading activity and display high intrinsic error rates. Besides, high recombination rates observed in retroviruses are promoted by poor processivity that causes template switching, a hallmark of reverse transcription. HIV-1 RT inhibitors acting on its polymerase activity constitute the backbone of current antiretroviral therapies, although novel drugs, including ribonuclease H inhibitors, are still necessary to fight HIV infections. In Hepadnaviridae and Caulimoviridae, reverse transcription leads to the formation of nicked circular DNAs that will be converted into episomal DNA in the host cell nucleus. Structural and biochemical information on their polymerases is limited, although several drugs inhibiting HIV-1 RT are known to be effective against the human hepatitis B virus polymerase. In this review, we summarize current knowledge on reverse transcription in the five virus families and discuss available biochemical and structural information on RTs, including their biosynthesis, enzymatic activities, and potential inhibition.

Effect of Glycosylation on an Immunodominant Region in the V1V2 Variable Domain of the HIV-1 Envelope gp120 Protein

Heavy glycosylation of the envelope (Env) surface subunit, gp120, is a key adaptation of HIV-1; however, the precise effects of glycosylation on the folding, conformation and dynamics of this protein are poorly understood. Here we explore the patterns of HIV-1 Env gp120 glycosylation, and particularly the enrichment in glycosylation sites proximal to the disulfide linkages at the base of the surface-exposed variable domains. To dissect the influence of glycans on the conformation these regions, we focused on an antigenic peptide fragment from a disulfide bridge-bounded region spanning the V1 and V2 hyper-variable domains of HIV-1 gp120. We used replica exchange molecular dynamics (MD) simulations to investigate how glycosylation influences its conformation and stability. Simulations were performed with and without N-linked glycosylation at two sites that are highly conserved across HIV-1 isolates (N156 and N160); both are contacts for recognition by V1V2-targeted broadly neutralizing antibodies against HIV-1. Glycosylation stabilized the pre-existing conformations of this peptide construct, reduced its propensity to adopt other secondary structures, and provided resistance against thermal unfolding. Simulations performed in the context of the Env trimer also indicated that glycosylation reduces flexibility of the V1V2 region, and provided insight into glycan-glycan interactions in this region. These stabilizing effects were influenced by a combination of factors, including the presence of a disulfide bond between the Cysteines at 131 and 157, which increased the formation of beta-strands. Together, these results provide a mechanism for conservation of disulfide linkage proximal glycosylation adjacent to the variable domains of gp120 and begin to explain how this could be exploited to enhance the immunogenicity of those regions. These studies suggest that glycopeptide immunogens can be designed to stabilize the most relevant Env conformations to focus the immune response on key neutralizing epitopes.

Heavy glycosylation of the envelope surface subunit, gp120, is a key adaptation of HIV-1, however, the precise effects of glycosylation on the folding, conformation and dynamics of this protein are poorly understood. The network of glycans on gp120 is of particular interest with regards to vaccine design, because the glycans both serve as targets for many classes of broadly neutralizing antibodies, and contribute to patterns of immune evasion and escape during HIV-1 infection. In this manuscript, we report on how glycosylation influences an immunogenic but disordered region of gp120. Glycosylation stabilizes the pre-existing conformation, and reduces its propensity to form other secondary structures. It also stabilizes preformed conformation against thermal unfolding. These complementary effects originate from a combination of multiple factors, including the observation that having a glycosylation site adjacent to the disulfide bond further promotes the formation of beta-strand structure in this peptide.

Evaluating Clonal Expansion of HIV-Infected Cells: Optimization of PCR Strategies to Predict Clonality

In HIV-infected individuals receiving suppressive antiretroviral therapy, the virus persists indefinitely in a reservoir of latently infected cells. The proliferation of these cells may contribute to the stability of the reservoir and thus to the lifelong persistence of HIV-1 in infected individuals. Because the HIV-1 replication process is highly error-prone, the detection of identical viral genomes in distinct host cells provides evidence for the clonal expansion of infected cells. We evaluated alignments of unique, near-full-length HIV-1 sequences to determine the relationship between clonality in a short region and clonality in the full genome. Although it is common to amplify and sequence short, subgenomic regions of the viral genome for phylogenetic analysis, we show that sequence identity of these amplicons does not guarantee clonality across the full viral genome. We show that although longer amplicons capture more diversity, no subgenomic region can recapitulate the diversity of full viral genomes. Consequently, some identical subgenomic amplicons should be expected even from the analysis of completely unique viral genomes, and the presence of identical amplicons alone is not proof of clonally expanded HIV-1. We present a method for evaluating evidence of clonal expansion in the context of these findings.

Although antiretroviral therapy effectively blocks HIV-1 replication, the virus persists indefinitely in a reservoir of latently infected cells. This reservoir is a major barrier to HIV-1 cure. Recent studies have identified the proliferation of latently infected cells as a mechanism that may contribute to the lifelong persistence of HIV-1. In contrast with cellular proliferation, viral replication is highly error-prone; therefore, the detection of identical viral genomes in distinct host cells provides evidence that those cells are the progeny of an infected cell that underwent clonal expansion. For this reason, the accurate and reliable identification of identical HIV-1 genomes derived from distinct cells is critical for understanding HIV-1 persistence and advancing cure research. Studies of HIV-1 clonal expansion often present clonality of short, subgenomic sequences as evidence for the clonality of full viral genomes. In this study, we quantified the relationship between sequence identity in short regions and sequence identity of near-full-length genomes to demonstrate that no subgenomic region completely captures the diversity of the full viral genome. Consequently, identical subgenomic sequences are not proof of identical full-length HIV-1 genomes. In the context of these findings, we developed a method to evaluate identical subgenomic sequences as possible evidence for clonal expansion.

HIV-1 Protease, Reverse Transcriptase, and Integrase Variation

HIV-1 protease (PR), reverse transcriptase (RT), and integrase (IN) variability presents a challenge to laboratories performing genotypic resistance testing. This challenge will grow with increased sequencing of samples enriched for proviral DNA such as dried blood spots and increased use of next-generation sequencing (NGS) to detect low-abundance HIV-1 variants. We analyzed PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to characterize variation at each amino acid position, identify mutations indicating APOBEC-mediated G-to-A editing, and identify mutations resulting from selective drug pressure. Forty-seven percent of PR, 37% of RT, and 34% of IN positions had one or more amino acid variants with a prevalence of ≥1%. Seventy percent of PR, 60% of RT, and 60% of IN positions had one or more variants with a prevalence of ≥0.1%. Overall 201 PR, 636 RT, and 346 IN variants had a prevalence of ≥0.1%. The median intersubtype prevalence ratios were 2.9-, 2.1-, and 1.9-fold for these PR, RT, and IN variants, respectively. Only 5.0% of PR, 3.7% of RT, and 2.0% of IN variants had a median intersubtype prevalence ratio of ≥10-fold. Variants at lower prevalences were more likely to differ biochemically and to be part of an electrophoretic mixture compared to high-prevalence variants. There were 209 mutations indicative of APOBEC-mediated G-to-A editing and 326 mutations nonpolymorphic treatment selected. Identification of viruses with a high number of APOBEC-associated mutations will facilitate the quality control of dried blood spot sequencing. Identifying sequences with a high proportion of rare mutations will facilitate the quality control of NGS.

IMPORTANCE Most antiretroviral drugs target three HIV-1 proteins: PR, RT, and IN. These proteins are highly variable: many different amino acids can be present at the same position in viruses from different individuals. Some of the amino acid variants cause drug resistance and occur mainly in individuals receiving antiretroviral drugs. Some variants result from a human cellular defense mechanism called APOBEC-mediated hypermutation. Many variants result from naturally occurring mutation. Some variants may represent technical artifacts. We studied PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to quantify variation at each amino acid position in these three HIV-1 proteins. We performed analyses to determine which amino acid variants resulted from antiretroviral drug selection pressure, APOBEC-mediated editing, and naturally occurring variation. Our results provide information essential to clinical, research, and public health laboratories performing genotypic resistance testing by sequencing HIV-1 PR, RT, and IN.

Minimal Contribution of APOBEC3-Induced G-to-A Hypermutation to HIV-1 Recombination and Genetic Variation

Although the predominant effect of host restriction APOBEC3 proteins on HIV-1 infection is to block viral replication, they might inadvertently increase retroviral genetic variation by inducing G-to-A hypermutation. Numerous studies have disagreed on the contribution of hypermutation to viral genetic diversity and evolution. Confounding factors contributing to the debate include the extent of lethal (stop codon) and sublethal hypermutation induced by different APOBEC3 proteins, the inability to distinguish between G-to-A mutations induced by APOBEC3 proteins and error-prone viral replication, the potential impact of hypermutation on the frequency of retroviral recombination, and the extent to which viral recombination occurs in vivo, which can reassort mutations in hypermutated genomes. Here, we determined the effects of hypermutation on the HIV-1 recombination rate and its contribution to genetic variation through recombination to generate progeny genomes containing portions of hypermutated genomes without lethal mutations. We found that hypermutation did not significantly affect the rate of recombination, and recombination between hypermutated and wild-type genomes only increased the viral mutation rate by 3.9 × 10-5 mutations/bp/replication cycle in heterozygous virions, which is similar to the HIV-1 mutation rate. Since copackaging of hypermutated and wild-type genomes occurs very rarely in vivo, recombination between hypermutated and wild-type genomes does not significantly contribute to the genetic variation of replicating HIV-1. We also analyzed previously reported hypermutated sequences from infected patients and determined that the frequency of sublethal mutagenesis for A3G and A3F is negligible (4 × 10-21 and1 × 10-11, respectively) and its contribution to viral mutations is far below mutations generated during error-prone reverse transcription. Taken together, we conclude that the contribution of APOBEC3-induced hypermutation to HIV-1 genetic variation is substantially lower than that from mutations during error-prone replication.