Next Article in Journal
Mexico’s Laboratory-Confirmed Human Case of Infection with the Influenza A(H5N2) Virus
Previous Article in Journal
Genetic Diversity and Spatiotemporal Distribution of SARS-CoV-2 Variants in Guinea: A Meta-Analysis of Sequence Data (2020–2023)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viral Suppression and HIV Drug Resistance Among Patients on Second-Line Antiretroviral Therapy in Selected Health Facility in Ethiopia

1
Infectious Diseases Research Directorate, Ethiopian Public Health Institute, Addis Ababa 1242, Ethiopia
2
Department of Translational Medicine, Lund University, 22185 Lund, Sweden
*
Author to whom correspondence should be addressed.
Viruses 2025, 17(2), 206; https://doi.org/10.3390/v17020206
Submission received: 27 December 2024 / Revised: 22 January 2025 / Accepted: 28 January 2025 / Published: 31 January 2025
(This article belongs to the Section Human Virology and Viral Diseases)

Abstract

:
HIV drug resistance (HIVDR) presents a significant challenge to antiretroviral therapy (ART) success, particularly in resource-limited settings like Ethiopia. This cross-sectional study investigated viral suppression rates and resistance patterns among patients on second-line ART across 28 Ethiopian health facilities. Blood samples collected from 586 participants were analyzed to measure CD4 count and viral load and assess HIVDR in patients experiencing virological failure (VF) (viral load ≥ 1000 copies/mL). Demographic and clinical data were analyzed using logistic regression to identify factors associated with VF. Results showed that 13.82% of participants experienced VF, with 67.57% of genotyped samples exhibiting at least one drug resistance mutation. Resistance to nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs) was detected in 48.64%, 64.86%, and 18.92% of cases, respectively. Dual-class resistance was identified in 48.64% of patients, while triple-class resistance was detected in 18.92%. VF was more likely among students and those with CD4 counts below 200 cells/mm³, but less likely in patients on second-line treatment for 12 months or more. Our findings highlight a substantial HIVDR burden among patients on second-line ART with VF, emphasizing the need for comprehensive HIV care, including adherence support, regular viral load monitoring, and HIVDR testing.

1. Introduction

Antiretroviral therapy (ART) roll-out has significantly improved the prognosis of HIV-infected individuals, transforming HIV from a fatal disease into a manageable chronic condition and potentially preventable condition [1,2]. However, the effectiveness of ART is increasingly threatened by the emergence and transmission of HIV drug resistance (HIVDR) [3,4,5]. This challenge is particularly significant in resource-limited settings where routine viral load (VL) and HIVDR monitoring are not consistently implemented due to financial, infrastructural, and logistical constraints. The lack of regular VL testing can result in undetected treatment failure for prolonged periods, leading to the accumulation of drug resistance mutations (DRMs) and undermining the effectiveness of both the current and future treatment options [3,5]. The absence of resistance testing in HIV treatment management can lead to premature switching to more costly and complex second-line regimens when first-line treatments could still be effective with improved adherence support. It may also result in the selection of suboptimal regimens that foster further resistance evolution, compromising the efficacy of second-line therapy and leaving patients with limited alternatives for long-term HIV management [4,6]. At the population level, insufficient monitoring can lead to the increased transmission of HIVDR strains, which complicates treatment strategies for newly infected individuals and potentially undermines the overall effectiveness of ART programs [5,6].
Ethiopia is one of the Sub-Saharan African countries most affected by the HIV-1 epidemic. In Ethiopia, ART is delivered through a public health approach that utilizes standardized first- and second-line treatment regimens along with simplified laboratory monitoring, which includes at least one viral load test annually. At the time of study enrollment, Ethiopia’s national ART guidelines recommended first-line regimens consisting of two nucleoside reverse transcriptase inhibitors (NRTIs)—tenofovir (TDF) with lamivudine (3TC)—combined with a non-nucleoside reverse transcriptase inhibitor (NNRTI) such as efavirenz (EFV) and nevirapine (NVP). However, the guidelines have evolved over time to incorporate newer, more effective drug combinations. Notably, dolutegravir (DTG), an integrase strand transfer inhibitor (INSTI), has been adopted as a preferred option in first-line regimens [7]. For second-line regimens, the guidelines advised the use of two NRTIs (such as zidovudine (ZDV) or abacavir (ABC) with 3TC) along with a ritonavir-boosted protease inhibitor, either lopinavir/ritonavir (LPV/r) or atazanavir/ritonavir (ATV/r) [7].
At the initiation of this study in 2017, Ethiopia had approximately 414,854 adults receiving ART out of an estimated 610,000 adults living with HIV, according to data from the Ethiopian Federal Ministry of Health. Of the patients receiving ART, approximately 6552 were on protease inhibitor/ritonavir (PI/r)-based second-line regimens [7,8]. Despite significant progress in scaling up ART in Ethiopia, routine virological monitoring and HIVDR testing are not yet standard practices in HIV care. As the number of individuals initiating first-line ART continues to grow, there is an anticipated increase in HIVDR and treatment failure, leading to a significant number of patients transitioning to second-line regimens [3,4]. Patients on second-line ART represent a unique population with histories of ART failure, longer durations of HIV infection, and a higher potential for accumulating HIVDRMs [9]. However, there is a notable lack of data on long-term outcomes and HIVDR patterns among patients on second-line regimens in Ethiopia. This study aims to assess the rate of virological failure (VF), the prevalence of HIVDR, and identify the factors associated with VF among patients on second-line ART in Ethiopia.

2. Methods

2.1. Study Design

In 2017, a cross-sectional study was conducted as part of a national survey to assess the prevalence of HIVDR in Ethiopia. This study included 28 ART sites, selected using probability proportional to size (PPS) sampling to ensure representation across various geographic regions and clinic sizes [10]. After selecting the clinics, we included all consecutive adults (≥18 years) who had been receiving second-line ART for at least six months during the study period based on a predetermined sample size that factored in the proportion of patients on second-line regimens at each facility. Written informed consent was obtained from all participants. Blood samples (10 mL) were collected via venipuncture from each participant for CD4+ T-cell counts, viral load measurements, and HIVDR genotyping. Basic demographic and clinical information were also collected using a standardized questionnaire. Plasma samples were transported on dry ice to the Ethiopian Public Health Institute (EPHI) for viral load testing and long-term storage at −80 °C. VL was quantified using the Abbott RealTime HIV-1 assay (Abbott Molecular Inc., Des Plaines, IL, USA). Samples with VL ≥ 1000 copies/mL were shipped to the National Institute of Respiratory Diseases (INER) laboratory in Mexico for HIVDR testing.

2.2. HIV-1 Sequencing

The HIVDR testing was performed at a WHO-designated laboratory in Mexico City, according to the WHO/HIV ResNet Laboratory Operational Framework [11]. Briefly, HIV RNA was extracted from 1 mL of plasma using the QIAamp Viral RNA Kit (QIAGEN, Valencia, CA, USA). Subsequently, the protease-reverse transcriptase (PR-RT) regions were amplified using validated in-house protocols [12]. The purified PCR fragments were then sequenced using the BigDye technology on the ABI Prism 3730 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Sequence assembly and editing were performed using the RECallV 2.0 HIV-1 sequencing analysis tool (University of British Columbia, Vancouver, BC, Canada) [13]. Sequence quality control was performed using the WHO tool https://sequenceqc-dev.bccfe.ca/who_qc (accessed on 16 March 2024) and the Quality Control program of the Los Alamos HIV sequence database https://www.hiv.lanl.gov/ (accessed on 16 March 2024).

2.3. HIV-1 Drug Resistance Analysis

The HIVDR analysis was conducted using the Stanford HIV Drug Resistance Database (HIVdB v9.1) (https://hivdb.stanford.edu/hivdb/by-mutations, accessed on 29 June 2023). The prevalence of mutation patterns and their frequencies were calculated for each of the following drug classes: NRTIs, NNRTIs, and PIs. Additionally, an HIVDR susceptibility analysis was performed using the HIVdb algorithm, which assigns a drug penalty score for each antiretroviral drug. The cumulative scores were then used to classify the level of resistance into four categories: susceptible (score < 15), low-level resistance (score 15–29), intermediate resistance (score 30–59), or high-level resistance (score ≥ 60).

2.4. HIV-1 Subtype Determination

The HIV-1 subtyping was carried out using the online automated subtyping tools REGA v3.0 [14], COMET [15], and the jumping profile Hidden Markov Model (jpHMM) [16,17]. Subtyping was further validated by maximum likelihood (ML) phylogenetic tree analysis with the reference sequences from HIV-1 subtypes (A–K) and recombinant viruses, which were downloaded from the Los Alamos database (http://www.hiv.lanl.gov, accessed on 3 June 2023). Multiple sequence alignment was performed using MAFFT version 7 [18] and then manually edited using BioEdit V7.0.9.0 to achieve a perfect codon alignment [19,20]. The ML tree topology was constructed using the online version of PhyML v3.0 [21] with the GTR + I + Γ nucleotide-substitution model (utilizing the estimated proportion of invariable sites and four gamma categories) [16]. A heuristic tree search was conducted using the SPR branch-swapping algorithm. Branch support was evaluated with the approximate likelihood ratio test (aLRT-SH, Shimodaira Hasegawa-like) [22]. Clusters were defined as monophyletic clades with aLRT-SH support ≥ 0.9 [16]. The subtype-resolved ML phylogeny trees were visualized using the FigTree v1.4.0 program [23]. Sequences forming clusters with the reference sequences of the same subtype were assigned to that subtype.

2.5. Data Quality Management

Data quality assurance was rigorously maintained through a multi-step process. Every tenth entry of the retrospectively collected data was cross-checked, and all case and laboratory report forms were verified against the source data to ensure accuracy. Following this initial verification, the data were independently entered into Epi-Info V6.04 by two separate data clerks. This software was equipped with check programs to prevent the entry of erroneous data, adding an additional layer of quality control. Furthermore, manual consistency checks were performed to further enhance data integrity. Once these quality control measures were completed, the data were exported to SPSS v.26 for comprehensive analysis.

2.6. Statistical Analysis

Descriptive statistics were used to summarize demographic and clinical characteristics. The prevalence of VF (defined as having a viral load ≥ 1000 copies/mL) and acquired HIVDR was calculated, and associations between clinical characteristics and virological outcomes were assessed using logistic regression models. Different patient-related characteristics, including socio-demographic, clinical, immunological, and virological characteristics, were included in the analysis. Initially, variables were examined individually, and those demonstrating an association with VF at a p-value of < 0.2 were included in a multivariable logistic regression model. Prior to conducting the multivariable analysis, correlation analysis was performed to evaluate potential multicollinearity among variables, and no multicollinearity was detected. In the final multivariable model, predictors with p-values < 0.05 were considered statistically significant. Crude and adjusted odds ratios (ORs) with 95% confidence intervals (CIs) were reported to quantify the strength of the associations.

2.7. Ethical Considerations

The study protocol received approval from the Scientific and Ethical Research Office of the Ethiopian Public Health Institute (EPHI) and the Institutional Review Board of the Ethiopian Ministry of Science and Technology.

3. Results

3.1. Demographic Characteristics of Study Participants

A total of 586 participants were included in this study. Females comprised a greater portion of the sample at 54.44%, compared to males at 45.56%. The largest age group was individuals aged between 35 and 44 years (41.81%), followed by those aged 25 to 34 years (27.30%), over 45 years (24.91%), and 18 to 24 years (5.97%). Most participants resided in urban areas (84.47%), while 15.53% lived in rural regions. Regarding marital status, 39.93% were in monogamous marriages, 21.84% were unmarried, 19.80% were divorced, 9.39% were widowed, and 9.04% were in polygamous marriages. Educational attainment varied as follows: 44.20% had primary education, 23.21% had secondary education, 12.63% had college or university education, and 19.97% had no formal education. Occupational statuses were diverse, with housewives constituting 18.09%, government employees 17.06%, merchants 16.38%, unemployed individuals 15.36%, and students 4.27%. Other occupations, including drivers and various other jobs, accounted for 24.74% of the participants (Table 1).

3.2. ART Regimen and Clinical Characteristics of Patients on Second-Line Therapy

At the initiation of first-line ART, the most frequently prescribed regimen was AZT + 3TC + NVP, accounting for 26.79% of patients. This was followed by TDF + 3TC + EFV at 18.09%. Other combinations included D4T + 3TC + NVP (15.19%) and TDF + 3TC + NVP (14.33%). Less commonly used regimens were AZT + 3TC + EFV (12.46%), D4T + 3TC + EFV (5.12%), ABC + 3TC + NVP (6.14%), and ABC + 3TC + EFV (1.88%). The majority of patients (69.97%) switched after more than 91 months on first-line ARVs, while smaller proportions switched earlier as follows: 1.71% before 30 months, 7.17% between 31 and 60 months, and 21.16% between 61 and 90 months.
The most commonly prescribed second-line regimens were TDF + 3TC + ATV/r, used by 44.20% of participants, followed by ABC + 3TC + ATV/r (19.45%) and AZT + 3TC + ATV/r (17.58%). Smaller proportions were on TDF + 3TC + LPV/r (9.73%), AZT + 3TC + LPV/r (5.63%), and ABC + 3TC + LPV/r (3.41%). The duration of therapy on the second-line regimen varied as follows: the majority (55.80%) had been on second-line therapy for 16–48 months, 29.69% for 6–15 months, and 14.51% for over 49 months. At the initiation of ART, all participants had a CD4 count below 200 cells/mm3. By the time they switched to second-line treatment, 70.30% still had a CD4 count below 200 cells/mm3. However, at the time of this study, 71.00% of participants had achieved CD4 counts above 200 cells/mm3, while 29.00% still had CD4 counts below this threshold. Viral suppression, defined as an RNA viral load below 1000 copies/mL, was observed in 86.18% of participants, as shown in Table 2.

3.3. HIV Drug Resistance Mutation Profile

Among the 81 (13.82%) participants experiencing VF, genotyping was successfully performed for 37 individuals. Of these, 25 participants (67.56%) harbored at least one DRM. Dual-class DRMs, affecting both NRTI and NNRTI and/or PI classes, were observed in 18 participants (48.64%), while triple-class DRMs (involving PI, NRTI, and NNRTI) were identified in 7 participants (18.92%) (see Figure 1).
NRTI-associated DRMs were detected in 18 participants (48.64%), with M184V being the most frequent NRTI mutation, occurring in 22.86% of patients. Thymidine analogue mutations (TAMs) were detected in 38.57% of patients, with specific mutations, including T215Y/V (10.00%), K70E (8.57%), D67N (8.57%), and K219E (7.14%), while K65R and Y115F were found in 7.14% and 8.57%, respectively. NNRTI-associated DRMs were detected in 18 patients (48.64%). The most predominant NNRTI-associated mutations detected included K103N (18.87%), Y181C (13.21%), and G190A (13.21%). Additional mutations such as H221Y (9.43%), V106I/M (5.66%), V108I (5.66%), A98G (7.55%), and L100I (5.66%) were also detected. PI-associated DRMs were identified in 7 patients (18.92%), with the V82A mutation being the most prevalent, observed in 29.41% of cases. The I54V and M46I mutations were each detected in 23.53% of participants, while I50L (11.76%), L90M (5.88%), and N88S (5.88%) mutations were also observed (Table 3).

3.4. Factors Associated with Virological Failure

Although various factors were examined, occupation, duration on second-line ART, and current CD4 count were identified as factors associated with second-line ART VF. Among occupational groups, students demonstrated a markedly higher likelihood of experiencing VF, with an adjusted odds ratio (AOR) of 6.14 (95% CI: 1.93–19.48) compared to government employees, the reference group. Other occupations, including drivers (AOR = 2.10, 95% CI: 0.61–7.20), housewives (AOR = 1.02, 95% CI: 0.33–3.19), unemployed individuals (AOR = 1.67, 95% CI: 0.65–4.31), and those in the “other” category (AOR = 1.57, 95% CI: 0.67–3.68), did not show statistically significant associations. The duration of second-line ART was another critical factor associated with VF. Participants who had been on second-line ART for over 12 months had a significantly lower likelihood of virological failure (AOR = 0.52, 95% CI: 0.28–0.97) compared to those on ART for less than 12 months. Additionally, current CD4 count emerged as a significant determinant of outcomes. Participants with a current CD4 count of less than 200 cells/mm³ were more likely to experience VF (AOR = 2.55, 95% CI: 1.52–4.28) compared to those with a CD4 count at or above 200 cells/mm3 (Table 4).

3.5. HIV Drug Resistance Susceptibility Testing

Our antiretroviral drug resistance analysis reveals varying degrees of susceptibility across the three primary ARV classes (NRTIs, NNRTIs, and PIs). Among NRTIs, drugs like AZT and TDF showed higher susceptibility rates at 78.38% and 62.16%, respectively, while drugs like emtricitabine (FTC) and 3TC had significant resistance, with 43.24% exhibiting high levels of resistance. Drugs like D4T and DDI had 59.46% susceptibility, with a notable portion (24.32% and 27.03%, respectively) having high resistance. Among NNRTIs, DOR and ETR had 54.05% susceptibility, with high-level resistance in 8.11% of samples. EFV and NVP showed lower susceptibility rates of 35.14%, accompanied by high-level resistance in 48.65% and 59.46% of samples, respectively. RPV exhibited 51.35% susceptibility, while high-level resistance was observed in 37.84% of samples. PIs generally maintained high levels of susceptibility. DRV was fully susceptible (100%), while ATV demonstrated susceptibility in 81.08% of samples and high-level resistance in 16.22%. FPV, IDV, LPV, NFV, and SQV exhibited susceptibility rates of 83.78% to 86.49%, with high-level resistance detected in 5.41% to 13.51% of samples. TPV showed 86.49% susceptibility, with no case of high-level resistance (Table 5).

3.6. Maximum-likelihood Phylogenetic Tree

The phylogenetic tree in Figure 2 comprises 459 sequences, including 37 Ethiopian sequences and 422 reference sequences for HIV-1 subtypes (A–K) and circulating recombinant forms, which were downloaded from the HIV-1 LANL database. An ML tree was constructed using the online version of PhyML v 3.0. In the figure, the reference sequences from the Los Alamos National Laboratory are shown in black. Ethiopian sequences clustering with the HIV-1 subtype C reference sequence are depicted in blue, while non-subtype C Ethiopian sequences are shown in pink. Phylogenetic analysis of the sequences showed that 94.59% (35 out of 37) clustered with HIV-1 subtype C. The remaining two sequences, accounting for 2.70% each, were identified as non-subtype C, specifically subtype B, and a recombinant form of subtypes D and A1 (see Figure 2).

4. Discussions

This study provides the first national data on viral suppression and HIVDR among patients on second-line ART in selected health facilities in Ethiopia. Our findings showed that 13.82% of patients on second-line ART experienced VF, defined as a VL ≥ 1000 copies/mL. Notably, 32.43% of study participants did not harbor any DRMs, suggesting that poor adherence to treatment may be the primary cause of VF in these cases. Among patients with VF, 67.57% exhibited at least one class of DRMs, with 48.64% showing dual-class DRMs and 18.921% displaying triple-class DRMs. The low PI-associated DRMs (18.92%) and high susceptibility to drugs like darunavir, atazanavir, and lopinavir suggest that PI-based regimens remain an effective backbone for second-line treatments in Ethiopia if adherence issues are addressed. We also identified risk factors for VF, including being a student (AOR = 6.14; 95% CI: 1.93–19.48), shorter duration on second-line treatment (less than 12 months, AOR = 0.52; 95% CI: 0.28–0.97), and a CD4 count below 200 cells/mm³ (AOR = 2.55, 95% CI: 1.52–4.28).
Our study showed that 13.82% of patients on second-line antiretroviral therapy (ART) experienced VF (viral load ≥1000 copies/mL). This finding is consistent with VF rates reported in several other African countries, including Tanzania, Rwanda, Uganda, South Africa, and Ethiopia, where rates ranged from 12% to 41% [24,25,26,27,28,29,30,31]. A systematic review and meta-analysis focusing on patients undergoing second-line ART regimens in resource-limited settings further corroborated our results, demonstrating cumulative VF rates between 21.8% and 38.0% [25]. Our findings, along with those from other studies, underscore the persistent challenge of maintaining viral suppression in patients on second-line therapy across various resource-limited settings.
Among our study participants, 32.4% did not harbor any DRMs, suggesting that poor adherence to treatment may be the primary cause of VF. This is consistent with previous research in resource-limited settings, which revealed that 33–67% of patients on second-line PI-based regimens with VF had no major DRMs, suggesting adherence issues rather than drug resistance as the primary cause of failure [10,24,29,32,33]. For example, among patients with VF, 37% in Rwanda [29], 58.5% in Namibia [10], and 33% in South Africa [32] did not harbor any DRMs. Identified risk factors for second-line failure include suboptimal adherence, prior exposure to non-standard regimens, and pre-existing PI resistance [34,35]. However, it also raises concerns that some patients may have been prematurely switched to second-line regimens without resistance testing, underscoring the need for resistance testing to ensure optimal treatment decisions.
Among patients experiencing VF, a significant proportion demonstrated complex resistance patterns. Dual-class DRMs affecting both NRTI and NNRTI classes and/or protease inhibitors were identified in 48.64% of patients. Even more concerning was the presence of triple-class DRMs in 18.92% of cases, involving mutations across all the following three major drug classes: PIs, NRTIs, and NNRTIs. These results are consistent with previous studies in resource-limited settings, where dual-class resistance in patients on second-line ART was found in 40% to 60% of cases [29,36]. Furthermore, the observed rate of triple-class resistance (18.92%) in this cohort closely aligns with findings from other studies in sub-Saharan Africa, where 22% of participants experiencing failure on second-line therapy in routine care settings also showed triple-class resistance [36]. However, this rate is lower than that reported in a study from India, which showed that 53% of patients harbored triple-class resistance [36]. The findings from this study, along with evidence from prior research, highlight a persistent and widespread issue of HIVDR among patients on second-line ART across various resource-constrained environments. The high prevalence of dual- and triple-class DRMs observed highlights the significant challenges in managing second-line ART failures in these settings. Such multi-class resistance substantially limits future treatment options, potentially necessitating the use of third-line regimens that are typically more expensive and less readily available. These results emphasize the critical need for enhanced viral load monitoring, resistance testing, and improved access to third-line ART in Ethiopia and similar resource-limited contexts [37,38]. Additionally, our findings underscore the importance of assessing and supporting medication adherence, as poor adherence remains a major contributor to VF [29].
We observed a 48.64% prevalence of NRTI-associated DRMs, with the most common being M184V and TAMs, present in 22.86% and 38.57% of patients, respectively. M184V is known to emerge rapidly in patients receiving 3TC or FTC, which are common components of first-line regimens [39,40]. The prevalence of TAMs (T215Y/V, K70E, D67N, and K219E) is particularly concerning, as these mutations accumulate over time and can confer cross-resistance to multiple NRTIs, significantly impacting the effectiveness of NRTI-based regimens [39]. Additionally, the K65R mutation, detected in 7.14% of patients, raises additional concern due to its association with broad resistance, especially to TDF and ABC, which are critical components of second-line regimens in resource-limited settings like Ethiopia. Similar resistance patterns have been documented in studies conducted among patients on second-line ART in various African countries, including Rwanda, Uganda, Namibia, Gabon, South Africa, and Zambia [29,41,42,43,44]. These complex mutation patterns, characterized by high rates of M184V mutations and TAMs, are typically attributed to several common factors in these settings as follows: limited availability of resistance testing, prolonged periods of VF while on first-line therapy, the use of suboptimal treatment regimens, and challenges with medication adherence [45,46]. These circumstances collectively contribute to the development and accumulation of DRMs, potentially compromising the efficacy of second-line ART regimens and complicating long-term HIV management in these regions.
We detected NNRTI-associated DRMs in 65% of patients, even though NNRTIs are typically excluded from second-line ART regimens. This indicates persistence of NNRTI-resistant viral strains even after discontinuation of NNRTI-based therapies [39]. Similar results from other studies have shown that NNRTI-associated DRMs often remain detectable in patients transitioning to second-line regimens following failure of first-line treatments [43,44,45,47]. The most frequently observed NNRTI-associated mutations include K103N, Y181C, and G190A. These mutations are known to confer high-level resistance to first-generation NNRTIs such as efavirenz and nevirapine, which are commonly used in first-line regimens [46]. Furthermore, the additional mutations detected, such as H221Y, V106M, V108I, A98G, and L100I, can contribute to cross-resistance within the NNRTI class and may impact the effectiveness of newer NNRTIs like etravirine and rilpivirine [48]. Our findings are consistent with findings from several studies conducted in resource-limited settings [44,49,50,51].
Our study’s finding of an 18.92% prevalence of PI-associated DRMs among patients on second-line regimens is supported by previous research in similar settings [44,45]. A pooled analysis of PI-based treatment failures in sub-Saharan Africa reported that 17% of patients had at least one major PI-resistance mutation at the time of treatment failure [43]. Studies from Namibia, Uganda, and South Africa reported 13.1%, 19.4%, and 22% of PI-associated HIVDRMs, respectively, among people receiving PI/r-based second-line ART [29,41,44,52]. However, our results are lower than those reported in other studies. For instance, Chimukangara et al. found in South Africa that among 348 samples analyzed, 287 (82.5%) had at least one DRM and 114 (32.8%) had at least one major PI-resistance DRM [53]. Similarly, a systematic review conducted in Asia by Ross et al. revealed that 13 out of 39 patients (33%) had major PI DRMs [54].
The V82A mutation, detected in 29.41% of patients, is a major PI resistance mutation that can significantly reduce susceptibility to several PIs, including lopinavir and indinavir. The I54IV/V and M46I mutations, each found in 23.53% of patients, are also major PI resistance mutations that can contribute to reduced PI efficacy. The presence of I50L, L90M, and N88S mutations, albeit at lower frequencies, further underscores the complexity of PI resistance patterns in this cohort. The I50L mutation is particularly associated with resistance to atazanavir, while L90M can confer broad cross-resistance to multiple PIs [25,36]. Several studies among patients receiving second-line PI-based regimens have documented similar PI-associated mutations [42,43,44,45].
The relatively low prevalence of PI mutations (18.92%) and the high level of susceptibility to DRV, ATV, and LPV among the majority of participants in our study suggest that PI-based regimens may still be an effective backbone of most second-line regimens for many patients failing second-line therapy, provided adherence issues are addressed. However, the presence of PI mutations in some patients underscores the need for careful monitoring and resistance testing to ensure these patients are prescribed effective third-line ART regimens.
Our result highlights a significant association between the duration of second-line ART and viral load suppression. Participants who had been on second-line ART for over 12 months had significantly lower odds of viral load failure compared to those treated for less than 12 months (AOR: 0.52, 95% CI: 1.52–4.28). The improved outcomes associated with longer ART duration may be attributed to enhanced medication adherence over time [55]. As patients become more accustomed to their treatment regimens and learn to manage side effects, adherence likely improves, resulting in better viral suppression. Additionally, extended periods of viral suppression reduce the risk of DRMs, which are more likely to arise when viral replication is poorly controlled. This is especially important in the context of second-line therapy, where treatment options become more limited, and maintaining the effectiveness of available drugs is crucial. Early identification of at-risk individuals through viral load monitoring and adherence support can improve long-term outcomes. Consistent with previous research, our findings suggest that providing intensive adherence support and regular VL monitoring during the early stages of second-line therapy is critical to enhance the likelihood of achieving and maintaining viral suppression and treatment success over time [45,53]. This is more significant in resource-limited settings where third-line ART options are limited.
Our study also reveals significant associations between occupation and VF among participants on second-line ART. Students emerged as the group with the highest risk (AOR: 6.14, 95% CI: 1.93–19.48) compared to government employees. Students may face unique adherence difficulties, such as balancing academic responsibilities, social life, and health management [56]. Furthermore, as young adults, they may experience greater stigma or lack of access to consistent healthcare, making adherence more challenging [56]. Additionally, stress, depression, and anxiety are prevalent among students and can disrupt consistent ART adherence. Our findings underscore the importance of designing interventions tailored to students’ unique circumstances, including flexible healthcare delivery models (e.g., mobile health units and telemedicine), enhanced adherence support, and psychosocial counseling to address their specific challenges [57].
Although not statistically significant, individuals in occupations such as drivers, unemployed persons, and those classified under “other” categories demonstrated increased odds of VF. For example, drivers often contend with irregular work schedules and extended periods away from home, which can disrupt medication adherence and hinder regular clinic attendance. To address these challenges, adherence interventions for drivers could involve strategies such as providing multi-month medication supplies or offering flexible clinic appointment times. These measures would help ensure that drivers can maintain their treatment regimens despite the demands and unpredictability of their work environments [58]. Overall, these findings highlight the critical importance of ensuring comprehensive adherence support for all patients on second-line ART, regardless of their occupation.
Our results showed a relationship between CD4 count and the likelihood of viral load failure among individuals on second-line ART. Specifically, participants with a CD4 count below 200 cells/mm3 were found to have 2.55 times higher odds (AOR: 2.55; 95% CI: 1.52–4.28) of experiencing VF compared to those with CD4 counts at and above this threshold. This suggests that individuals with lower CD4 counts, indicative of more pronounced immunosuppression, are at heightened risk of treatment failure. Our findings align with other studies that similarly report an increased risk of VF associated with low baseline CD4 counts, which is often accompanied by HIVDR [27,28,30,59,60].
The genetic diversity of HIV-1 in the study group aligns with the typical diversity observed in different studies in Ethiopia. It highlights the clear dominance of subtype C in the country, with 94.59% (35 out of 37) of the sequences clustering within this subtype [17,61,62,63].

Limitation of This Study

This study marks a significant milestone as the first national survey to encompass multiple health facilities, yet it is not without its limitations. The relatively small sample size poses a challenge, potentially failing to accurately represent the broader population of patients on second-line therapy, which in turn restricts the generalizability of the findings. We defined VF as a single viral load measurement of ≥1000 copies/mL, without the requirement for a second confirmatory test. This approach could potentially overestimate the incidence of VF. The overall genotyping success rate was 45.68%, which might have influenced the study results. Moreover, the cross-sectional nature of the study design only provides a snapshot of resistance mutations at a single point in time, preventing any assessment of longitudinal changes or the progression of these mutations over an extended period. To address these limitations and gain a more comprehensive understanding of the evolution of resistance mutations and their clinical implications among patients on second-line regimens, future research should focus on incorporating larger, more diverse cohorts and implementing longitudinal follow-up strategies.

5. Conclusions

This study, the first national HIVDR survey conducted in selected health facilities in Ethiopia, indicates that 13.82% of patients on second-line ART did not achieve viral suppression. This study identified both dual-class and triple-class resistance among those experiencing virological failure, indicating a substantial challenge in managing HIVDR. This underscores the need for comprehensive HIV care in Ethiopia, which should include adherence support, regular viral load monitoring, and HIVDR testing for patients on second-line ART. Notably, the relatively low prevalence of PIs mutations and the high level of susceptibility to PIs suggest that PI-based regimens could remain effective as the backbone of second-line regimens, provided that adherence issues are adequately addressed.

Author Contributions

Conceptualization, K.Z., A.G., P.M. and D.A.A.; methodology, K.Z., A.G. and D.A.A.; data curation, K.Z. and D.A.A.; formal analysis, D.A.A., K.Z. and P.M.; funding acquisition, K.Z., A.G. and D.A.A.; investigation, K.Z., D.A.A., A.G., S.S., Y.G., E.K., B.W., A.Y. (Ajanaw Yizengaw), G.G., S.A., M.Y., A.Y. (Amelework Yilma) and S.T.; project administration, K.Z., A.G. and D.A.A.; supervision, K.Z., A.G. and D.A.A.; visualization, K.Z. and D.A.A.; writing—original draft, K.Z. and D.A.A.; writing—review and editing, K.Z., S.S., P.M. and D.A.A. All authors revised the manuscript and provided important intellectual content. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the WHO Ethiopia Country Office and the Ethiopian Public Health Institute (EPHI) under WHO-HIVDR-5385 and by the Swedish Research Council (grant number 2020-02344), the Österlund Foundation, and a donation through the Medical Faculty of Lund.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Research and Ethical Clearance Committee of EPHI and the National Health Research Ethics Review Committee of the Ministry of Science and Technology of Ethiopia (SERO-60-2016; Version 001, on 11 May 2016).

Informed Consent Statement

Written informed consent was obtained from all participants for the national HIV drug-resistance survey.

Data Availability Statement

All the sequences from this study will be deposited in GenBank.

Acknowledgments

We extend our sincere gratitude to the study participants and all the staff who contributed to the survey. We also acknowledge the support provided by the WHO and the National Institute of Respiratory Diseases-Mexico (INER) laboratory for HIVDR genotyping.

Conflicts of Interest

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. McMahon, R. UNAIDS Issues New Fast-Track Strategy to End AIDS by 2030; Elizabeth Glaser Pediatric AIDS Foundation: Washington, DC, USA, 2014. [Google Scholar]
  2. WHO. Policy Brief: Update of Recommendations on First- and Second-Line Antiretroviral Regimens. 2019. Available online: https://apps.who.int/iris/handle/10665/325892 (accessed on 20 September 2023).
  3. Gregson, J.; Tang, M.; Ndembi, N.; Hamers, R.L.; Rhee, S.-Y.; Marconi, V.C.; Diero, L.; Brooks, K.A.; Theys, K.; Rinke de Wit, T.; et al. Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: A multicentre retrospective cohort study. Lancet Infect. Dis. 2016, 16, 565–575. [Google Scholar] [CrossRef] [PubMed]
  4. World Health Organization. Global Action Plan on HIV Drug Resistance 2017–2021; World Health Organization: Geneva, Switzerland, 2017. [Google Scholar]
  5. Bertagnolio, S.; Beanland, R.L.; Jordan, M.R.; Doherty, M.; Hirnschall, G. The World Health Organization’s response to emerging human immunodeficiency virus drug resistance and a call for global action. J. Infect. Dis. 2017, 216, S801–S804. [Google Scholar] [CrossRef] [PubMed]
  6. Bernabé, K.J.; Siedner, M.; Tsai, A.C.; Marconi, V.C.; Murphy, R.A. Detection of HIV Virologic Failure and Switch to Second-Line Therapy: A Systematic Review and Meta-analysis of Data from Sub-Saharan Africa. Open Forum Infect. Dis. 2022, 9, ofac121. [Google Scholar] [CrossRef] [PubMed]
  7. Federal Ministry of Health of Ethiopia. National Consolidated Guidelines for Comprehensive HIV Prevention, Care and Treatment; Fmoh: Addis Ababa, Ethiopia, 2018; pp. 1–238. [Google Scholar]
  8. Ethiopian Public Health Institute (EPHI). Ethiopia Population-based HIV Impact Assessment (EPHIA) 2017–2018. Final Report; EPHI: Addis Ababa, Ethiopia, 2020. [Google Scholar]
  9. Sawadogo, S.; Shiningavamwe, A.N.; Roscoe, C.; Baughman, A.L.; Negussie, T.; Mutandi, G.; Yang, C.; Hamunime, N.; Agolory, S. Human Immunodeficiency Virus-1 Drug Resistance Patterns Among Adult Patients Failing Second-Line Protease Inhibitor-Containing Regimens in Namibia, 2010–2015. Open Forum Infect. Dis. 2018, 5, ofy014. [Google Scholar] [CrossRef] [PubMed]
  10. World Health Organization. Surveillance of HIV Drug Resistance in Populations Initiating Antiretroviral Therapy (Pre-Treatment HIV Drug Resistance): Concept Note; World Health Organization: Geneva, Switzerland, 2014. [Google Scholar]
  11. World Health Organization. WHO HIVResNet HIV Drug Resistance Laboratory Operational Framework; World Health Organization: Geneva, Switzerland, 2020. [Google Scholar]
  12. Ávila-Ríos, S.; García-Morales, C.; Valenzuela-Lara, M.; Chaillon, A.; Tapia-Trejo, D.; Pérez-García, M.; López-Sánchez, D.M.; Maza-Sánchez, L.; Del Arenal-Sánchez, S.J.; Paz-Juárez, H.E.; et al. HIV-1 drug resistance before initiation or re-initiation of first-line ART in eight regions of Mexico: A sub-nationally representative survey. J. Antimicrob. Chemother. 2019, 74, 1044–1055. [Google Scholar] [CrossRef]
  13. Woods, C.K.; Brumme, C.J.; Liu, T.F.; Chui, C.K.; Chu, A.L.; Wynhoven, B.; Hall, T.A.; Trevino, C.; Shafer, R.W.; Harrigan, P.R. Automating HIV drug resistance genotyping with RECall, a freely accessible sequence analysis tool. J. Clin. Microbiol. 2012, 50, 1936–1942. [Google Scholar] [CrossRef]
  14. Pineda-Peña, A.C.; Faria, N.R.; Imbrechts, S.; Libin, P.; Abecasis, A.B.; Deforche, K.; Gómez-López, A.; Camacho, R.J.; de Oliveira, T.; Vandamme, A.M. Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: Performance evaluation of the new REGA version 3 and seven other tools. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2013, 19, 337–348. [Google Scholar] [CrossRef]
  15. Struck, D.; Lawyer, G.; Ternes, A.M.; Schmit, J.C.; Bercoff, D.P. COMET: Adaptive context-based modeling for ultrafast HIV-1 subtype identification. Nucleic Acids Res. 2014, 42, e144. [Google Scholar] [CrossRef]
  16. Arimide, D.A.; Esquivel-Gómez, L.R.; Kebede, Y.; Sasinovich, S.; Balcha, T.; Björkman, P.; Kühnert, D.; Medstrand, P. Molecular Epidemiology and Transmission Dynamics of the HIV-1 Epidemic in Ethiopia: Epidemic Decline Coincided with Behavioral Interventions Before ART Scale-Up. Front. Microbiol. 2022, 13, 821006. [Google Scholar] [CrossRef]
  17. Arimide, D.A.; Abebe, A.; Kebede, Y.; Adugna, F.; Tilahun, T.; Kassa, D.; Assefa, Y.; Balcha, T.T.; Björkman, P.; Medstrand, P. HIV-genetic diversity and drug resistance transmission clusters in Gondar, Northern Ethiopia, 2003–2013. PLoS ONE 2018, 13, e0205446. [Google Scholar] [CrossRef]
  18. Katoh, K.; Kuma, K.-I.; Miyata, T.; Toh, H. Improvement in the accuracy of multiple sequence alignment program MAFFT. Genome Inform. 2005, 16, 22–33. [Google Scholar] [PubMed]
  19. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In Nucleic Acids Symposium Series; Academia: Oxford, UK, 1999; pp. 95–98. [Google Scholar]
  20. Arimide, D.A.; Szojka, Z.I.; Zealiyas, K.; Gebreegziabxier, A.; Adugna, F.; Sasinovich, S.; Björkman, P.; Medstrand, P. Pre-Treatment Integrase Inhibitor Resistance and Natural Polymorphisms among HIV-1 Subtype C Infected Patients in Ethiopia. Viruses 2022, 14, 729. [Google Scholar] [CrossRef] [PubMed]
  21. Guindon, S.; Dufayard, J.-F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed]
  22. Guindon, S.; Lethiec, F.; Duroux, P.; Gascuel, O. PHYML Online—A web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33, W557–W559. [Google Scholar] [CrossRef]
  23. Rambaut, A. FigTree 1.4. 3. 2016. Available online: http://tree.bio.ed.ac.uk/software/figtree (accessed on 16 March 2023).
  24. Hosseinipour, M.C.; Gupta, R.K.; Van Zyl, G.; Eron, J.J.; Nachega, J.B. Emergence of HIV drug resistance during first- and second-line antiretroviral therapy in resource-limited settings. J. Infect. Dis. 2013, 207 (Suppl. 2), S49–S56. [Google Scholar] [CrossRef]
  25. Ajose, O.; Mookerjee, S.; Mills, E.J.; Boulle, A.; Ford, N. Treatment outcomes of patients on second-line antiretroviral therapy in resource-limited settings: A systematic review and meta-analysis. Aids 2012, 26, 929–938. [Google Scholar] [CrossRef]
  26. Murphy, R.A.; Sunpath, H.; Castilla, C.; Ebrahim, S.; Court, R.; Nguyen, H.; Kuritzkes, D.; Marconi, V.C.; Nachega, J.B. Second-line antiretroviral therapy: Long-term outcomes in South Africa. J. Acquir Immune Defic. Syndr. 2012, 61, 158–163. [Google Scholar] [CrossRef]
  27. Mwavika, E.T.; Kunambi, P.P.; Masasi, S.J.; Lema, N.; Kamori, D.; Matee, M. Prevalence, rate, and predictors of virologic failure among adult HIV-Infected clients on second-line antiretroviral therapy (ART) in Tanzania (2018–2020): A retrospective cohort study. Bull. Natl. Res. Cent. 2024, 48, 96. [Google Scholar] [CrossRef]
  28. Gunda, D.W.; Kilonzo, S.B.; Mtaki, T.; Bernard, D.M.; Kalluvya, S.E.; Shao, E.R. Magnitude and correlates of virological failure among adult HIV patients receiving PI based second line ART regimens in north western Tanzania; A case control study. BMC Infect. Dis. 2019, 19, 235. [Google Scholar] [CrossRef]
  29. Ndahimana, J.d.A.; Riedel, D.J.; Muhayimpundu, R.; Nsanzimana, S.; Niyibizi, G.; Mutaganzwa, E.; Mulindabigwi, A.; Baribwira, C.; Kiromera, A.; Jagodzinski, L.L.; et al. HIV drug Resistance Mutations among Patients Failing Second-Line Antiretroviral Therapy in Rwanda. Antivir. Ther. 2016, 21, 253–259. [Google Scholar] [CrossRef]
  30. Zakaria, H.F.; Raru, T.B.; Hassen, F.A.; Ayana, G.M.; Merga, B.T.; Debele, G.R.; Kiflemariam, G.; Kebede, S.A.; Ayele, T.A. Incidence and Predictors of Virological Failure Among Adult HIV/AIDS Patients on Second-Line Anti-Retroviral Therapy, in Selected Public Hospital of Addis Ababa, Ethiopia: Retrospective Follow-Up Study. HIV AIDS (Auckl.) 2022, 14, 319–329. [Google Scholar] [CrossRef] [PubMed]
  31. Levison, J.H.; Orrell, C.; Losina, E.; Lu, Z.; Freedberg, K.A.; Wood, R. Early Outcomes and the Virological Effect of Delayed Treatment Switching to Second-Line Therapy in An Antiretroviral Roll-Out Programme in South Africa. Antivir. Ther. 2011, 16, 853–861. [Google Scholar] [CrossRef] [PubMed]
  32. Levison, J.H.; Orrell, C.; Gallien, S.; Kuritzkes, D.R.; Fu, N.; Losina, E.; Freedberg, K.A.; Wood, R. Virologic Failure of Protease Inhibitor-Based Second-Line Antiretroviral Therapy without Resistance in a Large HIV Treatment Program in South Africa. PLoS ONE 2012, 7, e32144. [Google Scholar] [CrossRef] [PubMed]
  33. Moore, D.M.; Nanfuka, M.; Zhang, W.; Okoboi, S.; Kaleebu, P.; Kapaata, A.; Zhu, J.; Tibengana, B.; Birungi, J. Limited Impact of First-Line Drug Resistance Mutations on Virologic Response Among Patients Receiving Second-Line Antiretroviral Therapy in Rural Uganda. JAIDS J. Acquir. Immune Defic. Syndr. 2023, 92, 84–88. [Google Scholar] [CrossRef] [PubMed]
  34. Boender, T.S.; Hamers, R.L.; Ondoa, P.; Wellington, M.; Chimbetete, C.; Siwale, M.; Labib Maksimos, E.E.F.; Balinda, S.N.; Kityo, C.M.; Adeyemo, T.A.; et al. Protease Inhibitor Resistance in the First 3 Years of Second-Line Antiretroviral Therapy for HIV-1 in Sub-Saharan Africa. J. Infect. Dis. 2016, 214, 873–883. [Google Scholar] [CrossRef]
  35. Musana, H.; Ssensamba, J.T.; Nakafeero, M.; Mugerwa, H.; Kiweewa, F.M.; Serwadda, D.; Ssali, F. Predictors of failure on second-line antiretroviral therapy with protease inhibitor mutations in Uganda. AIDS Res. Ther. 2021, 18, 17. [Google Scholar] [CrossRef]
  36. Saravanan, S.; Vidya, M.; Balakrishnan, P.; Kantor, R.; Solomon, S.S.; Katzenstein, D.; Kumarasamy, N.; Yeptomi, T.; Sivamalar, S.; Rifkin, S.; et al. Viremia and HIV-1 Drug Resistance Mutations Among Patients Receiving Second-Line Highly Active Antiretroviral Therapy in Chennai, Southern India. Clin. Infect. Dis. 2012, 54, 995–1000. [Google Scholar] [CrossRef]
  37. Van Zyl, G.U.; Liu, T.F.; Claassen, M.; Engelbrecht, S.; de Oliveira, T.; Preiser, W.; Wood, N.T.; Travers, S.; Shafer, R.W. Trends in Genotypic HIV-1 Antiretroviral Resistance between 2006 and 2012 in South African Patients Receiving First- and Second-Line Antiretroviral Treatment Regimens. PLoS ONE 2013, 8, e67188. [Google Scholar] [CrossRef]
  38. Wallis, C.L.; Mellors, J.W.; Venter, W.D.; Sanne, I.; Stevens, W. Protease Inhibitor Resistance Is Uncommon in HIV-1 Subtype C Infected Patients on Failing Second-Line Lopinavir/r-Containing Antiretroviral Therapy in South Africa. AIDS Res. Treat 2011, 2011, 769627. [Google Scholar] [CrossRef]
  39. Sivamalar, S.; Dinesha, T.R.; Gomathi, S.; Pradeep, A.; Boobalan, J.; Solomon, S.S.; Poongulali, S.; Solomon, S.; Balakrishnan, P.; Saravanan, S. Accumulation of HIV-1 Drug Resistance Mutations After First-Line Immunological Failure to Evaluate the Options of Recycling NRTI Drugs in Second-Line Treatment: A Study from South India. AIDS Res. Hum. Retroviruses 2016, 33, 271–274. [Google Scholar] [CrossRef]
  40. Arimide, D.A.; Amogne, M.D.; Kebede, Y.; Balcha, T.T.; Adugna, F.; Ramos, A.; DeVos, J.; Zeh, C.; Agardh, A.; Chih-Wei Chang, J.; et al. High level of HIV drug resistance and virological non-suppression among female sex workers in Ethiopia: A nation-wide cross-sectional study. J. Acquir. Immune Defic. Syndr. 2021, 89, 566–574. [Google Scholar] [CrossRef] [PubMed]
  41. Namakoola, I.; Kasamba, I.; Mayanja, B.N.; Kazooba, P.; Lutaakome, J.; Lyagoba, F.; Kapaata, A.A.; Kaleebu, P.; Munderi, P.; On behalf of the Co, L.s.t. From antiretroviral therapy access to provision of third line regimens: Evidence of HIV Drug resistance mutations to first and second line regimens among Ugandan adults. BMC Res. Notes 2016, 9, 515. [Google Scholar] [CrossRef] [PubMed]
  42. Jordan, M.R.; Hamunime, N.; Bikinesi, L.; Sawadogo, S.; Agolory, S.; Shiningavamwe, A.N.; Negussie, T.; Fisher-Walker, C.L.; Raizes, E.G.; Mutenda, N.; et al. High levels of HIV drug resistance among adults failing second-line antiretroviral therapy in Namibia. Medicine 2020, 99, e21661. [Google Scholar] [CrossRef] [PubMed]
  43. Nzengui-Nzengui, G.F.; Mourembou, G.; M’boyis-Kamdem, H.; Kombila-Koumavor, A.C.; Ndjoyi-Mbiguino, A. HIV protease resistance mutations in patients receiving second-line antiretroviral therapy in Libreville, Gabon. BMC Infect. Dis. 2024, 24, 316. [Google Scholar] [CrossRef]
  44. Obasa, A.E.; Mikasi, S.G.; Brado, D.; Cloete, R.; Singh, K.; Neogi, U.; Jacobs, G.B. Drug Resistance Mutations Against Protease, Reverse Transcriptase and Integrase Inhibitors in People Living With HIV-1 Receiving Boosted Protease Inhibitors in South Africa. Front. Microbiol. 2020, 11, 438. [Google Scholar] [CrossRef]
  45. Makwaga, O.; Adhiambo, M.; Mulama, D.H.; Muoma, J.; Adungo, F.; Wanjiku, H.; Ongaya, A.; Maitha, G.M.; Mwau, M. Prevalence of human immunodeficiency virus-1 drug-resistant mutations among adults on first- and second-line antiretroviral therapy in a resource-limited health facility in Busia County, Kenya. Pan Afr. Med. J. 2020, 37, 311. [Google Scholar] [CrossRef]
  46. Hamers, R.L.; Sigaloff, K.C.E.; Wensing, A.M.; Wallis, C.L.; Kityo, C.; Siwale, M.; Mandaliya, K.; Ive, P.; Botes, M.E.; Wellington, M.; et al. Patterns of HIV-1 Drug Resistance After First-Line Antiretroviral Therapy (ART) Failure in 6 Sub-Saharan African Countries: Implications for Second-Line ART Strategies. Clin. Infect. Dis. 2012, 54, 1660–1669. [Google Scholar] [CrossRef]
  47. Wang, J.; Wang, Z.; Liu, J.; Yue, Y.; Yang, S.; Huang, H.; He, C.; Liao, L.; Xing, H.; Ruan, Y.; et al. Efficacy and HIV drug resistance profile of second-line ART among patients having received long-term first-line regimens in rural China. Sci. Rep. 2015, 5, 14823. [Google Scholar] [CrossRef]
  48. Melikian, G.L.; Rhee, S.-Y.; Varghese, V.; Porter, D.P.; White, K.L.; Taylor, J.; Towner, W.J.; Troia, P.D.; Burack, J.H.; Dejesus, E.; et al. Non-nucleoside reverse transcriptase inhibitor (NNRTI) cross-resistance: Implications for preclinical evaluation of novel NNRTIs and clinical genotypic resistance testing. J. Antimicrob. Chemother. 2014, 69, 12–20. [Google Scholar] [CrossRef]
  49. Ismael, N.; Gemusse, H.; Mahumane, I.; Laurindo, O.; Magul, C.; Baxter, C.; Wilkinson, E.; Hofstra, L.M.; Wagar, N.; Bila, D.; et al. HIV-1 pretreatment and acquired antiretroviral drug resistance before tenofovir/ lamivudine /dolutegravir (TLD) roll-out in Mozambique. BMC Infect. Dis. 2024, 24, 748. [Google Scholar] [CrossRef]
  50. Derache, A.; Wallis, C.L.; Vardhanabhuti, S.; Bartlett, J.; Kumarasamy, N.; Katzenstein, D. Phenotype, Genotype, and Drug Resistance in Subtype C HIV-1 Infection. J. Infect. Dis. 2015, 213, 250–256. [Google Scholar] [CrossRef] [PubMed]
  51. Chimukangara, B.; Lessells, R.J.; Sartorius, B.; Gounder, L.; Manyana, S.; Pillay, M.; Singh, L.; Giandhari, J.; Govender, K.; Samuel, R.; et al. HIV-1 drug resistance in adults and adolescents on protease inhibitor-based antiretroviral therapy in KwaZulu-Natal Province, South Africa. J. Glob. Antimicrob. Resist. 2022, 29, 468–475. [Google Scholar] [CrossRef] [PubMed]
  52. Obasa, A.E.; Ambikan, A.T.; Gupta, S.; Neogi, U.; Jacobs, G.B. Increased acquired protease inhibitor drug resistance mutations in minor HIV-1 quasispecies from infected patients suspected of failing on national second-line therapy in South Africa. BMC Infect. Dis. 2021, 21, 214. [Google Scholar] [CrossRef] [PubMed]
  53. Fox, M.P.; Berhanu, R.; Steegen, K.; Firnhaber, C.; Ive, P.; Spencer, D.; Mashamaite, S.; Sheik, S.; Jonker, I.; Howell, P.; et al. Intensive adherence counselling for HIV-infected individuals failing second-line antiretroviral therapy in Johannesburg, South Africa. Trop. Med. Int. Health 2016, 21, 1131–1137. [Google Scholar] [CrossRef]
  54. Ross, J.; Jiamsakul, A.; Kumarasamy, N.; Azwa, I.; Merati, T.P.; Do, C.D.; Lee, M.P.; Ly, P.S.; Yunihastuti, E.; Nguyen, K.V.; et al. Virological failure and HIV drug resistance among adults living with HIV on second-line antiretroviral therapy in the Asia-Pacific. HIV Med. 2021, 22, 201–211. [Google Scholar] [CrossRef]
  55. Chen, J.; Zhang, M.; Shang, M.; Yang, W.; Wang, Z.; Shang, H. Research on the treatment effects and drug resistances of long-term second-line antiretroviral therapy among HIV-infected patients from Henan Province in China. BMC Infect. Dis. 2018, 18, 571. [Google Scholar] [CrossRef]
  56. Audi, C.; Jahanpour, O.; Antelman, G.; Guay, L.; Rutaihwa, M.; van de Ven, R.; Woelk, G.; Baird, S.J. Facilitators and barriers to antiretroviral therapy adherence among HIV-positive adolescents living in Tanzania. BMC Public Health 2021, 21, 2274. [Google Scholar] [CrossRef]
  57. Zurbachew, Y.; Hiko, D.; Bacha, G.; Merga, H. Adolescent’s and youth’s adherence to antiretroviral therapy for better treatment outcome and its determinants: Multi-center study in public health facilities. AIDS Res. Ther. 2023, 20, 91. [Google Scholar] [CrossRef]
  58. Moosa, A.; Gengiah, T.N.; Lewis, L.; Naidoo, K. Long-term adherence to antiretroviral therapy in a South African adult patient cohort: A retrospective study. BMC Infect. Dis. 2019, 19, 775. [Google Scholar] [CrossRef]
  59. Gumede, S.B.; Venter, F.; de Wit, J.; Wensing, A.; Lalla-Edward, S.T. Antiretroviral therapy uptake and predictors of virological failure in patients with HIV receiving first-line and second-line regimens in Johannesburg, South Africa: A retrospective cohort data analysis. BMJ Open 2022, 12, e054019. [Google Scholar] [CrossRef]
  60. Gidey, K.; Mache, A.; Hailu, B.Y.; Asgedom, S.W.; Tassew, S.G.; Nirayo, Y.L. Second-Line Antiretroviral Treatment Outcomes and Predictors in Tigray Region, Ethiopia. Infect. Drug Resist 2023, 16, 4903–4912. [Google Scholar] [CrossRef]
  61. Amogne, W.; Bontell, I.; Grossmann, S.; Aderaye, G.; Lindquist, L.; Sönnerborg, A.; Neogi, U. Phylogenetic analysis of Ethiopian HIV-1 subtype C near full-length genomes reveals high intrasubtype diversity and a strong geographical cluster. AIDS Res. Hum. Retroviruses 2016, 32, 471–474. [Google Scholar] [CrossRef] [PubMed]
  62. Abebe, A.; Pollakis, G.; Fontanet, A.L.; Fisseha, B.; Tegbaru, B.; Kliphuis, A.; Tesfaye, G.; Negassa, H.; Cornelissen, M.; Goudsmit, J. Identification of a genetic subcluster of HIV type 1 subtype C (C′) widespread in Ethiopia. AIDS Res. Hum. Retroviruses 2000, 16, 1909–1914. [Google Scholar] [CrossRef] [PubMed]
  63. Pollakis, G.; Abebe, A.; Kliphuis, A.; De Wit, T.F.R.; Fisseha, B.; Tegbaru, B.; Tesfaye, G.; Negassa, H.; Mengistu, Y.; Fontanet, A.L. Recombination of HIV type 1C (C′/C″) in Ethiopia: Possible link of EthHIV-1C′ to subtype C sequences from the high-prevalence epidemics in India and Southern Africa. AIDS Res. Hum. Retroviruses 2003, 19, 999–1008. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The proportion of patients with HIV drug resistance mutation.
Figure 1. The proportion of patients with HIV drug resistance mutation.
Viruses 17 00206 g001
Figure 2. The maximum-likelihood phylogenetics tree of HIV-1 subtypes and circulating recombinant forms circulating in Ethiopia. Subtype C reference sequences are depicted in blue, while non-subtype C Ethiopian sequences are shown in pink.
Figure 2. The maximum-likelihood phylogenetics tree of HIV-1 subtypes and circulating recombinant forms circulating in Ethiopia. Subtype C reference sequences are depicted in blue, while non-subtype C Ethiopian sequences are shown in pink.
Viruses 17 00206 g002
Table 1. Demographic characteristics of patients on second-line treatment in Ethiopia (n = 586).
Table 1. Demographic characteristics of patients on second-line treatment in Ethiopia (n = 586).
Characteristics N (n = ) Percent (%)
Age (Year)
  18–24355.97
  25–3416027.30
  35–4424541.81
  ≥4514624.91
Sex
  Male 26745.56
  Female 31954.44
Residence
  Urban 49584.47
  Rural 9115.53
Marital status
  Unmarried 12821.84
  Married, monogamous 23439.93
  Married, polygamous 539.04
  Divorced11619.80
  Widowed 559.39
Educational level
  No formal education 11719.97
  Primary education 25944.20
  Secondary education 13623.21
  College or university 7412.63
Occupational status
  Government employee 10017.06
  Merchant 9616.38
  Driver 244.10
  Housewife10618.09
  Student254.27
  Unemployed 9015.36
  Other 14524.74
Table 2. Clinical characteristics and laboratory parameters of participants with HIV on second-line treatment.
Table 2. Clinical characteristics and laboratory parameters of participants with HIV on second-line treatment.
CharacteristicsN (n = )Percent (%)
Second-line ARV regimensn = 586
  ABC + 3TC + ATV/r11419.45
  ABC + 3TC + LPV/r203.41
  AZT + 3TC + ATV/r10317.58
  AZT + 3TC + LPV/r335.63
  TDF + 3TC + ATV/r25944.20
  TDF + 3TC + LPV/r579.73
Duration on second-line regimen (months)n = 586
  6–1517429.69
  16–4832755.80
  ≥498514.51
First-line ARV regimensn = 586
  ABC + 3TC + EFV111.88
  ABC + 3TC + NVP366.14
  AZT + 3TC + EFV7312.46
  AZT + 3TC + NVP15726.79
  D4T + 3TC + EFV305.12
  D4T + 3TC + NVP8915.19
  TDF + 3TC + EFV10618.09
  TDF + 3TC + NVP8414.33
Duration on first-line regimens (months)n = 586
  0–30101.71
  31–60427.17
  61–9012421.16
  ≥9141069.97
Baseline CD4 at ART initiation (cell/mm3)n = 548
  <20000
  ≥200548100
CD4 at second-line initiation (cell/mm3)n = 541
  <200541100
  ≥20000
Current CD4 (during study (cell/mm3)n = 555
  <20016129.00
  ≥20039471.0
Viral load (RNA copies/mL)n = 586
  <100050586.18
  ≥10008113.82
Abbreviations: antiretroviral (ARV), antiretroviral therapy (ART), abacavir (ABC), lamivudine (3TC), efavirenz (EFV), zidovudine (AZT), atazanavir/ritonavir (ATV/r), lopinavir/ritonavir (LPV/r), nevirapine (NVP), tenofovir (TDF), stavudine (D4T).
Table 3. Drug resistance mutation profiles among patients on second-line regimens with virological failure.
Table 3. Drug resistance mutation profiles among patients on second-line regimens with virological failure.
Mutation TypeN (n = )Percent (%)
NRTIs
  A62V22.86
  D67N/G68.57
  E44ED45.71
  K65R57.14
  K70E/G/R/T68.57
  L210W34.29
  L74I11.43
  M184V1622.86
  K219E57.14
  M41L22.86
  T215Y/V710.00
  T69SADN22.86
  Y115F68.57
  S68G57.14
NNRTIs
  A98G47.55
  E138A23.77
  G190A/E713.21
  H221Y59.43
  K101E/A47.55
  K103N1018.87
  K238T11.89
  L100I35.66
  P225H11.89
  V106I/M35.66
  V108I35.66
  V179T/D23.77
  Y181C713.21
  Y188L11.89
PIs
  I50L211.76
  I54V423.53
  L90M15.88
  M46I423.53
  N88S15.88
  V82A529.41
Table 4. Multivariate analysis of factors associated with virological failure among patients on second-line ART.
Table 4. Multivariate analysis of factors associated with virological failure among patients on second-line ART.
VariablesCOR (95% CI)p-ValueAOR (95% CI)p-Value
Sex
  MaleRef. 0.07
  Female0.53 (0.33–0.85)0.0080.58 (0.32–1.05)
Occupation
  Government employee Ref.
  Merchant 0.95 (0.40–2.30)0.9061.15 (0.44–3.06)0.772
  Driver 1.93 (0.61–6.13)0.2652.10 (0.61–7.20)0.243
  Housewife0.60 (0.23–1.53)0.2851.02 (0.33–3.19)0.973
  Student4.13 (1.50–11.40)0.0066.14 (1.93–19.48) 0.002 *
  Unemployed 1.35 (0.59–3.10)0.4781.67 (0.65–4.31)0.29
  Other 1.31 (0.62–2.80)0.4811.57 (0.67–3.68)0.298
Duration on ART at second-line (months)
  <12Ref.
  >120.53 (0.30–0.93)0.0280.52 (0.28–0.97)0.039 *
Current CD4 (cell/mm3)
  <2002.55 (1.54–4.21)<0.0012.55 (1.52–4.28)<0.001 *
  ≥200Ref.
Abbreviations: AOR, adjusted odds ratios; COR, crude odds ratios; CI, confidence interval. (*) Statistically significant at p ≤ 0.05.
Table 5. HIV drug resistance susceptibility test among the study participants (n = 37).
Table 5. HIV drug resistance susceptibility test among the study participants (n = 37).
Resistance Level (n (%))
Antiretroviral DrugSusceptibleLow Level of
Resistance
Intermediate Level of ResistanceHigh Level of
Resistance
NRTIs
  Abacavir (ABC) 21 (56.78%) 2 (5.41%)2 (5.41%)12 (32.43%)
  Zidovudine (AZT)29 (78.38%)1 (2.70%)07 (18.92%)
  Stavudine (D4T)22 (59.46%)3 (8.11%)3 (8.11%)9 (24.32%)
  Didanosine (DDI)22 (59.46%)2 (5.41%)3 (8.11%)10 (27.03%)
  Emtricitabine (FTC) 21 (56.76%)0016 (43.24%)
  Lamivudine (3TC) 21 (56.76%)0016 (43.24%)
  Tenofovir (TDF) 23 (62.16%)2 (5.41%)6 (16.22%)6 (16.22%)
NNRTIs
  Doravirine (DOR)20 (54.05%)6 (16.12%)8 (21.62%)3 (8.11%)
  Efavirenz (EFV) 13 (35.14%)2 (5.41%)4 (10.81%)18 (48.65%)
  Etravirine (ETR)20 (54.05%)6 (16.21%)8 (21.62%)3 (8.11%)
  Nevirapine (NVP) 13 (35.14%)02 (5.41%)22 (59.46%)
  Rilpivirine (RPV)19 (51.35%)2 (5.41%)2 (5.41%)14 (37.84%)
PIs
  Atazanavir (ATV)30 (81.08%)1 (2.70%)06 (16.22%)
  Darunavir (DRV)37 (100%)000
  Fosamprenavir (FPV) 32 (86.49%)1 (2.70%)2 (5.41%)2 (5.41%)
  Indinavir (IDV) 31 (83.78%)11 (2.70%)4 (10.81%)
  Lopinavir (LPV) 31 (83.78%)1 (2.70%)05 (13.51%)
  Nelfinavir (NFV) 31 (83.78%)01 (2.70%)5 (13.51%)
  Saquinavir (SQV) 31 (83.78%)1 (2.70%)1 (2.70%)4 (10.81%)
  Tipranavir (TPV)32 (86.49%)3 (8.11%)2 (5.41%)0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zealiyas, K.; Gebreegziabxier, A.; Getaneh, Y.; Kidane, E.; Woldesemayat, B.; Yizengaw, A.; Gutema, G.; Adane, S.; Yimer, M.; Yilma, A.; et al. Viral Suppression and HIV Drug Resistance Among Patients on Second-Line Antiretroviral Therapy in Selected Health Facility in Ethiopia. Viruses 2025, 17, 206. https://doi.org/10.3390/v17020206

AMA Style

Zealiyas K, Gebreegziabxier A, Getaneh Y, Kidane E, Woldesemayat B, Yizengaw A, Gutema G, Adane S, Yimer M, Yilma A, et al. Viral Suppression and HIV Drug Resistance Among Patients on Second-Line Antiretroviral Therapy in Selected Health Facility in Ethiopia. Viruses. 2025; 17(2):206. https://doi.org/10.3390/v17020206

Chicago/Turabian Style

Zealiyas, Kidist, Atsbeha Gebreegziabxier, Yimam Getaneh, Eleni Kidane, Belete Woldesemayat, Ajanaw Yizengaw, Gadisa Gutema, Sisay Adane, Mengistu Yimer, Amelework Yilma, and et al. 2025. "Viral Suppression and HIV Drug Resistance Among Patients on Second-Line Antiretroviral Therapy in Selected Health Facility in Ethiopia" Viruses 17, no. 2: 206. https://doi.org/10.3390/v17020206

APA Style

Zealiyas, K., Gebreegziabxier, A., Getaneh, Y., Kidane, E., Woldesemayat, B., Yizengaw, A., Gutema, G., Adane, S., Yimer, M., Yilma, A., Tadele, S., Sasinovich, S., Medstrand, P., & Arimide, D. A. (2025). Viral Suppression and HIV Drug Resistance Among Patients on Second-Line Antiretroviral Therapy in Selected Health Facility in Ethiopia. Viruses, 17(2), 206. https://doi.org/10.3390/v17020206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop