Phylogenetic analysis of HIV-1 archived DNA in blood and gut-associated lymphoid tissue in two patients under antiretroviral therapy
Gut Pathogens volume 13, Article number: 20 (2021)
One of the approaches to cure human immunodeficiency virus (HIV) is the use of therapeutic vaccination. We have launched the Provir/Latitude 45 study to identify conserved CTL epitopes in archived HIV-1 DNA according to the HLA class I alleles in aviremic patients under antiretroviral therapy (ART). A HIV-1 polypeptidic therapeutic vaccine based on viral sequence data obtained from circulating blood was proposed; here, our aim was to compare the proviral DNA in blood and gut-associated lymphoid tissue (GALT). Peripheral blood mononuclear cells and gut biopsies were obtained from two HIV-1 infected patients under successful antiretroviral therapy. Total DNA was extracted including the proviral DNA. The HIV-1 reverse transcriptase was sequenced in both compartments using next generation sequencing followed by single genome sequencing; phylogenetic trees were established and compared. The proviral sequences of both compartments intra-patient exhibited a very low genetic divergence while it was possible to differentiate the sequences inter-patients; single genome sequencing analysis of two couples of samples confirmed that there was no compartmentalization of the sequences intra-patient. We conclude that, considering these two cases, the proviral DNA sequences in blood and GALT are similar and that the epitope analysis of HIV-1 provirus in blood should be considered as relevant to that observed in the GALT, a hard-to-reach major compartment, and can therefore be used for therapeutic vaccine approaches.
HIV-1 infection can be managed by ART, leading to the control of viral replication and improving the health of people living with HIV. However, ART cannot be interrupted since this would lead to a rebound of viral replication [1, 2] as virus establishes cellular (latently infected resting CD4 + memory T cells) and anatomical reservoirs very early during infection [3,4,5,6,7,8]. Gut Associated Lymphoid Tissue (GALT) is considered to be one of the main reservoirs of Simian Immunodeficiency Virus (SIV) and HIV infection [9,10,11]. Cure strategies for HIV-1 include therapeutic vaccination , although immune response observed was not able to control viral replication after ART discontinuation . In this context, we launched the Provir/Latitude 45 project to identify conserved CTL epitopes in the proviral HIV-1 DNA of patients with long-term ART . The study involves in silico modeling based on the HIV-1 proviral DNA sequences, the HLA alleles and the HIV-1 CTL epitopes following sequencing of the archived DNA from peripheral blood mononuclear cells (PBMC), i.e., from circulating blood. Since our initial work was based on proviral DNA in PBMC, we assessed whether our observations would be the same in another compartment, namely GALT. Herein, we present a phylogenetic comparison of the archived HIV-1 DNA in PBMC and GALT from two HIV-1 infected patients at success of ART.
Materials and methods
The study participants (referenced as patients 5A and 10) were recruited at Centre Hospitalier Universitaire de Montreal (CHUM). They were Caucasian individuals infected with an HIV-1 subtype B and under successful antiretroviral therapy (ART) for more than 4 years. For each patient, a sigmoid biopsy was collected during colonoscopy 4 years after initiation of successful ART and processed as previously described . Matched peripheral blood was collected on the same day and immediately processed with Ficoll for PBMC isolation, performed in parallel with cell extraction from biopsy tissue. DNA was extracted from both compartments and used for next generation sequencing (NGS) analysis of HIV-1 provirus; we used the method already published [16, 17] to amplify fragment B, i.e., polymerase (Pol) region including RT and integrase. The PCR products were purified and quantitated, the library was prepared using the Nextera XT DNA Sample Preparation kit; each individual library was then sequenced on a MiSeq Illumina platform. Raw data (FASTQ files) were submitted to the SmartGene® NGS HIV-1 module to generate a BAM file for each patient and each sample was processed for further analysis . The study was carried out using only the Pol RT part region of the sequences obtained. Using Galaxy and Clustal software, RT gene sequences from the two compartments (GALT and PBMC) were selected for neighbor-joining analysis from matrix distances calculated after gapstripping of alignments, with a Kimura two-parameter algorithm and bootstrap analysis. To do so, an alignment was generated that included only reads with lengths > 400 bp corresponding to a given region of RT (variable according to the patient, from amino-acids 39 to 202). This length limitation explains the small number of reads used for this analysis compared to the total number of reads covering this region. Phylogenetic trees were visualized using Interactive Tree of Life (ITOL) software . Single genome sequencing (SGS) was carried out according to the method of Palmer et al. . The total extracted DNA of both compartments was diluted in TE buffer at a dilution yielding a PCR product in three out of 10 PCRs. In this case, according to Poisson’s distribution, the dilution contains one copy of cDNA per positive PCR at about 80% of the time. Two rounds of PCR for RT amplification were followed by visualization of the PCR products. The 1:9 dilution was found to be optimal for Sanger sequencing and the sequences (assuming that there was no mixture of population) of PBMC and GALT obtained were aligned by Clustal to obtain a neighbor-joining tree. For evaluation of evolutionary divergence, the median, mean and range of the number of base substitutions per site between RT sequences were calculated. Analyses were conducted using the Maximum Composite Likelihood model . The analysis involved 4100 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. There was a total of 376 positions in the final dataset. Evolutionary analyses were conducted in Molecular Evolutionary Genetics Analysis (MEGA) 7 software .
The phylogenetic trees are presented in Fig. 1.
NGS analysis of patient 5A shows that the sequences from blood and GALT compartments exhibit a low genetic divergence and are located on the same branch; however, there is a slight divergence between these sequences and the SGS analysis allows to demonstrate that there is a true intermingling of the sequences, therefore evidencing a lack of compartmentalization.
NGS analysis of patient 10 first indicates, compared to the NGS data from patient 5A, that we can fully differentiate the HIV-1 isolates from both patients although they are of the same subtype (B subtype); focusing again on patient 10, we can draw the same conclusion, as for patient 5A, on the low genetic divergence of the sequences between both compartments; when analysing the SGS data, the clonal PBMC sequences are located at the origin of the GALT NGS part of the tree, then found in the PBMC part of the NGS tree while GALT clonal sequences are located at the end of the GALT tree and the origin of the PBMC tree; as for patient 5A we can conclude that there is no evidence of compartmentalization.
To confirm that all the sequences were clustered by patient, we estimated the evolutionary divergence between sequences (Table 1) considering patients 5A and 10. HIV-1 clusters were identified at maximum genetic distances between 4.5 and 7.5% and bootstrap support threshold varied between 70 and 99% . As sequences from patients 5A and 10 are assembled with a median divergence of 5.3% and 2.2% respectively with bootstrap values of 80% and 93%, we confirm that these sequences from GALT and PBMC formed a specific cluster per patient.
Archived viral DNA is found in intestinal tissue at a higher concentration than in PBMCs in ART patients , although the distribution of DNA in CD4 + CCR7 + , transitional memory and effector memory CD4 + T cells is different in blood and intestinal compartments . GALT is therefore a compartment of major importance in the pathophysiology of HIV infection [26, 27]. The data presented here are not related to the quantification of archived DNA but rather to a comparison of the RT sequences obtained by NGS in GALT and PBMC. In that regard, Van Marle et al.  have studied biopsies from infected untreated individuals and sequenced the nef and RT genes of the viral RNA from blood (PBMC) and different parts of the gut by cloning and Sanger sequencing; they concluded that there is a compartmentalization of the virus in the gut reservoir. On the other hand, Lerner et al.  found a low diversity of the GALT and PBMC viruses in patients having experienced a voluntary treatment interruption while Imamichi et al.  did not demonstrate any difference between RNA and DNA sequences from gut and blood of patients chronically infected with HIV-1. Studying HIV-1 infected patients at early and chronic infection stages, Rozera et al.  found a more pronounced compartmentalization of proviral quasispecies in gut compared with PBMC samples in patients with early infection compared with chronic patients. The loss of gut/PBMC compartmentalization in more advanced stages of HIV infection was confirmed by longitudinal observation.
Regarding ART treated patients, Evering et al.  have studied the variability of the proviral DNA in the gut and blood compartments by SGS of the env part of the virus. They showed absence of evolution of the env sequences in the GALT and in PBMC; the authors mention that they cannot rule out the possibility of evolution in other viral genomic regions of HIV-1 such as pol which have not been investigated. Josefsson et al.  have compared the HIV DNA in PBMC and GALT from patients being on successful ART; they have used SGS technology and showed that there is no significant difference between the sequences from both compartments. They concluded that the HIV reservoir is stable on long-term suppressive ART and raise the hypothesis that the population of infected cells exhibiting a low variability of the virus could be maintained by homeostatic cell proliferation.
The patients of our study are similar to those of the two publications mentioned above i.e., ART- treated patients with controlled viral load and therefore a stable viral reservoir; the phylogenetic inferences obtained after NGS evidenced a very low genetic distance between the GALT and PBMC compartments intra-patient. On the other hand, it is possible to differentiate the GALT/PBMC sequences inter-patients; the SGS analysis performed plus the genetic divergence values after NGS and SGS are concordant with a high similarity between proviruses intra-patient. It must be underlined that the SGS technique decreases taq-induced recombination and nucleotide mis-incorporation, providing therefore a more reliable conclusion than conventional cloning .
Among the limitations of our study, we must note the fact that only the RT part of the proviral DNA has been considered and also that we have analyzed global archived DNA molecules without differentiating noninfectious and replication competent genomes ; however, more recent data show that even defective proviral DNA molecules can be expressed and yield viral proteins recognized by CTL T CD8 + lymphocytes .
In conclusion, our results confirm that the proviruses in GALT and PBMC are very similar in these patients under ART and who could be the target population of choice for a therapeutic vaccine and indicate that the analysis of the blood compartment can provide results that can extrapolated to the gut compartment, a major reservoir of HIV.
The strategies used in the future, whether they will be associated with a shock-and-kill approach (and a therapeutic vaccine is a part of this approach) or block-and-lock effect based on HIV silencing , will have to consider the proviral DNA of the archived virus not only in the blood but also in different tissue reservoirs including the gut.
Availability of data and materials
The NGS sequences are available in GenBank under accession number PRJN557560; the Sanger sequences of the SGS study are available in GenBank under accession numbers MN250226 to MN250287.
Human immunodeficiency virus
Simian immunodeficiency virus
Gut associated lymphoid tissue
Next generation sequencing
Peripheral blood mononuclear cells
Single genome sequencing
Interactive tree of life
Molecular evolutionary genetics analysis
Chun T-W, Stuyver L, Mizell SB, Ehler LA, Mican JA, Baseler M, Lloyd AL, Nowak MA, Fauci AS. Presence of an inducible HIV-1 latent reservoir during highly active antiretroviral therapy. Proc Natl Acad Sci USA. 1997;94:13193–7.
Chun T-W, Justement JS, Murray D, Hallahan CW, Maenza J, Collier AC, Sheth PM, Kaul R, Ostrowski M, Moi S, Kovacs C, Fauci AS. Rebound of plasma viremia following cessation of antiretroviral therapy despite profoundly low levels of HIV reservoir: implications for eradication. AIDS. 2010;24:2803–8.
Finzi D, Hermankova M, Pierson T, Carruth LM, Buck C, Chaisson RE, Quinn TC, Chadwick K, Margolick J, Brookmeyer R, Gallant J, Markowitz M, Ho DD, Richman DD, Siliciano RF. Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. Science. 1997;278:1295–300.
Chun T-W, Fauci AS. Viral persistence in HIV infection: much known, much to learn. J Infect Dis. 2013;208:1356–8.
Chun T-W, Moir S, Fauci AS. HIV reservoirs as obstacles and opportunities for an HIV cure. Nat Immunol. 2015;16:584–9.
Wong JK, Yukl SA. Tissue reservoirs of HIV. Curr Opin HIV AIDS. 2016;11:362–70.
Barton K, Winckelmann A, Palmer S. HIV-1 reservoirs during suppressive therapy. Trends Microbiol. 2016;24:345–55.
Miller RL, Ponte R, Jones BR, Kinloch NN, Omondi FH, Jenabian M-A, Dupuy FP, Fromentin R, Brassard P, Mehraj V, Chomont N, Poon AFY, Brumme ZL, Routy J-P, The ORCHID Study Group. HIV diversity and genetic compartmentalization in blood and testes during suppressive antiretroviral therapy. J Virol. 2019. https://doi.org/10.1128/JVI.00755-19.
Mattapallil JJ, Douek DC, Hill B, Nishimura Y, Martin M, Roederer M. Massive infection and loss of memory CD4 + T cells in multiple tissues during acute SIV infection. Nature. 2005;434:1093–7.
Chun T-W, Nickle DC, Justement JS, Meyers JH, Roby G, Hallahan CW, Kottilil S, Moir S, Mican JM, Mullins JI, Ward DJ, Kovacs JA, Mannon PJ, Fauci AS. Persistence of HIV in gut-associated lymphoid tissue despite long-term antiretroviral therapy. J Infect Dis. 2008;197:714–20.
Estes JD, Kityo C, Ssali F, Swainson L, Makamdop KN, Del Prete GQ, Deeks SG, Luciw PA, Chipman JG, Beilman GJ, Hoskuldsson T, Khoruts A, Anderson J, Deleage C, Jasurda J, Schmidt TE, Hafertepe M, Callisto SP, Pearson H, Reimann T, Schuster J, Schoephoerster J, Southern P, Perkey K, Shang L, Wietgrefe SW, Fletcher CV, Lifson JD, Douek DC, McCune JM, Haase AT, Schacker TW. Defining total-body AIDS-virus burden with implications for curative strategies. Nat Med. 2017;23:1271–6.
Deeks SG, Lewin SR, Ross AL, Ananworanich J, Benkirane M, Cannon P, et al. International AIDS society global scientific strategy: towards an HIV cure. Nat Med. 2016;22:839–50.
Fleury H, Tumiotto C, Bellecave P, Recordon-Pinson P. Therapeutic vaccine against HIV, viral variability, cytotoxic T lymphocyte epitopes, and genetics of patients. AIDS Res Hum Retrovir. 2018;34:27–30.
Tumiotto C, Alves BM, Recordon-Pinson P, Jourdain M, Bellecave P, Guidicelli G-L, Visentin J, Bonnet F, Hessamfar M, Neau D, Sanchez J, Brander C, Sajadi M, Eyzaguirre L, Soares EA, Routy JP, Soares MA, Fleury H. Provir/latitude 45 study: a step towards a multi-epitopic CTL vaccine designed on archived HIV-1 DNA and according to dominant HLA I alleles. PLoS ONE. 2019;14:e0212347.
Gosselin A, Wiche Salinas TR, Planas D, Wacleche VS, Zhang Y, Fromentin R, Chomont N, Cohen EA, Shacklett B, Mehraj V, Ghali MP, Routy JP, Ancuta P. HIV persists in CCR6 + CD4 + T cells from colon and blood during antiretroviral therapy. AIDS. 2017;31:35–48.
Ode H, Matsuda M, Matsuoka K, Hachiya A, Hattori J, Kito Y, Yokomaku Y, Iwatani Y, Sugiura W. Quasispecies analyses of the hiv-1 near-full-length genome with illumina MiSeq. Front Microbiol. 2015;6:1258.
Alves BM, Siqueira JD, Garrido MM, Botelh OM, Prellwitz IM, Ribeiro SR, Soares EA, Soares MA. Characterization of HIV-1 near full-length proviral genome quasispecies from patients with undetectable viral load undergoing first-line HAART therapy. Viruses. 2017;9:E392.
Bellecave P, Recordon-Pinson P, Fleury H. Evaluation of automatic analysis of ultradeep pyrosequencing raw data to determine percentages of HIV resistance mutations in patients followed-up in hospital. AIDS Res Hum Retrovir. 2016;32:85–92.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242-5.
Palmer S, Kearney M, Maldarelli F, Halvas EK, Bixby CJ, Bazmi H, Rock D, Falloon J, Davey RTJ, Dewar RL, Metcalf JA, Hammer S, Mellors JW, Coffin JM. Multiple linked human immunodeficiency virus type 1 drug resistance mutations in treatment-experienced patients are missed by standard genotype analysis. J Clin Microbiol. 2005;43:406–13.
Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. 2004;101:11030–5.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Ragonnet-Cronin M, Hodcroft E, Hué S, Fearnhill E, Delpech V, Brown AJ, Lycett S. UK HIV drug resistance database. Automated analysis of phylogenetic clusters. BMC Bioinform. 2013;14:317.
Moron-Lopez S, Puertas MC, Galvez C, Navarro J, Carrasco A, Esteve M, Manye J, Crespo M, Salgado M, Martinez-Picado J. Sensitive quantification of the HIV-1 reservoir in gut-associated lymphoid tissue. PLoS ONE. 2017;12:e0175899.
Yukl SA, Shergill AK, Ho T, Killian M, Girling V, Epling L, Li P, Wong LK, Crouch P, Deeks SG, Havlir DV, McQuaid K, Sinclair E, Wong JK. The distribution of HIV DNA and RNA in cell subsets differs in gut and blood of. J Infect Dis. 2013;208:1212–20.
Thompson CG, Gay CL, Kashuba ADM. HIV persistence in gut-associated lymphoid tissues: pharmacological challenges and opportunities. AIDS Res Hum Retrovir. 2017;33:513–23.
Augustin M, Horn C, Ercanoglu MS, Bondet V, Suarez I, Chon S-H, Nierhoff D, Knops E, Heger E, Oette M, Klein F, Fätkenheuer G, Duffy D, Müller-Trutwin M, Lehmann C. Treasure hunt: HIV DNA in distinct subsets in the blood, the terminal ileum, and rectum. Conference on Retroviruses and opportunistic infections. Seattle, USA. 2019. Abstract 335.
van Marle G, Gill MJ, Kolodka D, McManus L, Grant T, Church DL. Compartmentalization of the gut viral reservoir in HIV-1 infected patients. Retrovirology. 2007;4:87.
Lerner P, Guadalupe M, Donovan R, Hung J, Flamm J, Prindiville T, Sankaran-Walters S, Syvanen M, Wong JK, George MD, Dandekar S. The gut mucosal viral reservoir in HIV-infected patients is not the major source of rebound plasma viremia following interruption of highly active antiretroviral therapy. J Virol. 2011;85:4772–82.
Imamichi H, DeGray G, Dewar RL, Mannon P, Yao M, Chairez C, Sereti I, Kovacs JA. Lack of compartmentalization of HIV-1 quasispecies between the gut and peripheral blood compartments. J Infect Dis. 2011;204:309–14.
Rozera G, Abbate I, Vlassi C, Giombini E, Lionetti R, Selleri M, Zaccaro P, Bartolini B, Corpolongo A, D’Offizi G, Baiocchini A, Del Nonno F, Ippolito G, Capobianchi MR. Quasispecies tropism and compartmentalization in gut and peripheral blood during early and chronic phases of HIV-1 infection: possible correlation with immune activation markers. Clin Microbiol Infect. 2014;20:157–66.
Evering TH, Mehandru S, Racz P, Tenner-Racz K, Poles MA, Figueroa A, Mohri H, Markowitz M. Absence of HIV-1 evolution in the gut-associated lymphoid tissue from patients on combination antiviral therapy initiated during primary infection. PLoS Pathog. 2012;8:e1002506.
Josefsson L, von Stockenstrom S, Faria NR, Sinclair E, Bacchetti P, Killian M, Epling L, Tan A, Ho T, Lemey P, Shao W, Hunt PW, Somsouk M, Wylie W, Douek DC, Loeb L, Custer J, Hoh R, Poole L, Deeks SG, Hecht F, Palmer S. 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 USA. 2013;110:E4987-96.
Bruner KM, Murray AJ, Pollack RA, Soliman MG, Laskey SB, Capoferri AA, Lai J, Strain MJ, Lada SM, Hoh R, Ho Y-C, Richman DD, Deeks SG, Siliciano JD, Siliciano RF. Defective proviruses rapidly accumulate during acute HIV-1 infection. Nat Med. 2016;22:1043–9.
Pollack RA, Jones RB, Pertea M, Bruner KM, Martin AR, Thomas AS, Capoferri AA, Beg SA, Huang S-H, Karandish S, Hao H, Halper-Stromberg E, Yong PC, Kovacs C, Benko E, Siliciano RF, Ho Y-C. Defective HIV-1 proviruses are expressed and can be recognized by cytotoxic T lymphocytes which shapes the proviral landscape. Cell Host Microbe. 2017;21:494–50636.
Vansant G, Bruggemans A, Janssens J, Debyser Z. Block-and-lock strategies to cure HIV infection. Viruses. 2020;12:84.
We thank Lucie Richard, an engineering student at the university of Bordeaux, for her excellent technical collaboration. The valuable advices and comments by Z Brumme (British Columbia Centre for Excellence in HIV/AIDS, Vancouver, BC, Canada) were greatly appreciated. We would like to thank Maria Fraraccio for her administrative support in Montreal and Fredéric Perry for his administrative assistance at the DRCI CHU in Bordeaux. Ray Cooke is also thanked for copyediting.
The study was funded by a grant from MSD Avenir (DS-2016-005) and a financial contribution from GERMATAN (Groupe d'Etude et de Recherche sur les Maladies Animales Transmissibles et les Anthropozoonoses). Sequencing of proviral DNA by Illumina technology was performed at the Genome Transcriptome Facility in Bordeaux (grants from the Conseil Régional d’Aquitaine n 20030304002FA and 20040305003FA, from the European Union FEDER n 2003227 and from Investissements d'Avenir ANR-10-EQPX-16-01).
Ethics approval and consent to participate
This study using colon biopsies and blood from patients living with HIV (PLWH) was approved by the Research Ethics Boards of the Research Institute of McGill University Health Centre (MUHC), Montreal, QC, Canada and the Centre Hospitalier de l’Université de Montréal (CHUM) Research Centre, Montréal, QC, Canada. The study also received approval from the Comité d’Ethique du Sud Ouest et Outremer (DC 2012/48) and was conducted in accordance with the Declaration of Helsinki. Each participant provided written informed consent before any study procedure.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Recordon-Pinson, P., Gosselin, A., Ancuta, P. et al. Phylogenetic analysis of HIV-1 archived DNA in blood and gut-associated lymphoid tissue in two patients under antiretroviral therapy. Gut Pathog 13, 20 (2021). https://doi.org/10.1186/s13099-021-00416-6
- Proviral DNA
- Gut associated lymphoid tissue (GALT)
- Phylogenetic comparison of archived viral sequences