Figures
Abstract
Background
We have previously reported in Xishuangbanna (Banna) Dai Autonomous Prefecture, a well-developed tourist destination in the southwest border of China, that HIV-1 transmitted dominantly through heterosexual contact with less divergent genotypes and few drug resistant mutations [1]. Due to the rapid increase of newly diagnosed HIV-1 cases per year in Banna in recent years, it’s important to evaluate the evolution of HIV-1 molecular epidemiology for the better understanding of ongoing HIV-1 outbreak in this region.
Methodology/Principal Findings
By sequencing of HIV-1 pol genes and phylogenetic analysis, we conducted a molecular epidemiologic study in 352 HIV-1-seropositive highly active antiretroviral treatment (HAART)-naïve individuals newly diagnosed at the Banna Center for Disease Control and Prevention between 2009 and 2011. Of 283 samples (84.1% taken from heterosexually acquired adults, 10.6% from needle-sharing drug users, 2.8% from men who have sex with men, 0.4% from children born from HIV-1-infected mothers, and 2.1% remained unknown) with successful sequencing for pol gene, we identified 108 (38.2%) HIV-1 subtype CRF08_BC, 101 (35.7%) CRF01_AE, 49 (17.3%) CRF07_BC, 5 (1.8%) C/CRF57_BC, 3 (1.1%) B’, 1 (0.4%) B/CRF51_01B, and 16 (5.7%) unique recombinants forms. Among these infected individuals, 104 (36.7%) cases showed drug resistant or resistance-relevant mutations, and 4 of them conferring high-level resistance to 3TC/FTC, EFV/NVP or NFV. Phylogenetic analysis revealed 21 clusters (2–7 sequences) with only 21.2% (60/283) sequences involved.
Citation: Chen Y, Chen S, Kang J, Fang H, Dao H, Guo W, et al. (2014) Evolving Molecular Epidemiological Profile of Human Immunodeficiency Virus 1 in the Southwest Border of China. PLoS ONE 9(9): e107578. https://doi.org/10.1371/journal.pone.0107578
Editor: Dimitrios Paraskevis, University of Athens, Medical School, Greece
Received: April 16, 2014; Accepted: August 19, 2014; Published: September 10, 2014
Copyright: © 2014 Chen 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.
Funding: This work was supported by the Institut de Recherche pour le Développement (IRD) and the Institut de Recherche sur les Vaccins et l’Immunothérapie des Cancers et du Sida (IRVICS). 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
The Xishuangbanna (Banna) Dai Autonomous Prefecture of China’s Yunnan province, geographically located at the southwestern border of China, has two counties (Menghai and Mengla) and one city (Jinghong) with an estimated population of 1.13 million populations (including 13 ethnics) in 2012. Due to its distinctive biogeographic characteristics and diversified ethnic cultures, Banna prefecture has become a well-developed tourist destination in China since 1990s.
Before 2005, HIV-1 transmission in Banna prefecture had been characterized as mainly through heterosexual contact (77.3%) with less genotypic diversity and uncommon drug resistance mutations (DRMs) [1]. Although multiple HIV-1 genotypes (including B, B’, C, CRF07_BC, CRF08_BC, and CRF01_AE) had been identified in Yunnan province [2], only CRF01_AE and CRF08_BC had been found to be dominant genotypes (62.2% and 33.3% respectively). Moreover, only one case with subtype CRF01_AE had been identified to carry mutations conferring high-level resistance to antiretroviral drugs NRTI and NNRTI in our previous study [1].
Due to the rapidly increasing numbers of newly diagnosed HIV-1 cases per year in the past years (75 cases in 2005 and nearly 200 cases in 2011), it is important to evaluate the evolution of HIV-1 molecular epidemiology for the better understanding of ongoing HIV-1 outbreak in this region. In addition, the free HAART drugs have been delivered to registered HIV-1-infected individuals with a CD4 count less than 350 cells/µl according to the national “free antiretroviral treatment program” (NFATP) since 2005. Therefore, the transmission of HIV-1 with drug-related mutations might be emerged, as we observed recently in Guangdong province [3]. In the present study, we conducted a molecular epidemiological investigation in 357 cases of newly diagnosed HIV-1 infection from July 2009 to June 2011 at the Banna Center for Disease Control and Prevention (CDC).
Materials and Methods
Participants and specimens
From July 2009 to June 2011, a total of 352 individuals who had been newly diagnosed with HIV-1 infection at the Banna CDC had been enrolled into this study. All individuals were required to complete standardized questionnaires (describing sex, age, risk factors, mode of transmission, occupation, geographic location). An enzyme immunoassay for screening HIV-1/HIV-2/O antibodies (Vironostika Uni-form II plus O kit; Organon Teknika BV, Turnhout, Belgium) in sera was performed in the Banna CDC. The positive sera were confirmed by western blot (HIVBlot2.2 kit; Genelabs Diagnostics, Singapore).
After giving their written informed consents, 322 HIV-1 infected individuals were recruited into the molecular epidemiological and drug resistance survey. Briefly, a total of 5 ml of EDTA-treated whole blood was taken from each patient, and plasma samples were sent to our laboratory of Tropical Medicine Institute, Guangzhou University of Chinese Medicine (TMI, GUCM) by flight within 4 hours. The blood samples were used immediately for routine blood count and CD4 T-cell count measurements as well as to separate plasma. Plasma samples were stored at −80°C until use.
Ethics statement
The institutional ethics committees of TMI, GUCM had approved the study protocol (No. 2009C015).
Viral load measurement
All plasma samples were thawed at the same time and were used for viral RNA measurement using the Food and Drug Administration (FDA) -approved Amplicor HIV-1 Monitor Test kit (version 1.5) (Roche Molecular Systems, Inc., Branchburg, New Jersey, USA) according to the manufacturer’s instructions.
Sequencing of HIV-1 pol gene
The pol gene of HIV-1 encodes the viral enzymes including protease, reverse transcriptase, and integrase, and holds sufficient variability to permit the phylogenetic reconstruction of transmissions [4]. Although it is hard to identify the difference between subtypes CRF01_AE and CRF15_01B by sequencing pol region, CRF15_01B so far has not yet been identified in Yunnan province [5] and few cases with CRF15_01B have been reported only in Hebei province and Beijing city [6]. In addition, the newly identified circulating recombinant forms, CRF57_BC [7]/CRF65_cpx [8] and CRF51_01B [9], could not be distinguished from subtypes B and C by sequencing pol region. However, considering that subtypes B and C had not-yet been identified in this area, we conducted the HIV-1 genotyping by sequencing the pol region as described previously [1]. Briefly, viral RNA was extracted from the HIV-1-infected individual’s plasma (150 µl) using the QIAamp Viral RNA Mini kit (Qiagen, Valencia, California, United States) according to the manufacturer’s instructions initially. The viral RNA was then subjected to a one-step reverse transcription polymerase chain reaction (RT-PCR) to generate a fragment of pol gene (1864 base pairs) spanning protease and reverse transcriptase regions as previously described [1], [3], [10]. The PCR products were purified (Qiagen, Valencia, Spain) and directly sequenced. The sequences generated were edited using the SeqMan II software program from the DNAStar package v.5.08 (Lasergene, Madison, WI). To eliminate potential contamination, all of the sequences obtained were first subjected to an HIV-1 Blast search to compare with related reference sequences in the HIV Databases, funded by the Division of AIDS of the National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH) (http://hiv-web.lanl.gov/content/index). Finally, a total of 283 sequences were successfully obtained from the 322 blood samples.
Drug resistant mutations
Drug resistance mutations in protease and reverse transcriptase genes were identified by using the last updated online Stanford Resistance Database tool: HIVdb program-Genotypic Resistance Interpretation Algorithm (version 7.0, http://sierra2.stanford.edu/sierra/servlet/JSierra?action=sequenceInput). In additional, sequence data were also submitted and evaluated by using the last updated list for surveillance of transmitted drug-resistant strains in untreated patients (ver.6.0 http://cpr.stanford.edu/cpr.cgi). Resistance mutations identified with HIVdb algorithm were also interpreted referring to the “Update of the drug resistance mutations in HIV-1: March 2013” published by IAS–USA Drug Resistance Mutations Group.
HIV-1 genotyping and phylogenetic analysis
HIV-1 subtype and CRF designations were determined by uploading sequences into the REGA HIV-1 automated Subtyping Tool version 2.0 (http://www.bioafrica.net/rega- genotype/html/subtypinghiv.html), and confirmed by in-house phylogenetic analysis on nucleotide acid sequences of pol as previously described [1], [3], [10]. For phylogenetic analysis, reference sequences representing overall HIV-1 group M genetic variability were obtained from the National Institutes of Health/National Institute of Allergy and Infectious Diseases (NIH/NIAID)-funded HIV database, including all subtypes, sub-subtypes, and circulating recombinant forms (CRFs) references. CRFs identified most recently, such as CRF62_BC [11], CRF64_BC [12], and CRF65_CPX [8], etc., were also included. All reference sequences and our newly obtained nucleotide sequences were aligned using Muscle [13] and followed by manual editing. Neighbor-joining (N-J) tree was drawn under the HKY model of evolution with 1000 bootstrap replicates in SeaView v4.3 [14]. Sequences assigned to specific HIV-1 subtypes were finally confirmed by constructing maximum likelihood (ML) phylogenetic sub-trees as described below. To determine a recombinant virus, similarity analysis and bootscanning were performed with the Simplot version 3.5.1 software [15].
Phylogenetic inter-relationships among viral sequences were estimated using maximum likelihood (ML) phylogenies with PHYML 3.0 [16] with the minimal number of reference sequences. The whole sequence alignment was split into 6 sub-datasets based on the neighbor-joining tree, followed by gap-stripping and ML phylogenetic reconstruction. An approximate likelihood ratio test (SH-like) was used to assess confidence in topology. Phylogenies were inferred using a general-time reversible model of nucleotide substitution, an estimated proportion of invariant sites, and gamma distributed rates among sites. The best of SPR and NNI heuristic options was selected to search the tree space, and bootstrap values with 1000 replicates were used to assess confidence in topology. The existence of transmission clusters was determined using the statistical robustness of the maximum likelihood topologies assessed by high bootstrap values (>98%) with 1000 re-samplings and short branch lengths (genetic distances<0.015) of HIV-1 pol gene sequence [3], [4]. The phylogenetic tree was drawn with FigTree v.1.42 (tree.bio.ed.ac.uk/software/figtree/).
Results
Characteristics of HIV-1 participants in Banna prefecture
A total of 352 newly HIV-1 infected individuals had been diagnosed in Banna CDC from July 2009 to June 2011, including 345 natives inhabited in Jinghong city (221 cases, 62.8%), Menghai county (108 cases, 30.7%), Mengla county (16 cases, 4.5%) and immigrants (7 cases, 2.0%) (Figure 1A). The majority of them were men (224 cases, 63.6%), 20–49 years old (296 cases, 84.1%), poorly educated (illiterate or primary school) (155 cases, 44.3%), engaged in farmer (97 cases, 27.6%), unemployment or uncertain (284 cases, 80.7%), and infected mainly by heterosexual contact (294 cases, 83.5%) (Figure 1B). The means (±SD) of CD4 counts and of viral loads were 430±266 cells/ml and 3.76±1.35 log10 copies/ml respectively. All patients had not been previously treated with any antiretroviral drug.
The geographical location of Banna prefecture and geographical distribution of infections were presented (parentheses showed cases and percentage of local infected individuals enrolled in this region (A). The characteristics (including their age, occupation, risk factors and education) of all enrolled individuals and individuals with successful sequencing (B).
HIV-1 Genotyping and Unique Recombinant Forms (URF) identification
Among 322 infected individuals who had given their written consent to participate in this study, we obtained the complete sequencing from 283 samples, while 39 samples failed in amplification. However, the characteristics of the 283 individuals with a successful pol sequencing were comparable with that of the whole study group (352 infected individuals) (Figure 1B). Their means (±SD) of CD4 counts and of viral loads were 426±257 cells/ml and 3.84±1.32 log10 copies/ml respectively. Both values were similar to that of the whole study group (see above).
Of 283 samples with genotyping analysis, 108 (38.2%) CRF_08BC, 101 (35.7%) CRF_01AE, 49 (17.3%) CRF_07BC, 5 (1.8%) C/CRF57_BC, 3 (1.1%) B’, 1 (0.4%) B/CRF51_01B, and 16 (5.7%) URFs had been identified (Table 1 and Figure 2). As showed in Figure 2, sequences with ambiguous genotyping were highlighted in blue and red. Blue lines represented 3 sequences designated as B/CRF51_01B or C/CRF57_BC respectively. Red lines indicated 16 URFs (including 12 B’/C, 2 01_AE/B’ and 2 01_AE/C recombinants). Bootscan analysis revealed that these URFs possessed breakpoints differed apparently from any known reference sequences (Figure 3), even though we included all novel B’/C recombinants circulated in western China (such as CRF62_BC [11], CRF64_BC [12] and CRF65_CPX [8]). As the representative data showed in Figures 3A and 3B, 4 B’/C URF samples (including YN022, YN114, YN151 and YN153) within a cluster of phylogenetic tree, possessed a distinct multiple-mixed B’/C mosaic structure in pol region. The other 8 B’/C URFs carried similar structural pol region with CRF62_BC. Specifically, they presented a mosaic pol structure of C/B’/C with an about 600 bp B’ fragment in the middle (Figures 3C and 3D). Nonetheless, in our sequence alignment, the breakpoints of CRF62_BC located at positions 1034 and 1610, about 200 bp or 300 bp behind than breakpoints of YN139 or YN178 respectively. Similarly, the remaining URFs (2 01_AE/B’ and 2 CRF01_AE/C) also could not be defined as any known CRF based on their bootscan analysis (Figure 3E and 3F). Full genome sequencing is required for confirming any novel circulating recombinant form in these samples.
The sample sequences with ambiguous genotype were highlighted in blue and red, and labeled with numbers from 1 to 6 simultaneously. 1 and 2 indicated sequences can’t be well distinguished by this sequencing region, and were designated as B/CRF51_01B and C/CRF57_BC respectively. 3 to 6 indicated 16 clustered unique recombinant forms. Other sample sequences and reference sequences were showed in light brown and black respectively.
Simplot bootscan analyses were implemented based on sequences of pol region. YN151 (A) and YN022 (B) represented 4 mosaic B’/C recombinants close to CRF07_BC in phylogenetic tree. YN178 (C) and YN139 (D) served as examples of 8 novel B’/C recombinants carried breakpoints similar to those of CRF62_BC. YN527 (E) represented 2 undefined CRF01_AE/C recombinants and YN186 (F) represented 2 undefined CRF01_AE/B’ recombinants.
The whole alignment was split into 6 sub-datasets based on their genotypes, and phylogenetic reconstructions were implemented on each sub-dataset after gag stripping. A total of 21 transmission clusters were identified according to the criteria (bootstrap value>98% and branch lengths<0.015). Sequences composed these clusters were highlighted in blue (male) or red (female), and behind the sequence names, their transmission routes except for heterosexual contact were indicated with asterisk (*, injection drug use), pound sign (#, men who had sex with men) or question mark (?, uncertain) respectively. Unrelated sequences and reference sequences were showed in blown or black respectively.
Drug Resistance Mutation Analysis (DRMs)
As showed in Table 2, a total of 33 DRMs sites in 104/283 (36.7%) infected individuals were identified with HIVdb Algorithm in our study, including 6 DRMs identified by CPR ver.6.0 simultaneously. The 6 DRMs identified by both HIVdb Algorithm and CPR included M184V/I, T215S, K219N (NRTI), K103N (NNRTI), and D30N, M46I/L (PI). M184V is selected by 3TC/FTC and reduces susceptibility to these drugs more than 100-fold. M184I usually emerges before M184V and has similar resistance profiles. T215S is a known revertant of the resistance mutation T215Y/F. K103N is selected by NVP and EFV and reduces susceptibility to them by about 50 and 20-fold, respectively. D30N is a NFV-selected substrate-cleft mutation that causes high-level resistance to NFV. M46I/L are selected primarily by NFV, ATV and LPV etc., and reduce susceptibility to these drugs. Other mutations such as A98G, V179D/T and Q58E are also need to be considered. A98G and V179D reduces NVP and EFV susceptibility. Q58E is a major mutation for Tipranavir.
Identification of Transmission Clusters
To evaluate the global profile of the local HIV-1 transmission, we split the sequence alignment into 6 sub-datasets, and implemented maximum-likelihood phylogenetic reconstruction on each sub-dataset to obtain the genetic distances. These sub-datasets were designate with the data of a representative references, such as CRF08_BC (1467 base pairs remained after gap stripping), 07_BC (1666 base pairs), B/CRF51 (1711 base pairs), CRF01_AE (1194 base pairs), C/CRF57 (1752 base pairs) and B’/C URF (1708 base pairs). As screened by our criteria (bootstrap value>98% and branch lengths<0.015), 21.2% (60/283) of all sequenced samples analyzed segregated into 21 transmission clusters (Figure 4). Of those 21 transmission clusters identified, 17 chains had been supported by contact tracing (including the sex, age, habitation, occupation, risk infectors in the same cluster), while four clusters (clusters 3, 6, 8 and 9) remained ambiguous. Clusters 6 and 8 were both consisted of 2 women lived in the same city, who had reported to get HIV-1 from heterosexual contact. Both clusters 3 and 9 derived from 2 men lived in the same city, who claimed to contract HIV-1 either by heterosexual contact or by uncertain route.
Among the 21 clusters, most of them were composed of 2 sequences (12 clusters), while 5 clusters were composed of 3 sequences and 4 clusters were composed of multiple sequences (cluster 15 consisted of 4 sequences; both clusters 5 and cluster 13 consisted of 5 sequences; and the biggest cluster 1, consisted of 7 sequences). The new B’/C recombinant virus spreading in cluster 1 represented a potential regional outbreak of novel B’/C recombinant, as reported elsewhere in Yunnan province [17], [18]. The characteristics of patients in the 21 transmission clusters were summarized in Table 3.
We further investigated whether the transmission of drug resistant strains occurs among these 21 clusters. There were 6/21 (28.6%) clusters consisted of individuals who were infected with strains carrying the same DRMs. These cases included V179T in cluster 4, V90I in cluster 5, E138A in cluster 9, Q58E/T74S in cluster 14, V90I/V179D in cluster 15, and A98G in cluster 18.
Discussion
In Yunnan province, there were a total of 104981 cases of HIV-1 infection and 7671 newly diagnosed HIV-1 infected individuals in 2012 (http://yn.yunnan.cn/html/2012-11/27/content_2510358.htm). Although nearly 200 newly diagnosed HIV-1 infected individuals per year in Banna prefecture represents only 2% of new infections per year in Yunnan, our findings in the present study did raise serious concerns on the severity of regional HIV-1 transmission. Firstly, newly diagnosed infected individuals increased rapidly (averaged 18% per year deduced from cases in 2005 and 2011). Heterosexual transmission was consistently dominant in Banna. In addition, the unsafe sexual behavior of injecting drug users could further increased the risk of HIV-1 transmission in the general population [19].
Secondly, our results suggested that HIV-1 population evolved rapidly in Banna and the evidence of new sources of infection further increased the genetic complexity. The profiles of HIV-1 evolution in Yunnan province (including the emergence of numerous novel circulating recombinant forms) had been described by previous reports [5], [20], [21], [22]. Compared with our previous study in 2008 [1], CRF08_BC had replaced CRF01_AE to become the most prevalent genotype in this region. A shift of dominant HIV-1 subtype had also been described in other regions of Yunnan province [5], [20], [21]. Moreover, we identified in the present study 6 common genotypes and 3 URF genotypes (undefined B’/C, CRF01_AE/C, and CRF01_B’), as compared to 3 common genotypes and only 1 URF reported in our previous study in the same region [1]. These results indicated the local HIV-1 genotypes divergence did increase rapidly in this region.
A series of novel circulating recombinant forms have recently been identified in China, especially in Yunnan province, including CRF55_01B (identified from MSM in China) [23], CRF57_BC [7], CRF59_01B (from MSM in northeast China) [24], CRF61_BC (CRF found among the heterosexual population in different regions in China) [25], CRF62_BC [11], CRF64_BC [12], CRF65_CPX (first novel HIV-1 second-generation inter-CRF in China) [8]. In the present study, we included available sequences of CRFs in HIV-1 sequence database into our reference pool. The breakpoints of our URFs differed apparently from those of any known CRF, implying the unidentified novel CRFs emerged in this region.
The evolving profile of HIV-1 molecular epidemiology raises also a challenge for the antiretroviral therapy. The diagnosis of DRMs associated with free HAART drugs (supplied under the NFATP program) had been increased rapidly in the past 5 years. The rapid genotyping evolution and transmission of HIV-1 with DRMs in Banna might increase the risk of treatment failure in this region. Of 283 sequences analysed, 102 (36.0%) showed DRMs, even if we excluded 6 DRMs from the IAS-USA 2013 mutation list (due to the coexistence of multiple DRMs). This proportion of untreated HIV-1 patients with DRMs was apparently higher than those reported previously (13.3%) [1]. Moreover, not only the rapidly increasing prevalence of resistant strains, but also the emerging of mutations conferring high level resistance to the primary antiretroviral drugs represent a serious challenge for the control of epidemic in this region. As in this study, there were 4 (1.4%) cases carrying mutations M184V (2 cases), K103N (1 case) and D30N (1case), which should confer high level resistance to 3TC/FTC, EFV/NVP and NFV, the primary drugs composed the first-line antiretroviral therapy regimen recommended in China. Moreover, a higher proportion (28.6%) of transmission clusters were also found to be associated with DRMs transmission. Poor adherence and compliance may be the most important contributors related to the increasing prevalence of DRMs in this region. The “HIV-1 Education and Training Programs” need to be enhanced in Banna to guarantee the infected individuals taking medicine as advised. Since as high as 54.8% of infected individuals were farmer and 25.9% were unemployed (Figure. 1b), such a poor education status would frustrate our efforts on prevention and control of HIV-1 dissemination. Thus, the efforts to include people in low social-status into the NFATP program will be the key to control HIV-1 outbreak in this region.
Acknowledgments
We acknowledge Liang Yu, Wengdi Deng, Véronique Jagot, Ju Yuan, Jun Zheng, and Yingyun Zhou for technical assistance. We thank the Laboratory of National Key Discipline of Clinical Foundation of Chinese Medicine and the “Program 211” of Guangdong province at the Guangzhou University of Chinese Medicine for providing some laboratory facilities used in this study. Thank the TH-1/GZ, GuangZhou Supercomputing Center to supply support on PhyML phylogenetic reconstruction computing.
Author Contributions
Conceived and designed the experiments: YC SC JK HF HD WG CL ML JF LF JMA WL. Performed the experiments: YC SC JK HF HD WG CL ML. Analyzed the data: YC SC WL. Contributed reagents/materials/analysis tools: YC SC WL. Wrote the paper: YC SC WL
References
- 1. Bao L, Vidal N, Fang H, Deng W, Chen S, et al. (2008) Molecular tracing of sexual HIV Type 1 transmission in the southwest border of China. AIDS Res Hum Retroviruses 24: 733–742.
- 2. Yu XF, Wang X, Mao P, Wang S, Li Z, et al. (2003) Characterization of HIV type 1 heterosexual transmission in Yunnan, China. AIDS Res Hum Retroviruses 19: 1051–1055.
- 3. Chen S, Cai W, He J, Vidal N, Lai C, et al. (2012) Molecular epidemiology of human immunodeficiency virus type 1 in Guangdong province of southern China. PLoS One 7: e48747.
- 4. Hue S, Clewley JP, Cane PA, Pillay D (2004) HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy. AIDS 18: 719–728.
- 5. Chen M, Yang L, Ma Y, Su Y, Yang C, et al. (2013) Emerging variability in HIV-1 genetics among recently infected individuals in Yunnan, China. PLoS One 8: e60101.
- 6. Chu XG, Zhan FX, Peng GP, Chen HP, Peng TH, et al. (2012) Subtype and sequence analysis of the gag genes for HIV-1 strains isolated in Hubei province (in Chinese). Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 26: 460–463.
- 7.
Wei H, Liu Y, Feng Y, Hsi JH, Xing H, et al.. (2013) Genome Sequence of a Novel HIV-1 Circulating Recombinant Form (CRF57_BC) Identified from Yunnan, China. AIDS Res Hum Retroviruses.
- 8.
Feng Y, Wei H, Hsi JH, Xing H, He X, et al.. (2013) Identification of a Novel HIV-1 Circulating Recombinant Form (CRF65_cpx) comprised of CRF01_AE and subtypes B and C in western Yunnan, China. AIDS Res Hum Retroviruses.
- 9. Ng OT, Eyzaguirre LM, Carr JK, Chew KK, Lin L, et al. (2012) Identification of new CRF51_01B in Singapore using full genome analysis of three HIV type 1 isolates. AIDS Res Hum Retroviruses 28: 527–530.
- 10. Deng W, Fu P, Bao L, Vidal N, He Q, et al. (2009) Molecular epidemiological tracing of HIV-1 outbreaks in Hainan island of southern China. AIDS 23: 977–985.
- 11.
Wei H, His J, Feng Y, Xing H, He X, et al.. (2013) Identification of a Novel HIV-1 Circulating Recombinant Form (CRF62_BC) in Western Yunnan of China. AIDS Res Hum Retroviruses.
- 12.
Hsi JH, Wei H, Xing H, Feng Y, He X, et al.. (2013) Genome Sequence of a Novel HIV-1 Circulating Recombinant Form (CRF64_BC) Identified from Yunnan, China. AIDS Res Hum Retroviruses.
- 13. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.
- 14. Gouy M, Guindon S, Gascuel O (2010) SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 27: 221–224.
- 15. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, et al. (1999) Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 73: 152–160.
- 16. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, et al. (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59: 307–321.
- 17. Chen L, Yang S, Li J, Liu Y, Wang Z, et al. (2011) Near full-length genomic characterization of a novel HIV type 1 CRF07_ BC/CRF08_ BC recombinant strain from Yunnan, China. AIDS Res Hum Retroviruses 27: 693–699.
- 18. Liu Y, Li L, Yang S, Bao Z, Li H, et al. (2011) Identification and characterization of two new HIV type 1 Unique (B/C) recombinant forms in China. AIDS Res Hum Retroviruses 27: 445–451.
- 19. Yao Y, Wang N, Chu J, Ding G, Jin X, et al. (2009) Sexual behavior and risks for HIV infection and transmission among male injecting drug users in Yunnan, China. Int J Infect Dis 13: 154–161.
- 20. Li L, Chen L, Yang S, Li T, Li J, et al. (2012) Recombination form and epidemiology of HIV-1 unique recombinant strains identified in Yunnan, China. PLoS One 7: e46777.
- 21. Chen M, Ma Y, Duan S, Xing H, Yao S, et al. (2012) Genetic diversity and drug resistance among newly diagnosed and antiretroviral treatment-naive HIV-infected individuals in western Yunnan: a hot area of viral recombination in China. BMC Infect Dis 12: 382.
- 22.
Su YZ, Ma YL, Jia MH, He X, Yang L, et al.. (2013) Update on diversity and distribution of HIV-1 subtypes in Yunnan province. Epidemiol Infect: 1–10.
- 23.
Han X, An M, Zhang W, Cai W, Chen X, et al.. (2013) Genome Sequences of a Novel HIV-1 Circulating Recombinant Form, CRF55_01B, Identified in China. Genome Announc 1.
- 24.
Han X, An M, Zhang W, Zhao B, Chu Z, et al.. (2013) Genome Sequences of a Novel HIV-1 Circulating Recombinant Form (CRF59_01B) Identified among Men Who Have Sex with Men in Northeastern China. Genome Announc 1.
- 25.
Li X, Ning C, He X, Yang Y, Li F, et al.. (2013) Genome Sequences of a Novel HIV-1 Circulating Recombinant Form (CRF61_BC) Identified among Heterosexuals in China. Genome Announc 1.