- Open Access
Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China
BMC Infectious Diseases volume 22, Article number: 602 (2022)
Sichuan province is one of the highest AIDS epidemic provinces in China, with a large number of floating population. The annual number of cases of HIV/AIDS reported in Sichuan has been the highest province in China for several successive years. There is a lack of widespread and representative data on the distribution of HIV genotypes in Sichuan. We aim to investigate the characteristics of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China.
Archived plasma samples (n = 1524) from HIV-1 newly-diagnosed individuals in April 2019 were selected by cross-sectional investigation from all 21 cities in Sichuan province. Phylogenetic relationship, transmission cluster, and genotypic drug resistance analyses were performed using HIV-1 polymerase (pol) gene sequences. We also analysed the association of demographic and virological factors with transmitted drug-resistance (TDR) and transmission clusters.
Partial pol gene sequences were obtained from 1297 cases. HIV-1 epidemic strains in Sichuan province: the majority of genotypes were circulating recombinant form (CRF) 07_BC (675, 52.04%), CRF01_AE (343, 26.45%), CRF08_BC (115, 8.87%), CRF85_BC (67, 5.17%), subtype B (33, 2.54%), the other genotypes only accounted for 4.93%, and unique recombinant forms (URFs) (23, 1.77%) were observed in the study, and the difference of age, ethnicity, education, occupation, region and transmission pathway of different genotypes were statistically significant. According to WHO HIVDR surveillance threshold, the level of TDR has reached a medium level, with 72 of 1297 (5.55%) cases carrying drug-resistance mutation sites, TDR mutation frequency to nonnucleoside reverse transcriptase inhibitors (NNRTIs, 3.85%) was much higher than nucleoside reverse transcriptase inhibitors (NRTIs, 0.31%) and protease inhibitors (PIs, 1.70%), and CRF08_BC was a risk factor for TDR (odds ratio, 8.32; 95% CI 4.38–15.80 for CRF07_BC, P < 0.05). The most common drug resistance HIV-1 mutation pattern for NNRTI was V106 (1.31%, 17/1297) and E138 (1.16%, 15/1297), and for PI was M46 (0.69%, 9/1297). A total of 205 (15.8%) pol sequences were involved in the genetic transmission network clusters, CRF01_AE (odds ratio, 2.369; 95% CI 1.659–3.382; P < 0.05), subtype B (odds ratio, 13.723; 95% CI 6.338–29.71; P < 0.05), drug resistance (odds ratio, 0.306; 95% CI 0.106–0.881; P < 0.05) and different levels of education (P < 0.05) were significantly associated to be in clusters.
The distribution of HIV-1 genotypes in Sichuan is more diverse and complex, and the Men who have sex with men (MSM) is underrated, arguing for behavior scaling up intervention in this specific population besides the elderly people with heterosexual transmission risk groups. The risk of TDR mutation frequency increased in newly diagnosed patients highlights the significance of genotypic drug resistance monitoring and molecular surveillance of pretreatment HIV-1 drug resistance. The regimen composed of TDF, 3TC and EFV was still currently the preferred solution used free first-line therapy.
The prevalence of human immunodeficiency virus type 1 (HIV-1) remains a major public health burden in China. Sichuan province is in the interior of Southwest China, adjacent to Yunnan and Guizhou in the South and Tibet in the West. The epidemic of HIV/AIDS has kept rising in Sichuan province since reported the first case of imported AIDS in 1991, the number of newly reported cases per year has increased by about 20% since 2006, by 2019, Sichuan province reported more than 170 thousand cases of HIV/AIDS, with the number of surviving cases ranking first in the country .
Among the four phylogeny groups that constitute HIV-1 (M, N, O, P), M group is the virulence factor of AIDS pandemic . There were nine subtypes (A, B, C, D, F, G, H, J and K) and at least 118 circulating recombinant forms (CRFs) (https://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html) within HIV-1 group M. With the genetic variation and cross recombination of genes between different genotypes, numerous new circulating recombinant forms (CRFs)/unique recombinant forms (URFs) are emerging, and the global epidemic situations of HIV has changed significantly. HIV enormous genetic variability and rapid evolution have led to its epidemic and therapeutic challenges. Worldwide, the proportion of subtype C was the highest, accounting for 48%, followed by subtype A (12%), subtype B (11%) and CRF02_AG (8%), CRF01_AE (5%), G subtype (5%) and D subtype (2%) . A national study showed that 8 subtypes and 21 CRFs have been identified in China since 2013 . CRF01_AE, CRF07_BC, CRF08_BC, subtype B were the four main subtypes, accounting for 89.0% of all HIV-1 infections . CRF01_AE mainly transmitted by sexual route spread from the southeast coast and southwest border to the whole country . With the wide prevalence of subtypes B and C among IDUs in Yunnan, it was provided sufficient conditions to the hybridization and recombination of HIV-1 strain of and resulting in the formation of recombinant subtypes CRF07_BC and CRF08_BC . Additionally, B subtype mainly originated in Thailand, and this subtype also became the main HIV-1 subtype transmitted by blood in Central China in the 1990s .
Since Chinese government officially launched the National Free Antiretroviral Treatment Program (NFATP) in 2003, the number of patients enrolled in this program has increased rapidly year by year, with the increasing coverage of antiviral therapy; the occurrence of drug resistance is also increasing. The emergence of HIV drug resistance (HIVDR) threatens the global scale-up of antiretroviral treatment (ART) for treating HIV infection, which could increase the risk of ART failure. It is important to conduct HIVDR surveys to estimate the transmitted drug-resistance (TDR) rate of ART-naive people living with HIV (PLWH) to provide baseline information for effective ART programs, delaying HIVDR occurrence and developing a rational public health strategy to control HIV prevalence. However, HIVDR testing for ART-naive PLWH is not routinely performed in China as a developing country with limited conditions. The rates of TDR vary throughout the world. In China, the overall prevalence of TDR among ART-naïve individuals in 2004 and 2005 was 3.8% , and the rate of TDR was 3.6% in 2015 .
Sichuan province has a population of approximately 85 million people, and is a developing area, Sichuan province is the most severely HIV affected area in China . In a previous study, we performed a comprehensive investigation of the HIV epidemic in Sichuan province, 2014  and some cities in Sichuan [12, 13]. With the development of HIV-1 epidemics. With the HIV epidemic growing and treatment scaling up, it is essential to investigate the changing trend of HIV-1 genetics in the province, and to conduct a province-wide transmitted drug resistance (TDR) survey to understand the frequency of transmitted drug-resistant viruses. Therefore, we conducted a large cross-sectional study of the recently infected population identified in April 2019 in the entire Sichuan province. We have used molecular phylogenetic analyses to complement HIV surveillance tasks of Sichuan province. Specifically, we have used this information to infer the distribution and characteristics of HIV-1 subtypes, to analyse which risk groups are currently more vulnerable to TDR HIV infection and in molecular transmitted cluster. The results obtained from this work may be useful in establishing and reinforcing preventive measures in specific target groups.
Material and methods
A total of 1524 plasma samples were obtained from individuals newly diagnosed as HIV-1 infected patients in April 2019 were selected by cross-sectional investigation from all 21 cities in Sichuan province. The epidemiological data of demographic characteristics (i.e., ethnic, sex, age, and marital status) were acquired from China Information System for Disease Control and Prevention. After eliminating duplicate samples, pol sequence (covering 1, 060 base pairs, HXB2: 2, 254–3, 313) information was successfully exported from 1297 PLWH in this study.
RNA extraction, amplification, and sequencing
About 5 ml of venous blood was extracted. The EDTA-K2 anticoagulated plasma samples were isolated from each participant and preserved in a – 80 °C freezer before sending the cold chain to the Sichuan CDC.
The viral nucleic acid was obtained from 200 μl plasma of PLWH by extraction machines (MagNA Pure LC system, Roche, Branchburg, NJ). Sequences were generated from the HIV-1 pol. The Reverse Transcription-Polymerase Chain Reaction (RT-PCR) was used to amplify the full-length protease gene in the pol region and the first 300 codons of the reverse transcriptase gene. The PCR products were dealt with electrophoresis with 1% agarose gel, and the amplified positive products were purified and sequenced by Beijing Genomics Research Center Ltd. The detailed amplification and sequencing performed as previously described .
Sequence analysis and genotype determination
The obtained sequence was spliced by using the analysis software Sequencher 5.1. The BioEdit Sequence Alignment Editor was used to edit and correct the sequence. MEGA 7.0 software was used to determine the genotype. All assembled sequences were aligned together with the reference sequences (from the HIV sequence database of Los Alamos National Laboratory in the United States) using Clustal W program in MEGA 7.0, and then further checked manually in Bioedit . FastTree was used to estimate an approximately maximum likelihood phylogenetic tree for pol sequences using the GTR + G + I nucleotide substitution model , using the ShimodairaHasegawa (SH) test embedded in the software to calculated the node (branch point) value of the evolutionary tree, defined clades with SH-like support ≥ 0.70 [17, 18]. The genotype were preliminarily determined by clustering the sample sequences with international reference strains, then the online analysis tool HIV Databases BLAST (https://www.hiv.lanl.gov/content/index) was used to review the results. Presented phylogenetic tree by using FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/). If the sequence whose genotype cannot be confirmed by phylogenetic tree and HIV Databases BLAST were considered as URFs. For example, the sequence between the CRF_01AE and B subtype reference strains in the phylogenetic tree, and the similar sequences in HIV Databases BLAST were 01B genotype, we judged this URF as URF_01B.
Transmission cluster identification
The gene distance between the clusters was calculated using the software hyphy 2.2.4 in TN93 model, and the relationship between the two sequences was determined by the minimum gene distance method. The transmission network was constructed with a threshold of pairwise genetic distance of 0.5%, which is more appropriate for identifying rapidly growing clusters, and was based on study findings showing that within monoinfected individuals pol sequences do not diverge more than 1% during the first 3 years of infection . The molecular transmission clusters were deduced by cytoscape 3.7.2 .
Drug resistance mutation analysis
Drug-resistance were analyzed based on genotypic interpretations defined by the Surveillance Drug Resistance Mutation list recommended by the World Health Organization. The Stanford University HIV DRUG RESISTANCE DATABASE (https://hivdb.stanford.edu/) was used to screen the pol sequences for surveillance drug resistance mutations (DRMs), it was evaluated relevant resistance for nucleoside reverse transcriptase inhibitors (NRTIs), nonnucleoside reverse transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs) [21, 22]. If HIV-1 strains carrying at least one TDR mutation were defined as resistant.
Data were entered in a spreadsheet (Micorsoft Office Excel 2013) and analyzed using SPSS 25 software (IBM, Chicago, IL, USA). The potential demographic characteristics (sex, age, ethnicity, marital status, and transmission route) were expressed in percentage, differences between the variables group for genotypes were evaluated using the χ2 test. Univariate logistic regression model was used to test the independent variables (i.e. sociodemographic) and transmission drug resistance individuals (1 = drug-resistance, 0 = non drug-resistance), multivariate logistic regression model was used to further test the significance of all micro significant (P < 0.1), variables in univariate logistic regression model. Potential demographic differences between patients in cluster and those out of cluster were evaluated by uni-multivariate logistic regression analysis. P value < 0.05 was considered statistically significant.
Demographic features and distribution of HIV-1 genetic forms
1297 individuals with successfully obtained pol sequences from 1524 infected patients with HIV-1, 15 genotypes were identified (Fig. 1). In the phylogenetic tree, sequences belonging to different genotypes clustered separately. Furthermore, several clusters supported by high bootstrap values were also found among different genotypes, suggesting that various founder viruses were introduced into populations separately. Major HIV-1 epidemic strains in Sichuan province were: CRF07_BC (675, 52.04%), CRF01_AE (343, 26.45%), CRF08_BC (115, 8.87%), CRF85_BC (67, 5.17%), subtype B (33, 2.54%). The other genotypes only accounted for 4.93%, of which CRF55_01B (13, 1.00%), CRF105_0108 (11, 0.85%), CRF79_0107 (5, 0.39%), CRF59_01B (3, 0.23%), CRF86_BC (2, 0.15%), CRF57_BC (2, 0.15%), subtype C (2, 0.15%), CRF65_cpx (1, 0.08%), CRF88_BC (1, 0.08%), subtype D (1, 0.08%), and 23 URFs that could not be classified (1.77%). URF including URF_BC (10, 0.77%), URF_01/BC (7, 0.54%), URF_B (4, 0.31%), URF_01/B (2, 0.15%).
It can be seen in Fig. 2 that the genotypes in the central Sichuan Basin (Chengdu plain) are more abundant, the CRF07_BC (86/102, 84.31%) in the Western Sichuan Basin account for a higher proportion, and the CRF01_AE (26/46, 56.52%) in the Northeast Sichuan (Bazhong, Dazhou and Guangyuan) Basin account for more. CRF07_BC in Panzhihua-Xichang region accounts for 87.65% (71/81); most of the CRF85_BC was distributed in southern Sichuan (i.e., Yibin, 57/146, 39.04%); CRF08_BC is mainly distributed in Luzhou and Yibin (60%, 69/115); CRF07_BC (40/175, 22.29%), CRF08_BC (50/175, 28.57%), CRF01_AE (41/175, 23.43%) have similar composition ratios in Luzhou city, and CRF105_0108 of aggregation were found in there.
Demographic information is summarized in Table 1, among the six major genotypes, the difference of age, ethnicity, education, occupation, region and transmission pathway of different genotypes were statistically significant (P < 0.05). The majority groups of CRF85_BC (95.53%, 64/67) and CRF105_0108 (90.90%, 10/11) with age ≥ 50 years old were found in the heterosexual transmission; 89.39% of ethnic Yi infected with CRF07_BC; 49.75% of ethnic Han was infected with CRF07_BC, which was found their infection rate was higher than that of other age groups in the 25 year old (48/71, 67.60%) and 25–40 years old (129/183, 70.49%) age group. Among the transmission route of IDUs, mother to child transmission (MCT), and MSM, CRF07_BC was the dominant strain and the composition ratio of was 100%, 77.78% and 54.43%. Furthermore, 13/13 CRF55_01B, 5/5 CRF79_0107, 2/3 CRF59_01B was detected in men.
Characteristics of transmitted drug-resistant
Of the 1297 subjects, 72 had drug-resistant mutations, among which 12 had high drug resistance, 21 had moderate drug resistance and 55 had low drug resistance, potential drug resistance were 66 cases. 50 cases (3.85%, 50/1297) were found to have non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance mutation, the mutation probability of E138 (30.0%, 15/50), K103 (20.0%, 10/50) and V106 (34.0%, 17/50) was significantly higher than that of other sites. The mutation probability of drug resistance site in CRF07_BC, CRF08_BC, CRF01_AE was higher than that in other genotype, which were resistant to all NNRTIs in varying degrees, and the high resistance was mostly caused by K103N mutation. Four cases (0.31%, 4/1297) had NRTI resistance mutation, the mutation rate was relatively low, but d4T, ddI and TDF caused by K65KR/R mutation were highly resistant. 22 subjects (1.70%, 22/1297) had protease inhibitor (PI) resistance mutation, M46 (36%, 9/22) and Q58 (28%, 7/22) had the highest mutation probability, most of them were low-grade resistance. Additionally, the most frequent mutation was V106, which was observed in 23.6% (17/72) of patients, HIV-1 strains with this mutation exhibited the degree of potential low-level resistance, except for V106VIM resulted in NVF and EFV were highly resistant and dot was moderately resistant. Followed by E138 (20.8%, 15/72), V179 (15.3%, 11/72) and K103 (13.9%, 10/72), all of them were NNRTI related mutations. Moreover, the most resistant drugs were RPV (30.5%, 22/72), NVP (29.2%, 21/72), DOR (25.0%, 18/72). The distribution of drug resistance and mutation sites was shown in Fig. 3.
The drug resistance level is shown in Fig. 4. The level of HIV-transmitted drug resistance is 5.55% (72/1297) in Sichuan province. The level of transmissible drug-resistance in city of Liangshan, Luzhou, Deyang and Guang’an the high proportion, were 9.80%, 8.00%, 8.11% and 8.57% (Fig. 5). The overall prevalence of HIV-1 TDR to NNRTI (3.85%, 50/1297) was higher than that to NRTI (0.31%, 4/1297) and PI (1.70%, 22/1297). The majority (95.8%, 69/72) of HIV-1 drug-resistance variants displayed a single drug class resistance mutation, three cases (0.23%) contained mutations associated with drug resistance to both PIs and NNRTIs. One case (0.08%) contained mutations conferred to NRTIs sand NNRTIs. No HIV-1 strain with TDR mutations to triple classes of drugs was found in this study.
Analysis of risk factors of transmitted drug resistance
Taking drug resistance as an independent variable, 8 variables including gender, age at the time of HIV diagnosis, nationality, education level, marital status, occupation, route of infection and genotype as dependent variables, logistic stepwise regression model were used to analyze the risk factors of TDR (Table 2). The results showed that the HIV-1 genotype was a potential influencing factor associated with TDR. CRF07_BC is one of the three main epidemic strains in Sichuan province, which the resistance rate was 3.41% (23/675), the drug resistance rate of CRF01_AE was 4.08% (14/343), CRF08_BC (87%, 24/115) had a higher risk of TDR (odds ratio, 8.32; 95% CI 4.38–15.80; P < 0.05), while other variables had no significant effect on drug resistance (P > 0.05).
Transmission cluster analysis identified 76 clusters containing 205 sequences (205/1297, 15.8%), the cluster sizes ranged between 2 and 24, there were 72 (94.7%) clusters with size < 5, and 4 (5.3%) clusters with size ≥ 5. With regard to the characteristics of individuals involved in clusters in comparison with those out of clusters, CRF01_AE (odds ratio, 2.369; 95% CI 1.659–3.382; P < 0.05), subtype B (odds ratio, 13.723; 95% CI 6.338–29.71; P < 0.05), drug resistance(odds ratio, 0.306; 95% CI 0.106–0.881; P < 0.05) and different levels of education(P < 0.05) were significantly associated to be in clusters (Table 3). We also observed that only 4 drug-resistance cases were included in the transmission cluster, and which were concentrated in two clusters that contained only two individuals. It was worth noting that the clustering rate of Leshan, Zigong and Bazhong were much higher than that of other cities (Fig. 5). The largest cluster was composed of individuals infected with subtype B (n = 24), and Luzhou accounted for the majority of cases (15/24, 62.5%) in the largest molecular cluster, but, most clusters were composed of CRF07_BC and CRF01_AE cases. The other three relatively large clusters were two CRF01_AE and one CRF85_BC. We analyzed the infection routes and found that 98.5% of the putative transmission cluster might occurred in heterosexual individuals, among them 101 (49.3%) cases said they have commercial sexual.
In this study, we conducted a cross-sectional HIV-1 molecular epidemiological study to track the characteristics and distribution of HIV-1 genotypes and TDR in newly diagnosed infections in Sichuan. For the first time, the TDR of HIV-1 infection in Sichuan, which was one of the highest AIDS epidemic provinces, was analyzed.
This results were inconsistent with the national monitoring in 2016  and the survey results in some other areas [24, 25], while were similar to that of the study conducted in Sichuan 2014 . Compared with the 2014 survey, it is found that the PLWH in this study were older and less educated, shown that elderly people with low educational background in our province are the people who need to be paid attention to. Heterosexual transmission is likely to be the main route of infection, in which the proportion of commercial heterosexual is the highest, moreover, transmission cluster analysis showed similar results. It showed that commercial sexual activity is an infection factor worthy to be attached in Sichuan. 67.65% of HIV-infected men were over 50 year-old, who infection came mainly from commercial sexual activity, as such groups often exist with spouse separation and widowhood, their sexual needs are hard to be met, and their cognition of AIDS is insufficient. Note that although men younger than 25 years old account for a low proportion of infected people, most of their infection routes are male transmission, which increases with the decrease in age. There was little difference in transmission routes among women of different ages, mainly noncommercial heterosexual transmission, suggesting that further attention should be paid to the source of infection of such populations for more accurate intervention.
A total of 15 HIV-1 genotypes were detected, and 4 URFs were identified. It showed a much more gene diversity of HIV in Sichuan , which may reflect the active mobility of people across the province . Similar to the survey in 2014 , the main prevailing HIV-1 genotypes in Sichuan remain CRF07_BC and CRF01_AE, but the proportion decreased slightly, which was obviously different from other regions that CRF01_AE (i.e., Anhui, Liaoning and Guangxi), B (i.e., Henan), CRF08_BC (i.e., Yunnan) was the dominant strain [25, 27,28,29]. CRF01_AE was identified in the 1990s as being imported from Thailand to Southwestern China (i.e., Yunnan and Guangxi) in commercial sex worker (CSW) [30, 31]. The CRFs of CRF07_BC and CRF08_BC have common origin which were first circulating in intravenous drug users (IDUs) in Yunnan [7, 32], which was introduced from and from Liangshan to Sichuan . The proportion of CRF08_BC and CRF85_BC among HIV-1 PLWH in Sichuan increased from 4.96% and 3.39% to 8.87% and 5.17% in 2019, and eight different genotypes were identified more than in 2014 . The proportion of URFs also increased significantly, which may be due to the complexity of HIV-1 gene pool caused by long-term epidemic. The genetic diversity of HIV-1 is abundant in Sichuan, in addition to the known CRFs, some URFs were also detected in Sichuan . URFs contributed to the formation of novel CRFs, recently, new CRFs were identified in Sichuan [13, 14, 35].
In our study, most of the infected people entering the transmission network presumed be heterosexual transmission, but, research in Guangxi showed that most HIV-1 infection clusters were MSM , and Liu suggested that factors such as sex, mode of transmission, education level and ethnicity were not significantly correlated with access to the genetic transmission network in Liangshan . It showed that the situation of transmission networks in different regions was various, and targeted prevention and control measures need to be put forward for different regions, especially in cities with high clustering rate. In addition, we would focus on a high clustering rate of subtype B (66.67%, 24/36), which was much higher than other genotypes (10–25%), it is different from the local epidemic tendency. This result showed that timely genotype monitoring is conducive to more accurate prevention and control of HIV transmission.
The results also showed that there were significant differences in the distribution of genotypes in different regions. The frequent reconstitution of the HIV genome will accelerate the evolution of the HIV, which may lead to the emergence of a highly adaptive virus [38, 39]. It was worth noting that the CRFs in Leshan, Chengdu and Luzhou are more complex. The existence of multiple genotypes increases the probability of mutual recombination to form new CRFs/URFs, therefore, it is necessary to further strengthen the monitoring of HIV genotype in Sichuan province, timely grasp the epidemic trends and reduce the generation of CRFs.
Recently, the trend of virus strain diversification in Sichuan province was gradually obvious, CRF55_01B , CRF79_0107 , CRF59_01B [42, 43] were found for the first time in China’s MSM population. In this study, these three CRFs were detected in men, only one female was transmitted by spouse, but only 38.46% (5/13) CRF55_01B, 40% (2/5) CRF79_0107 PLWH were transmitted by MSM, and the rest were heterosexual sexual transmission. There may be a concealed sexual orientation because of the social homosexual cultural identity and discrimination . Some MSM will have sex with women inside and outside the marriage, and increase the difficulty of AIDS prevention and treatment , and suggest that the proportion of homosexual transmission in this province may be underestimated. CRF105_0108 was the genotype found for the first time among heterosexual people in Liangshan Prefecture, Sichuan province, which was found that there was an aggregated epidemic in Luzhou, and one case is also found in Meishan, suggesting that this CRF may have spread in Sichuan province, and 54.54% (6/11) PLWH infection with this CRF were transmitted by commercial heterosexual, indicating that commercial sex workers (CSWs) need to be further found in Luzhou to reduce virus transmission.
Another serious consequence of the high variability of HIV-1 is drug resistance, which is a new threat to epidemic control and can lead to treatment failure and further transmission of resistant HIV. Our study in Sichuan showed that the overall prevalence of TDR was 5.55% among the 1297 participants, belonging to the moderate drug resistance level (5–15%) according to WHO HIVDR surveillance threshold . The result was higher than the currently reported national total prevalence rate and transmissible drug resistance rate of 3–5% [5, 29, 47, 48]. Su  found through meta-analysis that the rate of transmissible drug resistance in Beijing, Henan and Hubei has reached the level of moderate drug resistance. The studies of Zhejiang  and Shanghai [51, 52] showed that the rate of transmissible drug resistance has been greater than 10%, which may be due to the early start of antiretroviral treatment in some areas and large treatment coverage. The pretreatment drug resistance rate was 9.9% in Liangshan Prefecture from 2017 to 2018 , the results were in keeping with our findings that the drug resistance rate of Liangshan was 9.8%. The high rate of TDR in some cities indicated that the monitoring of TDR rate should be carried out in these cities. It could be seen that the longer the antiviral time, the higher the proportion of drug resistance. With the increase in the number of people receiving antiviral treatment, the risk of HIV drug resistance also increases, which may lead to the increase the transmissible drug resistance rate, it directly affects the effect of antiviral therapy. The drug resistance survey results of infected people receiving antiviral treatment showed that the national acquired drug resistance rate was 8.6% , and the drug resistance rate of infected people who fail to inhibit the virus exceeds 50% [55, 56]. TDR was less frequent among individuals involved in clusters (1.95%) compared with those out of clusters (6.23%), which might be explained by several factors. First, we should consider that our analysis is only including sequences of newly-diagnosed individuals, we can’t rule out several infector infected with individuals who failed drug treatment. Secondly, due to the lack of CD4 and new infection detection data of these cases, we cannot judge that all individuals are recently infected. Some drug-resistant individuals may have been infected for a long time, resulting in drug-resistant mutations of virus strains in vivo.
According to our study, TDR mutation frequency to NNRTIs was much higher than NRTIs and PIs, due to NNRTI had a low resistance barrier, and it was more prone to drug resistance . Because of the limited availability of drugs in China, the regimen composed of TDF, 3TC and EFV was currently the most commonly used free first-line therapy. The above three kinds of drugs exhibited the degree of resistance was mainly at a potential low level, therefore, the first-line treatment drugs can still be used continuously in Sichuan province. The most frequent NNRTI-associated DRMs were V106 and E138, which were mainly resulted in low-resistance to RPV and DOR, whereas were K103N in whole country [5, 25] and V179D/E in Shanghai . There were 21 cases of DOR resistance, accounting for a high proportion, DOR is a new NNRTI, which has not provided through the NFATPI in China, suggesting that the use of DOR in the future needs to pay close attention to whether it will have a certain impact on the treatment effect. The proportion of K103N mutation in our study was also high, and it is highly resulted in resistant to EFV and NVP, these two drugs were free NNRTI drugs used in China, therefore, we need to be vigilant about the mutation of K103N. The NRTI-associated DRMs were K65KR/R, K70KR and M41L, among which K65KR/R mutation causes high resistance to d4T, DDI and TDF, but M184V , the most common mutation site of NRTI resistance, is not found in our study, which may be related to the current first-line treatment scheme and the low rate of NRTI resistance in Sichuan province. The main DRMs of PI were M46 and Q58, which produced low or intermediate-resistance to NFV, and low-resistance to TPV/r. The most frequent DRMs M46 was consistent in previous reports [5, 48, 51], and the main DRMs Q58 could be related to the high prevalence of CRF07_BC in Sichuan. Thus, it was considered that, the mutation of drug resistance discovered in the study could be induced by the mutation of exogenous population or non-drug selection pressure.
However, our study has limitations. First, a potential sampling bias, on the one hand, in our study only one cross-sectional study was conducted, and the results could not be used to observe the dynamics of the local HIV epidemic, on the other hand, we could analyze only the samples that had been diagnosed, but those that had been infected but not diagnosed could not be included in the analysis. Second, our analysis of TDR concentrated only on NRTI, NNRTI and PI, not containing integrase inhibitors, which are increasingly infected people began using this drug.
In summary, our study provides a comprehensive molecular epidemiologic dataset to understand the diversity and distribution of HIV genotypes in newly diagnosed HIV-infected patients, and supplements in transmitted drug-resistance in Sichuan province. The distribution of HIV-1 genotypes in Sichuan is more diverse and complex than 2014, with CRF07_BC and CRF01_AE still as the predominant genotypes. The MSM is underrated, arguing for behavior scaling up intervention in this specific population besides the elderly people with heterosexual transmission risk groups. The risk of TDR mutation frequency increased in newly diagnosed patients highlights the significance of genotypic drug resistance monitoring and molecular surveillance of pretreatment HIV-1 drug resistance, the regimen composed of TDF, 3TC and EFV was still currently the preferred solution used free first-line therapy.
Availability of data and materials
The international reference strain sequences determined in this study are available in the [Los Alamos National Laboratory HIV Database] repository, [www.hiv.lanl.gov]. Other datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
NCAIDS NCC. National AIDS venereal disease. 2019.
Hahn B, Shaw G, De Cock KM, Sharp P. AIDS as a zoonosis: scientific and public health implications. Science. 2000;287(5453):607–14.
Hemelaar J, Gouws E, Ghys PD, Osmanov S. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS (Lond, Engl). 2011;25(5):679–89.
Li X, Li W, Zhong P, Fang K, Zhu K, Musa TH, Song Y, Du G, Gao R. Nationwide trends in molecular epidemiology of HIV-1 in China. AIDS Res Hum Retrovir. 2016;32:851–9.
Zhao S, Feng Y, Hu J, Li Y, Zuo Z, Yan J, Zhang J, Cao P, Xu W, Li F et al. Prevalence of Transmitted HIV drug resistance in antiretroviral treatment naïve newly diagnosed individuals in China. SCI REP-UK 2018, 8(1).
Xing H, Liang H, Wan ZY, Chen X, Wei M, PengFei M, Guan Q, Quan Y, Hong KX, Shao YM. Distribution of recombinant human immunodeficiency virus type-1 CRF01_AE strains in China and its sequence variations in the env V3–C3 region. Chin J Prev Med. 2004;5(38):300–4.
Su L, Graf M, Zhang Y, Briesen HV, Xing H, Kostler J, Melzl H, Wolf H, Shao Y, Wagner R. Characterization of a virtually full-length human immunodeficiency virus type 1 genome of a prevalent intersubtype (C/B′) recombinant strain in China. J Virol. 2000;74(23):11367–76.
Zhang L, Chen Z, Cao Y, Yu J, Li G, Yu W, Yin N, Mei S, Li L, Balfe P, et al. Molecular characterization of human immunodeficiency virus type 1 and hepatitis C virus in paid blood donors and injection drug users in China. J Virol. 2004;78(24):13591–9.
Liao L, Xing H, Shang H, Li J, Zhong P. The prevalence of transmitted antiretroviral drug resistance in treatment-naive HIV-infected individuals in China. J Acquir Immune Defic Syndr. 2010;53:S10-14.
Zhao S, Yi F, Hu J, Yang L, Zuo Z, Yan J, Zhang J, Cao P, Xu W, Li F. Prevalence of transmitted HIV drug resistance in antiretroviral treatment nave newly diagnosed individuals in China. Sci Rep-UK. 2018;8:12273.
Su L, Yang H, Hu Y, Ye L, Yuan D, Wei D, Zhou G, Yu J, Liu H, Qin G. Molecular epidemiological study on HIV-1 in Sichuan province in 2014. Chin J AIDS STD. 2016;22(12):946–9.
Su L, Zhu J, Zeng Y, Hou X, Yan T, Li Y, Yang H, Ye L, Yuan D, Wei D, et al. Epidemiological characteristics of HIV-1 strain in Guangyuan city. Chin J AIDS STD. 2020;26(8):831–5.
Dong A, Liu L, Xiao L, Liang S, Li K, Hu J, Kang R, Xing J, Song C, Yin J, et al. First detection of a circulating recombinant form of HIV-1 CRF01_AE/08_BC (CRF105_0108) with drug-resistant mutations in Sichuan, China. AIDS Res Hum Retrovir. 2020;36(7):625–30.
Su L, Wei D, Yang H, Zeng Y, Hu Y, Yuan D, Feng L, Ruan Y, Qin G, Liang S. Identification of a Novel HIV-1 circulating recombinant form (CRF85_BC) in Sichuan, China. AIDS Res Hum Retrovir. 2016;32(9):895–9.
Sudhir K, Glen S, Koichiro T. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33(7):1870.
Liu K, Randal LC, Tandy W, Wu R. RAxML and FastTree: comparing two methods for large-scale maximum likelihood phylogeny estimation. PLoS ONE. 2011;6(11):e27731.
Patiño-Galindo JÁ, Torres-Puente M, Alma Bracho M. The molecular epidemiology of HIV-1 in the Comunidad Valenciana (Spain): analysis of transmission clusters. Sci Rep-UK. 2017;7(1):11584.
Anisimova M, Gascuel AO. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55(4):539–52.
CDC. Guidance: detecting, investigating, and responding to HIV transmission clusters 2018. Prevention TRCF, editor. 2018.
Yuan D, Liu M, Li Y, Jia P, Yang S. Genetic transmission networks of HIV-1 CRF07_BC strain among HIV-1 infections with virologic failure of ART in a minority area of China: a population-based study. BMC Infect Dis. 2020;20(1):631.
Shafer RW, Rhee SY, Pillay D, Miller V, Sandstrom P, Schapiro JM, Kuritzkes DR, Bennett D. HIV-1 protease and reverse transcriptase mutations for drug resistance surveillance. AIDS. 2007;21(2):215–23.
Liu TF, Shafer ARW. Web resources for HIV Type 1 genotypic-resistance test interpretation. Clin Infect Dis. 2006;42(11):1608–18.
Qiao Y, Xu Y, Jiang D, Wang X, Wang F, Yang J, Wei Y. Epidemiological analyses of regional and age differences of HIV/AIDS prevalence in China, 2004–2016. Int J Infect Dis. 2019;81:215–20.
Zhou P, Yu G, Kuang Y, Huang X, Li Y, Fu X, Lin P, Yan J, He X. Rapid and complicated HIV genotype expansion among high-risk groups in Guangdong Province, China. BMC Infect Dis. 2019;19(1):185.
Li L, Sun G, Liang S, Li J, Li T, Wang Z, Liu W, Yang S, Liu Y, Wang X, et al. Different distribution of HIV-1 subtype and drug resistance were found among treatment naive individuals in Henan, Guangxi, and Yunnan province of China. PLoS ONE. 2013;8(10):e75777.
Li X, Fang XD, Mao R, Wang J, Cottrell L, Harris C, Stanton B. HIV/STD risk behaviors and perceptions among rural-to-urban migrants in China. AIDS Educ Prev. 2004;16(6):538–56.
Zhang Z, Dai L, Jiang Y, Feng K, Liu L, Xia W, Yu F, Yao J, Xing W, Sun L, et al. Transmission network characteristics based on env and gag sequences from MSM during acute HIV-1 infection in Beijing, China. Arch Virol. 2017;162(11):3329–38.
Ma N, Chen X, Zhao Y, Kang X, Pan S, Yao W. HIV-1 molecular transmission network among sexually transmitted populations in Liaoning Province, China. Medicine. 2021;100(28):e26640.
Zhang D, Wu J, Zhang Y, Shen Y, Dai S, Wang X, Xing H, Lin J, Han J, Li J, et al. Genetic characterization of HIV-1 epidemic in Anhui Province, China. Virol J. 2020;17(1):17.
Feng Y, He X, Hsi JH, Li F, Li X, Wang Q, Ruan Y, Xing H, Lam TT, Pybus OG, et al. The rapidly expanding CRF01_AE epidemic in China is driven by multiple lineages of HIV-1 viruses introduced in the 1990s. AIDS. 2013;27(11):1793–802.
Cheng H, Zhang J, Capizzi J, Young NL, Mastro TD. HIV-1 subtype E in Yunnan, China. Lancet. 1994;344(8927):953–4.
Yang R, Kusagawa S, Zhang C, Xia X, Ben K, Takebe Y. Identification and characterization of a new class of human immunodeficiency virus type 1 recombinants comprised of two circulating recombinant forms, CRF07_BC and CRF08_BC, China. J Virol. 2003;77(1):685–95.
Li L, Wei D, Hsu WL, Li T, Gui T, Wood C, Liu Y, Li H, Bao Z, Liu S, et al. CRF07_BC strain dominates the HIV-1 epidemic in injection drug users in Liangshan Prefecture of Sichuan, China. AIDS Res Hum Retrovir. 2015;5(31):479–87.
Wang Y, Ma L, Xu W, Li F, Feng Y, Liang S, Shao Y, Zhang F. Near full-length genome identification of a novel HIV-1 unique (B/C) recombinant isolate in Sichuan, China. AIDS Res Hum Retrovir. 2016;32:1229.
Yan W, Wei H, Xu W, Yin Q, Ma L. Identification of a Novel CRF01_AE/CRF07_BC recombinant form in Men Who Have Sex with Men in Sichuan, China. AIDS Res Hum Retrovir. 2016;32(7):718.
Pang X, Tang K, He Q, Huang J, Fang N, Zhou X, Zhu Q, Wu X, Shen Z, Liang S. HIV drug resistance and HIV transmission risk factors among newly diagnosed individuals in Southwest China. BMC Infect Dis. 2021;21(1):160.
Liu L, Dong A, Liao L, Feng Y, Shao Y, Liang S, Ruan Y, Xing H. Survey of pretreatment HIV drug resistance and genetic transmission network analysis among HIV patients in a high drug-use area of Southwest China. Curr HIV Res. 2019;6(17):441–51.
Lau KA, Wong J. Current trends of HIV recombination worldwide. Infect Dis Rep. 2013;5(Suppl 1):e4.
Hemelaar J. The origin and diversity of the HIV-1 pandemic. Trends Mol Med. 2012;18(3):182–92.
Han X, An M, Zhang W, Cai W, Chen X, Takebe Y, Shang H. Genome sequences of a novel HIV-1 circulating recombinant form, CRF55_01B, identified in China. Genome Announc. 2013;1(1):e00050-12.
Li Y, Feng Y, Li F, Xue Z, Hu J, Xing H, Ruan Y, Shao Y. Genome sequence of a novel HIV-1 circulating recombinant form (CRF79_0107) Identified from Shanxi, China. AIDS Res Hum Retrovir. 2017;33(10):1056–60.
Li X, Ning C, Xiang H, Yao Y, Shao Y. Genome sequences of a novel HIV-1 circulating recombinant form (CRF59_01B) identified among MSM in China. AIDS Res Hum Retrovir. 2013;29(10):1395.
Han X, An M, Zhang W, Zhao B, Chu Z, Take Be Y, Shang H. 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):1395.
Liu H, Yang H, Li X, Wang N, Liu H, Wang B, Zhang L, Wang Q, Stanton B. Men who have sex with men and human immunodeficiency virus/sexually transmitted disease control in China. Sex Transm Dis. 2006;33(2):68–76.
Zhang D, Bi P, Lv F, Zhang J, Hiller JE. Changes in HIV prevalence and sexual behavior among men who have sex with men in a northern Chinese city: 2002–2006. J Infect. 2007;55(5):456–63.
WHO. Guidelines for surveillance of HIV drug resistance. Geneva: World Health Organization; 2003.
Li H, Chang S, Han Y, Zhuang D, Li L, Liu Y, Liu S, Bao Z, Zhang W, Song H. The prevalence of drug resistance among treatment-naïve HIV-1-infected individuals in China during pre- and post-2004. BMC Infect Dis. 2016;16(1):605.
Liao L, Xing H, Shang H, Li J, Zhong P, Kang L, Cheng H, Si X, Jiang S, Li X. The prevalence of transmitted antiretroviral drug resistance in treatment-naive HIV-infected individuals in China. J Acquir Immune Defic Syndr. 2010;53(Suppl. 1):10–4.
Su Y, Zhang F, Liu H, Smith MK, Zhu L, Wu J, Wang N. The prevalence of HIV-1 drug resistance among antiretroviral treatment naïve individuals in Mainland China: a meta-analysis. PLoS ONE. 2014;9(10):e110652.
Xu Y, Peng X, Xiuming J, Shujing C, Bin W. Characterization of HIV-1 subtypes and transmitted drug resistance among treatment-naive HIV-infected individuals in Zhejiang, China, 2014–2017. Arch Virol. 2018;163:2233–7.
Wang Z, Zhang M, Li R, Yinzhong L, Jiangrong S. Diversity of HIV-1 genotypes and high prevalence of pretreatment drug resistance in newly diagnosed HIV-infected patients in Shanghai, China. BMC Infect Dis. 2019;19:313.
Zhang F, Li L, Sun M, Sun J, Lu H. An analysis of drug resistance among people living with HIV/AIDS in Shanghai, China. PLoS ONE. 2017;12(2):e165110.
Dong AB, Xiao L, Liang S, Liu L, Xing H. HIV-1 drug resistance and influencing factors among people living with HIV/AIDS before antiretroviral therapy in Liangshan Yi Autonomous Prefecture. Chin J Epidemiol. 2019;40(6):648–53.
Liu P, Xing H, Liao L, Feng Y, Leng X, Wang J, Kan W, Yan J, Li Y, Zuo Z, et al. HIV drug resistance in patients in China’s national HIV treatment programme who have been on first-line ART for at least 9 months. AIDS Res Ther. 2020;17(1):9.
Zuo Z, Liang S, Sun X, Bussell S, Yan J, Kan W, Leng X, Liao L, Ruan Y, Shao Y, et al. Drug resistance and virological failure among HIV-infected patients after a decade of antiretroviral treatment expansion in eight provinces of China. PLoS ONE. 2016;11(12):e166661.
Wu J, Norris J, Liu HX, Li Z, Su YY, Zhu L, Wang N. The prevalence of HIV Drug resistance among treatment-failure individuals and treatment-nave individuals in China: a meta-analysis. Biomed Environ Sci. 2014;27(011):858–71.
Silvia B, Lucas H, Jordan MR, Santiago AR, Collins I, Anne D, Eric D, AnnEmArie W, Theresa A, Borhan A. Clinical impact of pretreatment human immunodeficiency virus drug resistance in people initiating nonnucleoside reverse transcriptase inhibitor containing antiretroviral therapy: a systematic review and meta-analysis. J Infect Dis. 2020;224(3):377–88.
We acknowledge the work of staffs involved in collecting patient blood sample the staff of 21 municipal CDCs. We appreciate all the volunteer patients enrolled in this study.
All methods in this study were operated according to relevant guidelines and regulations, and relevant references were cited in “Material and methods” section.
This study was supported by the Sichuan Science and Technology Project (2017JY0302) and Science and Technology Project of Sichuan Provincial Health Committee (20PJ121).
Ethics approval and consent to participate
This study was reviewed and approved by the ethics committees of Sichuan Center for Disease Prevention and Control. All participants signed written informed consents prior to sample donations.
Consent for publication
The authors declare no competing financial interests exist.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhou, C., Liang, S., Li, Y. et al. Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China. BMC Infect Dis 22, 602 (2022). https://doi.org/10.1186/s12879-022-07576-z
- Molecular epidemiology
- Transmitted drug resistance