HIV-1 Drug Resistance Profiles of Low-Level Viremia Patients and Factors Associated With the Treatment Effect of ART-Treated Patients: A Cross-Sectional Study in Jiangsu, China

Objectives: Evaluating the drug resistance (DR) profiles of LLV patients and the influencing factors of treatment effects in Jiangsu Province.

Method: The Pol gene (Reverse transcriptase and protease) was amplified and sequenced to identify the genotypes and DR profiles among LLV patients in 2021. Questionnaire survey was conducted among HIV/AIDS patients to investigate the potential influence factors of treatment effects.

Results: 242 Pol genes were amplified from 345 specimens, and ten genotypes were detected. The DR rate was 40.5%, with 66, 86, and 14 being resistant to NRTIs, NNRTIs, and PIs, respectively. Patients treated with the 2NRTIs+PIs regimen were detected with more DR; and drug resistance was less detected in married or cohabiting patients than unmarried patients. Non-smokers were less likely to develop LLV at follow-up than smokers; patients with stage II clinical stage at diagnosis and using 2NRTIs+PIs regimen were also more likely to develop LLV at follow-up.

Conclusion: Drug resistance profiles in LLV patients are severe and differ in treatment regimens and marital statuses. Meanwhile, smoking history, clinical stage, and treatment regimen may influence the therapeutic effect. It is necessary to include LLV people in the free drug resistance testing program.

Introduction

In 2020, about 1.5 million HIV cases were reported worldwide (1), and 62 167 AIDS cases were reported in China (2), There were 3 755 newly reported HIV/AIDS patients in Jiangsu Province, an increase of nearly 9.7% compared with last year. By the end of October, there were 35 284 living HIV/AIDS patients in Jiangsu province, among which the proportion of ART recipients had reached 99.1% (3).

ART can limit viral replication to <50 copies/ml and effectively prolong the life of HIV/AIDS patients. However, in ~10–30% of patients, standard ART does not fully control viral replication, and these patients show low levels of viral replication in their plasma, which is known as low-level viremia (LLV). According to WHO guidelines, LLV is defined as viral load (VL) between 50 and 1,000 copies/ml (4). However, LLV patients are not covered by free drug resistance monitoring in Jiangsu Province.

In previous studies that used 1,000 copies/ml as the Virological Failure (VF) criterion, Acquired Drug Resistance (ADR) was reported from 10.96 to 91.67% (5, 6), and as of 2017, the total ADR in China was 44.7% (7). Many studies pointed out that the mutation rate of drug resistance (DR) is 30–46% in LLV patients, and the highest can even reach 86% (8). However, there are few studies on drug resistance in LLV patients in China, and the drug resistance profiles of this population are still unclear.

The influencing factors of LLV and its impact on clinical prognosis remain controversial (9–15). In addition, there are few relevant research results in China, Only one long-term cohort study was conducted in Shenyang in 2020, and the results showed that baseline VL, CD4 cell counts, ethnicity, treatment duration, treatment regime, and virus subtype were all related to the occurrence of LLV, while different LLV groups and maintenance time were related to the occurrence of treatment failure (16). Nevertheless, there is still a lack of relevant research in Jiangsu Province.

In order to understand the drug resistance profiles of LLV patients in Jiangsu Province and explore possible factors related to the occurrence of LLV, this study tested the drug resistance of LLV patients who were followed up in Jiangsu province in 2021, and a questionnaire survey was also conducted among LLV patients and matched Viral Suppression (VS) patients, hoping to provide reference information for subsequent treatment effect evaluation, treatment regimen adjustment, life quality improvement, and DR strains' transmission reduction.

Materials and Methods

Study Design and Study Object

Molecular Epidemiological Investigation

A cross-sectional study design was used to study the DR profiles in LLV patients who were followed up in Jiangsu Province in 2021. Inclusion criteria included: (1) Management (treatment and follow-up) in Jiangsu Province; (2) The subjects were older than 18 years at follow-up; (3) ART time at follow-up was more than 6 months; (4) VL was in the range of 50–1,000 copies/ml during follow-up; (5) Informed consent was given to collect and use biological samples during follow-up.

Questionnaire Investigation

A matched case-control study design was adopted, LLV patients were selected as the case group, and matched sampling was conducted among VS patients in a 1:1 ratio according to the characteristics of included LLV patients (region, gender, and age) to carry out epidemiological case investigation. Inclusion criteria were: (1) age 18 and above; (2) HIV positive patients confirmed in Jiangsu Province; (3) The follow-up period was from January 2021 to December 2021; (4) the treatment duration was more than half a year (≥180 days) at the time of follow-up; (5) Patients in the case group were 50–1,000 copies/ml (including 50 and 1,000), while patients in the control group were VS patients (VL<50 copies/ml); (6) Voluntarily participate in the program and give informed consent.

Sample Size Calculation of Questionnaire Survey

The results of the follow-up treatment database showed that about 20% of the patients used the treatment regimen containing Protease Inhibitors (PIs), with the treatment regimen as the primary exposure factor to be observed, setting OR = 2.7, Phi = 0.2, α = 0.05, β = 0.10. The sample size calculated by PASS (Version: 2021) software is N = 272. Assuming that the non-response rate is 10% and the effective rate of the questionnaire is 90%, the sample size is N = 272/(1–0.1)/0.9 = 336 cases (168 LLV and 168 controls, respectively), and the volume of questionnaires should be appropriately increased according to the actual situation.

Sample Collection and Experiment Operation

The samples of this study were collected by the drug resistance monitoring program for HIV/AIDS patients in Jiangsu Province. The collected specimens were extracted by the Lizhu automatic nucleic acid extractor (4-channel). Target genes in the Pol region (Reverse transcriptase and protease) were amplified by the in-house method. Reaction conditions, amplification system, and primers are shown in Table 1. Gel electrophoresis was performed on the amplified products, and a gel imager (Tanon-3500B) was used to determine the electrophoresis results. PCR products were sent to Sangon Biotech for sequencing.

TABLE 1

Table 1. Amplification system of Pol gene.

Baseline Data Collection and Questionnaire Survey

The Ethics Review Committee of Jiangsu Provincial Center for Disease Control and Prevention has approved this study (Approval No: JSJK2021-B017-01). All the included patients' baseline data, follow-up data, and laboratory test data during treatment were obtained from the baseline information database and treatment information database of HIV patients in Jiangsu Province, and a supplementary information survey was conducted by online questionnaire and offline interview.

Data Management and Analysis

Sequence Processing and Subtype Determination

ChromasPro (Version: 1.6) was used for sequence spliced, website (https://comet.Lih.lu/) and HIV database (https://www.hiv.lanl.gov/content/sequence/BASIC_BLAST/basic_blast.html) were used for subtype determination. FastTree (Version: 1.4.3) was used to construct the phylogenetic tree to verify the online subtype determination results.

Drug Resistance Analysis

Pol region sequences were uploaded to the HIV Drug Resistance Database (https://hivdb.Stanford.edu/) for drug resistance analysis. The drug resistance degree of each drug was divided into five grades: sensitive (1 score), potential resistance (2 score), low-level resistance (3 score), moderate-level resistance (4 score), and high-level resistance (5 score). In this study, samples with sensitive and potential resistance were regarded as not resistant, while those with low or above level resistance grade (resistance score ≥3) were regarded as resistant. Drug resistance rate = number of drug resistance ÷ (number of drug resistance + number of non-drug resistance) * 100.

Statistical Analysis

Excel was used to manage the collected information, and SPSS (Version: 25) was used for statistical analysis. Quantitative data were described by mean ± standard deviation (X ± S), and comparison between groups was performed by two-sample independent T-test or analysis of variance. Qualitative data were described by ratio or constituent ratio. χ2 test, Fisher's exact probability method, or Kruskal-Wallis method were used for statistical analysis. Logistic regression model was used for univariate and multivariate analysis, and the α = 0.05.

Results of Molecular Epidemiological Investigation

Characteristics of LLV Patients

In this study, a total of 345 samples were collected throughout 2021. After nucleic acid extraction, PCR amplification and sequencing were performed on the plasma of 345 samples, Pol genes of 255 patients were obtained. Thirteen cases of sequences with poor quality or unmatched information were removed, and 242 cases of sequences meeting the research requirements were finally included. The success rate of amplification was 70.14% (VL 50–199: 59.3%; VL 200–499: 79.4%; VL: 500–1,000: 91.1%). The characteristics of the successfully amplified subjects are shown in Table 2.

TABLE 2

Table 2. Characteristics of LLV patients.

Chi-square test results showed that there were differences in drug resistance detection rates among patients with different marital statuses, viral load, and treatment regimens (Table 2). The detection rate of drug resistance was lower (33.1%) in married or cohabiting patients and lower in patients with VL 50–199 copies/ml (33.6%), while more patients with drug resistance were found in the VL 500–1,000 copies/ml group (52.9%). Resistance detection rates were also lower in patients treated with 2NRTIs+NNRTIs and 2NRTIs+INSTIs (35.2 and 35.7%, respectively) than in patients treated with 2NRTIs+PIs (63.6%).

Distribution of HIV-1 Genotypes in LLV Patients

Ten HIV-1 genotypes were detected based on the Pol gene fragments from 242 cases (Figure 1). A total of 10 HIV-1 genotypes were detected, the proportion from high to low were as follows: Circulating Recombinant Forms (CRF) 01_AE (94, 38.8%), CRF07_BC (85, 35.1%), B (21, 8.7%), CRF67_01B (14, 5.8%), CRF08_BC (13, 5.4%), CRF55_01B (5, 2.1%), CRF68_01B (4, 1.7%), 01BC (4, 1.7%), 0107 (1, 0.4%) and CRF65_cpx (1, 0.4%). The genotypes detected were mainly CRF with a small amount of Unique Recombinant Forms (URFs) (01BC and 0107). Figure 2 shows the phylogenetic tree constructed with Pol genes.

FIGURE 1

Figure 1. HIV-1 genotypes distribution among 242 LLV patients in Jiangsu Province, China.

FIGURE 2

Figure 2. Phylogenetic tree. Phylogenetic tree was constructed from 242 pol gene sequences of LLV patients in Jiangsu province and reference sequences by FastTree software, with parameters as GTR + G + I replacement model and Bootstrap value ≥0.90.

Drug Resistance Characteristics in LLV Patients

Frequency Analysis of Different Types of Drug Resistance

There are 98 drug resistance sequences among 242 pol genes (low or above level resistance to any class of drugs), and the drug resistance rate was 40.5% (98/242). Among them, 66 (67.3%), 86 (87.8%), and 14 (14.3%) were resistant to NRTIs, NNRTIs, and PIs, respectively. A total of 63 samples showed low or above level drug resistance to more than one class of drugs, including 55 cases of simultaneous resistance to NRTIs and NNRTIs, 1 case of simultaneous resistance to NRTIs and PIs, 2 cases of simultaneous resistance to NNRTIs and PIs, and 5 cases of drug resistance to all three types of drugs (Figure 3).

FIGURE 3

Figure 3. Frequency of different drug resistance types in 98 Drug resistance LLV patients in Jiangsu Province, China.

Influencing Factors of Drug Resistance in LLV Patients

Using drug resistance profiles as the dependent variable and multiple factors as independent variables, univariate regression analysis showed statistically significant differences among marital status (Married or cohabiting vs. Unmarried: OR = 0.537; 95%CI: 0.299–0.963; P = 0.037) and treatment regimens (2NRTIs+PIs vs. 2NRTIs+NNRITs: OR = 3.219; 95%CI: 1.606–6.450; P = 0.001), multivariate regression analysis was further performed with marital status and treatment regimen as independent variables (Table 3), and the results showed that compared with unmarried patients, drug resistance was less detected in married or cohabiting patients (OR = 0.509, 95%CI: 0.279–0.928). Compared with the 2NRTIs+NNRTIs regimen, more drug resistance was detected in patients treated with the 2NRTIs+PIs regimen (OR = 3.240; 95%CI: 1.595–6.584).

TABLE 3

Table 3. Factors associated with drug resistance among LLV patients.

Characteristics of Drug Resistance Mutations

According to the drug resistance analysis results, 66 drug resistance mutation (DRM) sites were detected in 98 patients with drug resistance. Fifty-four DRM sites were detected in the reverse transcription (RT) region, including 25 NRTIs related sites with 125 mutations. M184V mutation frequency was the most frequent (40.8%, 51/125), followed by K65R (12.0%, 15/125). There were 29 DRM sites associated with NNRTIs, with 184 mutations occurring in total, mainly K103N (24.5%, 45/184) and Y181C (8.2%, 15/184), followed by K101E (6.0%, 11/184) and V106M (6.0%, 11/184).

A total of 12 DRM sites were detected in the protease inhibitors (PIs) region, among which 5 were the major mutation sites, resulting in 11 mutations, with M46I (45.5%, 5/11) and M46L (27.3%, 3/11) having the highest frequency. The frequency and percentage of other DRM are shown in Figure 4.

FIGURE 4

Figure 4. Distribution and frequency of drug resistance mutation sites to NNRTIs, NRTIs and PIs (Major and Accessory).

The Influence of DRM Sites on Drug Sensitivity

The DRM sites and corresponding drug sensitivity of 98 patients with drug resistance are shown in the Supplementary Materials. Sixty-six of the 98 patients had NRTIs-related DRM sites, and all caused low or above level drug resistance, 37 of them carried a single DRM site, accounting for 56.1%, and the rest all carried two or more NRTIs DRM sites, which mainly c aused medium or high-level resistance to Abacavir (ABC), Emtricitabine (FTC), and Lamivudine (3TC).

Ninety patients developed NNRTIs-related DRM sites, 86 patients developed low or above level drug resistance to NNRTIs drugs, and 56 patients carried two or more DRM sites (62.2%). Except for four patients with only V179D/E mutation sites who did not develop DR to NNRTIs drugs, the rest had medium or high-level resistance to Efavirenz (EFV) and Nevirapine (NVP). Meanwhile, more than half of the patients had medium or high-level resistance to Doravirine (DOR) and Rilpivirine (RPV) (64.0%, 55/86).

Seventeen patients developed DRM sites associated with PIs drugs, and 14 patients had low or above level resistance, of which 71.4% carried only a single mutation site (10/14), mainly producing low or medium-level resistance to Nelfinavir (NFV).

Drug Resistance Profiles in Patients With Different Treatment Regimens

The treatment regimens were divided into four categories (Table 4). Patients using the 2NRTIs+NNRTIs regimen accounted for the highest proportion (57.1%, 56/98), with 39 patients resistant to NRTIs, ABC, FTC, and 3TC having the highest drug resistance frequency; 49 patients resistant to NNRTIs, and all were resistant to EFV and NVP. Patients using the 2NRTIs+PIs regimen accounted for 28.6% (28/98); among patients using this regimen, there were 20 and 25 patients resistant to NRTIs and NNRTIs, respectively, but only three patients were resistant to PIs.

TABLE 4

Table 4. Drug resistance profiles in patients with different treatment regimens.

Questionnaire Survey Results

Characteristics of Research Object

Four hundred seventy-two questionnaires were sent out, among which nine people refused the survey. Four hundred sixty-three questionnaires were returned with a recovery rate of 98.1%. The 17 questionnaires that did not fill in survey numbers and codes as required and whose information could not match were removed, and 446 questionnaires were available. Finally, 380 questionnaires (190 LLV patients and 190 matched VS patients) were analyzed.

Most VS and LLV patients were married or cohabiting (56.3 and 52.1%), and most VS and LLV patients did not smoke or drink alcohol. Most patients were in clinical stage I at the time of diagnosis (VS: 57.9%; LLV: 41.6%), and more than half of VS and LLV patients received ART for more than 3 years (59.5 and 56.8%). The proportion of the CD4 cell count group was from high to low: 200–499 cells/μL (VS: 48.4%; LLV: 51.6%), 500 cells/μL or more (VS: 38.9%; LLV: 32.6%) and 200 cells/μL or less (VS: 12.6% and 15.8%); In VS patients, 92.6% did not miss medication in the last month, and 85.3% in LLV patients; The majority of patients received 2NRTIs+NNRTIs at follow-up (VS: 83.7%; LLV: 68.9%), followed by 2NRTIs+PIs (VS: 10.0%; LLV: 20.5%) (Table 5).

TABLE 5

Table 5. Characteristics of surveyed HIV/AIDS patients.

At the same time, considering differences in usage and dosage among treatment regimens, which might affect patients' compliance, the treatment regimens were stratified, drug omissions in a recent month were compared, and no statistical difference was found (χ² = 5.133, P = 0.162).

Influencing Factors of Treatment Effect in HIV/AIDS Patients

With treatment effect (LLV or VS at follow-up) as the dependent variable, univariate logistic regression was performed with other independent variables (Table 6). There were statistically significant differences among self-assessment compliance (P = 0.033), treatment regimen (2NRTIs+PIs vs. 2NRTIs+NNRITs: P = 0.003), smoking history (No smoking vs. Smoking: P = 0.021), clinical stage at diagnosis (Stage II vs. Stage I: P = 0.032; Stage III vs. Stage I: P = 0.025; Stage IV vs. Stage I: P = 0.015), and medication omission in the last month (P = 0.024).

TABLE 6

Table 6. Factors associated with treatment effect among HIV/AIDS patients.

Independent variables with significant univariate analysis results were included in multivariate logistic regression analysis, and no correlation was found between LLV and self-assessment compliance or medication omission (P > 0.05). Smoking history (No smoking vs. Smoking: P = 0.026), clinical stage at diagnosis (Stage II vs. Stage I: P = 0.046), and treatment regimen (2NRTIs+PIs vs. 2NRTIs+NNRITs: P = 0.016) showed statistically significant differences in treatment effect among different groups. Non-smokers were less likely to develop LLV at follow-up than smokers (OR = 0.579, 95%CI: 0.358–0.936); patients with stage II clinical stage at diagnosis (OR = 1.700, 95%CI: 1.009–2.864) and using 2NRTIs+PIs regimen (OR = 2.171, 95%CI: 1.152–4.093) were also more likely to develop LLV at follow-up (Table 6).

Discussion

LLV patients are not included in the free drug resistance surveillance program in Jiangsu Province now. The study results showed that more than 40% of LLV patients had low or above level drug resistance to the three main ART drugs (NRTIs, NNRTIs, and PIs). Nearly 90% of the DR-patients were resistant to NNRTIs, and more than half were simultaneously resistant to NRTIs and NNRTIs. Suggesting that DR profiles in LLV patients are severe, so adjust the monitoring scope of DR and provide free DR testing for LLV patients, which can offer information for timely treatment regimen adjustment, treatment effect improving, disease progression deceleration, and patients' life quality improving, are necessary.

The information on LLV patients receiving ART shows that the elderly over 55 years old account for the highest proportion, followed by the young people between 25 and 34 years old. This may be due to the education lacking, insufficient disease understanding, poor compliance during ART, delayed diagnosis (17), other underlying diseases, and worse immune reconstitution effect is more prominent in the elderly (18). In recent years, the proportion of young people among all HIV/AIDS patients has increased, and the total number has also increased rapidly (19). In addition, most young patients have active thoughts, open sexual concepts, and unprotected homosexual behavior (20), which is at high risk of cross-infection and may also affect the treatment effect. Therefore, providing targeted consulting, testing services, and strengthening knowledge propaganda, warning education according to the characteristics of different populations can effectively improve the treatment effect and help promote prevention and control work.

Although there were no statistically significant differences in DR rates among LLV patients of different genders, DR rates were slightly higher among women. Previous studies on DR of LLV patients rarely provided detailed profiles of different genders. But previous studies on VF patients in China indicated that the DR rate of female VF patients was higher than that of males (21, 22), which was consistent with the results of this study. This may be due to higher PDR in the included female patients or because fewer female patients were collected and more resistant mutations occurred in these women. However, some studies also indicate that male VF patients have a higher drug resistance rate (23), so whether there are differences in DR rates between LLV patients of different genders needs further study.

HIV/AIDS patients in Jiangsu province mostly use the 2NRTI+NNRTIs regimen according to Chinese guidelines for diagnosis and treatment of HIV/AIDS (24). This study showed that NNRTIs resistance was the most in LLV patients, and the simultaneous resistance to NRTIs and NNRTIs was the main one among multiple-DR, same as previous results (25). More than half of DR patients developed simultaneous DR, but resistance to PIs was relatively rare. It is worth noting that this study was a cross-sectional study within 1 year; the treatment regimen of patients was identified as the one adopted at the follow-up time point, and the DR of patients may not have been generated during the use of the current regimen. However, this study did not obtain data on previous regimen changes and DR profiles of patients, so the statistical significance of the association between treatment regimens and DR in this study may not indicate the authenticity of the association. In order to verify this possibility based on existing data, we compared drug resistance profiles in patients with different regimens and found that among patients with 2NRTI+PIs regimens, only three patients developed low or medium-level DR to PIs but 25 patients developed medium or high-level resistance to NNRTIs. Although most patients are also resistant to NRTIs, it is not clear whether the resistance to NRTIs occurred before or after using the current regimen, so the difference in resistance risk between regimens still needs further study. It suggests that cross-drug resistance is severe in LLV patients; timely regimen adjustments should be made according to DR profiles, and using the PIs as the replaced drug is still a better choice, but the DR of PIs should also be monitored.

Multiple regression analysis also showed differences in the DR detection likelihood among patients with different marital statuses. DR detection likelihood in the married (Heterosexual marriage) or cohabitation (Both heterosexual and homosexual) patients was lower than that in unmarried patients; this may be due to the stable emotional state and their partner's encouragement, which promote patients' medication compliance and test timeliness, maintaining the better treatment effects, and reducing DR emergence.

Although no association between subtypes and DR was found in this study, the results showed that subtypes of drug-resistant LLV patients were mainly 01_AE and 07_BC, which had the same epidemic trend with subtypes in Jiangsu Province (26), and there were more drug-resistant mutation sites in 01_AE patients. There are few studies on the association between subtypes and DR in China; previous studies in Jiangsu province indicated that the proportion of 01_AE subtypes in the VF patients increased from 40.4% in 2009 to 56.8% in 2015 (25, 27). Meanwhile, higher X4 cell tropism (28), faster disease progression and CD4 cell counts decline (29) of 01_AE strain are also contribute to the higher DR likelihood (30, 31), which also indicates that subtype, especially 01_AE, may be related to the occurrence of drug resistance. It is suggested that DR monitoring should be strengthened in patients with these subtypes, but the differences in DR profiles among subtypes need further research.

The questionnaire survey found that patients with different smoking histories, clinical stages, and treatment regimens had different likelihoods of LLV at follow-up. Patients with smoking history were more likely to develop LLV than non-smokers at follow-up. Smoking is prevalent among HIV/AIDS patients (32) and is also considered the cause of reduced ART compliance and increased chance of complications (33, 34). Previous studies have pointed out that smoking has pro-inflammatory and immunosuppressive effects, which can change the activation of T cells (35) and affect the immune system response by changing the reactivity and differentiation of T cells (36, 37), thus damaging the immune system response to pathogens (including viruses). In addition, studies have also shown that cigarettes can induce oxidative stress and increase HIV replication (38), resulting in a higher VL in smokers than non-smokers (39). Smoking can promote HIV replication in cells has also been found in vitro studies (40). Although more research is needed to explore the specific cellular pathways by which cigarette smoke affects HIV replication, promoting smoking cessation and improving unhealthy lifestyles among HIV/AIDS patients may also be a supplemental measure to improving treatment effectiveness.

The WHO divides patients into four clinical stages (41). This study found that patients with different clinical stages at diagnosis have different likelihoods of LLV during subsequent treatment. Previous studies showed that patients' nutritional status recovery time in clinical stages I and II was shorter than in stages III and IV (4 vs. 8 months) (42, 43). When stages III and IV progress, the number of complications increases, the risk of infection and tumor increases (44), and the rate of weight gain is slower (45), leading to a more extended recovery period. The clinical-stage is the overall reaction of body status and influences therapeutic efficacy in various ways. Suggests that treatment regimens should combine with the baseline health condition, inhibit viral replication and cure the concomitant diseases simultaneously to improve the body's health and retard the disease process effectively.

Previous studies have pointed out that patients' compliance during treatment is related to the occurrence of LLV and is one of the main influencing factors. However, our study did not show a statistical association between compliance and treatment effect. Due to the small number of participants, the Power of negative conclusion was only 3%, so compliance cannot be considered irrelevant to the occurrence of LLV. In addition, the compliance in this study was surveyed by a self-filled questionnaire, which may be different from the actual situation, and lead to no statistical significance in the study.

Patients using PIs-based regimens were more likely to have LLV at follow-up. Cohort studies have indicated that PIs-based regimens will increase the risk of LLV in patients (46–48), and compared with PIs, patients treated with NNRTIs have a better performance in inhibiting plasma virus replication (49–51) and a lower incidence of LLV (52–55). While studies have suggested that differences between hospitals where patients recruit may influence the link between regimens and LLV, that study did not find significant differences in treatment regimens among patients from different sources (52). Leierer also noted that their results did not exclude the effect of compliance on this association (46), but Konstantopoulos still found that PIs-based regimens had nearly three times the risk of LLV compared with other regiments (56), which was consistent with the results of this study, indicating that patients using PIs-based regimens had more active virus replication and a higher possibility of virologic rebound. Viral replication is also related to different drugs' pharmacokinetics and tissue permeability. Previous studies have shown that insufficient drug concentration during ART will increase the risk of viral replication activation (57). The better pharmacokinetic characteristics of NNRTIs ensure that patients can maintain a high-level drug concentration even if poor compliance occurs (58), while the faster metabolic rate of PIs leads to a higher risk of LLV in drug omissions (59). However, further studies are needed to explore the mechanism of PIs leading to a higher risk of viral replication activation. The current drug development goals are better long-term sensitivity, longer half-life period, and fewer pharmacological interactions.

This questionnaire survey was carried out online and offline. The patients were surveyed after the management physician contacted and obtained informed consent. Hence, the cooperation degree of patients was high, the non-response rate was low, and the quality of the questionnaire was good. However, this study also has some limitations: Firstly, Among the patients included in this study, the number of patients in the VL 50–199 group was significantly higher than that in other groups, which may have biased the analysis to some extent. Secondly, due to the condition limitations, the compliance data obtained in this study were all from self-filled questionnaires (follow-up survey and supplementary survey), and lack of pharmacy re-fill records or Medication Event Monitoring System (MEMS) data, which might overestimate patients' compliance. Thirdly, the follow-up and treatment information were acquired in 2021, which was not combined with the multi-year VL variation trend and patients' PDR profiles, so the influence of factors on the treatment effect with the change of treatment duration in the whole treatment process could not be more accurately assessed. Further research and analysis are needed.

Conclusion

This study found that nearly 40% of LLV patients in Jiangsu province had low or above level drug resistance to the ART drugs, and drug resistance profiles differed in treatment regimens and marital statuses. Meanwhile, smoking history, clinical stage, and treatment regimen may influence the therapeutic effect. Based on the existing strategy, more detailed and targeted follow-up management standards should be formulated according to different patients' characters and virus strains' molecular characteristics. In the meantime, it is necessary to include LLV people in the free drug resistance testing program and strengthen the drug resistance monitoring of ART patients to improve the treatment effect.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Review Committee of Jiangsu Provincial Center for Disease Control and Prevention (Approval No: JSJK2021-B017-01). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

GF and BW guarantee the paper's project selection, overall conduct, and correction. DY is responsible for lab operation, questionnaire collection, and article writing. YZ and LS provided help for obtaining the specimen, questionnaire collection, constructing the framework of the article, and proofreading. YL assisted with data analysis and charting. At the same time, JC and JL assisted in writing, translating, and proofreading the paper. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by Molecular Network Analysis and Social Network Exploration of HIV-1 Infection Transmission among Young Students in Jiangsu Province (0701-184160070478).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.

Acknowledgments

We want to express our sincere appreciation to the Jiangsu Provincial CDC for providing the experiment condition and every doctor who assisted in sample collection and questionnaire survey.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.944990/full#supplementary-material

References

1. UNAIDS. UNAIDS data 2021 (2021). Available online at: https://www.unaids.org/sites/default/files/media_asset/JC3032_AIDS_Data_book_2021_En.pdf (accessed March 9, 2022).

Google Scholar

2. CDC C. General Situation of Notifiable Infectious Diseases in China in 2020. (2021). Available online at: http://www.nhc.gov.cn/jkj/s3578/202103/f1a448b7df7d4760976fea6d55834966.shtml (accessed March 8, 2022).

Google Scholar

3. Commission JPH. Progress of HIV/AIDS Prevention Treatment in 2021. (2021). Available online at: http://wjw.jiangsu.gov.cn/art/2021/11/30/art_7312_10137196.html (accessed March 9, 2022).

Google Scholar

4. Bennett DE, Bertagnolio S, Sutherland D, Gilks CF. The World Health Organization's global strategy for prevention and assessment of HIV drug resistance. Antivir Ther. (2008) 13(Suppl. 2):1–13. doi: 10.1177/135965350801302S03

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lu Y, Wang J, Chen Y, Li J. A retrospective study on drug resistance in HIV-1 infected patients failing firstline antiretroviral therapy in Chongqing. Chin J Infect Chemother. (2021) 21:427–31. doi: 10.16718/j.1009-7708.2021.04.009

CrossRef Full Text | Google Scholar

6. Qin C, Tan X, Liang X, Zhu G, Xu Z, Huang Y, et al. Analysis of genotype drug resistance of HIV-1 in commercial heterosexually transmitted population in Guigang city, Guangxi. Appl Prev Med. (2020) 26:374–8.

Google Scholar

7. Zuo L, Liu K, Liu H, Hu Y, Zhang Z, Qin J, et al. Trend of HIV-1 drug resistance in China: A systematic review and meta-analysis of data accumulated over 17 years (2001-2017). EClinicalMedicine. (2021) 33:100696. doi: 10.1016/j.eclinm.2020.100696

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Taramasso L, Magnasco L, Bruzzone B, Caligiuri P, Bozzi G, Mora S, et al. How relevant is the HIV low level viremia and how is its management changing in the era of modern ART? A large cohort analysis. J Clin Virol. (2020) 123:104255. doi: 10.1016/j.jcv.2019.104255

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Vancoillie L, Mortier V, Demecheleer E, Schauvliege M, Vandekerckhove L, Vogelaers D, et al. Drug resistance is rarely the cause or consequence of long-term persistent low-level viraemia in HIV-1-infected patients on ART. Antivir Ther. (2015) 20:789–94. doi: 10.3851/IMP2966

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Charpentier C, Landman R, Laouenan C, Joly V, Hamet G, Damond F, et al. Persistent low-level HIV-1 RNA between 20 and 50 copies/mL in antiretroviral-treated patients: associated factors and virological outcome. J Antimicrob Chemother. (2012) 67:2231–5. doi: 10.1093/jac/dks191

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Pernas B, Grandal M, Pertega S, Canizares A, Castro-Iglesias A, Mena A, et al. Any impact of blips and low-level viraemia episodes among HIV-infected patients with sustained virological suppression on ART? J Antimicrob Chemother. (2016) 71:1051–5. doi: 10.1093/jac/dkv433

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Antiretroviral Therapy Cohort C, Vandenhende MA, Ingle S, May M, Chene G, Zangerle R, et al. Impact of low-level viremia on clinical and virological outcomes in treated HIV-1-infected patients. AIDS. (2015) 29:373–83. doi: 10.1097/QAD.0000000000000544

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Paintsil E, Martin R, Goldenthal A, Bhandari S, Andiman W, Ghebremichael M. Frequent episodes of detectable viremia in HIV treatment-experienced children is associated with a decline in CD4+ T-cells over time. J AIDS Clin Res. (2016) 7:565. doi: 10.4172/2155-6113.1000565

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Leierer G, Grabmeier-Pfistershammer K, Steuer A, Sarcletti M, Geit M, Haas B, et al. A single quantifiable viral load is predictive of virological failure in human immunodeficiency virus (HIV)-infected patients on combination antiretroviral therapy: the Austrian HIV cohort study. Open Forum Infect Dis. (2016) 3:ofw089. doi: 10.1093/ofid/ofw089

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Doyle T, Geretti AM. Low-level viraemia on HAART: significance and management. Curr Opin Infect Dis. (2012) 25:17–25. doi: 10.1097/QCO.0b013e32834ef5d9

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang T, Ding H, An M, Wang X, Tian W, Zhao B, et al. Factors associated with high-risk low-level viremia leading to virologic failure: 16-year retrospective study of a Chinese antiretroviral therapy cohort. BMC Infect Dis. (2020) 20:147. doi: 10.1186/s12879-020-4837-y

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Mugavero MJ, Castellano C, Edelman D, Hicks C. Late diagnosis of HIV infection: the role of age and sex. Am J Med. (2007) 120:370–3. doi: 10.1016/j.amjmed.2006.05.050

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ding H, Liu J, Xu J, He Y, Guo X, Kang J, et al. The clinical outcomes and prognosis of older human immunodeficiency virus-infected patients under antiretroviral therapy. Chin J Infect Dis. (2018) 36:417–21. doi: 10.3760/cma.j.issn.1000-6680.2018.07.007

CrossRef Full Text | Google Scholar

19. Cai C, Tang H, Chen F, Li D, Lyu F. Characteristics and trends of newly reported HIV infection in young students in China, 2010-2019. Zhonghua Liu Xing Bing Xue Za Zhi. (2020) 41:1455–9. doi: 10.3760/cma.j.cn112338-20200417-00592

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chu J, Zhu Z, Wu S, Li X, Li W, Wei P. HIV-1 genotypes and drug resistance among newly diagnosed young students in Nanjing. AIDS Res Hum Retroviruses. (2020) 40:616–21. doi: 10.1089/AID.2019.0288

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Chen J, Liu Y, Liu S, Yuan D, Su L, Ye L, et al. HIV-1 Drug resistance, distribution of subtypes, and drug resistance-associated mutations in virologic failure individuals in Chengdu, Southwest China, 2014-2016. Biomed Res Int. (2020) 2020:1–8. doi: 10.1155/2020/5894124

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lu X, Zhao H, Zhang Y, Wang W, Zhao C, Li Y, et al. HIV-1 drug-resistant mutations and related risk factors among HIV-1-positive individuals experiencing treatment failure in Hebei Province, China. AIDS Res Ther. (2017) 14:4. doi: 10.1186/s12981-017-0133-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhou Y, Lu J, Wang J, Yan H, Li J, Xu X, et al. Prevalence of HIV antiretroviral drug resistance and its impacts on HIV-1 virological failures in Jiangsu, China: a cross-sectional study. Biomed Res Int. (2016) 2016:1752437. doi: 10.1155/2016/1752437

PubMed Abstract | CrossRef Full Text | Google Scholar

24. CDC C. Chinese guidelines for diagnosis and treatment of HIV/AIDS 2021. Chin J AIDS STD. (2021) 27:1182–201. doi: 10.13419/j.cnki.aids.2021.11.02

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhou Y, Xu X, Lu J, Zhang Z, Wang J, Fu G, et al. A survey of HIV drug resistance threshold in Jiangsu Province in 2015. J Nanjing Med Univ. (2017) 37:256–60. doi: 10.7655/NYDXBNS20170228

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yang W, An W, Ma Y, Li N, Zhu Q, Sun D, et al. Survival effect regarding earlier initiation of antiretroviral treatment on AIDS patients. Chin J Epidemiol. (2014) 35:1065–8. doi: 10.3760/cma.j.issn.0254-6450.2014.09.022

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Xu X, Guo H, Hu H, Fu G, Lu J, Zhou Y, et al. Investigation of HIV-1 drug resistance warning line and epidemic subtype in Jiangsu Province for three consecutive years. J Nanjing Med Univ. (2015) 35:1329–32. doi: 10.7655/NYDXBNS20150933

CrossRef Full Text | Google Scholar

28. Yao Z, Yang C, Zhang M, Yang B. Drug resistance mutations of human immunodeficiency virus-1 CRF08_BC, CRF07_BC and CRF01_AE subtypes. Chin J Infect Dis. (2017) 35:626–9. doi: 10.3760/cma.j.issn.1000-6680.2017.10.012

CrossRef Full Text | Google Scholar

29. Li Y, Han Y, Xie J, Gu L, Li W, Wang H, et al. CRF01_AE subtype is associated with X4 tropism and fast HIV progression in Chinese patients infected through sexual transmission. AIDS. (2014) 28:521–30. doi: 10.1097/QAD.0000000000000125

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zhou Y, Lu J, Zhang Z, Fu G. Risk factors associated with HIV drug resistance among ART virological failure patients taking first-line antiviral treatment from Jiangsu Province. Chin J Dis Control Prev. (2017) 21:1191–1194+1199. doi: 10.16462/j.cnki.zhjbkz.2017.12.002

CrossRef Full Text | Google Scholar

31. Yuan D, Su L, Liu H, Yu H, Yang H, Ye L, et al. Drug-resistance characteristics of CRF01_AE and CRF07_BC subtypes of HIV-1 strains in Sichuan province. Chin J Prev Med. (2015) 49:901–6. doi: 10.3760/cma.j.issn.0253-9624.2015.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Mdodo R, Frazier EL, Dube SR, Mattson CL, Sutton MY, Brooks JT, et al. Cigarette smoking prevalence among adults with HIV compared with the general adult population in the United States: cross-sectional surveys. Ann Intern Med. (2015) 162:335–44. doi: 10.7326/M14-0954

PubMed Abstract | CrossRef Full Text | Google Scholar

33. HIV.gov. Tobacco and HIV. (2020). Available online at: https://www.hiv.gov/hiv-basics/staying-in-hiv-care/other-related-health-issues/smoking (accessed March 8, 2022).

Google Scholar

34. CDC. Smoking and HIV. (2022). Available online at: https://www.cdc.gov/tobacco/campaign/tips/diseases/smoking-and-hiv.html (accessed March 9, 2022).

Google Scholar

35. Qiu F, Liang CL, Liu H, Zeng YQ, Hou S, Huang S, et al. Impacts of cigarette smoking on immune responsiveness: up and down or upside down? Oncotarget. (2017) 8:268–84. doi: 10.18632/oncotarget.13613

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Strzelak A, Ratajczak A, Adamiec A, Feleszko W. Tobacco smoke induces and alters immune responses in the lung triggering inflammation, allergy, asthma and other lung diseases: a mechanistic review. Int J Environ Res Public Health. (2018) 15:1033. doi: 10.3390/ijerph15051033

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Onor IO, Stirling DL, Williams SR, Bediako D, Borghol A, Harris MB, et al. Clinical effects of cigarette smoking: epidemiologic impact and review of pharmacotherapy options. Int J Environ Res Public Health. (2017) 14:1147. doi: 10.3390/ijerph14101147

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Rao P, Ande A, Sinha N, Kumar A, Kumar S. Effects of cigarette smoke condensate on oxidative stress, apoptotic cell death, and HIV replication in human monocytic cells. PLoS ONE. (2016) 11:e0155791. doi: 10.1371/journal.pone.0155791

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ande A, McArthur C, Ayuk L, Awasom C, Achu PN, Njinda A, et al. Effect of mild-to-moderate smoking on viral load, cytokines, oxidative stress, and cytochrome P450 enzymes in HIV-infected individuals. PLoS ONE. (2015) 10:e0122402. doi: 10.1371/journal.pone.0122402

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Zhao L, Li F, Zhang Y, Elbourkadi N, Wang Z, Yu C, et al. Mechanisms and genes involved in enhancement of HIV infectivity by tobacco smoke. Toxicology. (2010) 278:242–8. doi: 10.1016/j.tox.2010.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

41. WHO. Interim WHO Clinical Staging of HIV/AIDS and HIV/AIDS Case Definitions for Surveillance: African Region. Geneva: World Health Organization (2005).

Google Scholar

42. Woldie HG, Ketema DB, Ayana M, Alebel A. Predictors of recovery rate among undernourished HIV-positive adults treated with ready-to-use therapeutic food at debre markos comprehensive specialized hospital: a retrospective cohort study. PLoS ONE. (2021) 16:e0255465. doi: 10.1371/journal.pone.0255465

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Maldey B, Haile F, Shumye A. Outcome of ready to use food therapy among patients on HIV/AIDS care in Mekelle Hospital, Northern Ethiopia: retrospective cohort study. J AIDS Clin Res. (2014) 5:268. doi: 10.4172/2155-6113.1000268

CrossRef Full Text | Google Scholar

44. Siddiqui J, Phillips AL, Freedland ES, Sklar AR, Darkow T, Harley CR. Prevalence and cost of HIV-associated weight loss in a managed care population. Curr Med Res Opin. (2009) 25:1307–17. doi: 10.1185/03007990902902119

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Muzigaba M, Sartorius B, Puoane T, Van Wyk B, Sanders D. The impact of HIV infection and disease stage on the rate of weight gain and duration of refeeding and treatment in severely malnourished children in rural South African hospitals. S Afr J Child Health. (2017) 11:86–92 doi: 10.7196/SAJCH.2017.v11i2.1194

CrossRef Full Text | Google Scholar

46. Leierer G, Grabmeier-Pfistershammer K, Steuer A, Geit M, Sarcletti M, Haas B, et al. Factors associated with low-level viraemia and virological failure: results from the Austrian HIV cohort study. PLoS ONE. (2015) 10:e0142923. doi: 10.1371/journal.pone.0142923

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Grennan JT, Loutfy MR, Su D, Harrigan PR, Cooper C, Klein M, et al. Magnitude of virologic blips is associated with a higher risk for virologic rebound in HIV-infected individuals: a recurrent events analysis. J Infect Dis. (2012) 205:1230–8. doi: 10.1093/infdis/jis104

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Taiwo B, Gallien S, Aga E, Ribaudo H, Haubrich R, Kuritzkes DR, et al. Antiretroviral drug resistance in HIV-1-infected patients experiencing persistent low-level viremia during first-line therapy. J Infect Dis. (2011) 204:515–20. doi: 10.1093/infdis/jir353

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Martin-Blondel G, Saune K, Vu Hai V, Marchou B, Delobel P, Izopet J, et al. Factors associated with a strictly undetectable viral load in HIV-1-infected patients. HIV Med. (2012) 13:568–73. doi: 10.1111/j.1468-1293.2012.01012.x

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Haim-Boukobza S, Morand-Joubert L, Flandre P, Valin N, Fourati S, Sayon S, et al. Higher efficacy of nevirapine than efavirenz to achieve HIV-1 plasma viral load below 1 copy/ml. AIDS. (2011) 25:341–4. doi: 10.1097/QAD.0b013e3283427de3

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Tashima K, Staszewski S, Nelson M, Rachlis A, Skiest D, Stryker R, et al. Efficacy and tolerability of long-term efavirenz plus nucleoside reverse transcriptase inhibitors for HIV-1 infection. AIDS. (2008) 22:275–9. doi: 10.1097/QAD.0b013e3282f21b9d

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Vancoillie L, Demecheleer E, Callens S, Vogelaers D, Vandekerckhove L, Mortier V, et al. Markers associated with persisting low-level viraemia under antiretroviral therapy in HIV-1 infection. J Antimicrob Chemother. (2014) 69:1098–103. doi: 10.1093/jac/dkt484

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Geretti AM, Smith C, Haberl A, Garcia-Diaz A, Nebbia G, Johnson M, et al. Determinants of virological failure after successful viral load suppression in first-line highly active antiretroviral therapy. Antivir Ther. (2008) 13:927–36. doi: 10.1177/135965350801300707

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Maldarelli F, Palmer S, King MS, Wiegand A, Polis MA, Mican J, et al. ART suppresses plasma HIV-1 RNA to a stable set point predicted by pretherapy viremia. PLoS Pathog. (2007) 3:e46. doi: 10.1371/journal.ppat.0030046

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Palmisano L, Giuliano M, Nicastri E, Pirillo MF, Andreotti M, Galluzzo CM, et al. Residual viraemia in subjects with chronic HIV infection and viral load <50 copies/ml: the impact of highly active antiretroviral therapy. AIDS. (2005) 19:1843–7. doi: 10.1097/01.aids.0000188426.87538.ed

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Konstantopoulos C, Ribaudo H, Ragland K, Bangsberg DR, Li JZ. Antiretroviral regimen and suboptimal medication adherence are associated with low-level human immunodeficiency virus viremia. Open Forum Infect Dis. (2015) 2:ofu119. doi: 10.1093/ofid/ofu119

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Fletcher CV, Staskus K, Wietgrefe SW, Rothenberger M, Reilly C, Chipman JG, et al. Persistent HIV-1 replication is associated with lower antiretroviral drug concentrations in lymphatic tissues. Proc Natl Acad Sci U S A. (2014) 111:2307–12. doi: 10.1073/pnas.1318249111

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Kuritzkes DR, Ribaudo HJ, Squires KE, Koletar SL, Santana J, Riddler SA, et al. Plasma HIV-1 RNA dynamics in antiretroviral-naive subjects receiving either triple-nucleoside or efavirenz-containing regimens: ACTG A5166s. J Infect Dis. (2007) 195:1169–76. doi: 10.1086/512619

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Podsadecki TJ, Vrijens BC, Tousset EP, Rode RA, Hanna GJ. Decreased adherence to antiretroviral therapy observed prior to transient human immunodeficiency virus type 1 viremia. J Infect Dis. (2007) 196:1773–8. doi: 10.1086/523704

PubMed Abstract | CrossRef Full Text | Google Scholar