https://orcid.org/0000-0001-5333-8943
Figures
Abstract
Human immunodeficiency virus type 1 (HIV-1) recombinants in the world are believed to be generated through recombination between distinct HIV-1 strains among coinfection or superinfection cases. However, direct evidence to support transmission of HIV-1 recombinants from a coinfected/superinfected donor to putative recipient is lacking. Here, we report on the origin and evolutionary relationship between a set of recombinants from a CRF01_AE/CRF07_BC superinfected putative donor and diverse CRF01_AE/CRF07_BC recombinants from five putative recipients. Interviews on sociodemographic characteristics and sexual behaviors for these six HIV-1-infected men who have sex with men showed that they had similar ways of partner seeking: online dating sites and social circles. Phylogenetic and recombination analyses demonstrated that the near-full-length genome sequences from six patients formed a monophyletic cluster different from known HIV-1 genotypes in maximum likelihood phylogenetic trees, were all composed of CRF01_AE and CRF07_BC fragments with two common breakpoints on env, and shared 4–7 breakpoints with each other. Moreover, 3’ half-genomes of recombinant strains from five recipients had identical/similar recombinant structures with strains at longitudinal samples from the superinfected donor. Recombinants from the donor were paraphyletic, whereas five recipients were monophyletic or polyphyletic in the maximum clade credibility tree. Bayesian analyses confirmed that the estimated time to the most recent common ancestor (tMRCA) of CRF01_AE and CRF07_BC strains of the donor was 2009.2 and 2010.7, respectively, and all were earlier than the emergence of recombinants from five recipients. Our results demonstrated that the closely related unique recombinant forms of HIV-1 might be the descendent of a series of recombinants generated gradually in a superinfected patient. This finding highlights the importance of early initiation of antiretroviral therapy as well as tracing and testing of partners in patients with multiple HIV-1 infection.
Author summary
Recombination is a major mechanism for rapid evolution and diversification of human immunodeficiency virus type 1 (HIV-1). Multiple HIV-1 infection is believed to be a prerequisite for generating new recombinant strains. Recombinants sharing some breakpoints might have a common ancestor or direct parental/progenitor relationship or be irrelevant in evolution. However, evidence to support transmission of HIV-1 recombinants from coinfected/superinfected donor to putative recipients with recombinant descendants is lacking. Here, we report on a set of HIV-1 CRF01_AE/CRF07_BC unique recombinant forms (URFs) from six HIV-1-infected men who have sex with men (MSM) in China with similar recombination structures and homologous parental strains. Moreover, the 3’ half-genomes of all recombinants from five recipients had similar or identical forms with earlier continuous samples of a superinfected donor, and presented paraphyletic-monophyletic or paraphyletic-polyphyletic relationships in maximum clade credibility trees. These data suggest that HIV-1 superinfection might be a potential source of transmission of a series of closely related HIV-1 CRF01_AE/CRF07_BC URFs. Our study emphasizes the importance of continuous surveillance and early initiation of antiretroviral therapy in HIV-1-superinfected individuals with active high-risk behavior and substance abuse to prevent the generation and spread of recombinant strains among HIV-infected individuals.
Citation: Gao Y, He S, Tian W, Li D, An M, Zhao B, et al. (2021) First complete-genome documentation of HIV-1 intersubtype superinfection with transmissions of diverse recombinants over time to five recipients. PLoS Pathog 17(2): e1009258. https://doi.org/10.1371/journal.ppat.1009258
Editor: Ronald Swanstrom, University of North Carolina at Chapel Hill, UNITED STATES
Received: August 31, 2020; Accepted: December 22, 2020; Published: February 12, 2021
Copyright: © 2021 Gao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All sequences in this paper are available from the GenBank (accession numbers KX434794, KX434795, KX434797, KX434798 KX434799, MT857722, MW287665-MW287747 and MW344769-MW344807).
Funding: This work was supported by the Mega-projects of National Science Research for the 13th Five-Year Plan (2017ZX10201101) awarded to H.S. and Mega-projects of National Science Research for the 13th Five-Year Plan (2018ZX10721102006003) and the National Natural Science Foundation of China (81871637) awarded to X.H. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human immunodeficiency virus type 1 (HIV-1) is characterized by extensive genetic diversity. Recombination is a major mechanism for the rapid evolution and diversification of HIV-1 [1]. More than 100 circulating recombinant forms (CRFs), along with massive unique recombinant forms (URFs) in HIV-1 group M have been identified worldwide [2] (www.hiv.lanl.gov). Some recombinants have recombined further with other subtypes or CRFs to generate second-generation recombinants [3]. It has been estimated that HIV-1 recombinants, including CRFs (16.7%) and URFs (6.1%), accounted for 22.8% of epidemics globally between 2010 and 2015 [4]. Recombination of HIV-1 can potentially change biological characteristics, fitness, susceptibility to antiretroviral drugs, disease progression, as well as the diagnostic accuracy of serology- and molecular- based assays [5–8].
Identification of HIV-1 recombinants from primary infected individuals usually indicates the spreading of these recombinants among a population, but the origin and transmission history are incompletely understood [9]. Deciphering the origin of a group of recombinants with high genetic similarity using phylogenetic analyses is difficult. Studies have reported that recombinants sharing some breakpoints might have a common recombinant ancestor (e.g., CRF07_BC and CRF08_BC in China [10,11] and HIV-1 BF intersubtype recombinant viruses in Argentina [12]) or that there might be a direct parental/progenitor relationship between them (e.g., CRF48_01B and CRF74_01B were probably descended from CRF33_01B in Malaysia [13,14]) or that they may be irrelevant in evolution. Common breakpoints may be attributable to fragile sites or hairpin structure of genomic RNA, pause sites, high-pairing probability, or sequence similarity during reverse transcription. In this instance, the breakpoints often occur in well-conserved regions of viral genomes. These potential mechanisms of recombination have been supported by in vitro experiments and mathematical models, but they may be more complex in vivo [15–18]. Therefore, elucidation of the recombination mechanisms of HIV-1 may help to guide the surveillance and prevention of HIV spread.
Multiple HIV-1 strains infecting the same person concurrently (“coinfection”) or one after another (“superinfection”) is believed to be the prerequisite for the generation of recombinant strains [19–21], which is supported indirectly by the overlap of “hot areas” for HIV recombinants and multiple infections [22–24]. Moreover, recombinants composed of more than two parental viruses have also been found in HIV-1 superinfected cases [25,26]. However, there is no direct evidence to support the transmission of HIV-1 recombinants from coinfection/superinfection cases to a putative recipient with recombinant descendants, let alone the origin and transmission history of a group of HIV recombinant strains with genetic similarity.
We depicted the origin and evolutionary relationship among a group of closely related CRF01_AE/CRF07_BC URFs (0107 URFs) between an HIV-1 intersubtype superinfected donor and five putative CRF01_AE/CRF07_BC-infected recipients. This finding emphasized the importance of early initiation of antiretroviral therapy (ART) as well as tracing and testing of partners with multiple HIV-1 infection to prevent the spread of recombinant strains.
Results
Sociodemographic characteristics and sexual behaviors of six HIV-1-infected men who have sex with men (MSM)
Previously, we identified six patients infected with HIV-1 CRF01_AE/CRF07_BC strains in a newly diagnosed cohort of MSM in Liaoning, northeast China. Among them, the donor was diagnosed with recent HIV-1 infection on 3 March 2010 and started ART on 7 January 2014. His viral load was well-controlled to <100 copies/mL after that. He is a high-earning businessman and self-reported as seeking younger male partners via online dating sites and social circles. He reported sexual behaviors with one regular partner and ≥10 casual partners in the past 3 months before the diagnosis of HIV infection. Insertive and receptive positions were adopted when he had sex with other males without a condom. Moreover, methamphetamine and “rush poppers” were used constantly during sex.
All five recipients were diagnosed with HIV infection between 2013 and 2014. Among them (hereafter termed recipient), recipient 1, 2, 3, and 4 had been diagnosed with a recent HIV infection according to the results of limiting antigen (LAg)-avidity enzyme immunoassay (EIA). Recipient 5 was estimated to have become infected with HIV before the end of 2013. Five recipients were all ~10 years younger and had a lower income than that of the donor. They usually sought male partners via online dating sites and social circles. They self-reported to have had sex with 3–15 male partners within the last 3 months before the diagnosis of HIV infection except recipient 3, a “rent boy,” who had 80 commercial partners and one regular partner (recipient 5). When they had sex with other males, recipient 2, 3, and 4 adopted insertive and receptive positions; recipient 1 and 5 preferred insertive and receptive positions, respectively. Recipient 2 and 4 engaged in condom-less anal sex. Besides, recipient 1, 2, and 3 had a history of substance abuse. The sociodemographic characteristics and sexual behaviors of six HIV-1-infected MSM are outlined in Table 1.
Identification of a lineage of HIV-1 CRF01_AE/CRF07_BC URFs among MSM in Liaoning
To validate the lineage of HIV-1 new recombinants, the near-full-length genome (NFLG) (HXB2: 790–9601 bp) was used for phylogenetic analyses (Fig 1A). The NFLG of six patients formed a distinct monophyletic cluster that was separate from any other known subtypes, CRFs and URFs in the maximum likelihood (ML) phylogenetic tree with poster probability value 1, suggesting that this was a lineage of new HIV-1 URFs.
(A) The phylogenetic analysis of NFLG from six HIV-1-infected MSM (marked red) and reference sequences which were downloaded from HIV database (marked black). Maximum likelihood phylogenetic trees were constructed by Fast Tree and edited by Fig tree v1.4.2. The scale length represents 6% nucleotide sequence divergence. (B) Recombination analysis of six HIV-1 NFLG by using SimPlot (v3.5.1). CRF01_AE and CRF07_BC were used as putative parental sequences according to the online tool HIV Blast and subtype G (14003718) was used as an outgroup. The parameters of analysis were the default except the window size was 350 bp and step size was 50 bp, bootstrap replicate was 250. NFLG, near full-length genomes; MSM, men who have sex with men.
Six HIV-1 CRF01_AE/CRF07_BC URFs showed similar recombination forms and homologous parental strains
The recombination forms along the whole genomes of six strains were first screened with Recombinant Identification Program (RIP) and jumping profile hidden Markov model (jpHMM), and then validated with Simplot Six strains had similar (but not identical) recombination forms (Fig 1B). First, six strains were all CRF01_AE/ CRF07_BC recombinants with a CRF01_AE backbone and three or four CRF07_BC insertions in pol, vpr, tat/rev, env, and nef. Second, six strains shared two identical breakpoints: the fifth and sixth breakpoint in env; five strains shared the first and second breakpoints in pol except in the donor; four strains shared the seventh and eighth breakpoint in nef except in recipient 1 and the donor. Recipient 1, 2, and 3 and the donor shared the third and fourth breakpoint in vpr and first exon of tat/rev. Third, the main difference among recipient 3, 4, and 5 was the length of the CRF07_BC segment IV, which was 404 bp, 970 bp, and 644 bp, respectively.
To determine the evolutionary relationship and potential parental strains of the six strains, we undertook phylogenetic analyses on sub-regions between breakpoints (Fig 2). In general, six strains displayed monophyletic clustering in all the sub-region trees of CRF01_AE and CRF07_BC segments. In trees I, III, V, and VII+IX, the CRF01_AE segments from six strains belonged to CRF01_AE lineage 4 among the Chinese MSM population (poster probability = 1) (Fig 2A). In trees II, IV, VI, and VI+VIII, the CRF07_BC segments of six strains belonged to the CRF07_BC lineage predominant among the Chinese MSM population (poster probability ≥0.96) (Fig 2B). Taken together, this high genetic similarity and homologous parental strains suggested a close evolutionary relationship among the six CRF01_AE/CRF07_BC recombinant strains.
Six HIV-1 CRF01_AE/CRF07_BC recombinants were divided into nine segments (I-IX) according to the recombinant breakpoints shown in the HIV-1 structure diagram in Fig 1B. ML trees were constructed by Fast tree and edited by Fig tree v1.4.2, poster probability values over 0.90 were shown at each node. (A) The subtypes of segment I, III, V, VII, IX were CRF01_AE. (B) The subtypes of segment II, IV, VI, VIII were CRF07_BC. The segment VIII and IX were too short in some patients and were combined the segment VI and VII, respectively to construct ML trees. The scale bars were shown at the bottom of each sub-region tree. ML tree, Maximum likelihood phylogenetic trees.
Most likely origin of the series of HIV-1 CRF01_AE/CRF07_BC URFs
Among six patients, according to clinical records and results for LAg-avidity EIA, the donor was diagnosed to have been infected with an HIV-1 CRF01_AE strain on 3 March 2010 (Table 1). Moreover, the donor had been identified to have superinfected another CRF07_BC strain in our previous study [27]. To determine the recombination process between the primary infected CRF01_AE strain and superinfected CRF07_BC strain in the donor, the 3′ half-genome was obtained from longitudinal samples by a single-genome amplification (SGA) strategy and used for recombination analyses (Fig 3). The superinfected CRF07_BC strain was obtained first through the SGA strategy detected from the donor on 9 December 2010. The recombination between CRF01_AE and CRF07_BC was first detected ~3 months after superinfection. The predominant CRF01_AE/CRF07_BC recombinants were detected ~7 months and ~12 months after superinfection, respectively. Although the quasispecies and composition of the CRF01_AE/CRF07_BC recombinants among longitudinal samples of the donor changed continuously, some recombinants could be detected at ≥2 time points. Also, some identical/similar breakpoints were detected between distinct recombinants at different time points (S1A and S2 Figs), which suggested continuous evolution of CRF01_AE/CRF07_BC recombinants in vivo under immune selection from the host.
The first column represents the time interval from the first detection of superinfecting virus to this sampling date. The second and fifth columns represent the sampling dates of the donor and five recipients, respectively. The third and fourth columns described the recombination forms of the donor and five recipients, respectively. The sequence name is composed of patient’s ID+ sampling date + recombination forms + (the number of sequences with identical recombination forms at this time point). The different recombination forms at each sampling date are marked by Roman numerals. Initial CRF01_AE strains and superinfecting CRF07_BC strains were marked as light-coral and slate-blue, respectively. The arrow indicates similar or identical recombination structures were identified between the donor and 5 putative recipients.
We further compared the recombination forms of 3′ half-genomes of the donor at longitudinal samples and five recipients at baseline or after seroconversion (Fig 3). The viral quasispecies of the donor were more complex than those of the five recipients who had only one or two recombination forms. More importantly, each recombinant strain from five patients had at least one similar or identical recombination form with the strain from the donor sample at an earlier time point (Figs 3, and S1B and S2). For example, two forms of recombinants from recipient 1 (1 April 2013) were similar with recombinants from donor samples on 13 July 2011 and 29 November 2011, respectively. Similarly, two forms of recombinants from recipient 3 (12 May 2014) were similar or identical to recombinants from donor samples on 22 January 2013 and 23 September 2013, respectively. Moreover, the recombinants from recipients 2 (28 February 2013), recipient 4 (20 March 2014), and recipient 5 (12 May 2014) all had identical forms with those of earlier longitudinal samples from the donor. These data suggested that recombinants from five recipients were most likely being transmitted from the donor directly or indirectly.
Temporal evolutionary relationships among HIV-1 CRF01_AE/CRF07_BC URFs
To investigate the possible evolutionary relationship between these closely related CRF01_AE/CRF07_BC recombinants from six patients, the concatenated CRF01_AE segments and CRF07_BC segments were analyzed by Bayesian molecular clocks (Fig 4A and 4B). In the maximum clade credibility (MCC) tree of CRF01_AE and CRF07_BC, the sequences from six patients formed a monophyletic cluster within all reference sequences, respectively. Within each cluster, the initial infected CRF01_AE strains and the superinfected CRF07_BC strains from the donor were located at the root. The strains of five recipients formed internal branches following the CRF01_AE strains and CRF07_BC strains from the donor, respectively. Moreover, the strains from the donor were paraphyletic, whereas the strains from the five recipients were monophyletic or polyphyletic in the MCC tree. This result was further supported by the evolutionary ML trees constructed with concatenated CRF01_AE and CRF07_BC segments, respectively (S3 Fig).
Some representative sequences of the donor, all sequences of five recipients and references were selected to constructed MCC trees. Bayesian molecular clock analysis were performed using BEAST v2.5.1 for concatenated CRF01_AE segments (Regions V+VII + IX) (A) and CRF07_BC segments (Regions IV+VI) (B). The medians of the tMRCA and 95% HPD of nodes relevant to this study were shown as follows: African CRF01_AE (01.AF) ancestor, 1981.2 (1979.3–1983.0); Thai CRF01_AE (01.TH) ancestor, 1986.9 (1985.6–1988.1); Vietnam CRF01_AE (01.VN) ancestor, 1993.0 (1991.6–1994.5); CRF01_AE Lineage 5 (01.L5), 2001.3 (1999.7–2002.7); CRF01_AE Lineage 4 (01.L4), 1995.7 (1993.1–1998.1). C ancestor, 1987.8 (1985.6–1989.8); CRF07_BC ancestor 1995.9 (1995.2–1996.5). The tMRCA of concatenated CRF01_AE segments of the donor at sampling date 13-Jul-2011, 29-Nov-2011, 22-Jan-2013, 23-Sep-2013 are 2010.5 (2010.3–2010.7), 2010.5 (2010.3–2010.7), 2012.1 (2011.8–2012.6), 2011.7 (2011.4–2012.1). The tMRCA of concatenated CRF07_BC segments of the donor at sampling date 13-Jul-2011, 29-Nov-2011, 22-Jan-2013, 23-Sep-2013 are 2011.1 (2010.9–2011.3), 2011.2 (2011.1–2011.5), 2011.2 (2011.0–2011.5), 2012.3 (2012.0–2012.7). tMRCA, time of the most common ancestor; 95% HPD, 95% highest probability density.
The estimated time to the most recent common ancestor (tMRCA) for the concatenated initial infected CRF01_AE strain and superinfected CRF07_BC strain of the donor was dated to 2009.2 (95% highest probability density 2008.8–2009.7) and 2010.7 (2010.5–2010.9), respectively (Fig 4A and 4B). These observations were consistent with the results using LAg-avidity EIA and SGA on estimation and identification of initial infection and superinfection. The tMRCA of the concatenated CRF01_AE segments in the donor (13 July 2011, 29 November 2011, 22 January 2013, and 23 September 2013) was earlier than that of the corresponding similar/identical recombinants of the five recipients. The tMRCA of concatenated CRF01_AE segments of recipient 1, 2, 3, 4, and 5 was estimated to be 2012.0 (2011.8 to 2012.4), 2011.9 (2011.5 to 2012.2), 2013.1 (2012.6 to 2013.6), 2013.3 (2013.0 to 2013.7), and 2012.7 (2012.4 to 2013.2), respectively (Fig 4A). Similar results were found for the tMRCA of CRF07_BC segments of the six patients (Fig 4B).
Discussion
The rapid increase and high diversity of HIV-1 recombinants are a great challenge for the prevention and surveillance of HIV infection globally. The co-circulation of various strains and high rate of multiple infection are prerequisites for a generation of new HIV-1 recombinants. We identified a lineage of new 0107 URFs with homologous parental strains and similar recombination forms among six HIV-1-infected MSM in northeast China. Epidemiology, quasispecies diversity, timelines, and phylogenetics supported the notion that the potential origin of the ancestral virus of a series of closely related 0107 URFs could be traced back to an HIV-1-superinfected individual.
First, HIV-1-infected MSM carry a higher risk of multiple infection compared with that in heterosexuals [28]. Junjie Xu and colleagues reported that ~35% of MSM had >5 male sexual partners in the last 12 months in Shenyang, Liaoning [29]. In our study, six patients all sought multiple male sexual partners mainly through online dating sites and social circles in the same area. Moreover, most of them seldom used condoms, admitted substance abuse during anal/oral sex with their casual male sexual partners, and had a history of syphilis, which further increased the potential for HIV infection and acquisition of multiple viruses [30–34].
Second, the donor was an HIV-1 CRF01_AE/CRF07_BC-superinfected patient. This patient developed a series of CRF01_AE/CRF07_BC recombinants with 1–5 common breakpoints but non-identical recombination forms in longitudinal samples within a treatment-naïve period of ~4 years. This finding is consistent with the observation by McCutchan and colleagues that various recombinants from a heterosexual superinfected individual in Tanzania had common breakpoints in gag, env, and gp41/nef regions in serial samples within 30 months [26]. Compared with the method used by McCutchan and colleagues to amplify three regions of the HIV (multiple-region hybridization assays), we used the SGA method to amplify the relatively long fragments of the 3’ half-genome, which could fully reflect the complex structure of the recombinants in six patients [35]. Five recipients in our study had one or two genetically homogeneous CRF01_AE/CRF07 recombinants at baseline or after seroconversion, which showed less quasispecies complexity compared with that in the donor. Surprisingly, the recombination forms of all strains from five patients resembled those at different sampling dates from the donor. Studies have demonstrated the loss of genetic diversity of the HIV within donors to new hosts upon sexual transmission [36] and mother-to-child transmission [37], which may result from a “transmission bottleneck”. These data imply that this HIV-1-superinfected patient might have been the donor who transmitted a series of 0107 URFs to new hosts.
Third, according to clinical records and the results of LAg-avidity EIA, we found that the putative donor was diagnosed with HIV infection about 3–4 years earlier than the other five patients, and was the first to be detected with CRF01_AE/CRF07_BC recombinants. The tMRCA of the recombinants of the putative donor was also earlier than that of the other five patients. Furthermore, the putative donor did not start ART until 2014 (i.e., ~4 years after the diagnosis).
Finally, based on Bayesian analyses, the topology of the MCC tree provided compelling evidence of the source of these 0107 URFs. Strains from the putative donor and five recipients formed paraphyletic–polyphyletic or paraphyletic–monophyletic donor-recipient joint phylogeny. Studies have reported that phylogenetic methods might be used to infer the transmission history of epidemiologically linked hosts in an HIV-mono-infected population [38]. Paraphyletic–polyphyletic trees support direct transmission and are believed to exclude intervening transmission and a common source. Typically, paraphyletic–monophyletic trees result from direct or indirect transmission [39].
In recent years, several HIV-1 CRF01_AE and B-related second-generation recombinants have been identified in Asia [40–44], some of which (e.g., CRF55_01B and CRF59_01B) have spread widely around China [45]. More recently, massive new recombinants composed of CRF01_AE and CRF07_BC (the two predominant strains among MSM populations) have also been reported around China [46–49]. However, most of them were unrelated in terms of phylogenetics, and a few were clustered but the origin and evolutionary relationship among them were not clear. We demonstrated, for the first time, that a group of closely related HIV-1 new recombinants in MSM may derive from one superinfection case, which supported the hypothesis of a model of generation of HIV-1 BF intersubtype recombinants with coincident breakpoints from South America [50].
Recently, early initiation of ART has become a worldwide public-health prevention strategy for HIV transmission. In China, the standards of ART initiation have been updated several times in national treatment guidelines [51–54]. In 2014, it was suggested that people with HIV infection with CD4+ T-cell count <500 cells/μL should receive ART. In 2016, it was suggested that all HIV-infected people, regardless of the CD4+ T-cell count, should be treated. Therefore, the putative donor reported in the present study did not start ART until 2014 according to the policy in China at that time. During the treatment-naïve period (~4 years), he became superinfected with another HIV strain and developed a series of CRF01_AE/CRF07_BC recombinants, then infected the other five recipients directly or indirectly.
Treatment-as-prevention approach has been shown to reduce the risk of HIV transmission in serodiscordant couples [55], and has been proposed and implemented in many countries (including China). Therefore, the prevalence of superinfection (such as in the donor in the present study) might be reduced. However, detection of increased URFs in areas with multiple HIV strains suggest there are many undiagnosed multiple-infected cases. Hence, if a multiple-infected case is diagnosed, not only should ART be started early to reduce the transmission risk, strengthened tracing of partners should also be done immediately. In this way, persons infected with complicated HIV-1 strains and a higher risk of further transmission can be diagnosed rapidly.
Due to protection of personal privacy, we did not have sufficient epidemiological data to determine the transmission chain among the six patients. However, we provided evidence from different perspectives to suggest there might be a direct or indirect transmission relationship among our six patients. We also inferred that the patient with HIV-1 superinfection might be the source of a lineage of closely related 0107 URFs (Fig 5).
This hypothesis is based on the analysis of recombination forms and the estimated tMRCA of concatenated CRF01_AE and CRF07_BC segments of six HIV-1-infected patients. The superinfected patient is the putative donor and the other five patients are the putative recipients. The same color between donor and five recipients represents the similar or identical recombinants between them. The dotted line indicates that the recombinants might be transmitted directly or indirectly from the donor to putative five recipients. The tMRCA shown here is the mean tMRCA of concatenated CRF01_AE and CRF07_BC segments, 95% HPD is shown in parenthesis. tMRCA, time of the most common ancestor; 95% HPD, 95% highest probability density.
Our study suggests an important role of HIV-1 superinfection on the generation and transmission of new recombinants. Furthermore, recombinants with high genetic similarity (but distinct recombination forms) could share a common origin. This observation provides a new perspective to infer the evolutionary relationship between HIV-infected individuals harboring recombinants. The present study also calls for greater attention to the monitoring, early ART, and strengthened management of multiple-infected individuals, including the tracing and testing of partners.
Methods
Ethics statement
This study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University ([2018] 2015-140-5). Written informed consent to participate in this study was obtained from all patients before sample collection.
Study participants
The six study participants were newly diagnosed HIV-1-infected patients from a cohort of MSM in a voluntary HIV counseling and testing clinic of the First Affiliated Hospital of China Medical University. They were found to be infected with a lineage of CRF01_AE/CRF07_BC recombinant strains through phylogenetic analyses on pol sequences from routine genotypic testing for resistance to common anti-HIV drugs. The related laboratory testing has been described previously [56,57]. Donor (LNA819) was diagnosed with HIV-1 infection on 3 March 2010 and identified as having HIV-1 superinfection by next-generation sequencing [27]. Due to ART initiation or loss to follow-up, serial plasma samples between 2010 and 2013 were collected from the donor. One plasma sample at baseline was collected from recipient 3 (LN328575), recipient 4 (LN301538), and recipient 5 (LN328576), respectively. One plasma sample at the first or second HIV-positivity time-point after seroconversion was collected from recipient 1 (LN320639) and recipient 2 (LN320392), respectively.
HIV-1 limiting-antigen avidity enzyme immunoassay (LAg-avidity EIA)
Plasma samples were tested for recent HIV infections with LAg-avidity EIA (Maxim Biomedicals, Rockville, MD, USA) according to manufacturer instructions. Normalized optical density (ODn) of 2.0 was used as a threshold cutoff to distinguish long-term HIV infection from recent HIV infection. Plasma with ODn upon initial screening >2 was classified as “long-term infection”, whereas that with ODn ≤2 was retested in triplicate for confirmation. Plasma with median ODn >1.5 was classified as “long-term HIV infection”, whereas that with ODn >0.4 but ≤1.5 was classified as “recent seroconversion”; for ODn ≤0.4, a serology confirmation test was necessary to further ensure that the plasma sample was HIV-positive.
Amplification and sequencing of near-full-length genomes (NFLGs) and 3′ half-genomes
NFLGs and 3’ half-genomes were amplified and sequenced directly according to methods described previously [58]. In brief, HIV-1 RNA was extracted from 140-μL plasma sample using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA) and was transcribed into complementary DNA using Superscript III Reverse Transcriptase (Life Technologies, Carlsbad, CA) with primer 07Rev8 (5′-CCTARTGGGATGTGTACTTCTGAA CTT-3′, nt 5193–5219) and 1.R3.B3R (5′-ACTACTTGAAGCACTCAAGGCAAGCTTTATTG-3′, nt 9611–9642) for 5′ half-genome and 3′ half-genome, respectively. The sequences of the 5′ half-genome and 3′ half-genome were amplified by nested polymerase chain reaction (PCR) by single-genome amplification (SGA) using 1 unit of Platinum Taq DNA Polymerase High Fidelity (Life Technologies, Carlsbad, CA) according to manufacturer instructions. The primers used to amplify the 5’ half-genome in the first round of PCR was 172 (5′-ATCTCTAGCAGTGGCGCCCGAACAG-3′, nt 625–649) and 07Rev8, and in the second round of PCR was 174 (5′-CTCTCGACGCAGGACTCGGCTTGCT-3′, nt 683–707) and Rev11 (5′-ATCATCACCTGCCATCTGTTTTCCAT-3′, nt 5041–5066). The primers used to amplify the 3’ half-genome in the first round of PCR was 07For7 (5′-CAAATTAYAAAA ATTCAAAATTTTCGGGTTTATTACAG-3′, nt 4875–4912) and 2.R3.B6R (5′-TGA AGCACTCAAGGCAAGCTTTATTGAGGC-3′, nt 9607–9636), and in the second round of PCR was VIF1 (5′-GGGTTTATTACAGGGACAGCAGAG-3′, nt 4900–4923) and Low2C (5′-TGAGGCTTAAGCAGTGGGTTCC-3′, nt 9591–9612). The thermal cycling conditions for PCR were one cycle at 94°C for 2 min, followed by 35 cycles of 94°C for 15 s, 60°C for 30 s, 68°C for 5 min, followed by one cycle of extension at 68°C for 10 min. Positive amplification products were purified and sequenced directly using internal walking primers by BGI (Beijing, China).
Phylogenetic and recombination analyses
Sequences were assembled with Sequencher 5.4.6 (Gene Codes, Ann Arbor, MI, USA), and aligned using Gene Cutter within HIV databases (www.hiv.lanl.gov), then adjusted manually with BioEdit 7.0 (www.mbio.ncsu.edu/ BioEdit) [59]. All sequences obtained in this study were blasted in the local sequence library by our research team. We used BLAST within HIV databases (www.hiv.lanl.gov) to eliminate potential cross-contamination during the experiment. Reference sequences were downloaded from the Los Alamos HIV Database (www.hiv.lanl.gov). Maximum likelihood phylogenetic trees (ML trees) of the aligned NFLG, 3′ half-genomes, and concatenated CRF01_AE and CRF07_BC segments were constructed by Fast Tree [60] and edited by Fig Tree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree). Recombination analyses were first done with the Recombinant Identification Program (RIP) [61] and jumping profile Hidden Markov Model (jpHMM) within HIV databases (www.hiv.lanl.gov) [62]. Further confirmation was achieved with bootscanning in Simplot 3.5.1 [63] to define the recombination structures (window size: 350 nt; step size: 50 nt; bootstrap replicate: 250).
Bayesian Markov Chain Monte Carlo (MCMC) evolutionary analyses
We wished to explore the phylogenetic relationship and the time of the most recent common ancestor (tMRCA) of viruses from the six participants, Bayesian phylogenetic analyses were done using the MCMC inference implemented in BEAST v2.5.1 [64]. For concatenated CRF01_AE segments, strict molecular-clock analyses were undertaken under the model of general time reversible (GTR) + I+G nucleotide substitution. For concatenated CRF07_BC segments, relaxed molecular-clock analyses were undertaken under the model of Tamura Nei 93 (TN93) nucleotide substitution. The MCMC chains were run 200-million times and sampled every 20,000 steps. The output was tested for convergence using Tracer v1.6, and related parameters were estimated from an Effective Sample Size (ESS) more than 200. Phylogenetic trees were summarized using TreeAnnotator (with 10% burn-in) and then edited using Fig Tree v1.4.2.
Supporting information
S1 Fig. Highlighter plots of HIV-1 3′ half-genome diversity in six HIV-1-infected MSM.
The initial strain (CRF01_AE) and superinfected strain (CRF07_BC) from the donor were chosen as master sequences and are colored light-coral and slate-blue, respectively. The x-axis represents the base number. The y-axis represents the sampling dates of donor or recipients. The 3’ half-genome sequences obtained from the donor and five recipients are shown in panel A and B, respectively. Some recombination key sites were marked with ↑.
https://doi.org/10.1371/journal.ppat.1009258.s001
(TIF)
S2 Fig. Distribution of breakpoints across HIV-1 3’ half-genome sequences of all recombinant viruses (N = 88) from the donor and five recipients.
The initial strain (CRF01_AE) and superinfected strain (CRF07_BC) were chosen as parental strains. SimPlot v3.5 (window size = 300 nt; step size = 10 nt) was used to identify breakpoint locations. In general, 13 breakpoints were identified on the HIV-1 3’ half-genome and labeled by a to m, respectively. The x-axis represents the range of breakpoints. The y-axis represents the frequency of all recombinant viruses at the corresponding breakpoints. Distribution of 13 breakpoints on the HIV-1 3’ half-genome (HXB2: 4950–9589) (schematic) are below.
https://doi.org/10.1371/journal.ppat.1009258.s002
(TIF)
S3 Fig. Maximum likelihood (ML) trees of HIV-1 CRF01_AE/CRF07_BC recombinants from the donor and five recipients.
The same sequences as in the Bayesian Evolutionary Analysis Sampling Trees (BEAST) analysis (Fig 4) were chosen to construct ML trees. ML trees for concatenated CRF01_AE segments (regions V+VII + IX) (A) and CRF07_BC segments (regions IV+VI) (B) were constructed under the model of substitution of GTR+I+G sites.
https://doi.org/10.1371/journal.ppat.1009258.s003
(TIF)
Acknowledgments
We are grateful to the individuals who consented to participate and donate their blood samples for our study. We thank all the researchers for their contributions to the design and initiation of this study and manuscript preparation.
References
- 1. Ramirez BC, Simon-Loriere E, Galetto R, Negroni M. Implications of recombination for HIV diversity. Virus research. 2008;134(1–2):64–73. pmid:18308413.
- 2. Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, et al. HIV-1 nomenclature proposal. Science. 2000;288(5463):55–6. Epub 2000/04/15. pmid:10766634.
- 3. Hemelaar J. The origin and diversity of the HIV-1 pandemic. Trends Mol Med. 2012;18(3):182–92. Epub 2012/01/14. pmid:22240486.
- 4. Hemelaar J, Elangovan R, Yun J, Dickson-Tetteh L, Fleminger I, Kirtley S, et al. Global and regional molecular epidemiology of HIV-1, 1990–2015: a systematic review, global survey, and trend analysis. Lancet Infect Dis. 2019;19(2):143–55. pmid:30509777.
- 5. Gu Z, Quan Y, Li Z, Arts EJ, Wainberg MA. Effects of non-nucleoside inhibitors of human immunodeficiency virus type 1 in cell-free recombinant reverse transcriptase assays. J Biol Chem. 1995;270(52):31046–51. Epub 1995/12/29. pmid:8537362.
- 6. Konings FA, Burda ST, Urbanski MM, Zhong P, Nadas A, Nyambi PN. Human immunodeficiency virus type 1 (HIV-1) circulating recombinant form 02_AG (CRF02_AG) has a higher in vitro replicative capacity than its parental subtypes A and G. J Med Virol. 2006;78(5):523–34. Epub 2006/03/24. pmid:16555291.
- 7. Nora T, Charpentier C, Tenaillon O, Hoede C, Clavel F, Hance AJ. Contribution of recombination to the evolution of human immunodeficiency viruses expressing resistance to antiretroviral treatment. Journal of virology. 2007;81(14):7620–8. pmid:17494080.
- 8. Turk G, Carobene MG. Deciphering How HIV-1 Intersubtype Recombination Shapes Viral Fitness and Disease Progression. EBioMedicine. 2015;2(3):188–9. pmid:26137559.
- 9. Liu Y, Su B, Zhang Y, Jia L, Li H, Li Z, et al. Brief Report: Onward Transmission of Multiple HIV-1 Unique Recombinant Forms Among Men Who Have Sex With Men in Beijing, China. J Acquir Immune Defic Syndr. 2019;81(1):1–4. Epub 2019/02/16. pmid:30768488.
- 10. McClutchan FE, Carr JK, Murphy D, Piyasirisilp S, Gao F, Hahn B, et al. Precise mapping of recombination breakpoints suggests a common parent of two BC recombinant HIV type 1 strains circulating in China. AIDS research and human retroviruses. 2002;18(15):1135–40. Epub 2002/10/31. pmid:12402948.
- 11. Feng Y, Takebe Y, Wei H, He X, Hsi JH, Li Z, et al. Geographic origin and evolutionary history of China’s two predominant HIV-1 circulating recombinant forms, CRF07_BC and CRF08_BC. Sci Rep. 2016;6:19279. Epub 2016/01/15. pmid:26763952; PubMed Central PMCID: PMC4725877.
- 12. Thomson MM, Delgado E, Herrero I, Villahermosa ML, Vázquez-de Parga E, Cuevas MT, et al. Diversity of mosaic structures and common ancestry of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Argentina revealed by analysis of near full-length genome sequences. The Journal of general virology. 2002;83(Pt 1):107–19. pmid:11752707.
- 13. Li Y, Tee KK, Liao H, Hase S, Uenishi R, Li XJ, et al. Identification of a novel second-generation circulating recombinant form (CRF48_01B) in Malaysia: a descendant of the previously identified CRF33_01B. J Acquir Immune Defic Syndr. 2010;54(2):129–36. Epub 2010/04/14. pmid:20386110.
- 14. Cheong HT, Chow WZ, Takebe Y, Chook JB, Chan KG, Al-Darraji HA, et al. Genetic Characterization of a Novel HIV-1 Circulating Recombinant Form (CRF74_01B) Identified among Intravenous Drug Users in Malaysia: Recombination History and Phylogenetic Linkage with Previously Defined Recombinant Lineages. PLoS One. 2015;10(7):e0133883. Epub 2015/07/22. pmid:26196131; PubMed Central PMCID: PMC4510129.
- 15. Galetto R, Negroni M. Mechanistic features of recombination in HIV. Aids Rev. 2005;7(2). pmid:16092503.
- 16. Galetto R, Moumen A, Giacomoni V, Véron M, Charneau P, Negroni M. The structure of HIV-1 genomic RNA in the gp120 gene determines a recombination hot spot in vivo. The Journal of biological chemistry. 2004;279(35):36625–32. pmid:15218022.
- 17. Jia L, Li L, Gui T, Liu S, Li H, Han J, et al. Analysis of HIV-1 intersubtype recombination breakpoints suggests region with high pairing probability may be a more fundamental factor than sequence similarity affecting HIV-1 recombination. Virology journal. 2016;13(1):156. pmid:27655081.
- 18. Magiorkinis G, Paraskevis D, Vandamme AM, Magiorkinis E, Sypsa V, Hatzakis A. In vivo characteristics of human immunodeficiency virus type 1 intersubtype recombination: determination of hot spots and correlation with sequence similarity. J Gen Virol. 2003;84(Pt 10):2715–22. Epub 2003/09/19. pmid:13679605.
- 19. Delviks-Frankenberry K, Galli A, Nikolaitchik O, Mens H, Pathak VK, Hu W-S. Mechanisms and factors that influence high frequency retroviral recombination. Viruses. 2011;3(9):1650–80. pmid:21994801.
- 20. Fang G, Weiser B, Kuiken C, Philpott SM, Rowland-Jones S, Plummer F, et al. Recombination following superinfection by HIV-1. AIDS. 2004;18(2):153–9. Epub 2004/04/13. pmid:15075531.
- 21. Steain MC, Wang B, Dwyer DE, Saksena NK. HIV-1 co-infection, superinfection and recombination. Sex Health. 2004;1(4):239–50. pmid:16335754.
- 22. Han X, An M, Zhao B, Duan S, Yang S, Xu J, et al. High prevalence of HIV-1 intersubtype B’/C recombinants among injecting drug users in Dehong, China. PLoS One. 2013;8(5):e65337. Epub 2013/06/07. pmid:23741489; PubMed Central PMCID: PMC3669332.
- 23. Wei H, Xing H, Hsi JH, Jia M, Feng Y, Duan S, et al. The sexually driven epidemic in youths in China’s southwestern border region was caused by dynamic emerging multiple recombinant HIV-1 strains. Sci Rep. 2015;5:11323. pmid:26133091.
- 24. Vidal N, Diop H, Montavon C, Butel C, Bosch S, Ngole EM, et al. A novel multiregion hybridization assay reveals high frequency of dual inter-subtype infections among HIV-positive individuals in Cameroon, West Central Africa. Infect Genet Evol. 2013;14:73–82. Epub 2012/12/13. pmid:23232100.
- 25. Koning FA, Badhan A, Shaw S, Fisher M, Mbisa JL, Cane PA. Dynamics of HIV type 1 recombination following superinfection. AIDS research and human retroviruses. 2013;29(6):963–70. pmid:23495713.
- 26. McCutchan FE, Hoelscher M, Tovanabutra S, Piyasirisilp S, Sanders-Buell E, Ramos G, et al. In-depth analysis of a heterosexually acquired human immunodeficiency virus type 1 superinfection: evolution, temporal fluctuation, and intercompartment dynamics from the seronegative window period through 30 months postinfection. Journal of virology. 2005;79(18):11693–704. pmid:16140747.
- 27. Luan H, Han X, Yu X, An M, Zhang H, Zhao B, et al. Dual Infection Contributes to Rapid Disease Progression in Men Who Have Sex With Men in China. J Acquir Immune Defic Syndr. 2017;75(4):480–7. pmid:28490044.
- 28. Beyrer C, Baral SD, van Griensven F, Goodreau SM, Chariyalertsak S, Wirtz AL, et al. Global epidemiology of HIV infection in men who have sex with men. Lancet. 2012;380(9839):367–77. Epub 2012/07/24. pmid:22819660; PubMed Central PMCID: PMC3805037.
- 29. Xu J-J, Zhang M, Brown K, Reilly K, Wang H, Hu Q, et al. Syphilis and HIV seroconversion among a 12-month prospective cohort of men who have sex with men in Shenyang, China. Sex Transm Dis. 2010;37(7):432–9. pmid:20375928.
- 30. Cao B, Liu C, Stein G, Tang W, Best J, Zhang Y, et al. Faster and Riskier? Online Context of Sex Seeking Among Men Who Have Sex With Men in China. Sex Transm Dis. 2017;44(4):239–44. Epub 2017/03/11. pmid:28282651; PubMed Central PMCID: PMC5347461.
- 31. Hong H, Xu J, McGoogan J, Dong H, Xu G, Wu Z. Relationship between the use of gay mobile phone applications and HIV infection among men who have sex with men in Ningbo, China: a cross-sectional study. International journal of STD & AIDS. 2018;29(5):491–7. pmid:29099328.
- 32. Lewnard JA, Berrang-Ford L. Internet-based partner selection and risk for unprotected anal intercourse in sexual encounters among men who have sex with men: a meta-analysis of observational studies. Sex Transm Infect. 2014;90(4):290–6. pmid:24518249.
- 33. Liu Y, Wang J, Qian HZ, Liu H, Yin L, Lu H, et al. Seeking Male Sexual Partners via Internet and Traditional Venues among Chinese Men Who Have Sex with Men: Implications for HIV Risk Reduction Interventions. AIDS Behav. 2016;20(10):2222–30. Epub 2016/03/24. pmid:27000143.
- 34. Rajasingham R, Mimiaga MJ, White JM, Pinkston MM, Baden RP, Mitty JA. A systematic review of behavioral and treatment outcome studies among HIV-infected men who have sex with men who abuse crystal methamphetamine. AIDS Patient Care STDS pmid:22070609 2012;26(1):36–52.
- 35. Salazar-Gonzalez JF, Bailes E, Pham KT, Salazar MG, Guffey MB, Keele BF, et al. Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing. Journal of virology. 2008;82(8):3952–70. pmid:18256145.
- 36. Edwards CTT, Holmes EC, Wilson DJ, Viscidi RP, Abrams EJ, Phillips RE, et al. Population genetic estimation of the loss of genetic diversity during horizontal transmission of HIV-1. BMC Evol Biol. 2006;6:10. WOS:000236971800001. pmid:16464261
- 37. Wolinsky SM, Wike CM, Korber BT, Hutto C, Parks WP, Rosenblum LL, et al. Selective transmission of human immunodeficiency virus type-1 variants from mothers to infants. Science. 1992;255(5048):1134–7. pmid:1546316.
- 38. Leitner T. Phylogenetics in HIV transmission: taking within-host diversity into account. Current opinion in HIV and AIDS. 2019;14(3):181–7. pmid:30920395.
- 39. Romero-Severson EO, Bulla I, Leitner T. Phylogenetically resolving epidemiologic linkage. Proc Natl Acad Sci USA. 2016;113(10):2690–5. pmid:26903617.
- 40. Tovanabutra S, Watanaveeradej V, Viputtikul K, De Souza M, Razak MH, Suriyanon V, et al. A new circulating recombinant form, CRF15_01B, reinforces the linkage between IDU and heterosexual epidemics in Thailand. AIDS research and human retroviruses. 2003;19(7):561–7. pmid:12908933.
- 41. Han X, An M, Zhang W, Cai W, Chen X, Takebe Y, et al. Genome Sequences of a Novel HIV-1 Circulating Recombinant Form, CRF55_01B, Identified in China. Genome Announc. 2013;1(1). Epub 2013/02/14. pmid:23405298; PubMed Central PMCID: PMC3569284.
- 42. Zhang W, Han X, An M, Zhao B, Hu Q, Chu Z, et al. Identification and characterization of a novel HIV-1 circulating recombinant form (CRF59_01B) identified among men-who-have-sex-with-men in China. PLoS One. 2014;9(6):e99693. Epub 2014/07/01. pmid:24978029; PubMed Central PMCID: PMC4076182.
- 43. Wu J, Meng Z, Xu J, Lei Y, Jin L, Zhong P, et al. New emerging recombinant HIV-1 strains and close transmission linkage of HIV-1 strains in the Chinese MSM population indicate a new epidemic risk. PLoS One. 2013;8(1):e54322. Epub 2013/02/02. pmid:23372706; PubMed Central PMCID: PMC3553145.
- 44. Zhang C, Feng Y, Gao L, Zhang M, Miao J, Dong X, et al. Genetic characterization and recombinant history of a novel HIV-1 circulating recombinant form (CRF101_01B) identified in Yunnan, China. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases. 2019;73:109–12. pmid:31028881.
- 45. Han X, Takebe Y, Zhang W, An M, Zhao B, Hu Q, et al. A Large-scale Survey of CRF55_01B from Men-Who-Have-Sex-with-Men in China: implying the Evolutionary History and Public Health Impact. Sci Rep. 2015;5:18147. Epub 2015/12/17. pmid:26667846; PubMed Central PMCID: PMC4678862.
- 46. Han X, An M, Zhang M, Zhao B, Wu H, Liang S, et al. Identification of 3 distinct HIV-1 founding strains responsible for expanding epidemic among men who have sex with men in 9 Chinese cities. J Acquir Immune Defic Syndr. 2013;64(1):16–24. Epub 2013/04/02. pmid:23542640; PubMed Central PMCID: PMC3814940.
- 47. Guo H, Guo D, Wei J-F, Yang H, Huan X, Tsui SK-W, et al. First detection of a novel HIV Type 1 CRF01_AE/07_BC recombinant among an epidemiologically linked cohort of IDUs in Jiangsu, China. AIDS research and human retroviruses. 2009;25(4):463–7. pmid:19320602.
- 48. Li X, Ning C, He X, Yang Y, Xing H, Hong K, et al. Near full-length genome sequence of a novel HIV type 1 second-generation recombinant form (CRF01_AE/CRF07_BC) identified among men who have sex with men in Jilin, China. AIDS research and human retroviruses. 2013;29(12):1604–8. pmid:23809010.
- 49. Li F, Li Y, Feng Y, Hu J, Ruan Y, Xing H, et al. Four Closely Related HIV-1 CRF01_AE/CRF07_BC Recombinant Forms Identified in East China. AIDS research and human retroviruses. 2017;33(7):740–4. Epub 2017/03/17. pmid:28298138.
- 50. Sierra M, Thomson MM, Rios M, Casado G, Castro RO, Delgado E, et al. The analysis of near full-length genome sequences of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Chile, Venezuela and Spain reveals their relationship to diverse lineages of recombinant viruses related to CRF12_BF. Infect Genet Evol. 2005;5(3):209–17. Epub 2005/03/02. pmid:15737911.
- 51. AIDS-Professional-Group. AIDS Professional Group, Society of Infectious Diseases Chinese Medical Association: Guidelines for diagnosis and treatment of HIV/AIDS in China. ZHONGHUA CHUAN RAN BING ZA ZHI. 2011;29(10):629–40.
- 52. AIDS-Professional-Group. AIDS Professional Group, Society of Infectious Diseases, Chinese Medical Association: third edition of the guidelines for diagnosis and treatment of HIV/AIDS. Chin J Clin Infect Dis. 2015;8(5):385–401.
- 53. China-CDC. Chinese Medical Association, China CDC: Guidelines for diagnosis and treatment of HIV/AIDS in China. Chin Med J. 2005;119(19):1589–608.
- 54. China-CDC. AIDS and Hepatitis C Professional Group, Society of Infectious Diseases Chinese Medical Association, China CDC: [Chinese guidelines for diagnosis and treatment of HIV/AIDS]. Zhonghua Nei Ke Za Zhi. 2018;57(12):867–84. pmid:30486555
- 55. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365(6):493–505. Epub 2011/07/20. pmid:21767103; PubMed Central PMCID: PMC3200068.
- 56. Chen X, Lin M, Qian S, Zhang Z, Fu Y, Xu J, et al. The Early Antibody-Dependent Cell-Mediated Cytotoxicity Response Is Associated With Lower Viral Set Point in Individuals With Primary HIV Infection. Front Immunol. 2018;9:2322. pmid:30356637.
- 57. Xu J-J, Qian H-Z, Chu Z-X, Zhang J, Hu Q-H, Jiang Y-J, et al. Recreational drug use among Chinese men who have sex with men: a risky combination with unprotected sex for acquiring HIV infection. Biomed Res Int. 2014;2014:725361. pmid:24829916.
- 58. Sanchez AM, DeMarco CT, Hora B, Keinonen S, Chen Y, Brinkley C, et al. Development of a contemporary globally diverse HIV viral panel by the EQAPOL program. J Immunol Methods. 2014;409:117–30. pmid:24447533.
- 59. Tippmann H-F. Analysis for free: comparing programs for sequence analysis. Brief Bioinformatics. 2004;5(1):82–7. pmid:15153308.
- 60. Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PloS one. 2010;5(3):e9490. pmid:20224823.
- 61. Siepel AC, Halpern AL, Macken C, Korber BT. A computer program designed to screen rapidly for HIV type 1 intersubtype recombinant sequences. AIDS research and human retroviruses. 1995;11(11):1413–6. Epub 1995/11/01. pmid:8573400.
- 62. Zhang M, Schultz AK, Calef C, Kuiken C, Leitner T, Korber B, et al. jpHMM at GOBICS: a web server to detect genomic recombinations in HIV-1. Nucleic Acids Res. 2006;34(Web Server issue):W463–5. Epub 2006/07/18. pmid:16845050; PubMed Central PMCID: PMC1538796.
- 63. Salminen MO, Carr JK, Burke DS, McCutchan FE. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS research and human retroviruses. 1995;11(11):1423–5. pmid:8573403.
- 64. Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5. pmid:11524383.