Skip to main content

Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China



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.

Peer Review reports


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 [1].

Among the four phylogeny groups that constitute HIV-1 (M, N, O, P), M group is the virulence factor of AIDS pandemic [2]. There were nine subtypes (A, B, C, D, F, G, H, J and K) and at least 118 circulating recombinant forms (CRFs) ( 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%) [3]. A national study showed that 8 subtypes and 21 CRFs have been identified in China since 2013 [4]. CRF01_AE, CRF07_BC, CRF08_BC, subtype B were the four main subtypes, accounting for 89.0% of all HIV-1 infections [5]. CRF01_AE mainly transmitted by sexual route spread from the southeast coast and southwest border to the whole country [6]. 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 [7]. 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 [8].

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% [9], and the rate of TDR was 3.6% in 2015 [10].

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 [1]. In a previous study, we performed a comprehensive investigation of the HIV epidemic in Sichuan province, 2014 [11] 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

Study participants

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 [14].

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 [15]. FastTree was used to estimate an approximately maximum likelihood phylogenetic tree for pol sequences using the GTR + G + I nucleotide substitution model [16], 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 ( was used to review the results. Presented phylogenetic tree by using FigTree v1.4.3 ( 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 [19]. The molecular transmission clusters were deduced by cytoscape 3.7.2 [20].

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

Statistical analysis

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%).

Fig. 1
figure 1

The phylogenetic tree of HIV-1 pol sequences obtained from Sichuan. The maximum-likelihood tree obtained using 1062 bp corresponded to the protease region and part of the reverse transcriptase. A total of 675 HIV-1 sequences branched with CRF07_BC reference sequences (depicted in blue). 342 HIV-1 query sequences branched with CRF01_AE reference sequences (depicted in red). 115 HIV-1 CRF08_BC query sequences were identified (depicted in light blue). 67 HIV-1 CRF85_BC query sequences were identified (depicted in brown). 36 query sequences clustered with HIV-1 subtype B reference sequences (depicted in yellow). Other subtypes and URFs were not depicted color. The reference sequences from the Los Alamos HIV sequence database (

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.

Fig. 2
figure 2

Distribution of HIV-1 genotypes in different cities of Sichuan

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.

Table 1 General distribution characteristics of epidemic strains in Sichuan province

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.

Fig. 3
figure 3

Transmitted drug resistance mutations sites in treatment-naive patients

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.

Fig. 4
figure 4

The resistance of HIV-1 to antiretroviral drugs. 3TC: lamivudine; ABC: abacavir; AZT: zidovudine; d4T: stavudine; ddI: didanosine; FTC: emtricitabine; TDF: tenofovir disoproxil fumarate; EFV: efavirenz; DOR: doravirine; NVP: nevirapine; RPV: Rilpivirine; ETR: etravirine; ATV/r: atazanavir/ritonavir; FPV/r: Fosamprenavir/ritonavir; IDV/r: indinavir/ritonavir; LPV/r: lopinavir/ritonavir, NFV: nelfinavir, SQV/r: saquinavir/ritonavir; TPV/r: tipranavir/ritonavir

Fig. 5
figure 5

TDR and clustering rate of different cities in Sichuan

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

Table 2 Demographic characteristics and factors associated with drug resistance

Transmission clusters

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.

Table 3 Molecular transmission cluster characteristics of Sichuan patients


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 [23] and the survey results in some other areas [24, 25], while were similar to that of the study conducted in Sichuan 2014 [11]. 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 [11], which may reflect the active mobility of people across the province [26]. Similar to the survey in 2014 [11], 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 [33]. 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 [11]. 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 [34]. 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 [36], 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 [37]. 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 [40], CRF79_0107 [41], 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 [44]. Some MSM will have sex with women inside and outside the marriage, and increase the difficulty of AIDS prevention and treatment [45], 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 [46]. 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 [49] 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 [50] 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 [53], 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% [54], 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 [57]. 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 [51]. 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 [5], 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, []. Other datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.


  1. NCAIDS NCC. National AIDS venereal disease. 2019.

  2. Hahn B, Shaw G, De Cock KM, Sharp P. AIDS as a zoonosis: scientific and public health implications. Science. 2000;287(5453):607–14.

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

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

    Article  PubMed  CAS  Google Scholar 

  5. 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).

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

    Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Google Scholar 

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

    Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  15. Sudhir K, Glen S, Koichiro T. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33(7):1870.

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  Google Scholar 

  18. Anisimova M, Gascuel AO. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55(4):539–52.

    Article  PubMed  Google Scholar 

  19. CDC. Guidance: detecting, investigating, and responding to HIV transmission clusters 2018. Prevention TRCF, editor. 2018.

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  22. Liu TF, Shafer ARW. Web resources for HIV Type 1 genotypic-resistance test interpretation. Clin Infect Dis. 2006;42(11):1608–18.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  Google Scholar 

  31. 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  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Google Scholar 

  38. Lau KA, Wong J. Current trends of HIV recombination worldwide. Infect Dis Rep. 2013;5(Suppl 1):e4.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Hemelaar J. The origin and diversity of the HIV-1 pandemic. Trends Mol Med. 2012;18(3):182–92.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

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

    Article  Google Scholar 

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

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  46. WHO. Guidelines for surveillance of HIV drug resistance. Geneva: World Health Organization; 2003.

    Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Google Scholar 

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

    CAS  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  Google Scholar 

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

    PubMed  Google Scholar 

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

    Google Scholar 

Download references


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

Author information

Authors and Affiliations



CZ and LS designed the study; YL, YZ, LL, LY and DY participated in the process of sequences editing and phylogenetic analyses; CZ and LS performed statistical analysis; SL, DY and LS collected sample and the demographic data; CZ and LS participated in the writing process. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ling Su.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare no competing financial interests exist.

Additional information

Publisher's Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: