- Open Access
Highly divergent patterns of genetic diversity and evolution in proviral quasispecies from HIV controllers
Retrovirology volume 14, Article number: 29 (2017)
Ongoing intra-host HIV-1 evolution has been shown in individuals that naturally suppress the viremia to low levels (HIV controllers) by the analysis of the RNA in plasma compartment. Detection of evolution at the DNA proviral compartment in HIV controllers, however, has been more challenging and the precise correlation between the systemic viral suppression level and rate of reservoir’s reseeding in those individuals is not fully understood. In this sense, we examined the proviral DNA quasispecies by single genome amplification of the env gene in a cohort of 23 HIV controllers from Brazil, divided in three groups, according to the level of systemic viral suppression: (1) elite controllers with persistent undetectable viral load (PEC, n = 6); (2) elite controllers with occasional episodes of transient (51–400 copies/mL) viremia (EEC, n = 7); and (3) viremic controllers with persistent low-level (80–2000 copies/mL) viremia (VC, n = 10).
The HIV-1 diversity of the PBMC-associated proviral quasispecies in EC was significantly (P < 0.01) lower than in VC, but not significantly different between PEC and EEC groups. We detected a considerable variation in the average pairwise nucleotide distance and proportion of unique sequences in the HIV-1 proviral quasispecies of PEC and EEC. Some PEC and EEC displayed highly homogenous proviral populations with large clusters of identical sequences, while others exhibited relatively diverse proviral populations with a high proportion of unique sequences comparable to VC subjects. The long-term (10–15 years) follow-up of the HIV-1 proviral populations revealed a complete evolutionary stasis in one PEC and measurable divergence rates in one EEC [3.1 (1.2–5.6) × 10−3 substitutions/site/year and one VC [2.9 (0.7–5.1) × 10−3 substitutions/site/year].
There is no simple relationship between systemic viral suppression and intra-host proviral diversity or rate of reservoir’s reseeding in chronically infected HIV controllers. Our results demonstrate that very divergent patterns of intra-host viral diversity and divergence could be detected in the setting of natural suppression of HIV-1 replication and that ongoing evolution and reseeding of the PBMC proviral reservoir occurs in some elite controllers.
The natural history of human immunodeficiency virus type-1 (HIV-1) infections may display very divergent patterns among individuals. Most HIV-1 infected individuals, termed typical progressors (TP), display high plasma viral loads and progress to AIDS without treatment after 5–10 years of infection . Some individuals, termed long-term non-progressors (LTNPs), display longer asymptomatic periods (>10 years) and keep normal CD4+ T cell counts in the absence of treatment ; while others, termed HIV controllers, exhibit a durable control of viral replication maintaining at very low levels during chronic infection . Among HIV controllers, the viremic controllers (VC) suppress the viremia to levels <2000 HIV-1 RNA copies/mL and the elite controllers (EC) to levels <50-80 HIV-1 RNA copies/mL.
Intra-host HIV-1 evolution in TP follows a consistent pattern of temporal changes in viral diversity and divergence during the course of infection, that affect both proviral DNA populations in peripheral blood mononuclear cells (PBMC) and viral RNA populations in plasma . According to that pattern, infection is usually initiated by a relatively homogeneous viral population (with less than 1% envelope [env] diversity) that diversifies during the asymptomatic phase, reaching a peak of population diversity (up to 10% at the env gene) and divergence before leveling off or decrease towards the AIDS phase. A roughly similar pattern of intra-host HIV-1 evolution was described for LTNPs and HIV controllers in the plasma compartment . LTNPs display HIV-1 RNA populations that continuously evolve during chronic infection and reach an overall diversity comparable to that observed in TP . Several studies also demonstrate ongoing evolution and divergence of HIV-1 RNA sequences from most EC [5,6,7,8,9,10,11,12], although the mean diversity of plasma populations in EC is significantly lower than that observed for TP at chronic infection .
The HIV-1 diversity and divergence pattern of PBMC-associated proviral sequences from LTNPs and HIV controllers, however, differed strikingly from that observed in the plasma virus. In some LTNPs, DNA proviral populations are composed of a complex mixture of archival (dating close to the patient’s seroconversion time) and recent (dating close to the sampling time) variants  and displayed no temporal structure in the changes of diversity and divergence during chronic infection . In all chronically infected EC and some VC, DNA proviral populations are extremely homogenous (with less than 2% env diversity), mostly composed by ancestral sequences and with no measurable divergence over time [5, 9, 10, 12, 15,16,17,18,19]. A recent study demonstrates that most proviral sequences detected in PBMC from HIV controllers are largely representative of archival variants probably integrated during primary infection and propagated by clonal expansion of the memory CD4+ T cell latent reservoir, although rare proviral clones of recent origin could be detected in some patients .
These observations suggest that the virus is evolving in HIV controllers, but most evolving plasma viruses do not replenish the PBMC reservoir and the majority of PBMC-associated proviral sequences detected in chronically infected HIV controllers represent ancestral variants. The precise correlation between the systemic viral suppression level and the rate of reservoir’s reseeding in HIV controllers, however, is not fully understood. Furthermore, previous studies may have failed to detect proviral sequence replenishment and ongoing evolution in HIV controllers because of the narrow follow-up time (usually 2–6 years). To answer these questions, we performed a cross-sectional analysis of the DNA proviral quasispecies diversity at the env gene in 23 HIV controllers with different levels of systemic viral suppression and we also recover the long-term (10–15 years) pattern of changes of HIV-1 proviral populations in the setting of low/undetectable viremia.
A cohort of 23 HIV controllers, defined as subjects infected with HIV-1 for at least 5 years and maintaining RNA viral loads of <2000 copies/mL without antiretroviral therapy, has been followed-up at the Instituto Nacional de Infectologia Evandro Chagas (INI), Rio de Janeiro, Brazil. These subjects were classified in two categories according to the plasmatic viral load (VL) during follow-up : (1) elite controllers (EC) if most (≥70%) plasma VL determinations were below the limit of detection for the respective available assay (<50–80 copies/mL) (n = 13) and (2) viremic controllers (VC) if most (≥70%) VL determinations were between 80 and 2000 copies/mL (n = 10). The EC were further subdivided in two subgroups : persistent elite controllers (PEC) if 100% of VL measures were below the limit of detection (n = 6) and (2) ebbing elite controllers (EEC) if subjects had occasional (<30% of frequency) episodes of transient low-level (51–400 copies/mL) viremia (n = 7). Patients were followed at least once every 6–12 months to perform infection-monitoring tests such as RNA viral load quantification and CD4+ T lymphocyte count. In each visit, PBMC were obtained by Histopaque-1077 (Sigma, USA) density gradient and stored in liquid nitrogen until use. The present work was approved by the Brazilian National Human Research Ethics Committee (CONEP 14430/2011) and all subjects gave written informed consent.
CD4+ T cell counts and plasma HIV-1 RNA quantification
Absolute CD4+ T cell counts were obtained using the MultiTest TruCount-kit and the MultiSet software on a FACSCalibur flow cytometer (BD Biosciences San Jose, CA). Plasma VL were measured according to the Brazilian Ministry of Health guidelines, with methodologies being updated overtime to improve sensitivity: Nuclisens HIV-1 RNA QT assay (Organon Teknika, Durham, NC, limit of detection: 80 copies/mL) from 1999 to 2007; the Versant HIV-1 3.0 RNA assay (bDNA 3.0, Siemens, Tarrytown, NY, limit of detection: 50 copies/mL) from 2007 to 2013; and the Abbott RealTime HIV-1 assay (Abbott Laboratories, Wiesbaden, Germany, limit of detection: 40 copies/mL) from 2013 to 2016.
HIV-1 DNA extraction and single genome sequencing
Cryopreserved PBMC were thawed, washed and immediately after, the total genomic DNA was isolated with addition of the DNAzol® Reagent (Invitrogen, USA) as described . To limit template resampling, single genome amplification (SGA) was performed by limiting dilution nested PCR at a concentration of DNA that would produce less than 40% of positive PCR reactions, providing a >70% probability that a positive PCR originates from a single molecule . A fragment of nearly 600 bp of the HIV-1 env gene (including the C2–C4 regions of gp120) was amplified by PCR using AmpliTaq Gold® 360 DNA Polymerase (Applied Biosystems, USA) as described . The final PCR products were purified using the Illustra GFX PCR DNA purification kit (GE Healthcare, USA) and directly sequenced using the ABI BigDye Terminator v.3.1 reaction Kit (Applied Biosystems, Foster City, CA) in an ABI PRISM 3100 automate sequencer (Applied Biosystem). Chromatograms were assembled into contigs using the SeqMan 7.0 software (DNASTAR Inc., Madison, WI). Sequences resulting from low-quality chromatograms, from chromatograms with double peaks (indicative of more than one template per sequencing reaction), or showing APOBEC3G/F-mediated hypermutation as determined using Hypermut software  were discarded.
Env sequences from HIV controllers were aligned with HIV-1 subtype reference sequences using ClustalW and then manually edited, yielding a final alignment covering positions 7008–7650 relative to the HXB2 reference genome. Maximum-likelihood (ML) phylogenetic trees were reconstructed with the PhyML 3.0 program  using the most appropriate nucleotide substitution model selected using program jModeltest v. 3.7 , the SPR branch swapping heuristic tree search algorithm, and the approximate likelihood-ratio test (aLRT)  for branch support.
Prediction of coreceptor usage and CCR5 genotyping
The V3 region of env sequences was translated using MEGA7  and viral tropism was predicted using Geno2pheno (http://coreceptor.bioinf.mpisb.mpg.de/cgi-bin/coreceptor.pl) with a false positive rate (FPR) cutoff of 5% . The presence of the ∆32 variant in CCR5 was assessed by PCR amplification/agarose gel electrophoresis as previously described .
Analyses of viral diversity and divergence
The complexity of proviral quasispecies was characterized using two indices: the mean nucleotide diversity (π) and the normalized Shannon entropy (H SN). The π measures the average number of nucleotide differences between any two sequences of the quasispecies obtained at the same time point and was calculated using MEGA7  as described previously . The H SN provides a measure of haplotype (mutant) frequencies and was calculated by using the R package, Vegan , after rarefaction of samples to the small sample size (n = 10) for bias correction of sample size differences . The divergence rate of proviral env sequences was estimated for three patients (one PEC, one EEC and one VC) with available sequences sampled between 5 and 15 years ago . ML phylogenetic trees were reconstructed for each patient as described above and linear regression analysis of the root-to-tip distances against sampling time were performed using program Tempest  to verify the temporal structure of the datasets. The intra-host viral evolutionary (divergence) rate was then directly estimated from the sampling date of the sequences for those datasets with a good temporal structure using program BEAST v1.8 . Analyses were performed using the most appropriate nucleotide substitution model for each patient, a relaxed uncorrelated lognormal molecular clock model  with a CTMC rate reference prior  and a Bayesian coalescent tree prior . Three MCMC chains were run for 1 × 107 generations and then combined. Effective Sample Size (ESS) and 95% Highest Probability Density (HPD) values were inspected using Tracer v1.6 (http://tree.bio.ed.ac.uk/software/tracer/) to assess the convergence and uncertainty of parameter estimates.
Statistical analyses were performed using GraphPad v6 (Prism Software, USA). The Mann–Whitney test was use to compare the quasispecies diversity, the time since HIV-1 diagnosis and the CD4+ T cell counts between subjects groups. Tests were considered significant if the P value was ≤0.05.
Epidemiological, clinical and virological characteristics of HIV controllers
The main clinical and epidemiological characteristics of our HIV controllers’ cohort are shown in Table 1. Female gender (61%) was more frequent than male (39%), 70% of the patients identified themselves as heterosexual and 22% as men who have sex with men (MSM), while information regarding exposure behavior was not available for 9%. A higher proportion of females (77 vs 40%) and heterosexuals (69 vs 40%) was observed in the EC group than in the VC group. Participants had a median age of 49 years (IQR: 41–53 years old) and had documented HIV infection for a median of 11 years (IQR: 6–15 years). The EC and VC groups have a similar median age (52 vs 46 years, respectively) and median documented time of HIV infection (9 vs 10 years, respectively) at sampling time (Additional file 1: Figure S1A, B). None of the HIV controllers exhibited AIDS-related conditions and the CD4+ T cell counts were ≥500 cells/µL during follow-up (Fig. 1). Most of them (83%) also had documented HIV-infection for over 8 years, thus being classified as LTNPs. EC, however, displayed a higher median CD4+ T-cell count than VC at sampling time (1202 vs 735 cells/μL, respectively) (Additional file 1: Figure S1C). No significant differences in clinical and epidemiological characteristics were observed between PEC and EEC subgroups (data not shown).
A cross-sectional analysis of the HIV-1 proviral quasispecies in the 23 HIV controllers was performed by SGA of the env gene at between 5 and 20 years after HIV-diagnosis (Fig. 1). Similar median numbers of env clones per sample were obtained in EC (16, IQR: 14–18) and VC (18, IQR: 14–21) groups, as well as in PEC (15; IQR: 14–16) and EEC (17; IQR: 15–21) subgroups (Table 2). ML phylogenetic analysis revealed that env sequences from most individuals (n = 20, 87%) clustered by subject in highly supported (bootstrap >95%) monophyletic lineages (Fig. 2), thus supporting infection by a single variant. For three individuals (EEC09, VC06 and VC32), however, the env sequences branched in two independent monophyletic clades (Fig. 2), indicating dual infection. The subject VC06 was double infected by two HIV-1 subtype B variants, while individual EEC09 and VC32 were double infected by HIV-1 subtypes B and F1 variants. A second sample from these three individuals was analyzed confirming the previous result (data not shown). Overall, subtype B (74%) was the most frequent HIV-1 variant detected in our cohort, followed by subtypes F1 (13%), C (9%), and A1 (4%) (Fig. 2; Table 2). Prediction of coreceptor usage showed that most individuals (78%) presented only R5-tropic viral clones, two individuals (one EEC and one VC) presented a low frequency (5%) of X4-tropic clones, two individuals (one PEC and one VC) displayed a high frequency (30–45%) of X4-tropic viruses, and one PEC had only X4-tropic viral clones (Table 2). None of the subjects with high frequency of X4-tropic viruses is homozygous/heterozygous for the CCR5/Δ32 genotype (Table 2). Hypermutated proviral sequences were detected at a very low frequency (<5%) in only two individuals (Fig. 2).
Diversity of proviral quasispecies in HIV controllers
To address the potential relationship between systemic viral suppression level and reservoir’s reseeding among the 23 HIV controllers of our cohort, we calculated π and H SN indices that measure the average pairwise nucleotide distance and the mutant frequencies (proportion of unique sequences) in the set of aligned sequences of each individual, respectively. For double-infected patients, only sequences of the prevalent HIV-1 variant were considered. VC displayed quite diverse (π > 2%) and complex (H SN > 0.90) proviral quasipecies that were mostly (>70%) composed by unique sequences (Table 2; Additional file 2: Figure S2). The overall mean π and H SN estimated for HIV-1 quasispecies in the VC group were significantly higher than those estimated for the EC group (P < 0.01), despite the fact that the time since HIV-diagnosis was comparable among groups (Fig. 3). This supports that the PBMC reservoir of VC display higher rate of evolution and reseeding than that of EC.
A closer inspection of the EC group, however, reveals that both diversity and complexity of HIV-1 quasispecies extensively varied among subjects (Table 2; Fig. 3). The combined analysis of π and H SN allow us to detect two divergent patterns of intra-host viral diversity within the EC group (Fig. 4; Additional file 2: Figure S2). The first group (G1) comprises five EC (two PEC and three EEC) that present highly homogenous (π < 1%) proviral quasispecies with large clusters of identical sequences (H SN < 0.6). The second group (G2) comprises eight EC (four PEC and four EEC) showing relatively diverse (π ≥ 2%) proviral populations with high proportion of unique sequences (H SN > 0.7), comparable to those observed among VC subjects. Similar median values of π and H SN were estimated for EC with or without detection of occasional viremia above the limit of detection (Additional file 3: Figure S3). Quasispecies diversity and complexity were also not correlated with time since HIV diagnosis in EC (Additional file 4: Figure S4).
Rates of evolution of proviral quasispecies in HIV controllers
The pattern of intra-host viral diversity observed in EC-G1 is consistent with amplification of viral reservoir mostly by clonal expansion of infected memory CD4+ T cells; whereas the pattern observed in EC-G2 and VC patients supports a continuous reseeding of the proviral reservoir. To confirm that hypothesis, we investigated the long-term evolution of the PBMC proviral compartment in three individuals from groups EC-G1 (PEC52), EC-G2 (EEC42) and VC (VC05) by combining the env proviral sequences obtained in the present study with those obtained from the same patients 10–13 years ago and that were described previously . ML phylogenetic trees were reconstructed for each patient and the root-to-tip distances were plotted against sampling time. Despite the very long follow-up time (13 years), proviral env sequences of patient PEC52 were mostly identical and with no evidence of increasing root-to-tip distance over time (Fig. 5), thus confirming absence of reseeding and evolution of the PBMC reservoir in this patient. All proviral env sequences from patient EEC42 and most (70%) env sequences from patient VC05 sampled at the most recent time-point, by contrast, were different from those sampled 11–12 years earlier and with clear evidence of evolution (increasing root-to-tip distance over time) (Fig. 5). To estimate the intra-host HIV-1 evolutionary rate in subjects EEC42 and VC05, env sequences from different time points were analyzed using the BEAST program. For each cluster of identical sequences, only those env sequences sampled at the earliest time point were retained to reduce the impact of latency on intra-host evolutionary rate estimations. According to these analyses, the mean intra-host evolutionary rate of proviral env sequences estimated for patient EEC42 was 3.1 × 10−3 subst/site/year (95% HPD: 1.2–5.6 × 10−3 subst/site/year) and for patient VC was 2.9 × 10−3 subst/site/year (95% HPD: 0.7–5.1 × 10−3 subst/site/year).
In this study, we examined the DNA proviral quasispecies diversity at the env gene in 23 chronically infected HIV controllers with different levels of systemic viral suppression. Most HIV controllers included in our cohort were females (61%) and this percentage was higher for EC (77%) than for VC (40%). This may be a consequence of the greater frequency of women seeking health services for routine and preventive exams than men, enabling the diagnosis of HIV even in the absence of symptoms , and/or may reflect gender-specific differences in the plasma HIV-1 RNA levels [38,39,40]. Although no HIV controllers exhibited AIDS-related conditions and had CD4+ T cell counts ≥500 cells/µL during follow-up, EC (1202 cells/μL) displayed a significantly higher median CD4+ T cell counts than VC (735 cells/μL) at sampling time, supporting the relevance of persistent low-level viremia on the long-term CD4+ T cell decline [20, 41].
Analysis of proviral env sequences from HIV controllers revealed a diverse molecular epidemiologic profile with detection of HIV-1 subtypes B (74%), F1 (13%), C (9%) and A1 (4%). While subtypes B, F1 and C are common HIV-1 clades circulating in Brazil , subtype A1 has been only described in one case . Three individuals (one EC and two VC) were dually infected with strains of the same (B) or different (B and F1) subtypes, resulting in a prevalence of dual HIV-1 infection (13%) comparable to that previously estimated in a Spanish cohort of LTNP-EC (20%) . Prediction of coreceptor usage further revealed a significant frequency (30–100%) of X4-tropic clones in proviral quasispecies of two PEC and one VC. Reanalysis of proviral env sequences from EC and VC already published [12, 16, 18] showed that high frequency (>30%) of X4-tropic clones is a rare phenomenon, being detected in only one out of 25 subjects analyzed (data not shown). These results demonstrate that natural suppression of HIV-1 viremia below 2000 copies/mL can be achieved in the context of either single or dual HIV-1 infections, regardless of the subtype and coreceptor usage of infecting virus.
The HIV-1 proviral population continuously diversifies during untreated asymptomatic infection, although the rate of diversification greatly varies among individuals. In TP with RNA viral loads above 10,000 copies/mL, env gene diversity increases at a mean rate of 1%/year and reaches a peak (π = 6–10%) after 5–10 years post-infection . High levels of env proviral diversity (π = 4-8%) have been also described in samples taken 10–15 years after HIV diagnosis from LTNPs with plasma viremia between 2000 and 10,000 copies/mL [4, 16, 17, 19, 45]. Much lower levels of env proviral diversity (π = 0.1–6%), by contrast, were detected here in samples taken between 5 and 20 years after HIV diagnosis from HIV controllers (RNA viral load lower than 2000 copies/mL). This is consistent with previous studies [15,16,17,18,19] and with the notion that no viral diversification is expected when the host immune response greatly reduces the HIV-1 replication limiting the selection of escape mutants .
A closer inspection of the quasispecies diversity in different HIV controller groups here studied, however, revealed a more complex scenario. Particularly, the mean env diversity of proviral quasispecies in EC subjects varied over a large range (0.1–4.6%) and two distinct patterns of intra-host viral diversity were observed in that group. While some EC subjects (EC-G1) displayed highly homogeneous proviral populations (π < 1%) mainly composed by large clusters of identical sequences (H SN < 0.6), other EC subjects (EC-G2) showed more diverse (π ≥ 2%) proviral populations comprising high proportions of unique sequences (H SN > 0.7), comparable to those observed in VC subjects. Thus, contrary to initial expectations, the presence of a highly homogenous PBMC-associated HIV-1 proviral population is not a common characteristic of all EC subjects and no linear correlation could be observed between proviral quasispecies diversity and systemic viral suppression in HIV controllers.
Analysis of the long-term evolution of proviral populations revealed that the distinct patterns of intra-host viral diversity observed in HIV controllers might reflect different driving forces for the maintenance of the viral reservoir. Proviral env sequences of individual PEC52 (EC-G1 group) taken over a period of 13 years were mostly identical and displayed no evidence of divergence over time, demonstrating that most PBMC-associated proviral sequences detected in this chronically infected HIV controller represent ancestral variants that persist for >10 years of infection. This pattern supports the notion that the proviral reservoir, in some EC subjects, is mostly maintained by the clonal expansion of CD4+ T lymphocytes. Those cells were probably latently infected at the initial stage of infection, culminating in the absence of evolution and the preservation of a highly homogenous proviral population, similar to those observed in the majority of acutely infected patients [47,48,49,50,51,52,53,54].
In sharp contrast to patient PEC52, proviral populations of subjects EEC42 (EC-G2 group) and VC05 (VC group) displayed an increasing divergence and a partial or complete replacement of sequence variants over time. Although the mean env intra-host divergence rate here estimated for HIV controllers (~3 × 10−3 subst/site/year) was much lower than that previously estimated for TP (~10 × 10−3 subst/site/year) , the pattern observed is fully consistent with a continuous reseeding of the PBMC proviral reservoir in those HIV controllers. While several studies already demonstrate ongoing evolution and divergence of HIV-1 RNA sequences from the plasma compartment in VC and EC [5,6,7,8,9,10,11,12], this is the first study to quantify the intra-host divergence rate of DNA proviral sequences in the setting of undetectable viremia. These observations demonstrate that the HIV-1 in VC and in some EC is not only evolving, but also that the PBMC reservoir is continuously being resseded at a low, but measurable, rate leading to the partial or complete substitution of ancestral variants over time.
The divergent patterns of genetic diversity and evolution of proviral populations from EC here observed may be due to: (1) different levels of systemic suppression, (2) diverse mechanisms of natural control of HIV-1 replication, and/or (3) differences in the transmitted virus populations. Although a previous study conducted by our group demonstrated that rare episodes of detectable viremia in EC are associated to higher levels of systemic immune activation and a stronger HIV-1 specific immune response , pointing to lower levels of systemic viral supression in EEC than in PEC, we found no significant difference in the quasispecies diversity between both EC subgroups. It is possible that EC-G1 subjects display more efficient control mechanisms, capable of limiting new rounds of infection, particularly in the lymph nodes, than those present in EC-G2 subjects. Finally, it is also possible that the high proviral diversity detected in some EC was not due to intra-host evolution, but was present since the beginning of infection. Indeed, it was demonstrated that a substantial fraction of subjects (20–30%) displayed heterogeneous (2–5% env diversity) proviral populations in PBMC before seroconversion, most likely resulting from transmission of multiple HIV-1 variants [47,48,49,50,51,52,53,54]. If differences observed arise from multiple underlying mechanisms, definition of homogenous EC subgroups could become increasing challenging as more subjects are characterized.
These results reveal that very divergent patterns of intra-host viral diversity and divergence could be detected in the setting of natural suppression of HIV-1 replication, suggesting that HIV-1 may evolve differently in every patient. We found no simple relationship between systemic viral suppression and intra-host proviral diversity or rate of reservoir’s reseeding in chronically infected HIV controllers, although the influence of some potential confounding factors such as the transmission of multiple HIV-1 variants in some EC cannot be ruled out. Our study also demonstrates that ongoing evolution and reseeding of the PBMC proviral reservoir is possible in some EC. The long-term longitudinal follow-up of more EC patients will be important to elucidate the major driving forces of the different intra-host evolutionary patterns here detected as well as their impact on the long-lasting control of HIV-1 replication and disease progression.
Buchbinder S, Vittinghoff E. HIV-infected long-term nonprogressors: epidemiology, mechanisms of delayed progression, and clinical and research implications. Microbes Infect. 1999;1:1113–20.
Deeks SG, Walker BD. Human immunodeficiency virus controllers: mechanisms of durable virus control in the absence of antiretroviral therapy. Immunity. 2007;27:406–16.
Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, et al. Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999;73:10489–502.
Bagnarelli P, Mazzola F, Menzo S, Montroni M, Butini L, Clementi M. Host-specific modulation of the selective constraints driving human immunodeficiency virus type 1 env gene evolution. J Virol. 1999;73:3764–77.
Bailey JR, Williams TM, Siliciano RF, Blankson JN. Maintenance of viral suppression in HIV-1-infected HLA-B*57+ elite suppressors despite CTL escape mutations. J Exp Med. 2006;203:1357–69.
Miura T, Brumme CJ, Brockman MA, Brumme ZL, Pereyra F, Block BL, et al. HLA-associated viral mutations are common in human immunodeficiency virus type 1 elite controllers. J Virol. 2009;83:3407–12.
Miura T, Brockman MA, Schneidewind A, Lobritz M, Pereyra F, Rathod A, et al. HLA-B57/B*5801 human immunodeficiency virus type 1 elite controllers select for rare gag variants associated with reduced viral replication capacity and strong cytotoxic T-lymphotye recognition. J Virol. 2009;83:2743–55.
Mens H, Kearney M, Wiegand A, Shao W, Schønning K, Gerstoft J, et al. HIV-1 continues to replicate and evolve in patients with natural control of HIV infection. J Virol. 2010;84:12971–81.
Salgado M, Brennan TP, O’Connell KA, Bailey JR, Ray SC, Siliciano RF, et al. Evolution of the HIV-1 nef gene in HLA-B*57 positive elite suppressors. Retrovirology. 2010;7:94.
O’Connell KA, Brennan TP, Bailey JR, Ray SC, Siliciano RF, Blankson JN. Control of HIV-1 in elite suppressors despite ongoing replication and evolution in plasma virus. J Virol. 2010;84:7018–28.
O’Connell KA, Rabi SA, Siliciano RF, Blankson JN. CD4+ T cells from elite suppressors are more susceptible to HIV-1 but produce fewer virions than cells from chronic progressors. Proc Natl Acad Sci USA. 2011;108:E689–98.
Boritz EA, Darko S, Swaszek L, Wolf G, Wells D, Wu X, et al. Multiple origins of virus persistence during natural control of HIV infection. Cell. 2016;166:1004–15.
Bello G, Casado C, García S, Rodríguez C, del Romero J, López-Galíndez C. Co-existence of recent and ancestral nucleotide sequences in viral quasispecies of human immunodeficiency virus type 1 patients. J Gen Virol. 2004;85:399–407.
Bello G, Casado C, García S, Rodríguez C, del Romero J, Carvajal-Rodriguez A, et al. Lack of temporal structure in the short term HIV-1 evolution within asymptomatic naive patients. Virology. 2007;362:294–303.
Wang B, Mikhail M, Dyer WB, Zaunders JJ, Kelleher AD, Saksena NK. First demonstration of a lack of viral sequence evolution in a nonprogressor, defining replication-incompetent HIV-1 infection. Virology. 2003;312:135–50.
Bello G, Casado C, Sandonis V, Alonso-Nieto M, Vicario JL, García S, et al. A subset of human immunodeficiency virus type 1 long-term non-progressors is characterized by the unique presence of ancestral sequences in the viral population. J Gen Virol. 2005;86:355–64.
Bello G, Casado C, Sandonis V, Alvaro-Cifuentes T, Dos Santos CAR, García S, et al. Plasma viral load threshold for sustaining intrahost HIV type 1 evolution. AIDS Res Hum Retrovir. 2007;23:1242–50.
Bailey JR, Lassen KG, Yang H-C, Quinn TC, Ray SC, Blankson JN, et al. Neutralizing antibodies do not mediate suppression of human immunodeficiency virus type 1 in elite suppressors or selection of plasma virus variants in patients on highly active antiretroviral therapy. J Virol. 2006;80:4758–70.
Sandonís V, Casado C, Alvaro T, Pernas M, Olivares I, García S, et al. A combination of defective DNA and protective host factors are found in a set of HIV-1 ancestral LTNPs. Virology. 2009;391:73–82.
Okulicz JF, Marconi VC, Landrum ML, Wegner S, Weintrob A, Ganesan A, et al. Clinical outcomes of elite controllers, viremic controllers, and long-term nonprogressors in the US Department of Defense HIV natural history study. J Infect Dis. 2009;200:1714–23.
Côrtes FH, Passaes CPB, Bello G, Teixeira SLM, Vorsatz C, Babic D, et al. HIV controllers with different viral load cutoff levels have distinct virologic and immunologic profiles. J Acquir Immune Defic Syndr. 2015;68:377–85.
Sharkey M, Babic DZ, Greenough T, Gulick R, Kuritzkes DR, Stevenson M. Episomal viral cDNAs identify a reservoir that fuels viral rebound after treatment interruption and that contributes to treatment failure. PLoS Pathog. 2011;7:e1001303.
Rodrigo AG, Goracke PC, Rowhanian K, Mullins JI. Quantitation of target molecules from polymerase chain reaction-based limiting dilution assays. AIDS Res Hum Retrovir. 1997;13:737–42.
Rose PP, Korber BT. Detecting hypermutations in viral sequences with an emphasis on G-to-A hypermutation. Bioinformatics. 2000;16:400–1.
Guindon S, Lethiec F, Duroux P, Gascuel O. PHYML online—a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005;33:W557–9.
Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.
Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55:539–52.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Lengauer T, Sander O, Sierra S, Thielen A, Kaiser R. Bioinformatics prediction of HIV coreceptor usage. Nat Biotechnol. 2007;25:1407–10.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package [Internet]. R Packag. version 2.4-2. https://cran.r-project.org/package=vegan (2017).
Gregori J, Perales C, Rodriguez-Frias F, Esteban JI, Quer J, Domingo E. Viral quasispecies complexity measures. Virology. 2016;493:227–37.
Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2016;2:vew007.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.
Ferreira MAR, Suchard MA. Bayesian analysis of elapsed times in continuous-time Markov chains. Can J Stat. 2008;36:355–68.
Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.
Pinheiro RS, Viacava F, Travassos C, Brito AS. Gênero, morbidade, acesso e utilização de serviços de saúde no Brasil. Cien Saude Colet. 2002;7:687–707.
Grinsztejn B, Smeaton L, Barnett R, Klingman K, Hakim J, Flanigan T, et al. Sex-associated differences in pre-antiretroviral therapy plasma HIV-1 RNA in diverse areas of the world vary by CD4+ T-cell count. Antivir. Ther. 2011;16:1057–62.
Ruel TD, Zanoni BC, Ssewanyana I, Cao H, Havlir DV, Kamya M, et al. Sex differences in HIV RNA level and CD4 cell percentage during childhood. Clin Infect Dis. 2011;53:592–9.
Sterling TR, Lyles CM, Vlahov D, Astemborski J, Margolick JB, Quinn TC. Sex differences in longitudinal human immunodeficiency virus type 1 RNA levels among seroconverters. J Infect Dis. 1999;180:666–72.
Pereyra F, Palmer S, Miura T, Block BL, Wiegand A, Rothchild AC, et al. Persistent low-level viremia in HIV-1 elite controllers and relationship to immunologic parameters. J Infect Dis. 2009;200:984–90.
Avila-Rios S, Sued O, Rhee S-Y, Shafer RW, Reyes-Teran G, Ravasi G. Surveillance of HIV transmitted drug resistance in Latin America and the Caribbean: a systematic review and meta-analysis. PLoS ONE. 2016;11:e0158560.
Pimentel VF, Morgado MG, Bello G, Guimarães MDC, Castilho EA, Veloso VG, et al. Temporal trends and molecular epidemiology of HIV type 1 infection in Rio de Janeiro, Brazil. AIDS Res Hum Retrovir. 2013;29:1553–61.
Pernas M, Casado C, Sandonis V, Arcones C, Rodriguez C, Ruiz-Mateos E, et al. Prevalence of HIV-1 dual infection in long-term nonprogressor-elite controllers. J Acquir Immune Defic Syndr. 2013;64:225–31.
Wang B, Spira TJ, Owen S, Lal RB, Saksena NK. HIV-1 strains from a cohort of American subjects reveal the presence of a V2 region extension unique to slow progressors and non-progressors. AIDS. 2000;14:213–23.
Grenfell BT, Pybus OG, Gog JR, Wood JLN, Daly JM, Mumford JA, et al. Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004;303:327–32.
Sagar M, Kirkegaard E, Long EM, Celum C, Buchbinder S, Daar ES, et al. Human immunodeficiency virus type 1 (HIV-1) diversity at time of infection is not restricted to certain risk groups or specific HIV-1 subtypes. J Virol. 2004;78:7279–83.
Sagar M, Kirkegaard E, Lavreys L, Overbaugh J. Diversity in HIV-1 envelope V1–V3 sequences early in infection reflects sequence diversity throughout the HIV-1 genome but does not predict the extent of sequence diversity during chronic infection. AIDS Res Hum Retrovir. 2006;22:430–7.
Gottlieb GS, Heath L, Nickle DC, Wong KG, Leach SE, Jacobs B, et al. HIV-1 variation before seroconversion in men who have sex with men: analysis of acute/early HIV infection in the multicenter AIDS cohort study. J Infect Dis. 2008;197:1011–5.
Abrahams M-R, Anderson JA, Giorgi EE, Seoighe C, Mlisana K, Ping L-H, et al. Quantitating the multiplicity of infection with human immunodeficiency virus type 1 subtype C reveals a non-poisson distribution of transmitted variants. J Virol. 2009;83:3556–67.
Kearney M, Maldarelli F, Shao W, Margolick JB, Daar ES, Mellors JW, et al. Human immunodeficiency virus type 1 population genetics and adaptation in newly infected individuals. J Virol. 2009;83:2715–27.
Bar KJ, Li H, Chamberland A, Tremblay C, Routy JP, Grayson T, et al. Wide variation in the multiplicity of HIV-1 infection among injection drug users. J Virol. 2010;84:6241–7.
Li H, Bar KJ, Wang S, Decker JM, Chen Y, Sun C, et al. High multiplicity infection by HIV-1 in men who have sex with men. PLoS Pathog. 2010;6:e1000890.
Novitsky V, Wang R, Margolin L, Baca J, Rossenkhan R, Moyo S, et al. Transmission of single and multiple viral variants in primary HIV-1 subtype C infection. PLoS ONE. 2011;6:e16714.
SSDA conducted most of the experiments and analyzed the data. DGC participated in DNA extraction and single genome sequencing. FHC participated in sample processing and determination of CD4+ T-cell counts. SLMT and KSS conducted the CCR5 genotyping. BH, BG and VGV conducted the patient recruitment and follow-up. MGM participated in study design and supervised the project. GB conceived and designed the study and supervised the experiments. SSDA and GB wrote the first draft. All authors read and approved the final manuscript.
We thank Dr. Edson Delatorre for assisting in the Vegan analysis, Dr. Vera Bongertz for review of the manuscript and the Plataforma Genômica de Sequenciamento de DNA-RPT01A-PDTIS/FIOCRUZ for the nucleotide sequencing. We also thank all patients whose participation and commitment made this work possible.
The authors declare that they have no competing interests.
Availability of data and materials
HIV-1 sequences generated during the current study were deposited in Genbank under the Accession Numbers KY852518—KY852939.
Ethics approval and consent to participate
The present work was approved by the Brazilian National Human Research Ethics Committee (CONEP 14430/2011) and all subjects gave written informed consent.
SSDA is supported by funding from the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). This work was supported by Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro—FAPERJ (Grant Number E-26/110.123/2014).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.