- 1Nuffield Department of Population Health, Infectious Disease Epidemiology Unit, University of Oxford, Oxford, United Kingdom
- 2School of Public Health, Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 3Liangshan Prefecture Centre for Disease Control and Prevention, Xichang, Sichuan, China
- 4Sichuan Provincial Center for Disease Control and Prevention, Center for AIDS/STD Control and Prevention, Chengdu, Sichuan, China
Introduction: Despite expanded antiretroviral therapy (ART) in China, HIV transmission persists. Liangshan Prefecture is one of the areas in China most severely affected by HIV, with high levels of drug resistance. A deeper understanding of HIV-1 drug resistance can lead to improvements in current treatment policies.
Methods: We conducted an analysis of HIV drug resistance mutations (DRMs) among patients with treatment failure and people newly diagnosed with HIV in Liangshan Prefecture. 8,523 blood samples were collected from people living with HIV with treatment failure and newly diagnosed individuals in all 15 counties and two cities in Liangshan Prefecture between 2021 and 2023.
Results: 43.0% of patients with treatment failure acquired HIV through the heterosexual route, followed by injecting drug use (38.7%), while newly diagnosed individuals mainly acquired HIV through the heterosexual route (86.7%). 95.6% of patients with treatment failure were infected with HIV-1 variant CRF07_BC and 2.7% with CRF08_BC, and newly diagnosed individuals were also main infected with HIV-1 variant CRF07_BC (90.9), followed by CRF08_BC (4.0%) and CRF01_AE (2.5%). The overall prevalence of acquired drug resistance (ADR) among patients with treatment failure was 57.4%. The overall prevalence of pre-treatment drug resistance (PDR) among newly diagnosed individuals was 23.9%. A high prevalence of ADR and PDR (especially high-level resistance) to efavirenz (48.0% vs. 11.1%) and nevirapine (49.6% vs. 11.4%) was found. The main non-nucleoside reverse transcriptase inhibitor (NNRTI)-associated ADR and PDR mutations were K103, V106, and V179. Our findings highlight age <18 years, injecting drug use, and initiation on NNRTI-based regimen as independent risk factors for HIV ADR development. We found minor variants as a risk factor for PDR, and CRF01_AE was associated with a higher risk than CRF07_BC for nucleoside reverse transcriptase inhibitor (NRTI) PDR.
Discussion: Given the high levels of NNRTI ADR and PDR, future clinical treatment plans should minimize the use of NNRTI-based regimens and should instead adopt alternative ART regimens more frequently.
1 Introduction
HIV/AIDS remains a major health problem worldwide. According to the Joint United Nations Program on HIV/AIDS (UNAIDS), approximately 39.9 million people globally were living with HIV, and 30.7 million of them were accessing antiretroviral therapy (ART) in 2023 (1). In China, HIV/AIDS is a significant health issue, with around 1,280,700 people living with HIV (PLWH) at the end of December 2023 (2). The HIV epidemic in China is unevenly distributed, and the province of Sichuan has the highest number of PLWH (3, 4).
Treatment options for HIV/AIDS have expanded significantly in recent years, with ART proving highly effective at reducing disease progression and improving the quality of life for PLWH (5–7). However, drug resistance has been a key challenge in the management and control of HIV-1, impacting the success of antiretroviral treatment, and heightening HIV-related morbidity and mortality (8, 9). Drug resistance mutations (DRMs) are the major cause of antiretroviral treatment failure (10). Since the initiation of the National Free Antiretroviral Treatment Program (NFATP) by the Chinese government in 2003, drug resistance has risen with the escalation of treatment coverage (11). Despite expanded ART access in China, drug resistance could, therefore, threaten ART effectiveness. DRMs to nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs) are complex and common and were the most common cause of first-line antiretroviral treatment failure (12, 13).
Drug resistance can be divided into primary drug resistance and acquired drug resistance (ADR). Primary drug resistance can be divided into pre-treatment drug resistance (PDR) and transmitted drug resistance (TDR) (14). TDR refers to HIV drug resistance present when a new HIV infection is acquired, whether through horizontal or vertical transmission; PDR is a broader term encompassing TDR and any HIV drug resistance present at the initiation of ART and is more commonly used to describe infections of unknown duration and when TDR cannot be conclusively ascertained (15). PDR can result in the failure of first-line regimens, especially when not recognized at the time of ART initiation (16). The prevalence of PDR and ADR in China varies significantly across regions (17), with the reported PDR ranging from 5.5% in Sichuan (3) to 18.3% in Xian (18), and ADR ranging from approximately 46.6% in Chongqing (19) to 81.2% in Henan (20).
Liangshan Prefecture, an autonomous prefecture with the largest population of Yi people in Sichuan Province, is located on the key route through which drugs from the “Golden Triangle” flow into mainland China (21). Liangshan Prefecture is one of the areas most severely affected by HIV/AIDS in China and had a total of 41,623 PLWH in 2019 (22). Within Liangshan Prefecture, prevalence rates of HIV are over 1% in Butuo County, Zhaojue County, Meigu county, Yuexi County, and Jinyang County, which are the counties with the highest HIV prevalence in China (23). In Liangshan Prefecture, CRF07_BC was the main HIV variant of drug users in the past two decades (24). While the HIV epidemic in Liangshan was initially driven by injecting drug use, sexual transmission became the predominant mode of HIV transmission in Liangshan Prefecture from 2014 (25). The genotype of HIV-1 is related to the route of transmission, and there are differences in epidemic scale and distribution characteristics of HIV-1 variants (26, 27). Liangshan Prefecture is an important labor-exporting area, with socioeconomic development and the growth of population mobility. HIV-1 infected populations moving frequently may result in a more wide-spread and complex distribution of HIV-1 genotypes, which may increase the risk of virological failure of antiretroviral treatment, increase the spread of drug-resistant strains, and bring a huge challenge to controlling the HIV epidemic (28). To date, there has not been any large-scale study on HIV drug resistance in Liangshan Prefecture. We conducted an analysis of HIV DRMs among patients with treatment failure and newly diagnosed people in Liangshan, aiming to delineate the characteristics of DRM and their associated risk factors.
2 Materials and methods
2.1 Study population and sample collection
Liangshan Prefecture is located in the southwestern mountainous area of Sichuan Province (Figure 1). The area primarily depends on agriculture and transportation infrastructure is basic, resulting in a relatively low level of economic development. Blood samples were collected for testing from all 15 counties and 2 cities in Liangshan Prefecture between 2021 and 2023. The samples included all patients with treatment failure (defined as having viral load ≥1,000 copies/ml or being assessed by physicians as having poor treatment outcomes) reported to Center for Disease Control and Prevention (CDC) as well as newly diagnosed individuals.
Plasma samples were collected in accordance with standard procedures (29) by laboratory personnel of the local CDC and were transported to the Sichuan Provincial Center for Disease Control and Prevention for HIV-1 drug resistance testing (30).
Patients’ demographic information—including sex, age, marital status, ethnicity, education level, and transmission route—was collected from the National HIV/AIDS Comprehensive Response Information Management System, a web-based real-time database managed by the National Center for AIDS/STD Control and Prevention (NCAIDS) of the Chinese Center for Disease Control and Prevention (CDC) (31).
2.2 HIV-1 gene amplification and drug resistance analysis
Viral nucleic acid was obtained from 200 μL plasma of PLWH by extraction machines (MagNA Pure LC system, Roche, Branchburg, NJ). HIV-1 pol sequences were amplified and sequenced. Reverse Transcription-Polymerase Chain Reaction (RT-PCR) was used to amplify the full-length protease gene and the first 300 codons of the reverse transcriptase gene. Two rounds of PCR amplification were used, following the HIV-1 Genotype Drug Resistance Detection and Quality Assurance Guidelines (2013 Edition) (3, 32). The target bands were subjected to 1% agarose gel electrophoresis for validation, and the amplified product was sequenced. Sequences were then spliced and edited, and were submitted to the drug resistance database of Stanford University1 for resistance mutation analysis.
For HIV-1 subtyping, the HIV sequences were aligned with HIV-1 reference sequences available in the Los Alamos National Laboratory database.2 Multiple alignments were made automatically using Mega version 7.0 with minor manual adjustments.
Each drug resistance mutation is assigned a drug penalty score. The total score for a drug is derived by adding the scores for all mutations associated with resistance to that drug to infer 1 of 5 levels of resistance: susceptible, potential low-level resistance, low-level resistance, intermediate resistance, and high-level resistance (33–35).
2.3 Statistical analysis
Demographic characteristics of the study participants were summarized using frequencies and percentages. Multivariate logistic regression was employed to assess the risk factors for HIV PDR and ADR. Unadjusted odds ratios (OR), adjusted odds ratios (AOR), and their respective 95% confidence intervals (CI) were calculated. Significance was established at p < 0.05. IBM SPSS 22 software was utilized for statistical analysis.
3 Results
3.1 Participant characteristics
A total of 8,523 participants were enrolled in this study (Figure 2). After excluding samples with sequencing failure and missing information, the final analysis included 3,924 individuals with treatment failure for ADR analysis and 1,999 newly diagnosed individuals for PDR analysis.
Characteristics of patients with treatment failure are summarized in Table 1. Of the 3,924 individuals with treatment failure in this study, the majority were male (71.5%). 63.8% were aged 25–44 years and 92.3% were of Yi ethnicity. The most common mode of HIV-1 transmission was heterosexual contact (43.0%), followed by injecting drug use (38.7%). Regarding education, 57.2% were illiterate, and 33.6% attended primary schooling. The primary occupation was farming, comprising 75.1%. The most prevalent variant of HIV-1 was CRF07_BC (95.6%), followed by CRF08_BC (2.7%) and CRF01_AE (0.8%). ART regimens mainly consisted of TDF + 3TC + EFV (73.4%), AZT + 3TC+ NVP (11.0%) and AZT + 3TC+ EFV (9.3%).
Among the 1,999 newly diagnosed individuals in the study, 53.5% were male, and 50.6% were aged 25–44 years (Table 2). The majority of newly diagnosed individuals acquired HIV through heterosexual contact (86.7%). Fewer newly diagnosed patients were illiterate (45.9%) and more completed primary schooling (38.8%), compared to treatment failure patients. Most newly diagnosed patients were farmers (83.7%). Compared to treatment failure patients, the HIV-1 variants among newly diagnosed patients were more diverse, with a lower proportion of CRF07_BC (90.9%), and a higher prevalence of CRF08_BC (4.0%) and CRF01_AE (2.5%).
3.2 ADR and PDR prevalence
Of 3,924 pol sequences from individuals with treatment failure, 57.4% had one or more resistance mutations (Figure 3). The prevalence of ADR to NNRTI was 54.9%, followed by ADR to NRTI at 23.8%, and ADR to PI at 3.1%. Specifically, resistance to only NNRTI, NRTI, and PI were 31.8, 1.3, and 1.1%, respectively. The prevalence of dual-class resistance to NNRTI+NRTI was 21.2%, and the occurrence of triple-drug resistance (NNRTI+NRTI+PI) was relatively low at 1.2% (Figure 3).

Figure 3. Frequency of drug resistance to ART classes. ADR, Acquired drug resistance; PDR, Pre-treatment drug resistance; PI, Protease inhibitor; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor; DRM, Drug Resistance Mutation.
The prevalence of PDR among 1999 newly diagnosed individuals was 23.9%. The prevalence of PDR to NNRTI was 20.5%, followed by PDR to NRTI at 3.6%, and PDR to PI at 2.8% (Figure 3). Resistance to only NNRTI, NRTI, and PI were 17.7, 1.6, and 1.9%, respectively. The prevalence of PDR to only NNRTI (17.7%) was lower than of ADR, while the prevalence of PDR to only NRTI (1.6%) and PI (1.9%) was higher than of ADR. The most prevalent dual-class resistance was against NNRTI+NRTI, accounting for 2.0% of newly diagnosed individuals. The occurrence of triple-class (NNRTI+NRTI+PI) PDR was low at 0.1% (Figure 3).
3.3 Resistance to individual antiviral drugs
Figure 4 illustrates the predicted resistance of 3,924 individuals with treatment failure from Liangshan to 20 antiretroviral drugs. For NRTIs, the virus exhibited high/intermediate level resistance to FTC, 3TC, DDI, ABC, D4T, TDF and AZT at frequencies of 22.4, 21.0, 12.1, 10.5, 10.4, 8.0 and 4.0%, respectively. For NNRTIs, NVP, EFV, DOR, RPV and ETR had high/intermediate level frequency of 50.6, 50.1, 20.7, 11.8 and 8.6%, respectively. NVP (54.6%) and EFV (54.5%) were the most common NNRTI drugs with ADR. Conversely, frequencies of high/intermediate level to PI drugs remained low, all below 1%.

Figure 4. ART drug resistance among HIV patients with treatment failure in Liangshan Prefecture, China. ABC, Abacavir; AZT, Azidothymidine; D4T, Stavudine; DDI, Didanosine; FTC, Emtricitabine; 3TC, Lamivudine; TDF, Tenofovir disoproxil fumarate; DOR, Doravirine; EFV, Efavirenz; ETR, etravirine; NVP, Nevirapine; RPV, Rilpivirine; ATV, Atazanavir; DRV, Darunavir; FPV, Fosamprenavir; IDV, Indinavir; LPV/r, lopinavir/ritonavir; NFV, Nelfinavir; SQV, Saquinavir; TPV, Tipranavir; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor; PI, Protease inhibitor.
In contrast, newly diagnosed individuals showed lower resistance frequencies to 20 antiretroviral drugs than treatment failure patients (Figure 5). For NRTIs, the virus exhibited high/intermediate level resistance to FTC, 3TC, DDI, ABC, AZT, D4T and TDF at frequency of 1.6, 1.5, 1.0, 0.9, 0.8 and 0.6%, respectively. For NNRTIs, NVP, EFV, DOR, RPV and ETR had high/intermediate level frequency of 12.5, 11.9, 2.9, 1.7 and 0.8%, respectively. NVP (20.1%) and EFV (19.9%) were the most common NNRTI drugs with PDR. The frequencies of high/intermediate levels among newly diagnosed individuals to PIs drugs were all below 1%.

Figure 5. ART drug resistance among people newly diagnosed with HIV in Liangshan Prefecture, China. ABC, Abacavir; AZT, Azidothymidine; D4T, Stavudine; DDI, Didanosine; FTC, Emtricitabine; 3TC, Lamivudine; TDF, Tenofovir disoproxil fumarate; DOR, Doravirine; EFV, Efavirenz; ETR, etravirine; NVP, Nevirapine; RPV, Rilpivirine; ATV, Atazanavir; DRV, Darunavir; FPV, Fosamprenavir; IDV, Indinavir; LPV/r, lopinavir/ritonavir; NFV, Nelfinavir; SQV, Saquinavir; TPV, Tipranavir; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor; PI, Protease inhibitor.
3.4 Drug resistance mutations
Among the 3,924 successfully amplified samples from individuals with treatment failure, 80 NRTI-associated, 78 NNRTI-associated, and 43 PI-associated drug resistance mutations were identified (Figures 6, 7). The most prominent NRTI resistance mutations included M184I/V (89.9%), followed by K65E/R (31.2%), K70E/N/Q/R/T (26.9%), and K219E/N/Q/R (17.6%). The most common NNRTI-associated drug resistance mutations were K103N/S, accounting for 79.9%, followed by the V179D/E/F/L and V106A/M mutations, accounting for 30.9 and 23.7%, respectively. The most common PI-associated drug resistance mutations were Q58E (49.6%), followed by M46V (13.8%), I54T/V/M (13.8%) and L10F (13.8%; Figure 6).

Figure 6. DRMs to the main ART classes among HIV patients with treatment failure in Liangshan Prefecture, China. PI, Protease inhibitor; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor.

Figure 7. DRMs to the main ART classes among people newly diagnosed with HIV in Liangshan Prefecture, China. PI, Protease inhibitor; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor.
The number of drug-resistance-associated mutations among newly diagnosed individuals was fewer than treatment failure individuals. 34 NRTI-associated, 56 NNRTI-associated, and 20 PI-associated drug resistance mutations were identified. The most common NRTIs mutations were M184I/V (37.5%) and S68G/N/R (27.8%). The most prominent NNRTIs resistance mutations included K103N/S (46.2%) and V179D/E/F/L (40.6%). The most common PI mutations were Q58E (58.2%), followed by M46I/L (18.2%; Figure 7).
Analysis of NRTI mutations showed that S68-associated mutations were detected in only 5.7% of treatment failure individuals but in 27.8% of newly diagnosed individuals. M184-associated mutations were present in 89.9% of treatment failure participants, but in 37.5% of newly diagnosed individuals. For NNRTI mutations, the proportions of G190-associated and V106-associated mutations in newly diagnosed individuals were smaller than treatment failure participants, which only occurred at frequencies of <6%. For PI variants, I54-associated and V82-associated mutations were observed in both treatment failure and newly diagnosed participants, though the proportion of newly diagnosed individuals with these mutations was significantly smaller (1.8% vs. 13.8 and 1.8% vs. 13.0%) than in treatment failure participants (Figures 6, 7).
3.5 Risk factors for drug resistance mutations
Univariable and multivariable logistic regression analyses were performed to identify factors associated with DRM. The results showed that in our population, gender was not a predictor of ADR [OR (95% CI): 0.93 (0.81–1.06)], whereas there was a negative association between the presence of ADR and individuals treated with NRTI+PI containing ART regimens [AOR (95% CI): 0.64 (0.45–0.89)], compared with NRTI+NNRTI based ART regimen (Table 3). Compared with heterosexual transmission, injecting drug use was a risk factor for ADR [AOR (95% CI): 1.29 (1.11–1.49); Table 3]. In multivariable logistic regression, all of the variables were not significantly associated with PI ADR (Table 4). In the multivariable logistic regression model of NRTI ADR, age, transmission route, marital status, and ART regimen were significantly associated with DRM (Table 5). In the multivariate model, the transmission route and ART regimen were independently associated with NNRTI ADR (Table 6). Compared with heterosexual transmission, injecting drug use had a higher prevalence of NNRTI ADR [AOR (95% CI): 1.29 (1.12–1.50)]. NNRTI ADR was found less frequently in patients receiving ART regimens with NRTI+PI [AOR (95% CI): 0.51 (0.36–0.72)] than in those with NRTI+NNRTI (Table 6).

Table 3. Risk factors for ADR among HIV patients with treatment failure in Liangshan Prefecture, China.

Table 4. Risk factors for PI ADR among HIV patients with treatment failure in Liangshan Prefecture, China.

Table 5. Risk factors for NRTI ADR among HIV patients with treatment failure in Liangshan Prefecture, China.

Table 6. Risk factors for NNRTI ADR among HIV patients with treatment failure in Liangshan Prefecture, China.
In multivariate analysis of PDR (Table 7), factors that were significantly associated with PDR were age ≥ 65 years [AOR (95% CI): 0.55 (0.31–0.97)] compared to age 25–44 years, HIV-1 minor variants [AOR (95% CI): 4.35 (1.78–10.64)] compared CRF07_BC, recent CD4+ T lymphocyte count 200–350 pcs/μl [AOR (95% CI): 1.80 (1.31–2.49)] and <200 pcs/μl [AOR (95% CI): 1.34 (1.05–1.71)], compared with >350 pcs/μl. Gender, ethnicity, transmission route, education, marital status, occupation and HIV-1 variants were not significantly associated with PDR. Multivariate logistic regression analysis revealed that transmission route and occupation were the most significant factors associated with PI PDR (Table 8). Mother to child transmission [AOR (95% CI): 0.06 (0.01–0.98)] was associated with PI PDR compared to heterosexual transmission, while unemployment [AOR (95% CI): 3.69 (1.06–12.81)] and children and students [AOR (95% CI): 7.79 (1.16–52.61)] were associated with a higher risk of PI PDR, compared to farmers. In the multivariate model, the following four factors were independently correlated with NRTI PDR (Table 9): other ethnicity [AOR (95% CI): 3.92 (1.11–13.91)] had higher risk than Yi ethnicity, the rate of NRTI PDR among high school or technical secondary school [AOR (95% CI): 3.55 (1.34–9.43)] was higher than illiterate people, unmarried [AOR (95% CI): 1.95 (1.10–3.46)] had higher risk than married or cohabiting people, and HIV-1 variant CRF01_AE [AOR (95% CI): 3.99 (1.65–9.67)] was associated with a higher risk than CRF07_BC. Some of the factors that were independently associated with an increased risk of NNRTI PDR (Table 10), included education, HIV-1 variants and recent CD4+ T lymphocyte count. Primary school education [AOR (95% CI): 0.76 (0.59–0.98)] was associated with a lower risk than illiterate people, minor HIV-1 variants [AOR (95% CI): 5.13 (2.11–12.48)] with higher risk than CRF07_BC, and recent CD4+ T lymphocyte count<200 pcs/μl [AOR (95% CI): 1.53 (1.09–2.15)] with higher risk than>350 pcs/μl of NNRTI PDR.

Table 8. Risk factors for PI PDR among people newly diagnosed with HIV in Liangshan Prefecture, China.

Table 9. Risk factors for NRTI PDR among people newly diagnosed with HIV in Liangshan Prefecture, China.

Table 10. Risk factors for NNRTI PDR among people newly diagnosed with HIV in Liangshan Prefecture, China.
3.6 Trends of drug resistance mutations
The prevalence of ADR initially decreased from 60.9% in 2021 to 53.8% in 2022, but then increased to 57.7% in 2023, whereas PDR rose sharply from 21.1 to 32.1%, before declining to 25.9% (Figure 8; Table 11). Trends for NNRTI ADR and NNRTI PDR closely matched the overall trends of ADR and PDR. The prevalence of ADR and PDR associated with PI remained at low levels throughout the study period. PDR linked to NRTI increased from 2.1 to 14.1%, but then declined back to 4.1%. In contrast, ADR associated with NRTI displayed a continuous upward trend from 21.1% in 2021 to 29.3% in 2023.

Figure 8. The trend of prevalence of DRMs in Liangshan Prefecture, China. ADR, acquired drug resistance; PDR, pre-treatment drug resistance; PI, Protease inhibitor; NRTI, Nucleoside Reverse Transcriptase Inhibitor; NNRTI, Non-Nucleoside Reverse Transcriptase Inhibitor; DRM, Drug Resistance Mutation.

Table 11. The occurrence of drug resistance to ART classes over the years in Liangshan Prefecture, China.
4 Discussion
In this study, we used a large sample size and the latest data to investigate HIV-1 drug resistance in Liangshan Prefecture for the first time. The proportion of males in newly diagnosed individuals was lower than other newly diagnosed groups in China (36–38), which showed an increasing trend of HIV infection in females in Liangshan. For transmission routes, the proportion of sexual transmission in the newly diagnosed individuals is significantly higher than in the treatment failure patients. Compared to treatment failure patients, the proportion of CRF07_BC in newly diagnosed individuals was lower, while other HIV-1 variants were more common. Our study showed a more widespread distribution of HIV-1 genotypes and that the heterosexual route comprised the predominant transmission route in Liangshan Prefecture. The overall prevalence of ADR and PDR was significantly higher in Liangshan Prefecture, compared to the overall situation in China (34, 39). These findings provide important guidance for the adjustment of future HIV treatment regimens and contribute to the body of literature on drug resistance.
ADR prevalence (57.4%) was lower than other treatment failure cohorts in China and other countries (40–42). This may be because of the government’s increased focus on ART treatment adherence in Liangshan Prefecture in recent years, accompanied by free treatment provision strategies. The overall prevalence of PDR was higher than in other studies in China and other countries (20, 38, 43–45) and also higher than in a previous study in Liangshan (37). According to the WHO criteria, PDR prevalence of below 5%, 5–15%, and above 15% are considered low, moderate, and high levels of HIV-1 drug resistance (46). The PDR prevalence in Liangshan Prefecture was relatively high (23.9%), highlighting the need to pay more attention to treating newly diagnosed individuals. The proportion of NNRTI PDR, NRTI PDR, and PI PDR is higher than in other studies of newly diagnosed individuals in China (36, 38, 44). The results showed that DRMs against NNRTIs and NRTIs were the main DRMs with high-level drug resistance in both treatment failure patients and newly diagnosed individuals, signifying a substantial concern for dual-class resistance.
Consistent with studies conducted in other areas of China (40), the prevalence of ADR to NNRTIs was substantially higher than that of NRTIs and PIs among patients with ART failure in Liangshan Prefecture. In this study, the main regimens were TDF + 3TC + NVP/EFV or AZT + 3TC + NVP/EFV, which account for more than 95.0%. Under the pressure of drug selection, the prevalence of ADR and PDR in Liangshan Prefecture was mainly driven by NNRTI-associated DRMs. The top NNRTI-associated ADR and PDR mutations were K103, V106, and V179, which cause high-level resistance to NVP and EFV (47). This is different from previous analyses of HIV-1 drug resistance in newly diagnosed HIV-infected patients in Sichuan Province (V106 and E138) (3) and the overall situation in China (K103 and Y181) (34). The study also found that FTC and 3TC were the most critical NRTI drugs responsible for high-level drug resistance, which is slightly different from other studies in China (ABC, FTC, and 3TC) (40). This shows the specificity of HIV DRM in Liangshan Prefecture and needs more attention in future clinical decision-making regarding ART regimens. The prevalence of the M184 mutation was significant in treatment failure patients and present at a lower level in newly diagnosed individuals, consistent with other observations. This may also pose a risk to the future use of a long-acting regimen combination of the capsid maturation inhibitor lenacapavir with NRTTI islatravir (48, 49). Our report suggests a higher prevalence of K65 mutations, which confer intermediate to high-level TFV and ABC resistance (50). This is higher than observed in high-income settings (51) and South Africa (48). The main ADR and PDR for PI included Q58 and M64 mutations. Resistance to PI was rare, which correlates with the late introduction of PIs into southwest China, the short duration of clinical application, and the high resistance barrier (19).
Our findings highlight age <18 years, injecting drug use, and initiation on NNRTIs-based regimen as independent risk factors for HIV ADR development. A previous study showed that HIV-infected children and adolescents exhibit markedly reduced virological suppression compared to adults (52). For PLWH under the age of 18 in Liangshan Prefecture, it is crucial to focus on the effectiveness of ART while paying special attention to drug resistance. Medications should be adjusted promptly for HIV-infected children and adolescents who experience poor treatment outcomes. Among the current PLWH in Liangshan Prefecture, injecting drug users represent a significant proportion. Due to the high rates of drug resistance in this group, it is essential to enhance the monitoring of treatment effectiveness. Additionally, the risks associated with initial regimens based on NNRTIs suggest that ART regimens with higher genetic barriers to drug resistance should be selected as first-line therapy.
Different HIV-1 genotypes have different propensities to develop drug resistance and exhibit variable DRM profiles (38). We found minor variants as a risk factor for PDR, and CRF01_AE had a higher risk than CRF07_BC for NRTI PDR. Due to the trend of the proportion of CRF07_BC gradually decreasing and other HIV-1 variants increasing in newly diagnosed individuals, the findings emphasize the importance of continuous monitoring of HIV subtypes and resistance patterns, especially in the newly diagnosed group. Additionally, the incidence of NNRTI PDR was significantly higher in newly diagnosed individuals with recent CD4 T-cell counts <200/μl than in the >350/μl group. This highlights the importance of improving HIV test coverage to reduce late diagnoses. For future clinical treatment plans for newly diagnosed individuals, it is recommended to consider the results of genotypic and drug resistance testing.
The slight decrease in ADR over the last three years may be linked to the increased focus and strategic changes made recently, further highlighting that improving management and prioritizing drug resistance can lead to lower resistance rates. The increase in PDR indicates that our future efforts to address drug resistance should extend beyond just the treatment population. It is essential to focus more on newly diagnosed individuals and enhance pre-treatment drug resistance testing. Over the past three years (10), changes in resistance to NNRTIs, NRTIs, and PIs in Liangshan Prefecture have shown a relatively slow increase. However, there is a general upward trend that aligns with recent global patterns. This situation underscores the serious challenges posed by drug resistance and highlights the urgent need to bolster research and control efforts.
Our study has several key strengths. First, blood samples of all treatment failure patients and most newly diagnosed individuals were collected and analyzed from all counties and cities in Liangshan Prefecture. Second, the study was conducted from 2021 to 2023 and provides a recent and comprehensive overview of both ADR and PDR in Liangshan Prefecture. Third, our sample population’s characteristics and proportions closely align with the overall PLWH population in the region, making it representative of the situation in Liangshan Prefecture.
Our study has some limitations. First, as a cross-sectional survey, this study only captured resistance at a single time point rather than longitudinally during an infection course. Second, our study lacked information on ART treatment adherence, which prevented us from considering the impact of adherence on drug resistance. Third, the variety of treatment regimens included in our sample is too extensive, which limits our ability to perform direct comparisons between the regimens.
5 Conclusion
In conclusion, we reported the prevalence of ADR and PDR for the first time in Liangshan Prefecture. While current therapies remain largely effective, the high level of PDR requires ongoing surveillance and possibly adjustments in treatment protocols. It is necessary to pay attention to the obstacle caused by PDR in the prevention and control of HIV epidemics.
Data availability statement
The datasets presented in this article are not readily available due to privacy or ethical restrictions. Requests to access the datasets should be directed to cm9uZ2ZyeUAxNjMuY29t.
Ethics statement
The studies involving human participants were approved by Ethics Committee of the Chengdu University of TCM. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants.
Author contributions
RP: Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Resources, Software, Visualization, Writing – original draft. YuZ: Data curation, Writing – review & editing. CJ: Data curation, Project administration, Resources, Writing – review & editing. GY: Data curation, Project administration, Resources, Writing – review & editing. LS: Software, Writing – review & editing. JW: Data curation, Investigation, Writing – review & editing. LX: Data curation, Investigation, Writing – review & editing. YW: Data curation, Investigation, Writing – review & editing. MS: Data curation, Investigation, Writing – review & editing. JL: Data curation, Writing – review & editing. YiZ: Data curation, Writing – review & editing. JH: Conceptualization, Methodology, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Natural Science Foundations of China (no. 72304050), Liangshan state science and technology plan key research and development project (nos. 23ZDYF0112, 23ZDYF0025, and 22ZDYF0125), Xinglin scholars research promotion Program of Chengdu University of Traditional Chinese Medicine (no. 030058026), and the China Scholarship Council (no. 202308510244) for overseas research.
Acknowledgments
We sincerely appreciate all the volunteer patients who participated in this study. We also extend our gratitude to everyone involved for their valuable contributions and hard work.
Conflict of interest
The authors declare that this research was conducted without any commercial or financial relationships that could be interpreted as potential conflicts of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1. UNAIDS. Global HIV & AIDS statistics 2024. (2024). Available from: https://www.unaids.org/en/resources/fact-sheet (Accessed February 21, 2025).
3. Zhou, C, Liang, S, Li, Y, Zhang, Y, Li, L, Ye, L, et al. Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China. BMC Infect Dis. (2022) 22:602. doi: 10.1186/s12879-022-07576-z
4. Yang, H, Li, Y, Xu, M, Hu, Y, Yuan, F, Liu, L, et al. The Update of HIV-1 Prevalence and Incidence and Spatio-Temporal Analyses of HIV Recent Infection Among Four Sub-Groups in Sichuan, China During Surveillance Period Between 2016 and 2022. Infect Drug Resist. (2023) 16:6535–48. doi: 10.2147/IDR.S428744
5. Bekker, LG, Alleyne, G, Baral, S, Cepeda, J, Daskalakis, D, Dowdy, D, et al. Advancing global health and strengthening the HIV response in the era of the Sustainable Development Goals: the International AIDS Society-Lancet Commission. Lancet. (2018) 392:312–58. doi: 10.1016/S0140-6736(18)31070-5
6. Sulistina, DR, Martini, S, Prasetyo, B, Rahman, FS, Adji, AS, Li, CY, et al. A systematic review and meta-analysis of HIV transmission risk behaviors, genetic variations, and antiretroviral (ARV) resistance in LGBT populations. J Public Health Res. (2024) 13:1–11. doi: 10.1177/22799036241239464
7. Volkow, ND, and Montaner, J. The urgency of providing comprehensive and integrated treatment for substance abusers with HIV. Health Aff (Millwood). (2011) 30:1411–9. doi: 10.1377/hlthaff.2011.0663
8. Bertagnolio, S, Hermans, L, Jordan, MR, Avila-Rios, S, Iwuji, C, Derache, A, et al. 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. (2021) 224:377–88. doi: 10.1093/infdis/jiaa683
9. Collier, DA, Monit, C, and Gupta, RK. The Impact of HIV-1 Drug Escape on the Global Treatment Landscape. Cell Host Microbe. (2019) 26:48–60. doi: 10.1016/j.chom.2019.06.010
10. Chen, X, Chen, X, and Lai, Y. Development and emerging trends of drug resistance mutations in HIV: a bibliometric analysis based on CiteSpace. Front Microbiol. (2024) 15:1374582. doi: 10.3389/fmicb.2024.1374582
11. Hao, J, Zheng, S, Gan, M, Dong, A, Kang, R, Li, M, et al. Changing Proportions of HIV-1 Subtypes and Transmitted Drug Resistance Among Newly Diagnosed HIV/AIDS Individuals— China, 2015 and 2018. China CDC Weekly. (2021) 3:1133–8. doi: 10.46234/ccdcw2021.251
12. Boffito, M, Waters, L, Cahn, P, Paredes, R, Koteff, J, Van Wyk, J, et al. Perspectives on the Barrier to Resistance for Dolutegravir + Lamivudine, a Two-Drug Antiretroviral Therapy for HIV-1 Infection. AIDS Res Hum Retrovir. (2020) 36:13–8. doi: 10.1089/aid.2019.0171
13. Gupta, RK, Gregson, J, Parkin, N, Haile-Selassie, H, Tanuri, A, Andrade Forero, L, et al. HIV-1 drug resistance before initiation or re-initiation of first-line antiretroviral therapy in low-income and middle-income countries: a systematic review and meta-regression analysis. Lancet Infect Dis. (2018) 18:346–55. doi: 10.1016/S1473-3099(17)30702-8
14. Yuan, D, Liu, Y, Zhou, Y, Shi, L, Chen, J, Lu, J, et al. Men who have sex with men is the high-risk drug resistance population: A meta-analysis of HIV-1 drug resistance profiles and trends in China. J Clin Pharm Ther. (2022) 47:1729–37. doi: 10.1111/jcpt.13772
15. Trebelcock, WL, Lama, JR, Duerr, A, Sanchez, H, Cabello, R, Gilada, T, et al. HIV pretreatment drug resistance among cisgender MSM and transgender women from Lima, Peru. J Int AIDS Soc. (2019) 22:e25411. doi: 10.1002/jia2.25411
16. Milne, RS, Silverman, RA, Beck, IA, Jennifer Mckernan-Mullin, WD, Sibley, TR, Dross, S, et al. Minority and majority pre-treatment HIV-1 drug resistance associated with failure of 1st-line NNRTI ART in Kenyan women. AIDS. (2019) 33:941–51. doi: 10.1097/QAD.0000000000002134
17. Shi, YZ, Huang, HH, Wang, XH, Song, B, Jiang, TJ, Yu, MR, et al. Retrospective Study on Genetic Diversity and Drug Resistance among People Living with HIV at an AIDS Clinic in Beijing. Pharmaceuticals (Basel). (2024) 17:115. doi: 10.3390/ph17010115
18. Xia, H, Jin, J, Ba, H, Zhang, Y, Li, J, Guo, R, et al. Genetic Diversity and Characteristics of Drug Resistance Among Treatment-Naive People Living with HIV in Xi'an. China Drug Des Devel Ther. (2023) 17:1485–94. doi: 10.2147/DDDT.S406255
19. Gao, W, Zhou, G, Li, M, Wang, P, Li, J, and Deng, R. HIV drug resistance: analysis of viral genotypes and mutation loci in people living with HIV in Chongqing, China (2016-2023). AIDS Res Ther. (2024) 21:62. doi: 10.1186/s12981-024-00646-z
20. Liu, J, Chen, Z, Wei, S, Ma, J, Zhang, X, Zhao, S, et al. Drug resistance mutations among people living with HIV with treatment failure in Henan Province. China Chin Med J (Engl). (2023) 136:2744–6. doi: 10.1097/CM9.0000000000002881
21. Kanglin, RY, Hong, K, He, YMS, Liu, S, Zhou, F, and Qin, G. Community-Based Survey of HIV Transmission Modes Among Intravenous Drug Users in Sichuan, China. Sex Transm Dis. (2004) 31:623–7. doi: 10.1097/01.olq.0000140018.24262.4a
22. Liu, L, Dong, A, Liao, L, Feng, Y, Shao, Y, Liang, S, et al. 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. (2020) 17:441–51. doi: 10.2174/1570162X17666191128101426
23. Yuan, D, Liu, M, Li, Y, Jia, P, Su, L, Ye, L, et al. 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:631. doi: 10.1186/s12879-020-05347-2
24. Li, L, Wei, D, Hsu, WL, Li, T, Gui, T, Wood, C, et al. CRF07_BC Strain Dominates the HIV-1 Epidemic in Injection Drug Users in Liangshan Prefecture of Sichuan. China AIDS Res Hum Retroviruses. (2015) 31:479–87. doi: 10.1089/aid.2014.0120
25. Wang, QX, Zhang, SH, Yang, SJ, Liao, Q, Yu, G, Yao, YN, et al. Characteristics of transmission route among the reported HIV infections during 2011-2014 in Liangshan Prefecture. Chin J AIDS STD. (2017) 23:719–22.
26. Nchinda, N, Elangovan, R, Yun, J, Dickson-Tetteh, L, Kirtley, S, Hemelaar, J, et al. Global associations of key populations with HIV-1 recombinants: a systematic review, global survey, and individual participant data meta-analysis. Front Public Health. (2023) 11:1153638. doi: 10.3389/fpubh.2023.1153638
27. Hemelaar, J. Implications of HIV diversity for the HIV-1 pandemic. J Inf Secur. (2013) 66:391–400. doi: 10.1016/j.jinf.2012.10.026
28. Yuan, D, Liu, M, Jia, P, Li, Y, Huang, Y, Ye, L, et al. Prevalence and determinants of virological failure, genetic diversity and drug resistance among people living with HIV in a minority area in China: a population-based study. BMC Infect Dis. (2020) 20:443. doi: 10.1186/s12879-020-05124-1
29. Xing, H, Ruan, Y, Hsi, JH, Kan, W, Liao, L, Leng, X, et al. Reductions in virological failure and drug resistance in Chinese antiretroviral-treated patients due to lamivudine-based regimens, 2003-12. J Antimicrob Chemother. (2015) 70:2097–103. doi: 10.1093/jac/dkv078
30. Su, L, Liang, S, Hou, X, Zhong, P, Wei, D, Fu, Y, et al. Impact of worker emigration on HIV epidemics in labour export areas: a molecular epidemiology investigation in Guangyuan, China. Sci Rep. (2018) 8:16046. doi: 10.1038/s41598-018-33996-6
31. Mao, Y, Wu, Z, Poundstone, K, Wang, C, Qin, Q, Ma, Y, et al. Development of a unified web-based national HIV/AIDS information system in China. Int J Epidemiol. (2010) 39:ii79–89. doi: 10.1093/ije/dyq213
32. Su, L, Wei, D, Yang, H, Zeng, Y, Hu, Y, Yuan, D, et al. Identification of a Novel HIV-1 Circulating Recombinant Form (CRF85_BC) in Sichuan. AIDS Res Hum Retrovir. (2016) 32:895–9. doi: 10.1089/aid.2016.0053
33. Liu, TF, and Shafer, RW. Web Resources for HIV Type 1 Genotypic-Resistance Test Interpretation. Clin Infect Dis. (2006) 42:1608–18. doi: 10.1086/503914
34. Liu, P, Xing, H, Liao, L, Feng, Y, Leng, X, Wang, J, 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:9. doi: 10.1186/s12981-020-00264-5
35. Rhee, SY, Gonzales, MJ, Kantor, R, Betts, BJ, Ravela, J, and Shafer, RW. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. (2003) 31:298–303. doi: 10.1093/nar/gkg100
36. Hong, H, Tang, C, Liu, Y, Jiang, H, Fang, T, and Xu, G. HIV-1 drug resistance and genetic transmission network among newly diagnosed people living with HIV/AIDS in Ningbo, China between 2018 and 2021. Virol J. (2023) 20:233. doi: 10.1186/s12985-023-02193-x
37. Pang, X, Liang, S, Tang, K, Huang, J, He, Q, Fang, N, et al. Disparity of HIV-1 Pretreatment Drug Resistance in Men Who Have Sex With Men and the Heterosexual Population in Guangxi, China. Open Forum. Infect Dis Ther. (2023) 10:016. doi: 10.1093/ofid/ofad016
38. Zhang, M, Ma, Y, Wang, Z, Wang, G, Wang, Q, Li, X, et al. Prevalence and transmission of pretreatment drug resistance in people living with HIV-1 in Shanghai China, 2017-2021. Virulence. (2024) 15:2373105. doi: 10.1080/21505594.2024.2373105
39. Chen, H, Hao, J, Hu, J, Song, C, Zhou, Y, Li, M, et al. Pretreatment HIV Drug Resistance and the Molecular Transmission Network Among HIV-Positive Individuals in China in 2022: Multicenter Observational Study. JMIR Public Health Surveill. (2023) 9:e50894. doi: 10.2196/50894
40. Yu, E, Xu, YJ, Li, M, Yang, Y, Liang, HY, Zhong, SM, et al. HIV-1 Subtype Diversity and Factors Affecting Drug Resistance among Patients with Virologic Failure in Antiretroviral Therapy in Hainan Province, China, 2014-2020. Biomed Environ Sci. (2023) 36:800–13. doi: 10.3967/bes2023.077
41. Santos-Pereira, A, Triunfante, V, Araujo, PMM, Martins, J, Soares, H, Poveda, E, et al. Nationwide Study of Drug Resistance Mutations in HIV-1 Infected Individuals under Antiretroviral Therapy in Brazil. Int J Mol Sci. (2021) 22:304. doi: 10.3390/ijms22105304
42. Shchemelev, AN, Ostankova, YV, Zueva, EB, Semenov, AV, and Totolian, AA. Detection of Patient HIV-1 Drug Resistance Mutations in Russia's Northwestern Federal District in Patients with Treatment Failure. Diagnostics (Basel). (2022) 12:821. doi: 10.3390/diagnostics12081821
43. Zhang, J, Sun, B, Sheng, Z, Ding, X, Fan, Q, Huang, G, et al. Full-Spectrum Surveillance of Pre-Treatment HIV Drug Resistance in Southeastern China. Pharmaceuticals (Basel). (2024) 17:900. doi: 10.3390/ph17070900
44. Geremia, N, Basso, M, De Vito, A, Scaggiante, R, Giobbia, M, Battagin, G, et al. Patterns of Transmitted Drug Resistance Mutations and HIV-1 Subtype Dynamics in ART-Naive Individuals in Veneto, Italy, from 2017 to 2024. Viruses. (2024) 16:393. doi: 10.3390/v16091393
45. Fiebig, U, Altmann, B, Hauser, A, Koppe, U, Hanke, K, Gunsenheimer-Bartmeyer, B, et al. Transmitted drug resistance and subtype patterns of viruses from reported new HIV diagnoses in Germany, 2017-2020. BMC Infect Dis. (2023) 23:673. doi: 10.1186/s12879-023-08649-3
46. WHO. World Health Organization global strategy for the surveillance and monitoring of HIV drug resistance. (2012). Available from: https://iris.who.int/bitstream/handle/10665/77349/9789241504768_eng.pdf?sequence=1 (Accessed February 21, 2025).
47. Mazzuti, L, Melengu, T, Falasca, F, Calabretto, M, Cella, E, Ciccozzi, M, et al. Transmitted drug resistance mutations and trends of HIV-1 subtypes in treatment-naive patients: A single-centre experience. J Glob Antimicrob Resist. (2020) 20:298–303. doi: 10.1016/j.jgar.2019.08.024
48. Kemp, SA, Kamelian, K, Cuadros, DF, Consortium, P, Vukuzazi, T, Cheng, MTK, et al. HIV transmission dynamics and population-wide drug resistance in rural South Africa. Nat Commun. (2024) 15:3644. doi: 10.1038/s41467-024-47254-z
49. Paredes, R, Sagar, M, Marconi, VC, Hoh, R, Martin, JN, Parkin, NT, et al. In vivo fitness cost of the M184V mutation in multidrug-resistant human immunodeficiency virus type 1 in the absence of lamivudine. J Virol. (2009) 83:2038–43. doi: 10.1128/JVI.02154-08
50. Rhee, SY, Varghese, V, Holmes, SP, Van Zyl, GU, Steegen, K, Boyd, MA, et al. Mutational Correlates of Virological Failure in Individuals Receiving a WHO-Recommended Tenofovir-Containing First-Line Regimen: An International Collaboration. EBioMedicine. (2017) 18:225–35. doi: 10.1016/j.ebiom.2017.03.024
51. El Bouzidi, K, Datir, RP, Kwaghe, V, Roy, S, Frampton, D, Breuer, J, et al. Deep sequencing of HIV-1 reveals extensive subtype variation and drug resistance after failure of first-line antiretroviral regimens in Nigeria. J Antimicrob Chemother. (2022) 77:474–82. doi: 10.1093/jac/dkab385
Keywords: HIV-1, treatment failure, drug resistance, subtype, antiretroviral therapy, mutation
Citation: Pei R, Zhang Y, Jike C, Yu G, Su L, Wang J, Xiao L, Wang Y, Shen M, Liao J, Zheng Y and Hemelaar J (2025) Prevalence and transmission of HIV-1 drug resistance mutations among patients with treatment failure and newly diagnosed people in Liangshan Prefecture, China, in 2021–2023. Front. Public Health. 13:1550121. doi: 10.3389/fpubh.2025.1550121
Edited by:
Yoon-Seok Chung, Korea Disease Control and Prevention Agencyorea, Republic of KoreaReviewed by:
Bingyu Liang, Guangxi Medical University, ChinaNam Yee Kim, Incheon Institute of Health and Environment, Republic of Korea
Copyright © 2025 Pei, Zhang, Jike, Yu, Su, Wang, Xiao, Wang, Shen, Liao, Zheng and Hemelaar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joris Hemelaar, am9yaXMuaGVtZWxhYXJAbmRwaC5veC5hYy51aw==