Reduction of HIV-1 Reservoir Size and Diversity After 1 Year of cART Among Brazilian Individuals Starting Treatment During Early Stages of Acute Infection

The aim of early combined antiretroviral therapy (cART) of HIV is to limit the seeding of the viral reservoir during the initial phase of infection and, consequently, decrease intrahost viral diversity. Here, we assessed the effect of early cART on size and complexity of the proviral reservoir. Peripheral blood mononuclear cell (PBMC) and plasma samples were obtained from ten HIV-infected Brazilian individuals diagnosed at the acute phase of infection, before (PRE_ART) and 12 months (M12_ART) after suppressive cART. HIV proviral reservoir size was determined by quantitative real-time PCR; intrahost viral diversity of the env C2-V3 region was assessed by single genome amplification or next-generation sequencing in PBMC and plasma, respectively. Mean nucleotide diversity (π) and normalized Shannon entropy (H_SN) were used to infer the complexity of the viral population. Compared to PRE_ART, M12_ART saw an immunological recovery with a gain of ∼200 CD4⁺ T cells (P = 0.008) and a normalization of the CD4/CD8 ratio [1.0 (IQR: 0.88–1.18), P = 0.016], as well as a significant decrease in HIV-1 RNA (∼4 log, P = 0.004) and DNA (∼1 log, P = 0.002) levels. The median time to achieve viral suppression was 3 months (IQR: 2.8–5.8 months). The high intermixing between sequences from both visits suggests that the HIV-1 DNA reservoir remained remarkably stable under cART. After 1 year of cART, there was a minor reduction in proviral π (Pre_ART = 0.20 vs. M12_ART = 0.10; P = 0.156) but a significant decrease in H_SN (Pre_ART = 0.41 vs. M12_ART = 0.25; P = 0.019). We found no correlation between π or H_SN at Pre_ART and the rate of HIV DNA decay, T CD4⁺ counts, or CD4/CD8 ratio at M12_ART. Based on a small cohort of Brazilian infected individuals under early cART and analyses of the env region, 1 year of follow-up suggested a reservoir size reduction, allowed a significant decrease of HIV-1 complexity, and achieved immunological restoration regardless of the initial HIV-1 plasma viral load, CD4⁺ T cell counts, or HIV-1 subtype. However, further studies in the Brazilian setting aiming a longer follow-up and larger cohort are required in this field.

Introduction

Combined antiretroviral therapy (cART) suppresses HIV-1 replication and reduces morbidity and mortality, but does not eradicate HIV-1 infection, as a low but persistent level of HIV-1 can still be detected in plasma and cell reservoirs (Chun et al., 1997; Kiselinova et al., 2015; Ghosn et al., 2018). Latent infected resting memory CD4⁺ T lymphocytes are the best known HIV-1 reservoir, which is established already during early infection. This reservoir includes cells with an integrated copy of the HIV-1 genome that is not expressed while the cells remain in a resting state (Chun et al., 2002; Douek et al., 2002) and is maintained mainly by the cells’ clonal expansion (Chomont et al., 2009, 2011; von Stockenstrom et al., 2015). However, it has been suggested that persistent virus replication may be an important contributor to its maintenance (Buzón et al., 2010; Yukl et al., 2010), particularly in lymphoid tissue sanctuary sites (Lorenzo-Redondo et al., 2016). Elimination/reduction of this latent reservoir represents a great hope for curing HIV-1 infection. As the reservoir is established during the acute phase of infection, early cART has been proposed as a means to restrict reservoir size and genetic complexity (Lori et al., 1999; Strain et al., 2005; Chomont et al., 2009; Josefsson et al., 2013).

HIV-1 sequence diversity is limited by the “genetic bottleneck effect” during sexual transmission, which selects viruses with the highest overall fitness (Carlson et al., 2014) but results in a more homogeneous viral population. However, in cohorts of men who have sex with men (MSM), the selection for fitter variants appears less stringent, resulting in infection being established by multiple founding viruses (Gottlieb et al., 2008; Li et al., 2010). Accordingly, the complexity of sequence diversity during early HIV-1 infection may affect the efficacy of cART in decreasing the reservoir, mainly due to immune and therapeutic escape mutations. Another potential escape route is through viral recombinants (Abrahams et al., 2009; Kearney et al., 2009; Bar et al., 2010; Li et al., 2010; Batorsky et al., 2011; Novitsky et al., 2011). Therefore, investigating the diversity of viral variants during primary infection may help evaluate viral evolution and predict clinical outcomes.

Residual replication has been proposed as a mechanism that maintains the HIV-1 reservoir during cART. Failure to block the viral replication cycle, enables renewed cellular infections and continuous replenishment of the HIV-1 DNA reservoir (Josefsson et al., 2013). Individuals who initiated cART during early infection, may present low viral diversity in peripheral blood mononuclear cell (PBMC) and plasma reservoirs after years on suppressive therapy, and the slight sign of viral replication indicates that the reservoir was maintained by homeostatic cell proliferation (Josefsson et al., 2013). A recent study demonstrated a faster decline in HIV DNA levels in early cART-treated patients with homogeneous viral populations after 6 months of therapy (Wang et al., 2017).

Understanding the dynamics of long-lived cellular HIV-1 reservoirs in individuals treated during primary infection, can direct the choice of long-term treatment regimens to achieve post-treatment control of HIV-1 in Brazil. Here, we investigated the effect of early cART during the initial Fiebig stages of HIV-1 acute infection on the size, diversity, and complexity of HIV-1 total DNA in the PBMC reservoir by examining HIV-1 populations prior and 1 year after cART initiation.

Materials and Methods

Subject Characteristics

A cohort of Brazilian individuals presenting recent HIV-1 infection has been followed-up since August 2013 at the Instituto Nacional de Infectologia Evandro Chagas (INI-FIOCRUZ). These individuals initiated cART immediately after diagnosis and have been described elsewhere (Ferreira et al., 2017). For the present study, we selected only individuals, who started cART during the acute phase (Fiebig II–V) (Fiebig et al., 2003) of HIV-1 infection (n = 10). Participants were recruited between December 2014 and October 2015, and had at least 1 year of successful cART after that. PBMC and plasma samples were obtained at the baseline visit (PRE_ART) and 12 months after cART onset (M12_ART), and were stored until use. The processing of all HIV samples was performed in accordance with institutional standard biosecurity and safety procedures at biosafety level 2. The study was approved by the INI Ethical Review Board (approval number 36859614.8.0000.5262), and all subjects gave written informed consent in accordance with the Declaration of Helsinki.

CD4⁺ and CD8⁺ T Cell Counts and HIV-1 RNA Quantification

Peripheral blood CD4⁺ and CD8⁺ T cell counts were determined by flow cytometry using the MultiTest TruCount-Kit and MultiSet software on a FACSCalibur flow cytometer (BD Biosciences, United States). HIV-1 RNA in plasma was measured by the Abbot Real-Time HIV-1 Assay, whose lower limit of detection was 40 copies/mL (Abbott Laboratories, Germany).

HIV-1 Total DNA Measurement in PBMCs

Total cellular DNA was extracted from cryopreserved PBMCs (1 × 10⁷ cells) obtained at PRE_ART and M12_ART using the QIAamp DNA Mini Kit (Qiagen, Germany). Cell-associated HIV-1 DNA was quantified using the Generic HIV^® DNA Cell Kit (Biocentric, France), following the manufacturer’s instructions. The assay’s lower limit of detection was 40 HIV DNA copies/10⁶ cells.

HIV-1 DNA Single Genome Amplification (SGA)

Proviral DNA was extracted from PBMCs using the QIAamp DNA Blood Mini Kit (Qiagen, United States) according to the manufacturer’s instructions. HIV-1 quasispecies was obtained by SGA of a 552-bp fragment from the C2-V3 env region through nested PCR using Platinum Taq DNA polymerase (Invitrogen, United States) as described elsewhere (Delwart et al., 1993). Considering a Poisson distribution, at a dilution in which approximately 30% of amplicons are positive, a single amplifiable molecule is present about 80% of the time (Palmer et al., 2005). The PCR products were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, United Kingdom). Sequences were obtained using the ABI BigDye Terminator v.3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems, United States) on an ABI 3130 Genetic Analyzer (Applied Biosystems). Sequences were assembled and edited using SeqMan 7.0 software (DNASTAR Inc., United States). APOBEC3G/F-mediated hypermutations were revealed by Hypermut software (Rose and Korber, 2000) and sequences showing ambiguous bases were excluded.

HIV-1 env RNA Haplotypes Reconstruction From NGS Data

Viral RNA from plasma samples collected at PRE_ART (baseline) was extracted using the QIAamp Viral RNA Mini Kit (Qiagen, Germany). The cDNA was obtained by reverse-transcribed PCR using the SuperScript^TM III Reverse Transcriptase (Invitrogen, United States) and was then subjected to nested PCR for amplification of the env gene as described above. The resulting amplicons were made into a library using the Nextera^® XT DNA Library Prep Kit with unique barcodes from the Nextera^® XT Index Kit (Illumina, United States), following the manufacturer’s instructions. DNA sequencing was performed on a MiSeq instrument using MiSeq Reagent Nano Kit, v2 (500 cycles; Illumina, United States). Demultiplexed reads were trimmed to remove adaptors, low-quality bases (Q <25), and short reads (<100 bp), and then mapped against single-genome amplification consensus sequences from each patient using Geneious software v.9.1.8 (Kearse et al., 2012) with high mapping quality (MAQ ≥30). Alignment regions with at least 500× coverage were used for haplotype reconstruction with QuasiRecomb 1.2 (Marz et al., 2014), employing the flag “-conservative” to increase specificity. Only haplotypes with frequencies ≥1% were used for further analysis.

HIV-1 Subtyping

Sequences were aligned with HIV-1 reference sequences from the Los Alamos database¹ using ClustalW in MEGA 6 and were manually edited. The final env alignment covered positions 6,840–7,372 relative to the HXB2 genome. Maximum-likelihood (ML) phylogenetic trees were reconstructed with PhyML 3.0 (Guindon et al., 2010) using the most appropriate nucleotide substitution model selected with jModeltest v. 3.7 (Darriba et al., 2012). The approximate likelihood-ratio test (aLRT) was used to estimate the confidence of branching on the tree.

Analyses of Viral Diversity

Complexity of the intrahost viral population was assessed through two diversity measures. The mean nucleotide diversity (π), an abundance-based functional index representing the average number of nucleotide differences between any two representatives of the population, was calculated in MEGA7 (Kumar et al., 2016). The normalized Shannon entropy (H_SN), an abundance-based index that measures viral population diversity based on haplotype frequencies, was calculated in the R package vegan (Oksanen et al., 2017), after sample rarefaction to correct for bias in sample size (Gregori et al., 2016).

Phenotypic Prediction of Co-receptor Usage

The quasispecies viral tropism was predicted based on the V3 amino acid sequence through the Geno2pheno algorithm (available at http://coreceptor.geno2pheno.org). The false positive rate cut-off was 10% for DNA, whereas RNA sequences were classified as CCR5 or non-CCR5-using viruses as described previously (Lengauer et al., 2007; Hayashida et al., 2017).

Statistical Analyses

All statistical analyses were performed in GraphPad Prism v6 (GraphPad Software, United States). Variables between groups (unpaired) or from the same group from different visits (paired) were compared using the Mann–Whitney U-test or Wilcoxon test, respectively. Association between variables was evaluated using the Spearman’s rank correlation. Linear regression was used to calculate the rate of HIV-1 DNA variation between visits in each subject. P-values <0.05 were considered statistically significant.

Sequence Availability

SGA env sequences have been submitted to GenBank under accession numbers MH765045–MH765329. NGS data have been deposited to the NCBI BioProject database under accession number PRJNA487221.

Results

Clinical, Epidemiological, and Immunological Characteristics of Early-Treated Individuals

The study included ten HIV-1-infected Brazilian individuals diagnosed during the acute phase of infection, with a median age at diagnosis of 28 years [interquartile range (IQR): 26–42 years]. All participants were MSM and started cART immediately after HIV-1 diagnosis. The cART regimen varied among participants, however, all included two co-formulated nucleoside reverse transcriptase inhibitors (lamivudine plus tenofovir) in combination with a non-nucleoside reverse transcriptase inhibitor (efavirenz, n = 5), a protease inhibitor (fosamprenavir/atazanavir, n = 2 each), or an integrase inhibitor (raltegravir, n = 1).

Prior to cART initiation (Pre_ART visit), median CD4⁺ and CD8⁺ T cell counts of all individuals were 634 and 1473 cells/mm³ (IQR: 420–886 for CD4⁺ and 526–1840 for CD8⁺ cells), respectively, and the CD4/CD8 ratio was 0.52 (IQR: 0.42–0.83). After 12 months on cART (M12_ART), we observed a significant immunological recovery, with median CD4⁺ counts of 836 cells/mm³ (IQR: 719–1122) (Figures 1A, 2). This meant an average increase of ∼200 CD4⁺ T cells (P = 0.008) and a significant rise in the CD4/CD8 ratio [median of 1.0 (IQR: 0.88–1.18), P = 0.016] (Figures 1B, 2).

FIGURE 1

Figure 1. Immunological and virological measurements before and after cART initiation. CD4⁺ T cell counts (A), CD4⁺/CD8⁺ ratios (B), HIV-1 proviral load in PBMCs (C), and HIV-1 viral load in plasma (D) were measured at Pre_ART and M12_ART visits (pink circles and blue squares, respectively). P-values <0.05 were considered statistically significant.

FIGURE 2

Figure 2. Clinical follow-up of the ten early-treated individuals. T cell counts (cells/μL; CD4⁺, blue circles; CD8⁺, red squares) and plasma RNA viral loads (log copies/mL, green triangles) values over time after cART onset (months) are shown on the left and right Y-axis, respectively. The individuals’ identifications are indicated at the top of each graph.

The immunological improvement observed 12 months after cART initiation was accompanied by a significant decrease in the median levels of HIV-1 total DNA in PBMCs and HIV-1 RNA in plasma (Figures 1C,D). The plasma HIV-1 viral load decreased drastically (∼4 log) between PRE_ART and M12_ART visits (from 5.86 to <1.6 log copies/mL, P = 0.004), whereas HIV-1 DNA in the PBMC compartment displayed a modest yet significant reduction (∼1 log, from 3.04 to 2.18 log copies/10⁶ PBMC, P = 0.002). The median time to achieve viral suppression (<1.6 log HIV RNA copies/mL) after the onset of cART was approximately 3 months (IQR: 2.8–5.8 months), and only one individual (IVA41) exhibited a minor viral rebound (2.66 log HIV RNA copies/mL) within the 12 months of follow-up (Figure 2). It should be noted, however, that sustained viral loads were recovered after that (data not shown).

Sequence Analysis

To understand the effect of early cART initiation on HIV-reservoir diversity and complexity, we analyzed the intrahost viral population by SGA of the HIV-1 env gene before and after cART. Similar quantities of SGA HIV-1 proviral sequences were obtained from both visits’ samples: 15 (IQR: 15–16) in Pre_ART and 14.5 (IQR: 11–16) in M12_ART (P = 0.14). Additionally, we used NGS to evaluate HIV diversity in the individuals’ plasma samples at the Pre_ART visit. The median coverage per sample was 8,182 (IQR: 4,755–10,904) reads per base (Supplementary Figure S1), and after assembly between one to eight HIV-1 env haplotypes were reconstructed per sample.

All sequences obtained from proviral DNA and plasma RNA branched together in highly supported monophyletic clusters (aLRT = 1) by subject (Figure 3A); suggesting that in all subjects infection resulted from a single or limited number of closely related viral variants. It is noteworthy that in four subjects (IVA32, IVA21, IVA17, and IVA37), the HIV-1 RNA sequences were the closest to the common ancestor of the subjects’ HIV-1-infecting lineage and outside the main cluster comprising the proviral variants. In the remaining subjects, the HIV-1 RNA sequences were intermixed with proviral DNA sequences. HIV-1 subtype B was the most frequent variant detected (n = 6), followed by subtypes F1 and C (n = 2 in each). R5-tropic viruses dominated the intrahost HIV-1 population in plasma and PBMCs at both visits, with only two individuals (IVA21 and IVA40) presenting a low frequency (6%) of X4-tropic viral clones in PBMCs at Pre_ART (Table 1 and Figure 3A). APOBEC3G/F-induced hypermutations were identified at Pre_ART in proviral sequences from only three individuals (IVA09, IVA36, and IVA40), with frequencies of 6–12%.

FIGURE 3

Figure 3. Impact of early cART initiation on intrahost HIV-1 population diversity. (A) ML phylogenetic tree of the env sequences from Pre_ART and M12_ART visits. Tips’ shapes represent the viral compartment (proviral DNA, circles; plasma RNA, triangles) and are color-coded according to the visit (Pre_ART, pink; M12_ART, purple). The branches’ colors agree with the subtype assignment as indicated in the legend. Clusters from each individual are indicated by shaded gray boxes. Branch supports (aLRT-SH) are indicated at key nodes. Tips shapes marked with an “H” indicate the presence of APOBEC3G-mediated G to A hypermutations and X4 labels highlight X4-tropic sequences. Horizontal branch lengths are proportional to the bar at the bottom indicating nucleotide substitutions per site. (B) Mean nucleotide diversity (π) and normalized Shannon entropy (H_SN) indices were calculated from the proviral env sequences obtained at Pre_ART and M12_ART visits. Thick and thin lines represent the median and interquartile ranges, respectively. P-values <0.05 were considered statistically significant.

TABLE 1

Table 1. Virological characteristics of HIV-1 early-infected individuals.

To address the influence of early cART on reservoir complexity, we measured the average pairwise nucleotide differences of the viral population (π) and the uniformity of the haplotype distribution (H_SN) of viral DNA and RNA from each visit. Hypermutated sequences were excluded from these analyses. All individuals exhibited low π values (<1.5%) despite the Fiebig stage at cART initiation. We observed a minor decline in median proviral π after 1 year of cART (Pre_ART = 0.20, IQR: 0.10–0.58 vs. M12_ART = 0.10, IQR: 0.20–0.33; P = 0.156) (Figure 3B). There was, however, a strong positive correlation between proviral and plasma π values before cART onset (r = 0.77; P = 0.0157, Supplementary Figure S2), which was not observed at M12_ART nor when comparing π values from HIV DNA between visits. The overall median H_SN estimated for the HIV-1 proviral population decreased significantly between Pre_ART and M12_ART visits (0.41, IQR: 0.28–0.56 vs. 0.25, IQR: 0.14–0.30; P = 0.019) (Figure 3B). However, no correlation was found when comparing the H_SN values calculated from different compartments or visits.

In spite of relatively low π values at Pre_ART, the values varied among participants (range: 0.1–1.5%). A comparison revealed no significant correlation between π or H_SN at Pre_ART with immunological and virological parameters, such as rate of HIV DNA decay, T CD4⁺ cell change, or the CD4/CD8 ratio, at M12_ART (Figure 4).

FIGURE 4

Figure 4. Correlations between proviral HIV-1 diversity indices and immunological and virological measurements. Mean nucleotide diversity (π) of the proviral population at Pre_ART was compared with HIV decay in PBMCs (A), T CD4⁺ cell change (B), and CD4⁺/CD8⁺ ratio in M12_ART (C). Normalized Shannon entropy (H_SN) of the proviral population at Pre_ART was compared with HIV decay in PBMCs (D), T CD4⁺ cell change (E), and CD4⁺/CD8⁺ ratio in M12_ART (F). P-values <0.05 were considered statistically significant.

Discussion

A latent HIV-1 reservoir in resting memory CD4+ T cells is recognized as the major barrier to HIV-1 eradication, as its establishment and long-term persistence enables renewed viremia after treatment failure or interruption (Richman et al., 2009). Early cART has been proposed as a means for achieving long-term control of viral replication upon cART interruption, by delaying the viral rebound and potentially inducing a post-treatment controller status (Archin et al., 2012; Avettand-Fènoel et al., 2016). Here, we analyzed ten acutely HIV-infected patients starting cART during the early Fiebig stages, before and 1 year after treatment, to investigate the effect of early therapy on the size and complexity of HIV-1 total DNA in the PBMC reservoir.

Benefits of early cART include a reduction in residual viral replication, a restraint of viral diversity and reservoir development, and accelerated immune restoration (Yerly et al., 2000; Ngo-Giang-Huong et al., 2001; Delwart et al., 2002; Hecht et al., 2006). We observed significant immunological restoration in all subjects 1 year after cART onset during acute HIV infection, with increases of both CD4⁺ T cell counts and the CD4/CD8 ratio (∼1.0), in agreement with previous results (Hoenigl et al., 2016). The rise of the CD4/CD8 ratio was previously linked with reduced levels of HIV DNA in peripheral blood cells (Chun et al., 2002), mainly if cART was established during primary HIV-1 infection (Hocqueloux et al., 2013). Thus, initiating cART during the acute phase of HIV infection offers an opportunity to reduce HIV reservoirs and achieve optimal immune reconstitution. Nevertheless, a recent study (Colby et al., 2018) has demonstrated that a CD4/CD8 ratio <1 is associated with a rapid viral rebound following treatment interruption, in spite of early cART initiation and sustained undetectable viral load along treatment follow-up. The immunological recovery was accompanied here by a reduction in both HIV RNA and DNA levels. The median viral load decay time to undetectable levels was about 3 months after cART start, similar to previous studies (Pilcher et al., 2004; Hoenigl et al., 2016), indicating that a rapid and efficient virologic suppression was achieved in this cohort. The median HIV-1 DNA load of our subjects at baseline visit was ∼3 log copies/10⁶ PBMCs, similar to the one described in the ANRS PRIMO cohort (Schiffer et al., 2004). After 1 year of suppressive cART, HIV-1 DNA levels in PBMCs were reduced to ∼2 log copies/10⁶ PBMCs, reaching similar levels as those of HIV-1-infected patients, who had started intensive and standard cART during primary HIV-1 infection (Hocqueloux et al., 2013; Chéret et al., 2015).

A reduction of the HIV-1 latent reservoir may help HIV-1-infected patients reach at least a transient drug-free remission of their disease (Hocqueloux et al., 2013; Sáez-Cirión et al., 2013). Lower levels of HIV DNA in PBMCs at antiretroviral treatment interruption have been associated with a longer time to treatment resumption (Piketty et al., 2010). Despite the reduction in HIV-1 total DNA in the PBMC reservoir observed in our study, the values were still higher (median 2.18 log copies/10⁶ PBMCs) than in post-treatment controllers (PTC) from the VISCONTI cohort (median 1.71 log copies/10⁶ PBMCs), for whom long-term viral control was maintained after treatment interruption (Sáez-Cirión et al., 2013). The PTC’s median time to plasma viral load becoming undetectable was similar to the one calculated in this study; however, median cART duration was much longer (36.5 months in PTC vs. 12 months in this study). Moreover, a study conducted in Thailand showed that total HIV DNA mean values in PBMCs from early-treated acutely HIV-infected individuals decreased over time achieving very low levels after 3 years of cART (Ananworanich et al., 2016). Thus, with longer follow-up of the patients included in this study, total HIV DNA may reduce even further in the PBMC reservoir.

Here, analysis of HIV env sequences from plasma and PBMCs revealed very high population intermixing between the clones obtained before and after cART onset, indicating that the HIV-1 DNA reservoir remained remarkably stable. This result agrees with the findings of Josefsson et al. (2013), whereby viral populations in both pre-therapy plasma and cells isolated after suppressive cART initiated during early/acute infection were nearly monomorphic. The absence of genetic changes indicates that maintenance of the HIV reservoir during suppressive cART, as observed in this study, was the result of homeostatic cell proliferation rather than of ongoing viral replication. This is consistent with the results of other studies analyzing HIV reservoir sequences but from longer time intervals than those reported here (Chomont et al., 2009; Kearney et al., 2014; von Stockenstrom et al., 2015; Brodin et al., 2016; Van Zyl et al., 2017). However, the absence of divergence in the PBMC HIV DNA reservoir does not exclude the possibility that ongoing residual replication could be occurring in sanctuary sites, especially in lymphoid tissue (Lorenzo-Redondo et al., 2016). A more robust sampling scheme, assessing distinct compartments and over a longer time frame should be performed to increase the chances of identifying putatively evolving viruses.

Recent studies have demonstrated that patients, who were treated early during primary HIV-1 infection, had low viral diversity in both plasma or cells isolated after years on therapy (Josefsson et al., 2013; Kearney et al., 2014). Even though diversity values were low at the baseline visit, we detected a slight decline in HIV DNA diversity in PBMCs and a significant decline in proviral population complexity after 12 months on cART. As we found no evidence of ongoing viral replication, the reduction in intrahost proviral complexity following cART may be explained by the clearance of HIV-infected cells, probably long-lived T cells (Palmer et al., 2008), and coincided with a reduction in HIV DNA levels. Interestingly, the diversity indices exhibited by these subjects after 1 year of suppressive cART, were similar to those found in a subgroup of rare individuals capable of naturally controlling HIV-1 replication and maintaining it at low levels (named elite controllers). The diversity indices exhibited by elite controllers are associated with the presence of a putative more efficient mechanism to control HIV-1 replication and disease progression (de Azevedo et al., 2017). A previous study found an association between HIV genetic diversity during the early phase of infection and a faster HIV DNA decline following cART (Wang et al., 2017). In contrast, we did not find any significant correlation between immunological recovery or virologic control achieved after 1 year of suppressive cART and HIV diversity before cART initiation.

We detected three HIV-1 env subtypes (B, C, and F1) among the ten patients included in this study. This subtype distribution is in accordance with the molecular epidemiology found in Brazil (Bello et al., 2011). The different HIV-1 clades seem not to have influenced cART outcomes observed in the different subjects. However, it is noteworthy that IVA05 and IVA38, both infected with HIV-1 subtype C strains, presented the highest values of viral diversity and complexity in both PBMC and plasma compartments. This finding agrees with a previous study showing that subtype C displays higher median diversity than subtype B (Rieder et al., 2011). Additionally, the majority of HIV-1 variants were R5-tropic during the entire follow-up period. The only exceptions were two subjects infected with HIV-1 subtype B, who exhibited a low frequency of X4-tropic clones even at the very beginning of HIV-1 infection, in agreement with previous studies (Sheppard et al., 2002; Raymond et al., 2010; Chalmet et al., 2012).

Despite the intensification of offering rapid tests for HIV detection inside and outside the health services, covering key populations at risk, and the ongoing follow-up of cohorts of serodiscordant couples conducted by our institution (INI-FIOCRUZ), the identification of individuals in the early primary infection phases, according to the Fiebig classification, and their follow-up are still hard tasks, as well as, the obtaining of high quality cryopreserved PMBC samples from long term storage to perform the laboratory analyses. Consequently, most of the studies addressing this group, considering the complex strategies for characterizing the early stages of HIV-acute infection, management of early cART starting and the clinical follow-up at specific time points limits the scope, mainly in developing countries where resources to closely following clinical cohorts are more restricted. Thus, the effects of early treatment initiation on the viral evolution, size and distribution of the viral reservoir are still poorly explored and need to be better characterized and, although having a relatively small number of individuals included in our study, we still justified conducting the present analyzes. As of our knowledge, this is the first Brazilian study discussing the viral evolution and reservoir size overtime in a cohort of HIV infected individuals starting cART in the early Fiebig stages of HIV primary infection. Further studies with larger cohorts and longer follow-up periods should be performed to identify better predictive markers for individuals, who might reach the PTC status after cART interruption. Assessment of HIV DNA level and diversity in different T cell subsets and/or compartments would substantially improve our understanding of the effect early cART has in these subjects and would likely add to the growing body of evidence indicating the contribution of each cell population to maintaining the HIV reservoir (Buzon et al., 2014).

In summary, we demonstrate that based on this small cohort and analyses of only the envelope region, 1 year of suppressive cART initiated in the early stages of HIV infection suggested a reduction of the size and complexity of HIV-1 total DNA in the PBMC reservoir, as well as achieve immunological restoration. Immunological recovery and virologic suppression were not associated with proviral diversity before cART onset. The early initiation of cART in HIV-acutely infected individuals in current Brazilian setting may favor strategies to achieve post-treatment control of HIV and, ultimately, a functional cure by restricting the pool of variants and allowing a more focused targeting with therapeutic vaccines or other immune approaches. However, further studies are required to determine the follow-up time necessary to achieve the lowest viral endpoints for clinical management.

Author Contributions

MG, MM, VV, FC, and BG conceived and designed the study. TL, ED, and FC performed the experiments. AF, SC, and MA participated in patient recruitment. TL, ED, and MG analyzed the data and drafted the manuscript. All authors reviewed and approved the final manuscript.

Funding

ED was financed by a Postdoctoral fellowship from the “Programa Nacional de Pós-Doutorado (PNPD)” by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES) – Finance Code 001. MG and MM were recipients of a CNPq Fellowship, and TL was funded by a CAPES Ph.D. fellowship.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Soraia Moura and Sandro Nazer from Laboratório de Pesquisa Clínica em DST e AIDS for their assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00145/full#supplementary-material

Footnotes

1. ^ https://www.hiv.lanl.gov/

References

Abrahams, M.-R., Anderson, J. A., Giorgi, E. E., Seoighe, C., Mlisana, K., Ping, L.-H., et al. (2009). Quantitating the multiplicity of infection with human immunodeficiency virus type 1 subtype C reveals a non-poisson distribution of transmitted variants. J. Virol. 83, 3556–3567. doi: 10.1128/JVI.02132-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Ananworanich, J., Chomont, N., Eller, L. A., Kroon, E., Tovanabutra, S., Bose, M., et al. (2016). HIV DNA set point is rapidly established in acute HIV infection and dramatically reduced by early ART. EBioMedicine 11, 68–72. doi: 10.1016/j.ebiom.2016.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Archin, N. M., Vaidya, N. K., Kuruc, J. D., Liberty, A. L., Wiegand, A., Kearney, M. F., et al. (2012). Immediate antiviral therapy appears to restrict resting CD4+ cell HIV-1 infection without accelerating the decay of latent infection. Proc. Natl. Acad. Sci. U.S.A. 109, 9523–9528. doi: 10.1073/pnas.1120248109

PubMed Abstract | CrossRef Full Text | Google Scholar

Avettand-Fènoel, V., Hocqueloux, L., Ghosn, J., Cheret, A., Frange, P., Melard, A., et al. (2016). Total HIV-1 DNA, a marker of viral reservoir dynamics with clinical implications. Clin. Microbiol. Rev. 29, 859–880. doi: 10.1128/CMR.00015-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Bar, K. J., Li, H., Chamberland, A., Tremblay, C., Routy, J. P., Grayson, T., et al. (2010). Wide variation in the multiplicity of HIV-1 infection among injection drug users. J. Virol. 84, 6241–6247. doi: 10.1128/JVI.00077-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Batorsky, R., Kearney, M. F., Palmer, S. E., Maldarelli, F., Rouzine, I. M., and Coffin, J. M. (2011). Estimate of effective recombination rate and average selection coefficient for HIV in chronic infection. Proc. Natl. Acad. Sci. U.S.A. 108, 5661–5666. doi: 10.1073/pnas.1102036108

PubMed Abstract | CrossRef Full Text | Google Scholar

Bello, G., Soares, M. A., and Schrago, C. G. (2011). The Use of bioinformatics for studying HIV evolutionary and epidemiological history in South America. AIDS Res. Treat. 2011, 1–13. doi: 10.1155/2011/154945

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodin, J., Zanini, F., Thebo, L., Lanz, C., Bratt, G., Neher, R. A., et al. (2016). Establishment and stability of the latent HIV-1 DNA reservoir. eLife 5:e18889. doi: 10.7554/eLife.18889

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzon, M. J., Martin-Gayo, E., Pereyra, F., Ouyang, Z., Sun, H., Li, J. Z., et al. (2014). Long-term antiretroviral treatment initiated at primary HIV-1 infection affects the size, composition, and decay kinetics of the reservoir of HIV-1-infected CD4 T cells. J. Virol. 88, 10056–10065. doi: 10.1128/JVI.01046-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzón, M. J., Massanella, M., Llibre, J. M., Esteve, A., Dahl, V., Puertas, M. C., et al. (2010). HIV-1 replication and immune dynamics are affected by raltegravir intensification of HAART-suppressed subjects. Nat. Med. 16, 460–465. doi: 10.1038/nm.2111

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlson, J. M., Schaefer, M., Monaco, D. C., Batorsky, R., Claiborne, D. T., Prince, J., et al. (2014). Selection bias at the heterosexual HIV-1 transmission bottleneck. Science 345:1254031. doi: 10.1126/science.1254031

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalmet, K., Dauwe, K., Foquet, L., Baatz, F., Seguin-Devaux, C., Van Der Gucht, B., et al. (2012). Presence of CXCR4-using HIV-1 in patients with recently diagnosed infection: correlates and evidence for transmission. J. Infect. Dis. 205, 174–184. doi: 10.1093/infdis/jir714

PubMed Abstract | CrossRef Full Text | Google Scholar

Chéret, A., Nembot, G., Mélard, A., Lascoux, C., Slama, L., Miailhes, P., et al. (2015). Intensive five-drug antiretroviral therapy regimen versus standard triple-drug therapy during primary HIV-1 infection (OPTIPRIM-ANRS 147): a randomised, open-label, phase 3 trial. Lancet Infect. Dis. 15, 387–396. doi: 10.1016/S1473-3099(15)70021-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Chomont, N., DaFonsecaa, S., Vandergeetena, C., Ancuta, P., and Sékaly, R. P. (2011). Maintenance of CD4+ T-cell memory and HIV persistence: keeping memory, keeping HIV. Curr. Opin. HIV AIDS 6, 30–36. doi: 10.1097/COH.0b013e3283413775

PubMed Abstract | CrossRef Full Text | Google Scholar

Chomont, N., El-Far, M., Ancuta, P., Trautmann, L., Procopio, F. A., Yassine-Diab, B., et al. (2009). HIV reservoir size and persistence are driven by T cell survival and homeostatic proliferation. Nat. Med. 15, 893–900. doi: 10.1038/nm.1972

PubMed Abstract | CrossRef Full Text | Google Scholar

Chun, T.-W., Justement, J. S., Pandya, P., Hallahan, C. W., McLaughlin, M., Liu, S., et al. (2002). Relationship between the size of the human immunodeficiency virus type 1 (HIV-1) reservoir in peripheral blood CD4 + T cells and CD4 +:CD8 + T cell ratios in aviremic HIV-1–infected individuals receiving long-term highly active antiretroviral therapy. J. Infect. Dis. 185, 1672–1676. doi: 10.1086/340521

PubMed Abstract | CrossRef Full Text | Google Scholar

Chun, T.-W., Stuyver, L., Mizell, S. B., Ehler, L. A., Mican, J. A. M., Baseler, M., et al. (1997). Presence of an inducible HIV-1 latent reservoir during highly active antiretroviral therapy. Proc. Natl. Acad. Sci. U.S.A. 94, 13193–13197. doi: 10.1073/pnas.94.24.13193

CrossRef Full Text | Google Scholar

Colby, D. J., Trautmann, L., Pinyakorn, S., Leyre, L., Pagliuzza, A., Kroon, E., et al. (2018). Rapid HIV RNA rebound after antiretroviral treatment interruption in persons durably suppressed in Fiebig i acute HIV infection brief-communication. Nat. Med. 24, 923–926. doi: 10.1038/s41591-018-0026-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Darriba, D., Taboada, G. L., Doallo, R., and Posada, D. (2012). jModelTest2: more models, new heuristics and parallel computing. Nat. Methods 9:772. doi: 10.1038/nmeth.2109

PubMed Abstract | CrossRef Full Text | Google Scholar

de Azevedo, S. S. D., Caetano, D. G., Côrtes, F. H., Teixeira, S. L. M., dos Santos Silva, K., Hoagland, B., et al. (2017). Highly divergent patterns of genetic diversity and evolution in proviral quasispecies from HIV controllers. Retrovirology 14:29. doi: 10.1186/s12977-017-0354-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Delwart, E., Magierowska, M., Royz, M., Foley, B., Peddada, L., Smith, R., et al. (2002). Homogeneous quasispecies in 16 out of 17 individuals during very early HIV-1 primary infection. AIDS 16, 189–195. doi: 10.1097/00002030-200201250-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Delwart, E. L., Shpaer, E. G., Louwagie, J., McCutchan, F. E., Grez, M., Rübsamen-Waigmann, H., et al. (1993). Genetic relationships determined by a DNA heteroduplex mobility assay: analysis of HIV-1 env genes. Science 262, 1257–1261. doi: 10.1126/science.8235655

PubMed Abstract | CrossRef Full Text | Google Scholar

Douek, D. C., Brenchley, J. M., Betts, M. R., Ambrozak, D. R., Hill, B. J., Okamoto, Y., et al. (2002). HIV preferentially infects HIV-specific CD4+ T cells. Nature 417, 95–98. doi: 10.1038/417095a

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, A. C. G., Coelho, L. E., Grinsztejn, E., Jesus, C. S. D., Guimarães, M. L., Veloso, V. G., et al. (2017). Transmitted drug resistance in patients with acute/recent HIV infection in Brazil. Braz. J. Infect. Dis. 21, 396–401. doi: 10.1016/j.bjid.2017.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Fiebig, E. W., Wright, D. J., Rawal, B. D., Garrett, P. E., Schumacher, R. T., Peddada, L., et al. (2003). Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS 17, 1871–1879. doi: 10.1097/01.aids.0000076308.76477.b8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosn, J., Taiwo, B., Seedat, S., Autran, B., and Katlama, C. (2018). HIV. Lancet 392, 685–697. doi: 10.1016/S0140-6736(18)31311-4

CrossRef Full Text | Google Scholar

Gottlieb, G. S., Heath, L., Nickle, D. C., Wong, K. G., Leach, S. E., Jacobs, B., et al. (2008). 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. 197, 1011–1015. doi: 10.1086/529206

PubMed Abstract | CrossRef Full Text | Google Scholar

Gregori, J., Perales, C., Rodriguez-Frias, F., Esteban, J. I., Quer, J., and Domingo, E. (2016). Viral quasispecies complexity measures. Virology 493, 227–237. doi: 10.1016/j.virol.2016.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashida, T., Tsuchiya, K., Kikuchi, Y., Oka, S., and Gatanaga, H. (2017). Emergence of CXCR4-tropic HIV-1 variants followed by rapid disease progression in hemophiliac slow progressors. PLoS One 12:e0177033. doi: 10.1371/journal.pone.0177033

PubMed Abstract | CrossRef Full Text | Google Scholar

Hecht, F. M., Wang, L., Collier, A., Little, S., Markowitz, M., Margolick, J., et al. (2006). A multicenter observational study of the potential benefits of initiating combination antiretroviral therapy during acute HIV infection. J. Infect. Dis. 194, 725–733. doi: 10.1086/506616

PubMed Abstract | CrossRef Full Text | Google Scholar

Hocqueloux, L., Avettand-Fènoël, V., Jacquot, S., Prazuck, T., Legac, E., Mélard, A., et al. (2013). Long-term antiretroviral therapy initiated during primary HIV-1 infection is key to achieving both low HIV reservoirs and normal T cell counts. J. Antimicrob. Chemother. 68, 1169–1178. doi: 10.1093/jac/dks533

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoenigl, M., Chaillon, A., and Little, S. J. (2016). CD4/CD8 cell ratio in acute HIV infection and the impact of early antiretroviral therapy: table 1. Clin. Infect. Dis. 63, 425–426. doi: 10.1093/cid/ciw293

PubMed Abstract | CrossRef Full Text | Google Scholar

Josefsson, L., von Stockenstrom, S., Faria, N. R., Sinclair, E., Bacchetti, P., Killian, M., et al. (2013). The HIV-1 reservoir in eight patients on long-term suppressive antiretroviral therapy is stable with few genetic changes over time. Proc. Natl. Acad. Sci. U.S.A. 110, E4987–E4996. doi: 10.1073/pnas.1308313110

PubMed Abstract | CrossRef Full Text | Google Scholar

Kearney, M., Maldarelli, F., Shao, W., Margolick, J. B., Daar, E. S., Mellors, J. W., et al. (2009). Human immunodeficiency virus type 1 population genetics and adaptation in newly infected individuals. J. Virol. 83, 2715–2727. doi: 10.1128/JVI.01960-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Kearney, M. F., Spindler, J., Shao, W., Yu, S., Anderson, E. M., O’Shea, A., et al. (2014). Lack of detectable HIV-1 molecular evolution during suppressive antiretroviral therapy. PLoS Pathog. 10:e1004010. doi: 10.1371/journal.ppat.1004010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiselinova, M., Geretti, A. M., Malatinkova, E., Vervisch, K., Beloukas, A., Messiaen, P., et al. (2015). HIV-1 RNA and HIV-1 DNA persistence during suppressive ART with PI-based or nevirapine-based regimens. J. Antimicrob. Chemother. 70, 3311–3316. doi: 10.1093/jac/dkv250

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Lengauer, T., Sander, O., Sierra, S., Thielen, A., and Kaiser, R. (2007). Bioinformatics prediction of HIV coreceptor usage. Nat. Biotechnol. 25, 1407–1410. doi: 10.1038/nbt1371

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Bar, K. J., Wang, S., Decker, J. M., Chen, Y., Sun, C., et al. (2010). High multiplicity infection by HIV-1 in men who have sex with men. PLoS Pathog. 6:e1000890. doi: 10.1371/journal.ppat.1000890

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenzo-Redondo, R., Fryer, H. R., Bedford, T., Kim, E. Y., Archer, J., Kosakovsky Pond, S. L., et al. (2016). Persistent HIV-1 replication maintains the tissue reservoir during therapy. Nature 530, 51–56. doi: 10.1038/nature16933

PubMed Abstract | CrossRef Full Text | Google Scholar

Lori, F., Jessen, H., Lieberman, J., Finzi, D., Rosenberg, E., Tinelli, C., et al. (1999). Treatment of human immunodeficiency virus infection with Hydroxyurea, Didanosine, and a Protease inhibitor before seroconversion is associated with normalized immune parameters and limited viral reservoir. J. Infect. Dis. 180, 1827–1832. doi: 10.1086/315113

PubMed Abstract | CrossRef Full Text | Google Scholar

Marz, M., Beerenwinkel, N., Drosten, C., Fricke, M., Frishman, D., Hofacker, I. L., et al. (2014). Challenges in RNA virus bioinformatics. Bioinformatics 30, 1793–1799. doi: 10.1093/bioinformatics/btu105

PubMed Abstract | CrossRef Full Text | Google Scholar

Ngo-Giang-Huong, N., Deveau, C., Da Silva, I., Pellegrin, I., Venet, A., Harzic, M., et al. (2001). Proviral HIV-1 DNA in subjects followed since primary HIV-1 infection who suppress plasma viral load after one year of highly active antiretroviral therapy. AIDS 15, 665–673. doi: 10.1097/00002030-200104130-00001

PubMed Abstract | CrossRef Full Text | Google Scholar

Novitsky, V., Wang, R., Margolin, L., Baca, J., Rossenkhan, R., Moyo, S., et al. (2011). Transmission of single and multiple viral variants in primary HIV-1 subtype C infection. PLoS One 6:e16714. doi: 10.1371/journal.pone.0016714

PubMed Abstract | CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2017). Vegan: Community Ecology Package. R Package Version. 2.0-10.

Google Scholar

Palmer, S., Kearney, M., Maldarelli, F., Halvas, E. K., Bixby, C. J., Bazmi, H., et al. (2005). Multiple, linked human immunodeficiency virus type 1 drug resistance mutations in treatment-experienced patients are missed by standard genotype analysis. J. Clin. Microbiol. 43, 406–413. doi: 10.1128/JCM.43.1.406-413.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmer, S., Maldarelli, F., Wiegand, A., Bernstein, B., Hanna, G. J., Brun, S. C., et al. (2008). Low-level viremia persists for at least 7 years in patients on suppressive antiretroviral therapy. Proc. Natl. Acad. Sci. U.S.A. 105, 3879–3884. doi: 10.1073/pnas.0800050105

PubMed Abstract | CrossRef Full Text | Google Scholar

Piketty, C., Weiss, L., Assoumou, L., Burgard, M., Mélard, A., Ragnaud, J. M., et al. (2010). A high HIV DNA level in PBMCs at antiretroviral treatment interruption predicts a shorter time to treatment resumption, independently of the CD4 nadir. J. Med. Virol. 82, 1819–1828. doi: 10.1002/jmv.21907

PubMed Abstract | CrossRef Full Text | Google Scholar

Pilcher, C. D., Eron, J. J., Galvin, S., Gay, C., and Cohen, M. S. (2004). Acute HIV revisited: new opportunities for treatment and prevention. J. Clin. Invest. 113, 937–945. doi: 10.1172/JCI21540

PubMed Abstract | CrossRef Full Text | Google Scholar

Raymond, S., Delobel, P., Mavigner, M., Cazabat, M., Encinas, S., Souyris, C., et al. (2010). CXCR4-using viruses in plasma and peripheral blood mononuclear cells during primary HIV-1 infection and impact on disease progression. AIDS 24, 2305–2312. doi: 10.1097/QAD.0b013e32833e50bb

PubMed Abstract | CrossRef Full Text | Google Scholar

Richman, D. D., Margolis, D. M., Delaney, M., Greene, W. C., Hazuda, D., and Pomerantz, R. J. (2009). The challenge of finding a cure for HIV infection. Science 323, 1304–1307. doi: 10.1126/science.1165706

PubMed Abstract | CrossRef Full Text | Google Scholar

Rieder, P., Joos, B., Scherrer, A. U., Kuster, H., Braun, D., Grube, C., et al. (2011). Characterization of Human Immunodeficiency Virus Type 1 (HIV-1) diversity and tropism in 145 patients with primary HIV-1 infection. Clin. Infect. Dis. 53, 1271–1279. doi: 10.1093/cid/cir725

PubMed Abstract | CrossRef Full Text | Google Scholar

Rose, P. P., and Korber, B. T. (2000). Detecting hypermutations in viral sequences with an emphasis on G -&gt; a hypermutation. Bioinformatics 16, 400–401. doi: 10.1093/bioinformatics/16.4.400

PubMed Abstract | CrossRef Full Text | Google Scholar

Sáez-Cirión, A., Bacchus, C., Hocqueloux, L., Avettand-Fenoel, V., Girault, I., Lecuroux, C., et al. (2013). Post-treatment HIV-1 controllers with a long-term virological remission after the interruption of early initiated antiretroviral therapy ANRS VISCONTI study. PLoS Pathog. 9:e1003211. doi: 10.1371/journal.ppat.1003211

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiffer, V., Deveau, C., Meyer, L., Iraqui, I., Nguyen-Wartel, A., Chaix, M.-L., et al. (2004). Recent changes in the management of primary HIV-1 infection: results from the French PRIMO cohort. HIV Med. 5, 326–333. doi: 10.1111/j.1468-1293.2004.00231.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheppard, H. W., Celum, C., Michael, N. L., O’Brien, S., Dean, M., Carrington, M., et al. (2002). HIV-1 infection in individuals with the CCR5-Δ32/Δ32 genotype: acquisition of syncytium-inducing virus at seroconversion. JAIDS J. Acquir. Immune Defic. Syndr. 29, 307–313. doi: 10.1097/00126334-200203010-00013

CrossRef Full Text | Google Scholar

Strain, M. C., Little, S. J., Daar, E. S., Havlir, D. V., Günthard, H. F., Lam, R. Y., et al. (2005). Effect of treatment, during primary infection, on establishment and clearance of cellular reservoirs of HIV-1. J. Infect. Dis. 191, 1410–1418. doi: 10.1086/428777

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Zyl, G. U., Katusiime, M. G., Wiegand, A., McManus, W. R., Bale, M. J., Halvas, E. K., et al. (2017). No evidence of HIV replication in children on antiretroviral therapy. J. Clin. Invest. 127, 3827–3834. doi: 10.1172/JCI94582

PubMed Abstract | CrossRef Full Text | Google Scholar

von Stockenstrom, S., Odevall, L., Lee, E., Sinclair, E., Bacchetti, P., Killian, M., et al. (2015). Longitudinal genetic characterization reveals that cell proliferation maintains a persistent HIV type 1 DNA pool during effective HIV therapy. J. Infect. Dis. 212, 596–607. doi: 10.1093/infdis/jiv092

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Li, Y., Han, Y., Xie, J., and Li, T. (2017). HIV sequence diversity during the early phase of infection is associated with HIV DNA reductions during antiretroviral therapy. J. Med. Virol. 89, 982–988. doi: 10.1002/jmv.24723

PubMed Abstract | CrossRef Full Text | Google Scholar

Yerly, S., Kaiser, L., Perneger, T. V., Cone, R. W., Opravil, M., Chave, J. P., et al. (2000). Time of initiation of antiretroviral therapy: impact on HIV-1 viraemia. The Swiss HIV cohort study. AIDS 14, 243–249. doi: 10.1097/00002030-200002180-00006

CrossRef Full Text | Google Scholar

Yukl, S. A., Shergill, A. K., McQuaid, K., Gianella, S., Lampiris, H., Hare, C. B., et al. (2010). Effect of raltegravir-containing intensification on HIV burden and T-cell activation in multiple gut sites of HIV-positive adults on suppressive antiretroviral therapy. AIDS 24, 2451–2460. doi: 10.1097/QAD.0b013e32833ef7bb

PubMed Abstract | CrossRef Full Text | Google Scholar