Skip to main content

Spatial clusters of HIV-1 genotypes in a recently infected population in Yunnan, China

Abstract

Background

As a gateway for HIV-1 in China, Yunnan has experienced dramatic changes in HIV-1 epidemics, during which HIV-1 genotypes have become complex. To track dynamic changes in HIV-1 genotypes, an HIV-1 molecular epidemiological study was implemented in the recently infected population in Yunnan.

Methods

From 6,357 HIV-1-positive samples diagnosed during the first half of 2015 in Yunnan, 586 samples were identified as recent infections with BED-capture enzyme immunoassay (CEIA) and were subjected to phylogenetic analyses. Spatial scanning analyses for the main HIV-1 genotypes were also performed.

Results

Among the 439 specimens successfully genotyped, more than ten genotypes were detected, including CRF08_BC (45.3%), CRF07_BC (19.4%), unique recombinant forms (URFs) (18.2%), CRF01_AE (11.4%), subtype C (2.1%), CRF85_BC (1.1%), CRF55_01B (0.9%), subtype B (0.5%), CRF64_BC (0.5%), CRF59_01B (0.2%), CRF83_cpx (0.2%) and CRF87_cpx (0.2%). Females, Chinese, heterosexual contact and intravenous drug injection were significantly associated with CRF08_BC infection; homosexual contact was significantly associated with CRF01_AE and CRF07_BC infection; males and non-Chinese had a higher risk of URF infection than females. Among all HIV-1 genotypes, the geographic coverage of CRF08_BC was the largest. For CRF08_BC, CRF07_BC, URFs and CRF01_AE, spatial clusters were detected. The two CRF08_BC clusters and one URF cluster were associated with heterosexual transmission, and two of CRF01_AE clusters were associated with homosexual transmission. Transmitted drug resistance (TDR)-associated mutations were detected in 2.4% of individuals.

Conclusions

The diversity of HIV-1 genotypes increased in recent infections because of a long-term HIV-1 epidemic in Yunnan. The predominant HIV-1 strains showed distinct demographic characteristics and formed spatial clusters. These findings improved our understanding of the evolution of HIV-1 in Yunnan and provided information for further HIV-1 control and prevention.

Peer Review reports

Background

Since HIV-1 was first identified among intravenous drug users (IDUs) in western Yunnan in 1989, HIV-1 epidemics have dramatically increased in Yunnan [1,2,3]. Yunnan is considered the HIV-1 gateway in China. Originating in Yunnan, different HIV-1 genotypes have spread to other parts of China [1, 4]. Over the past nearly 30 years, unprotected sexual contact has replaced intravenous drug injection as the predominant route of HIV-1 transmission [2]. By the end of 2014, the number of people living with HIV/AIDS (PLWHA) in Yunnan was 80,610 (including 50,426 people living with HIV and 30,184 AIDS patients), which was the highest among all provinces in China. Among all PLWHA in Yunnan, 68.9% received antiretroviral therapy, which was higher than the national average (63.2%). Among the annual newly reported cases, the proportion of cases attributed to intravenous drug use decreased from 28.6% in 2009 to 9.3% in 2014; the proportion of heterosexually infected cases increased from 62.5% in 2009 to 85.1% in 2014, while the proportion of homosexually infected cases increased from 1.4% in 2009 to 4.3% in 2014.

Initially, subtype B and subtype C emerged among IDUs [5,6,7], and CRF01_AE was found among female sex workers (FSWs) [8]. In the 1990s, CRF07_BC and CRF08_BC formed through the recombination of subtype B and subtype C, which were circulating among IDUs in Yunnan and transmitted along drug trafficking routes [9, 10]. Through the populations with the dual risk factors of intravenous drug use and heterosexual contact [11], CRF08_BC and CRF07_BC emerged in the heterosexually transmitted population [12]. In the late 2000s, CRF08_BC became the chief HIV-1 genotype among both heterosexuals and IDUs, and the composition of HIV-1 genotypes was nearly the same in these two populations [12]. The first HIV-positive man who had sex with men (MSM) in Yunnan was reported in 2004. The surveillance showed that CRF01_AE and CRF07_BC were the main strains among MSM, accounting for more than 90% of cases [13]. Recently, CRF55_01B and CRF59_01B were also detected in MSM in Yunnan [14], which were first identified among MSM in southern China and northeastern China, respectively [15, 16]. Moreover, the spatial distribution of HIV-1 genotypes dynamically changed. The main HIV-1 genotypes spread from the original areas to other areas of Yunnan. For example, CRF01_AE was transmitted from western to eastern Yunnan, and CRF08_BC and CRF07_BC were transmitted from southeast of Yunnan to other areas.

The genetic diversity of HIV-1 is striking in Yunnan. In addition to the abovementioned subtypes and circulating recombinant forms (CRFs), a multitude of URFs were also detected in Yunnan [12, 17]. In the late 2000s, the proportion of URFs ranked second only to that of CRF08_BC in recent infections [12]. These URFs contributed to the formation of novel CRFs. In recent years, at least nine CRFs were identified in western Yunnan [18,19,20,21,22,23,24,25]. These CRFs consist of components from subtype B, subtype C and CRF01_AE. A respective phylogenetic study showed that some of these CRFs originated in the early 1990s [9].

With the development of HIV-1 epidemics, the demographic profile of those infected has changed, while the characteristics and distribution of HIV-1 genotypes have also changed. To track the changing trend of HIV-1 genotypes in Yunnan, we performed an HIV-1 molecular epidemiological study among a recently infected population in Yunnan in 2015.

Methods

Study participants and sample collection

From January 2015 to June 2015, 6357 HIV/AIDS cases were newly reported through different screening methods in Yunnan Province. Before carrying out the recent HIV-1 infection assay, 1755 long-term infected cases with CD4+ T lymphocytes < 200 cells/μl or AIDS-defining illnesses were excluded. The remaining 4602 cases were further tested for the status of recent HIV-1 infection. The adults’ written consents were provided by themselves. The juveniles’ consents for HIV testing were provided by their guardians, if they had HIV, written consents about this study for HIV control and prevention were further obtained from their guardians when informing testing results. The study was approved by the Biomedical Ethics Review Committee of Yunnan Province.

Recent HIV-1 infection tested with BED-capture enzyme immunoassay (CEIA)

BED-capture enzyme immunoassay (BED-CEIA) (Calypte Biomedical Corporation, Portland, OR, USA) was used to test for the recency of HIV-1 infection. According to the manufacturer’s instructions, plasma specimens were first tested individually. If the normalized OD (ODn) was > 1.2, the specimen was classified as having a long-term infection. Specimens with ODn < 1.2 were tested again in triplicate for confirmation. In confirmatory testing, specimens with ODn < 0.8 were classified as recent infections. According to the evaluation carried out in China, the duration of a recent HIV infection determined by BED-CEIA was 168 days after seroconversion.

Amplification of HIV-1 gene fragments

According to the manufacturer’s instructions, viral RNA was extracted from 140 μl of plasma using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA). The partial gag gene (HXB2: 781–1861), partial pol gene (HXB2: 2147–3462) and partial env gene (HXB2: 7002–7541) were amplified using nested polymerase chain reactions (PCR). The successfully amplified products were sent to SinoGenoMax Co. (Beijing, China) for sequencing. The primers and procedures for nested PCR and the primers for sequencing were described in a previous study [13].

Sequence analysis

The sequences were assembled using Sequencher 5.1 (Gene Codes, Ann Arbor, MI). The assembled sequences were aligned with Bio-Edit 7.0 software and further manually edited. HIV-1 reference sequences were selected and downloaded from the HIV databases of the Los Alamos National Laboratory (LANL) (https://www.hiv.lanl.gov). Neighbour-joining phylogenetic trees were construed with the Kimura 2-parameter model with 1000 bootstrap replicates with MEGA version 6.0 [26]. The sequences with possible intersubtype recombination were analysed with the Recombination Identification Programme (RIP, version 3.0; https://www.hiv.lanl.gov). HIV-1 genotyping was based on at least two segments of the gag, pol and env genes. Samples with only one gene segment were excluded. Finally, 260 samples had gag, pol and env segments, 108 had gag and pol segments, 45 had gag and env segments, and 26 had pol and env segments. The above 439 samples were considered to have achieved reliable HIV-1 genotyping.

Geographic distribution analysis and spatial scan statistics analysis of HIV-1 genotypes

The distribution density of each HIV-1 genotype at the county level was displayed using a dot density map, within which the number of dots in a county represented the proportion of cases with a given genotype in this county among all the cases successfully genotyped. One dot represented 0.025% of the population.

The spatial scanning statistic is used to detect possible spatial clusters by scanning across space with a radius varying circular window. The null hypothesis is that the risks of event occurrence are equal between the inside and outside windows, while the alternative hypothesis is that the risks are different between the inside and outside windows.

Under the Poisson assumption, the log likelihood ratio (LLR) is calculated as:

$$ LLR={\left(\frac{c}{E\left[c\right]}\right)}^c{\left(\frac{C-c}{C-E\left[c\right]}\right)}^{C-c}I\left(\right) $$

where C is the total number of cases; c is the observed number of cases within the window; E [c] is the covariate adjusted expected number of cases within the window under the null-hypothesis; I() is an indicator function.

For each location, by varying the size, a window with the maximum likelihood will be found. The p-value is obtained through Monte Carlo simulations, by comparing the rank of the maximum likelihood from the real data set with those from random simulated data sets. If this rank is R, then p = R / (1+ n), where n is the number of random simulations. Among the regions with statistical significance, the regions with the maximum likelihood are defined as the most likely cluster (primary cluster), and the others with a smaller likelihood are defined as secondary clusters.

The relative risk (RR) represents the aggregation risk of the cluster compared with the rest of the regions. RR is defined as the estimated risk within the cluster divided by the estimated risk outside the cluster. The mathematical notation is as follows:

$$ RR=\frac{c/E\left[c\right]}{\left(C-c\right)/\left(C-E\left[c\right]\right)} $$

where the meaning of the symbols are the same as above.

In this study, SaTScan 9.6 was used for the spatial scan statistics analysis [27]. A Poisson-based model combined with circular scan windows was used. For the analysis of each HIV-1 genotype, three types of files (case file, population file and coordinates file) were submitted into SaTScan. The case file contained the numbers of cases with a given HIV-1 genotype in each county. The population file contained the number of populations in each county. The coordinates file contained the coordinates of each county. The most commonly used maximum spatial cluster size is 50% of the population at risk [28,29,30,31,32,33,34]. To choose the maximum spatial cluster size, a sensitivity analysis was performed with 10, 30 and 50% of the population at risk in the spatial window (Additional file 1: Table S1). The significant clusters of CRF01_AE, CRF07_BC, CRF08_BC and URFs showed no difference when using these settings. Therefore, 10% of the population at risk was used as the maximum spatial cluster size. The spatial clusters were mapped using Quantum GIS [35].

Genotypic analysis of HIV-1 drug resistance

Surveillance drug resistance mutations (SDRMs) within the obtained pol sequences were screened using the Calibrated Population Resistance (CPR) Tool (Version 6.0) on Stanford HIV Drug Resistance Database (http://hivdb.stanford.edu) [36].

Statistical analysis

Statistical analyses were performed with the SPSS 21.0 statistical analysis software package (SPSS Inc. Chicago, IL). Demographic characteristics associated with HIV-1 genotypes were first analysed with univariate logistic regression. For each HIV-1 genotype, variables marginally significant with p < 0.10 in univariate logistic analysis were further analysed with multivariate logistic regression. The trend chi-square test was performed with the linear-by-linear association method under the chi-square test of SPSS. All tests were two-tailed and statistical significance was considered when p- value < 0.05.

Results

Demographic characteristics of the study participants

As shown in Fig. 1, during the first half of 2015, 6,357 HIV/AIDS cases were newly diagnosed through the different screening methods in Yunnan Province. Based on information regarding disease progression and a BED-CEIA assay, 586 samples were identified as recent infections and were subjected to HIV-1 gene amplification. Finally, 444 gag, 411 pol and 357 env sequences were obtained, which were further subjected to phylogenetic analyses (Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3). By combining the genotyping of at least two gene segments, 439 samples obtained explanative HIV-1 genotyping. The constituent of the 439 participants showed no significant differences with that of the total 586 recently infected individuals (Additional file 5: Table S2).

Fig. 1
figure 1

Flow chart of the study process. The flow chart shows the overall study process, including HIV-1 screening and confirmation, recent infection testing and HIV-1 genotyping

Of the 439 participants, 62.0% (272/439) were males, and 38.0% (167/439) were females; the median age was 35 years (range: 15–81 years); Han ethnicity accounted for 62.0% (272/439), while other ethnicities accounted for 38.0%, including Yi, Hani, Dai, Zhuang, Jingpo, Bai, Hui, Wa, Lahu, Lisu, Miao, Yao, Man, Achang, Buyi, Naxi and Nu. Regarding marital status, 33.7% (148/439) were single, 43.7% (192/439) were married, and 22.6% (99/439) were divorced or widowed. Of the 439 participants, 5.0% (22/439) were non-Chinese, most of which were Burmese living Dehong Prefecture (72.7%, 16/22). For the transmission routes, heterosexual contact, homosexual contact and intravenous drug injection accounted for 81.1% (356/439), 11.2% (49/439) and 7.7% (34/439), respectively.

HIV-1 genotypes in recent HIV infections

Among the 439 samples successfully genotyped, more than ten HIV-1 genotypes were identified, including two subtypes, nine CRFs and six discrete URFs. CRF08_BC was the most common (45.3%, 199/439), followed by CRF07_BC (19.4%, 85/439), URFs (18.2%, 80/439), CRF01_AE (11.4%, 50/439), subtype C (2.1%, 9/439), CRF85_BC (1.1%, 5/439), CRF55_01B (0.9%, 4/439), subtype B (0.5%, 2/439), CRF64_BC (0.5%, 2/439), CRF59_01B (0.2%, 1/439), CRF83_cpx (0.2%, 1/439) and CRF87_cpx (0.2%, 1/439). Among the six discrete URFs, BC recombinants were the predominant recombinant form (66.3%, 53/80), followed by BC/CRF01_AE (10.0%, 8/80), CRF07_BC/CRF01_AE (8.8%, 7/80), B/CRF01_AE (6.3%, 5/80), C/CRF01_AE (5.0%, 4/80) and CRF07_BC/CRF08_BC (3.8%, 3/80).

Demographic characteristics associated with HIV-1 genotypes

Demographic characteristics associated with HIV-1 genotypes were further analysed using logistic regression (Table 1). The results of the multivariate analyses showed that homosexual contact was significantly associated with CRF01_AE and CRF07_BC infection; females, Chinese, heterosexual contact and intravenous drug injection were significantly associated with CRF08_BC infection; males and non-Chinese had a higher risk of URF infection than females.

Table 1 Demographic characteristics associated with HIV-1 genotypes

Spatial distribution of HIV-1 genotypes

The spatial distribution of HIV-1 genotypes in recent infections was analysed at the county level (Fig. 2). CRF08_BC was the most widely distributed genotype, which was found in 54.3% (70/129) of the counties in Yunnan Province. However, more CRF08_BC cases were found in the six eastern prefectures (Zhaotong, Qujing, Kunming, Yuxi, Honghe and Wenshan). The other three main genotypes, CRF07_BC, URFs and CRF01_AE, were found in 38.0% (49/129), 34.1% (44/129) and 20.2% (26/129) of all counties, respectively. The distributions of subtype C, CRF85_BC, CRF55_01B, subtype B, CRF64_BC, CRF59_01B, CRF83_cpx and CRF87_cpx were limited.

Fig. 2
figure 2

Geographic distribution of HIV-1 genotypes in recent infections in Yunnan. The dot density map for subtype B, subtype C, CRF01_AE, CRF07_BC, CRF08_BC, CRF55_01B, CRF55_01B, CRF64_BC, CRF83_cpx, CRF85_BC, CRF87_cpx and URFs shows the percentage of each genotype in each county. One dot was defined as 0.025% of the population

Spatial clusters were further analysed using spatial scan statistics (Fig. 3). For CRF08_BC, two non-overlapping statistically significant clusters were detected. Among them, the most likely cluster (RR = 4.78) included 11 counties in Honghe Prefecture and Wenshan Prefecture, and the secondary clusters (RR = 4.24) included two counties in Zhaotong Prefecture. In these two spatial clusters, 95.2 and 90.0% of observed cases were infected through heterosexual contact (Table 2). For CRF07_BC, there was only one most likely cluster (RR = 3.68) detected in Kunming Prefecture (Table 2). For CRF01_AE, two non-overlapping statistically significant clusters were detected. Among them, the most likely cluster (RR = 17.41) was located in one county in Dali Prefecture, and the secondary cluster (RR = 5.36) included six counties in Kunming Prefecture and Yuxi Prefecture. In the most likely cluster and the secondary cluster, 80.0 and 63.6% of observed cases were infected through homosexual contact (Table 2). For URFs, two non-overlapping statistically significant clusters were detected. Among them, the most likely cluster (RR = 9.62) was located in the west, including six counties in Dehong Prefecture and Baoshan Prefecture, in which 60.9% of observed cases were infected through intravenous drug use; the secondary cluster (RR = 3.51) was located in the southeast and included nine counties in Honghe Prefecture and Wenshan Prefecture, in which 100.0% of observed cases were infected through heterosexual contact (Table 2).

Fig. 3
figure 3

Spatial clusters of HIV-1 genotypes in recent infections in Yunnan. The counties in the primary and secondary clusters of CRF01_AE, CRF07_BC, CRF08_BC and URFs are marked in the maps

Table 2 The characteristics of spatial clusters for CRF01_AE, CRF07_BC, CRF08_BC and URFs

Genotypic analysis of transmitted drug resistance (TDR)

Among the 411 participants with pol sequences, ten (2.4%) were identified as harbouring SDRMs. As shown in Table 3, 0.5% (2/411), 1.2% (5/411), and 0.7% (3/411) of the sequences harboured SDRMs to nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs), respectively. As shown in Table 4, the prevalence of SDRMs increased with age (trend χ2 = 4.739, p = 0.033).

Table 3 Demographic characteristics of ten individuals harbouring transmitted drug resistance associated mutations
Table 4 The prevalence of SDRMs in different age groups

Discussion

In this study, we performed a cross-sectional HIV-1 molecular epidemiological study to track the characteristics and distribution of HIV-1 genotypes in recent infections in Yunnan. For the first time, the spatial epidemiological characteristics of HIV-1 genotypes among recent infections were analysed.

As described in the previous studies [12], CRF08_BC, CRF07_BC, URFs and CRF01_AE were the main HIV-1 genotypes circulating in Yunnan. However, compared with the surveillance carried out in 2009, the ranks of these four genotypes changed. The proportion of CRF08_BC among recent infections remained the highest and showed no significant difference from 2009 to 2015 (40.8% in 2009, 47.0% in 2014 and 45.3% in 2015; χ2 = 1.605, p = 0.450) [12, 37]. The proportion of CRF07_BC increased from fourth to second place and showed a statistically increasing trend for 2009–2015 (9.2% in 2009, 14.4% in 2014 and 18.3% in 2015; trend χ2 = 8.557, p = 0.003). The proportion of URFs declined from second to third place and showed a decreasing trend (27.7% in 2009, 18.2% in 2014 and 18.2% in 2015; χ2 = 6.272, p = 0.045). The proportion of CRF01_AE declined from third to fourth place and showed a decreasing trend that was not statistically significant (18.5% in 2009, 15.8% in 2014 and 11.4% in 2015; χ2 = 5.437, p = 0.064). The overall increase in CRF07_BC may be partly due to the increase in CRF07_BC among homosexual contact. Among MSM, the proportion of CRF07_BC increased from 17.2% in 2014 to 42.9% in 2015 (χ2 = 5.379, p = 0.026).

In addition to the above four predominant genotypes, two subtypes and six CRFs were found in recent infections. The genetic diversity of HIV-1 constitutes a characteristic of the HIV-1 epidemic in Yunnan. HIV-1 recombination was active in the Yunnan-Myanmar border area, where URFs were highly prevalent [9, 38] and more than ten CRFs were identified recently, including CRF57_BC, CRF62_BC, CRF64_BC, CRF65_cpx, CRF78_cpx, CRF82_cpx, CRF83_cpx, CRF86_BC, CRF87_cpx, CRF88_BC and CRF96_cpx [18,19,20,21,22,23,24,25, 39]. In this study, CRF64_BC, CRF83_cpx and CRF87_cpx were found in western Yunnan, which suggested that they were circulating locally. CRF55_01B and CRF59_01B were originally identified in MSM outside Yunnan [15, 16]; CRF85_BC was originally identified in the neighbouring province, Sichuan Province [40]. In this study, these three CRFs were found in eastern Yunnan, where the floating population was relatively larger.

The demographic characteristics associated with HIV-1 genotypes were further analysed using multivariate logistic regression. Because heterosexual contact was the major transmission route (81.1%), the absolute numbers of the main genotypes in this transmission route were the largest. However, the proportions of CRF01_AE and CRF07_BC in homosexual contact were higher than those in heterosexual contact and intravenous drug injection, while the proportion of CRF08_BC in homosexual contact was lower than those in the other two transmission routes. These findings suggested that the homosexually transmitted population remained relatively unique from the other two populations. The proportion of CRF08_BC was higher in females. The reason for this finding could be that the predominant transmission route among the female cases was heterosexual contact (95.8%, 160/167), and in this transmission route, CRF08_BC was the predominant HIV-1 genotype. Furthermore, combined with gender and transmission routes, the subjects were divided into five subgroups. The proportion of CRF08_BC was the highest in the subgroup of female heterosexual contact and the lowest in the subgroup of male homosexual contact (Additional file 6: Table S3). We also found that the proportion of URFs among males was higher than that among females, that the proportion of URFs among non-Chinese was higher than that among Chinese. Usually, URFs were identified in dually or multiply infected individuals. This finding suggested that males and non-Chinese tended to have more complex risk factors. To fully understand HIV-1 genetic characteristics among non-Chinese, the further study should be carried out in the border areas.

For the first time, we described the geographic distribution of HIV-1 genotypes among recent infections at the county level, which was more elaborate than the previous study at the prefecture level [12]. Among the four predominant genotypes, the coverage of CRF08_BC was greater than those of CRF07_BC, URFs and CRF01_AE. However, the other genotypes were only found in one to six counties. In some areas, the four predominant genotypes displayed a clustering tendency. To reveal the spatial clustering characteristics, spatial scanning analysis was performed. There were two statistically significant spatial clusters for CRF08_BC, among which the most likely cluster was located in southeast Yunnan Province and bordered with Vietnam. Strikingly, more than 90% of observed cases in these two spatial clusters were infected through heterosexual contact, which suggested that more efforts for the prevention of heterosexual contact transmission should be taken in these areas.

Among the two statistically significant clusters for CRF01_AE, the most likely cluster in Dali Prefecture and the second secondary cluster in Kunming and Yuxi Prefecture were largely composed of the cases infected through homosexual contact. Among MSM identified as recent HIV infections in this study, the proportions of MSM identified in Kunming and Dali were 46.8% (24/59) and 16.9% (11/59), respectively, and ranked the top two (Additional file 7: Table S4). If the HIV incidence is assumed to be roughly the same among MSM, the population sizes of MSM in Kunming and Dali were the highest in Yunnan Province. Strikingly, the RR of the most likely cluster for CRF01_AE was the highest among all the detected spatial clusters. RR means how many times the estimated risk in cluster is as high as the estimated risk outside the cluster. The higher RR reflects the higher aggregation degree of the specific risk. Thus, this suggested that the risk of CRF01_AE infection displayed a high clustering tendency in this cluster.

As mentioned above, Dehong Prefecture, one of the prefectures bordering Myanmar, was a hot spot of active HIV-1 genetic recombination. In this study, the most likely cluster for URFs was detected in the China-Myanmar border area, including four counties in Dehong Prefecture and one county in Baoshan Prefecture. In this URF cluster, 65.2% (15/23) of the observed cases were Burmese residing in Dehong and Baoshan, and the transmission risks of heterosexual contact (60.9%, 14/23) and intravenous drug use (39.1%, 9/23) coexist. Recent studies also suggested that Burmese living in Yunnan contributed to cross-border transmission of HIV-1 in the China-Myanmar border areas [41, 42]. These results suggested that AIDS prevention and control in border areas should be further considered. The secondary cluster for URFs was located in the nine counties in Honghe Prefecture and Wenshan Prefecture, which almost overlapped with the most likely cluster for CRF08_BC. In this URF cluster, 100% of the observed cases were infected through heterosexual contact, which further suggested the high transmission risk of heterosexual contact in this area. Overall, the spatial clustering analysis for recent infections could provide clues for potential transmission risk, based on which the targeted measures for AIDS control and prevention should be developed.

The prevalence of SDRMs in recent infections was 2.4% in 2015, which was not significantly higher than that in 2014 (1.8%, χ2 = 0.263, p = 0.791) [37]. However, we found that the prevalence of TDR increased with age. In the group that was above 50 years of age, the prevalence of TDR was over 5% and reached a moderate level, which should be considered. To control HIV-1 TDR, a priority is to strengthen the standard management of patients and regularly perform drug resistance surveillance.

Conclusions

With the occurrence of HIV-1 epidemics in Yunnan, HIV-1 genotypes have become more diverse. Our study revealed two subtypes, nine CRFs and six discrete URFs circulating in recent infections, among which CRF08_BC, CRF07_BC, URFs, and CRF01_AE were still the predominant strains. The distribution of the four predominant genotypes in the population was associated with demographic characteristics. The geographic distribution of HIV-1 genotypes in recent infections was further analysed at the county level. The four predominant genotypes demonstrated a clustering tendency in some geographic areas, which was confirmed by the spatial scan statistics. Some spatial clusters were associated with a specific transmission route. TDR-associated mutations remained low in recent infections. Our findings provide valuable information to improve strategies to prevent new infections.

Availability of data and materials

The sequences obtained in this study were submitted to NCBI GenBank under accession numbers MH572689-MH573900. The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Abbreviations

BED-CEIA:

BED-capture enzyme immunoassay

CRF:

Circulating recombinant form

FSW:

Female sexual worker

HIV-1:

Human immunodeficiency virus-1

IDU:

Intravenous drug user

MSM:

Men who have sex with men

NNRTIs:

Non-nucleoside reverse transcriptase inhibitors

NRTIs:

Nucleoside reverse transcriptase inhibitors

PIs:

Protease inhibitors

SDRMs:

Surveillance drug resistance mutations

TDR:

Transmitted drug resistance

URF:

Unique recombinant form

References

  1. Lu L, Jia M, Ma Y, Yang L, Chen Z, Ho DD, et al. The changing face of HIV in China. Nature. 2008;455(7213):609–11.

    Article  CAS  Google Scholar 

  2. Jia M, Luo H, Ma Y, Wang N, Smith K, Mei J, et al. The HIV epidemic in Yunnan Province, China, 1989-2007. J Acquir Immune Defic Syndr. 2010;53(Suppl 1):S34–40.

    Article  Google Scholar 

  3. National Health and Family Planning Commission of the People’s Republic of China. 2015 China AIDS Response Progress Report. 2015.

    Google Scholar 

  4. Wu Z, Sullivan SG, Wang Y, Rotheram-Borus MJ, Detels R. Evolution of China’s response to HIV/AIDS. Lancet. 2007;369(9562):679–90.

    Article  Google Scholar 

  5. Graf M, Shao Y, Zhao Q, Seidl T, Kostler J, Wolf H, et al. Cloning and characterization of a virtually full-length HIV type 1 genome from a subtype B’-Thai strain representing the most prevalent B-clade isolate in China. AIDS Res Hum Retrovir. 1998;14(3):285–8.

    Article  CAS  Google Scholar 

  6. Luo CC, Tian C, Hu DJ, Kai M, Dondero T, Zheng X. HIV-1 subtype C in China. Lancet. 1995;345(8956):1051–2.

    Article  CAS  Google Scholar 

  7. Shao Y, Zhao Q, Wang B, Chen Z, Su L, Zeng Y, et al. Sequence analysis of HIV env gene among HIV infected IDUs in Yunnan epidemic area of China. Chin J Virol. 1994;10(4):291–9.

    CAS  Google Scholar 

  8. Cheng H, Zhang J, Capizzi J, Young NL, Mastro TD. HIV-1 subtype E in Yunnan, China. Lancet. 1994;344(8927):953–4.

    Article  CAS  Google Scholar 

  9. 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.

    Article  CAS  Google Scholar 

  10. Tee KK, Pybus OG, Li XJ, Han X, Shang H, Kamarulzaman A, et al. Temporal and spatial dynamics of human immunodeficiency virus type 1 circulating recombinant forms 08_BC and 07_BC in Asia. J Virol. 2008;82(18):9206–15.

    Article  CAS  Google Scholar 

  11. Yao Y, Wang N, Chu J, Ding G, Jin X, Sun Y, et al. Sexual behavior and risks for HIV infection and transmission among male injecting drug users in Yunnan, China. Int J Infect Dis. 2009;13(2):154–61.

    Article  Google Scholar 

  12. Chen M, Yang L, Ma Y, Su Y, Yang C, Luo H, et al. Emerging variability in HIV-1 genetics among recently infected individuals in Yunnan, China. PLoS One. 2013;8(3):e60101.

    Article  CAS  Google Scholar 

  13. Chen M, Ma Y, Su Y, Yang L, Zhang R, Yang C, et al. HIV-1 genetic characteristics and transmitted drug resistance among men who have sex with men in Kunming, China. PLoS One. 2014;9(1):e87033.

    Article  Google Scholar 

  14. Chen M, Ma Y, Chen H, Dai J, Dong L, Yang C, et al. HIV-1 genetic transmission networks among men who have sex with men in Kunming, China. PLoS One. 2018;13(4):e0196548.

    Article  Google Scholar 

  15. 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.

    Google Scholar 

  16. Han X, An M, Zhang W, Zhao B, Chu Z, Takebe Y, et al. 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. 2013;1:3.

    Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. Miao J, Ran J, Song Y, Liu Y, Gao L, Miao Z, et al. Characterization of a novel HIV-1 circulating recombinant form, CRF01_AE/B’/C (CRF96_cpx), in Yunnan, China. AIDS Res Hum Retrovir. 2018;34(4):393–7.

    Article  CAS  Google Scholar 

  19. Li Y, Miao J, Miao Z, Song Y, Wen M, Zhang Y, et al. Identification of a novel HIV type 1 circulating recombinant form (CRF86_BC) among heterosexuals in Yunnan, China. AIDS Res Hum Retrovir. 2017;33(3):279–83.

    Article  CAS  Google Scholar 

  20. Hu Y, Wan Z, Zhou YH, Smith D, Zheng YT, Zhang C. Identification of two new HIV-1 circulating recombinant forms (CRF87_cpx and CRF88_BC) from reported unique recombinant forms in Asia. AIDS Res Hum Retrovir. 2017;33(4):353–8.

    Article  CAS  Google Scholar 

  21. Song Y, Feng Y, Miao Z, Wang B, Yang M, Zhang AM, et al. Near-full-length genome sequences of a novel HIV-1 circulating recombinant form, CRF01_AE/B’/C (CRF78_cpx), in Yunnan, China. AIDS Res Hum Retrovir. 2016;32(6):601–6.

    Article  CAS  Google Scholar 

  22. Wei H, Liu Y, Feng Y, Hsi J, Xing H, He X, et al. Genome sequence of a novel HIV-1 circulating recombinant form (CRF57_BC) identified from Yunnan, China. AIDS Res Hum Retrovir. 2014;30(4):384–8.

    Article  CAS  Google Scholar 

  23. Wei H, His J, Feng Y, Xing H, He X, Liao L, et al. Identification of a novel HIV-1 circulating recombinant form (CRF62_BC) in western Yunnan of China. AIDS Res Hum Retrovir. 2014;30(4):380–3.

    Article  CAS  Google Scholar 

  24. Hsi J, Wei H, Xing H, Feng Y, He X, Liao L, et al. Genome sequence of a novel HIV-1 circulating recombinant form (CRF64_BC) identified from Yunnan, China. AIDS Res Hum Retrovir. 2014;30(4):389–93.

    Article  CAS  Google Scholar 

  25. Feng Y, Wei H, Hsi J, Xing H, He X, Liao L, et al. Identification of a novel HIV type 1 circulating recombinant form (CRF65_cpx) composed of CRF01_AE and subtypes B and C in Western Yunnan, China. AIDS Res Hum Retrovir. 2014;30(6):598–602.

    Article  CAS  Google Scholar 

  26. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  Google Scholar 

  27. Kulldorff M. and Information Management Services, Inc. SaTScanâ„¢ v9.6: Software for the spatial and space-time scan statistics, 2018 www.satscan.org.

    Google Scholar 

  28. Wand H, Ramjee G. Targeting the hotspots: investigating spatial and demographic variations in HIV infection in small communities in South Africa. J Int AIDS Soc. 2010;13:41.

    Article  Google Scholar 

  29. Gonzalez R, Augusto OJ, Munguambe K, Pierrat C, Pedro EN, Sacoor C, et al. HIV incidence and spatial clustering in a rural area of southern Mozambique. PLoS One. 2015;10(7):e0132053.

    Article  Google Scholar 

  30. Cuadros DF, Abu-Raddad LJ. Spatial variability in HIV prevalence declines in several countries in sub-Saharan Africa. Health Place. 2014;28:45–9.

    Article  Google Scholar 

  31. Hixson BA, Omer SB, del Rio C, Frew PM. Spatial clustering of HIV prevalence in Atlanta, Georgia and population characteristics associated with case concentrations. J Urban Health. 2011;88(1):129–41.

    Article  Google Scholar 

  32. Wylie JL, Cabral T, Jolly AM. Identification of networks of sexually transmitted infection: a molecular, geographic, and social network analysis. J Infect Dis. 2005;191(6):899–906.

    Article  Google Scholar 

  33. Jennings JM, Curriero FC, Celentano D, Ellen JM. Geographic identification of high gonorrhea transmission areas in Baltimore, Maryland. Am J Epidemiol. 2005;161(1):73–80.

    Article  Google Scholar 

  34. Acharya BK, Cao C, Lakes T, Chen W, Naeem S. Spatiotemporal analysis of dengue fever in Nepal from 2010 to 2014. BMC Public Health. 2016;16(1):849.

    Article  Google Scholar 

  35. Quantum GIS Development Team, Quantum GIS Geographic Information System. Open Source Geospatial Foundation Project, 2002, http://www.qgis.org.

  36. Gifford RJ, Liu TF, Rhee SY, Kiuchi M, Hue S, Pillay D, et al. The calibrated population resistance tool: standardized genotypic estimation of transmitted HIV-1 drug resistance. Bioinformatics. 2009;25(9):1197–8.

    Article  CAS  Google Scholar 

  37. Chen M, Jia MH, Ma YL, Luo HB, Chen HC, Yang CJ, et al. The changing HIV-1 genetic characteristics and transmitted drug resistance among recently infected population in Yunnan, China. Epidemiol Infect. 2018;146(6):775–81.

    Article  CAS  Google Scholar 

  38. Pang W, Zhang C, Duo L, Zhou YH, Yao ZH, Liu FL, et al. Extensive and complex HIV-1 recombination between B’, C and CRF01_AE among IDUs in south-East Asia. AIDS. 2012;26(9):1121–9.

    Article  Google Scholar 

  39. Chen X, Ye M, Duo L, Pang W, Smith D, Zhang C, et al. First description of two new HIV-1 recombinant forms CRF82_cpx and CRF83_cpx among drug users in northern Myanmar. Virulence. 2017;8(5):497–503.

    Article  Google Scholar 

  40. Su L, Wei D, Yang H, Zeng Y, Hu Y, Yuan D, et al. Identification of a novel HIV-1 circulating recombinant form (CRF85_BC) in Sichuan, China. AIDS Res Hum Retrovir. 2016;32(9):895–9.

    Article  CAS  Google Scholar 

  41. Chen X, Zhou YH, Ye M, Wang Y, Duo L, Pang W, et al. Burmese injecting drug users in Yunnan play a pivotal role in the cross-border transmission of HIV-1 in the China-Myanmar border region. Virulence. 2018;9(1):1195–204.

    Article  CAS  Google Scholar 

  42. Chen X, Duo L, Ye M, Zhang C, Zheng YT. Non-Chinese immigrants: challenge faced by Yunnan of China to achieve the 90-90-90 goals. Virol Sin. 2018;33(3):291–3.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge the staff at the local Centers for Disease Control and Prevention for their assistance in sample collection.

Funding

This work was supported by the National Natural Science Foundation of China (81560327) and the Research Project for Yunnan Science and Technology (2015FB200). The funders had no role in the study design, data collation and analysis, interpretation of data, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Conceived and designed the experiments: MC, YM, ZS and MJ. Collected the samples: MC, HC, JD, CY, LD, XJ, MY and LY. Performed the experiments: MC, HC, JD, CY, LD, XJ, MY and LY. Collected epidemiological data: YM, HL, LS, ZS and MJ. Analysed the data: MC. Wrote the paper: MC. Approved the final version of this manuscript: MC, YM, HC, JD, HL, CY, LD, XJ, MY, LY, LS, ZS and MJ.

Corresponding author

Correspondence to Manhong Jia.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Biomedical Ethics Review Committee of Yunnan Province. The adults’ written consents were provided by themselves. The juveniles’ consents for HIV testing were provided by their guardians, if they had HIV, written consents about this study for HIV control and prevention were further obtained from their guardians when informing testing results.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. The comparison of scan statistical results by using different percentages of the population at risk. (PDF 70 kb)

Additional file 2:

Figure S1. Neighbour-joining phylogenetic tree of the partial gag gene. The scale bar indicates 5% nucleotide sequence divergence. Values on the branches represent the percentage of 1000 bootstrap replicates. (PDF 680 kb)

Additional file 3:

Figure S2. Neighbour-joining phylogenetic tree of the partial pol gene. The scale bar indicates 5% nucleotide sequence divergence. Values on the branches represent the percentage of 1000 bootstrap replicates. (PDF 719 kb)

Additional file 4:

Figure S3. Neighbour-joining phylogenetic tree of the partial env gene. The scale bar indicates 10% nucleotide sequence divergence. Values on the branches represent the percentage of 1000 bootstrap replicates. (PDF 591 kb)

Additional file 5:

Table S2. The constituent of subjects successfully genotyped. (PDF 159 kb)

Additional file 6:

Table S3. The proportion of CRF08_BC in men and women with different transmission routes. (PDF 68 kb)

Additional file 7:

Table S4. The distribution of recent HIV infections was attributed to the main transmission routes in each prefecture. (PDF 67 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, M., Ma, Y., Chen, H. et al. Spatial clusters of HIV-1 genotypes in a recently infected population in Yunnan, China. BMC Infect Dis 19, 669 (2019). https://doi.org/10.1186/s12879-019-4276-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-019-4276-9

Keywords