Characterizing the phenotypic and genetic changes of pre-epidemic HIV-2 group F virus following serial passage in humanized mice.

HIV-2 Group F virus with an origin in NHPs was isolated from only two individuals. Two serial passages in hu-mice showed increased viral loads, CD4+ T cell decline and nonsynonymous genetic changes showing its capacity for further evolution, and spread in the human.

Evolution of SIVmac239 following serial passaging in humanized mice.

Serial passage of SIVmac239 allows for greater understanding of the genetic changes necessary for cross-species transmission of primate lentiviruses into humans. Using humanized mice, we show that adaptive mutations continue to accumulate in SIVmac239 during four serial passages, with persistent CD4+ T cell decline and increases in plasma viral loads.

Longitudinal within-host evolution of HIV Nef-mediated CD4, HLA and SERINC5 downregulation activity: a case study.

The HIV accessory protein Nef downregulates the viral entry receptor CD4, the Human Leukocyte Antigen (HLA)-A and -B molecules, the Serine incorporator 5 (SERINC5) protein and other molecules from the infected cell surface, thereby promoting viral infectivity, replication and immune evasion. The nef locus also represents one of the most genetically variable regions in the HIV genome, and nef sequences undergo substantial evolution within a single individual over the course of infection. Few studies however have simultaneously characterized the impact of within-host nef sequence evolution on Nef protein function over prolonged timescales. Here, we isolated 50 unique Nef clones by single-genome amplification over an 11-year period from the plasma of an individual who was largely naïve to antiretroviral treatment during this time. Together, these clones harbored nonsynonymous substitutions at 13% of nef's codons. We assessed their ability to downregulate cell-surface CD4, HLA and SERINC5 and observed that all three Nef functions declined modestly over time, where the reductions in CD4 and HLA downregulation (an average of 0.6% and 2.0% per year, respectively) achieved statistical significance. The results from this case study support all three Nef activities as being important to maintain throughout untreated HIV infection, but nevertheless suggest that, despite nef's mutational plasticity, within-host viral evolution can compromise Nef function, albeit modestly, over prolonged periods.

The PTAP sequence duplication in HIV-1 subtype C Gag p6 in drug-naive subjects of India and South Africa.

BACKGROUND: HIV-1 subtype C demonstrates several biological properties distinct from other viral subtypes. One such variation is the duplication of PTAP motif in p6 Gag. PTAP motif is a key player in viral budding. Here, we studied the prevalence of PTAP motif duplication in subtype C viral strains in a longitudinal study. METHODS: In a prospective follow-up study, 65 HIV-1 seropositive drug-naive subjects were monitored in two different clinical cohorts of India for 2 years with repeated sampling at 6-month intervals. The viral RNA was extracted from plasma, the gag segment was amplified and sequenced. From a subset of viral isolates the sequences of pol, env and LTR were sequenced. Using HIV-1 gag amino acid sequences available from public databases and additional sequences derived from the Indian and South-African cohorts, we examined the nature of PTAP motif duplication in subtype C. RESULTS: In 16% (8 of 50) of the primary viral strains of India, we identified a sequence duplication of the PTAP motif in Gag p6. The length of the sequence duplication varied from 6 to 14 amino acids in the viral isolates but remained fixed within a subject over a period of 24-36 month follow-up. In the duplicated motif, the core PTAP motif was invariable, but the flanking residues were highly variable. In an acute phase clinical cohort of South Africa, in a subset of 75 subjects, we found the presence of the PTAP duplication at a frequency of 29.3%. An analysis of the gag sequences from the extant databases showed that unlike other subtypes of HIV-1, subtype C has a natural propensity to generate the PTAP motif duplication at a significantly higher frequency and of greater length. Additionally, the global prevalence of PTAP duplication in subtype C appears to be increasing progressively over the past 30 years. CONCLUSION: We showed that in subtype C, the duplication of the PTAP motif in p6 Gag involves sequence stretches of greater length, and at a much higher frequency as compared to other HIV-1 subtypes. Given that subtype C naturally lacks the Alix binding motif, the acquisition of an additional PTAP motif may confer replication advantage on this HIV-1 subtype. Further investigation is warranted to examine the significance of PTAP motif duplication on the replicative fitness of HIV-1.

No Difference in the Prevalence of HIV-1 gag Cytotoxic T-Lymphocyte-Associated Escape Mutations in Viral Sequences from Early and Late Parts of the HIV-1 Subtype C Pandemic in Botswana.

HIV is known to accumulate escape mutations in the gag gene in response to the immune response from cytotoxic T lymphocytes (CTLs). These mutations can occur within an individual as well as at a population level. The population of Botswana exhibits a high prevalence of HLA*B57 and HLA*B58, which are associated with effective immune control of HIV. In this retrospective cross-sectional investigation, HIV-1 gag gene sequences were analyzed from recently infected participants across two time periods which were 10 years apart: the early time point (ETP) and late time point (LTP). The prevalence of CTL escape mutations was relatively similar between the two time points-ETP (10.6%) and LTP (9.7%). The P17 protein had the most mutations (9.4%) out of the 36 mutations that were identified. Three mutations (A83T, K18R, Y79H) in P17 and T190A in P24 were unique to the ETP sequences at a prevalence of 2.4%, 4.9%, 7.3%, and 5%, respectively. Mutations unique to the LTP sequences were all in the P24 protein, including T190V (3%), E177D (6%), R264K (3%), G248D (1%), and M228L (11%). Mutation K331R was statistically higher in the ETP (10%) compared to the LTP (1%) sequences (p < 0.01), while H219Q was higher in the LTP (21%) compared to the ETP (5%) (p < 0.01). Phylogenetically, the gag sequences clustered dependently on the time points. We observed a slower adaptation of HIV-1C to CTL immune pressure at a population level in Botswana. These insights into the genetic diversity and sequence clustering of HIV-1C can aid in the design of future vaccine strategies.

Early-Transmitted Variants and Their Evolution in a HIV-1 Positive Couple: NGS and Phylogenetic Analyses.

We had access to both components of a couple who became infected with human immunodeficiency virus (HIV)-1 through sexual behavior during the early initial phase of infection and before initiation of therapy. We analyzed blood samples obtained at the time of diagnosis and after six months of combined antiretroviral therapy. Next-generation sequencing (NGS) and phylogenetic analyses were used to investigate the transmission and evolution of HIV-1 quasispecies. Phylogenetic analyses were conducted using Bayesian inference methods. Both partners were infected with an HIV-1 B subtype. No evidence of viral recombination was observed. The lowest intrapersonal genetic distances were observed at baseline, before initiation of therapy, and in particular in the V1V2 fragment (distances ranging from 0.102 to 0.148). One HIV-1 single variant was concluded to be dominant in all of the HIV-1 regions analyzed, although some minor variants could be observed. The same tree structure was observed both at baseline and after six months of therapy. These are the first extended phylogenetic analyses performed on both members of a therapy-naïve couple within a few weeks of infection, and in which the effect of antiretroviral therapy on viral evolution was analyzed. Understanding which HIV-1 variants are most likely to be transmitted would allow a better understanding of viral evolution, possibly playing a role in vaccine design and prevention strategies.

Emergence of HIV-1 Unique DG Recombinant Form in Pakistan.

After the first case of HIV infection in 1987, the number of cases has continuously increased in Pakistan, turning from isolated incidents to outbreaks to the concentrated epidemic. The HIV epidemic in Pakistan is mainly driven by subtype A; however, the overlapping transmission chains facilitate recombination between subtypes and existing circulating recombinants forms (CRFs), leading to the emergence of unique recombinant forms (URFs). In this study, we report the first case of a URF (URF_DG) in a Pakistani HIV-infected patient. Phylogenetic and drug resistance analysis of the patient-derived sequence indicated that Pakistani URF_DG sequence was closely related to the URF_DG sequence reported from the United Kingdom, but had more drug resistance mutations than the U.K. sequence.

Modelling and in vitro testing of the HIV-1 Nef fitness landscape.

An effective vaccine is urgently required to curb the HIV-1 epidemic. We have previously described an approach to model the fitness landscape of several HIV-1 proteins, and have validated the results against experimental and clinical data. The fitness landscape may be used to identify mutation patterns harmful to virus viability, and consequently inform the design of immunogens that can target such regions for immunological control. Here we apply such an analysis and complementary experiments to HIV-1 Nef, a multifunctional protein which plays a key role in HIV-1 pathogenesis. We measured Nef-driven replication capacities as well as Nef-mediated CD4 and HLA-I down-modulation capacities of thirty-two different Nef mutants, and tested model predictions against these results. Furthermore, we evaluated the models using 448 patient-derived Nef sequences for which several Nef activities were previously measured. Model predictions correlated significantly with Nef-driven replication and CD4 down-modulation capacities, but not HLA-I down-modulation capacities, of the various Nef mutants. Similarly, in our analysis of patient-derived Nef sequences, CD4 down-modulation capacity correlated the most significantly with model predictions, suggesting that of the tested Nef functions, this is the most important in vivo. Overall, our results highlight how the fitness landscape inferred from patient-derived sequences captures, at least in part, the in vivo functional effects of mutations to Nef. However, the correlation between predictions of the fitness landscape and measured parameters of Nef function is not as accurate as the correlation observed in past studies for other proteins. This may be because of the additional complexity associated with inferring the cost of mutations on the diverse functions of Nef.

The Spleen Is an HIV-1 Sanctuary During Combined Antiretroviral Therapy.

Combined antiretroviral therapy (cART) does not eradicate HIV, which persists for years and can re-establish replication if treatment is stopped. The current challenge is identifying those tissues harboring virus through cART. Here, we used HIV env-nef single genome sequencing and HIV gag droplet digital PCR (ddPCR) to survey 50 tissues from five subjects on cART with no detectable plasma viral load at death. The spleen most consistently contained multiple proviral and expressed sequences (4/5 participants). Spleen-derived HIV demonstrated two distinct phylogenetic patterns: multiple identical sequences, often from different tissues, as well as diverse viral sequences on long terminal branches. Our results suggested that ddPCR may overestimate the size of the tissue-based viral reservoir. The spleen, a lymphatic organ at the intersection of the immune and circulatory systems, may play a key role in viral persistence.

Longitudinal sequencing of HIV-1 infected patients with low-level viremia for years while on ART shows no indications for genetic evolution of the virus.

HIV-infected patients on antiretroviral therapy (ART) may present low-level viremia (LLV) above the detection level of current viral load assays. In many cases LLV is persistent but does not result in overt treatment failure or selection of drug resistant viral variants. To elucidate whether LLV reflects active virus replication, we extensively sequenced pol and env genes of the viral populations present before and during LLV in 18 patients and searched for indications of genetic evolution. Maximum likelihood phylogenetic trees were inspected for temporal structure both visually and by linear regression analysis of root-to-tip and pairwise distances. Viral coreceptor tropism was assessed at different time points before and during LLV. In none of the patients consistent indications for genetic evolution were found over a median period of 4.8 years of LLV. As such these findings could not provide evidence that active virus replication is the main driver of LLV.

Ongoing HIV Replication During ART Reconsidered.

Lorenzo-Redondo et al. recently analyzed HIV RNA sequences in plasma virus and proviral DNA sequences in lymph nodes (LN) and peripheral blood mononuclear cells (PBMC) from samples collected over a 6-month period from 3 individuals following initiation of antiretroviral therapy (ART) and concluded that ongoing HIV replication occurred in LN despite ART and that this replication maintained the HIV reservoir. We analyzed the same sequences and found that the dataset was very limited (median of 5 unique RNA or DNA sequences per sample) after accounting for polymerase chain reaction resampling and hypermutation and that the few remaining DNA sequences after 3 and 6 months on ART were not more diverse or divergent from those in pre-ART in any of the individuals studied. These findings, and others, lead us to conclude that the claims of ongoing replication on ART made by Lorenzo-Redondo et al. are not justified from the dataset analyzed in their publication.

Coalescent Inference Using Serially Sampled, High-Throughput Sequencing Data from Intrahost HIV Infection.

Human immunodeficiency virus (HIV) is a rapidly evolving pathogen that causes chronic infections, so genetic diversity within a single infection can be very high. High-throughput "deep" sequencing can now measure this diversity in unprecedented detail, particularly since it can be performed at different time points during an infection, and this offers a potentially powerful way to infer the evolutionary dynamics of the intrahost viral population. However, population genomic inference from HIV sequence data is challenging because of high rates of mutation and recombination, rapid demographic changes, and ongoing selective pressures. In this article we develop a new method for inference using HIV deep sequencing data, using an approach based on importance sampling of ancestral recombination graphs under a multilocus coalescent model. The approach further extends recent progress in the approximation of so-called conditional sampling distributions, a quantity of key interest when approximating coalescent likelihoods. The chief novelties of our method are that it is able to infer rates of recombination and mutation, as well as the effective population size, while handling sampling over different time points and missing data without extra computational difficulty. We apply our method to a data set of HIV-1, in which several hundred sequences were obtained from an infected individual at seven time points over 2 years. We find mutation rate and effective population size estimates to be comparable to those produced by the software BEAST. Additionally, our method is able to produce local recombination rate estimates. The software underlying our method, Coalescenator, is freely available.

Cladograms with Path to Event (ClaPTE): a novel algorithm to detect associations between genotypes or phenotypes using phylogenies.

BACKGROUND: Associations between genotype and phenotype provide insight into the evolution of pathogenesis, drug resistance, and the spread of pathogens between hosts. However, common ancestry can lead to apparent associations between biologically unrelated features. The novel method Cladograms with Path to Event (ClaPTE) detects associations between character-pairs (either a pair of mutations or a mutation paired with a phenotype) while adjusting for common ancestry, using phylogenetic trees. METHODS: ClaPTE tests for character-pairs changing close together on the phylogenetic tree, consistent with an associated character-pair. ClaPTE is compared to three existing methods (independent contrasts, mixed model, and likelihood ratio) to detect character-pair associations adjusted for common ancestry. Comparisons utilize simulations on gene trees for: HIV Env, HIV promoter, and bacterial DnaJ and GuaB; and case studies for Oseltamavir resistance in Influenza, and for DnaJ and GuaB. Simulated data include both true-positive/associated character-pairs, and true-negative/not-associated character-pairs, used to assess type I (frequency of p-values in true-negatives) and type II (sensitivity to true-positives) error control. RESULTS AND CONCLUSIONS: ClaPTE has competitive sensitivity and better type I error control than existing methods. In the Influenza/Oseltamavir case study, ClaPTE reports no new permissive mutations but detects associations between adjacent (in primary sequence) amino acid positions which other methods miss. In the DnaJ and GuaB case study, ClaPTE reports more frequent associations between positions both from the same protein family than between positions from different families, in contrast to other methods. In both case studies, the results from ClaPTE are biologically plausible.

Inferring HIV Escape Rates from Multi-Locus Genotype Data.

Cytotoxic T-lymphocytes (CTLs) recognize viral protein fragments displayed by major histocompatibility complex molecules on the surface of virally infected cells and generate an anti-viral response that can kill the infected cells. Virus variants whose protein fragments are not efficiently presented on infected cells or whose fragments are presented but not recognized by CTLs therefore have a competitive advantage and spread rapidly through the population. We present a method that allows a more robust estimation of these escape rates from serially sampled sequence data. The proposed method accounts for competition between multiple escapes by explicitly modeling the accumulation of escape mutations and the stochastic effects of rare multiple mutants. Applying our method to serially sampled HIV sequence data, we estimate rates of HIV escape that are substantially larger than those previously reported. The method can be extended to complex escapes that require compensatory mutations. We expect our method to be applicable in other contexts such as cancer evolution where time series data is also available.