Quantifying the impact of cascade inequalities: a modelling study on the prevention impacts of antiretroviral therapy scale-up in Eswatini

1 Introduction

Early HIV treatment via antiretroviral therapy (ART) is a lifesaving intervention increasing both the quantity and quality of life [1]. A secondary benefit of early ART given Undetectable = Untransmittable (U=U) is that transmission risks are mitigated in serodifferent partnerships [2]. To realize these benefits, massive efforts are underway to achieve the UNAIDS 95-95-95 ART cascade targets [3] — i.e., to have: 95% diagnosed among people living with HIV, 95% on ART among those diagnosed, and 95% virally suppressed among those on ART. Botswana, Eswatini, Rwanda, Tanzania, and Zimbabwe have already surpassed 95-95-95 nationally [3], and and achieving these targets is expected to help reduce HIV incidence towards elimination.

Yet, there are growing concerns that inequalities in the ART cascade could undermine the population-level prevention impacts of ART anticipated from individual-level and model-based studies [4–6]. Specifically, available data suggest that cascade attainment is often lower among subpopulations at greater risk of HIV acquisition and/or transmission, including key populations, younger men and women, and highly mobile populations [5,7]. These inequalities can be driven by systemic barriers to engagement in care faced by marginalized populations, which intersect with individual, network, and structural determinants of HIV risk, such as economic insecurity, mobility, stigma, discrimination, and criminalization [8–11]. Moreover, cascade data may be lacking entirely for subpopulations experiencing the greatest barriers to care — i.e., the lowest ART cascades likely remain unmeasured [7].

Numerous transmission modelling studies have sought to estimate the prevention impacts of achieving 90-90-90+ across Sub-Saharan Africa [12,13]. Modelled populations are often stratified by risk, including key populations like female sex workers (FSW) and their clients, to capture important epidemic dynamics related to risk heterogeneity [14]. However, these studies have generally assumed that ART cascade attainment (i.e., proportions diagnosed, treated, and virally suppressed) or progression (i.e., rates of diagnosis, treatment initiation, and viral suppression) were equal across modelled subpopulations. For example, among the studies in [13] (see also Box 1), key populations were usually assumed to have equal cascade progression with the population overall, or greater in some scenarios, but never lesser. Thus, the potential influence of intersecting risk and cascade inequalities on ART prevention impacts has not been explored.

We therefore examined the following questions in an illustrative modelling analysis:

1. How are estimates of population-level ART prevention impacts influenced by inequalities in ART cascade across subpopulations?

2. Under which epidemic conditions do such inequalities have the largest influence?

We examined these questions using a deterministic compartmental model of heterosexual HIV transmission in Eswatini, focusing on inequalities related to sex work. Eswatini has the highest national HIV prevalence in the world, and recently surpassed 95-95-95 with minimal inequalities across subpopulations [3,15]. As such, we used observed ART cascade scale-up in Eswatini as a base case reflecting exemplary and evidently attainable scale-up, and explored counterfactual scenarios in which scale-up was slower, and where specific subpopulations could have been left behind.

2 Methods

We constructed a deterministic compartmental model of heterosexual HIV transmission, stratified by subpopulations defined by sex and sexual activity, and health states reflecting HIV natural history and ART cascade of care. The model includes eight subpopulations, including FSW at higher vs lower risk, and likewise for clients of FSW, and four partnership types, including main/spousal, casual, repeat sex work, and one-off sex work (Figure 1). We calibrated the model to reflect the HIV epidemic and ART scale-up in Eswatini (base case). We then explored counterfactual scenarios in which ART cascade was reduced among various combinations of subpopulations, and quantified ART prevention impacts by comparing base case and counterfactual scenarios.

• Download figure
• Open in new tab

Figure 1:

Model structure and transitions

Low: lowest activity; Med: medium activity; LR/HR: lower/higher risk; FSW: female sex workers; Clients: of FSW; CD4: CD4+ T-cell count per mm³; ART: antiretroviral therapy; rates — λ: force of infection; η: HIV progression; δ: diagnosis; τ: ART initiation; σ: viral suppression; σ^′: viral re-suppression; ς: ART failure / discontinuation; not shown: turnover amongst activity groups in (a).

3 Results

We first summarize modelled patterns of HIV transmission in the base case (1000 model fits), calibrated to reflect the Eswatini epidemic up to 2020. Our model suggests that transmission within repeat sex work partnerships was a dominant driver of early epidemic growth (Figures B.14 and B.15). However, from approximately 1994 onward, the majority of new yearly infections were transmitted within casual partnerships, including 50% (median) of infections in 2020 in the base case.³ By 2020, the majority of cumulative infections acquired were among women with lowest sexual activity (Figure B.12b), while the majority of cumulative infections transmitted were from clients (Figure B.12a). Overall HIV prevalence in 2020 was median (95% CI): 23.8 (22.4, 24.7) % (Figure B.4a), and overall incidence was 6.6 (5.3, 7.6) per 1000 person-years (Figure B.4c). The prevalence ratio between FSW and women overall was 1.78 (1.70, 1.87), and between clients and men overall it was 1.92 (1.49, 2.49) (Figure B.4b). Due to turnover and higher HIV incidence among FSW, achieving similar rates of diagnosis among FSW vs other women (Figure B.6a) required approximately twice the rate of testing. Sex work contributed a growing proportion of infections over 2020–2040: from 32% to 42% (Figure B.14).

3.1 Objective 1: Influence of cascade differences between subpopulations

Figure B.17 illustrates ART cascade attainment over time in each of the four counterfactual scenarios (80-80-90 overall by 2020), and the base case (95-95-95 overall by 2020). Figure B.16 illustrates overall HIV incidence in each scenario. Figure 2 then illustrates cumulative additional infections (CAI) and additional incidence rate (AIR) in each counterfactual scenario vs the base case. Leaving behind neither FSW and nor clients led to the fewest CAI: median (95% CI) 14.9 (10.4, 18.4) % more than the base case by 2030. By contrast, leaving behind both FSW and clients led to the most CAI: 26.3 (19.7, 33.0) % more than the base case by 2030 — a 73 (40, 149) % increase. Leaving behind either FSW or clients resulted in 19.7 (15.4, 23.6) % or 22.7 (17.6, 26.7) % cumulative additional infections vs the base case, respectively. Results were similar for additional incidence rate (AIR). In all counterfactual scenarios, the majority of additional infections were transmitted via casual partnerships (Figure B.18a), and acquired among non-FSW women (Figure B.18c).

• Download figure
• Open in new tab

Figure 2:

Relative additional infections under counterfactual scenarios vs the base case

Base case: 95-95-95 by 2020; counterfactual scenarios: 80-80-90 overall by 2020, with reduced cascade (60-40-80: left behind) among FSW, clients of FSW, both, or neither; whiskers, boxes, and midlines: 95% CI, 50% CI, median of model fits.

Patterns of onward transmission were also similar across scenarios (Figure B.18b), though subpopulation contributions increased if they were left behind.

3.2 Objective 2: Conditions that maximize the influence of cascade differences

The fitted regression models Eq. (2) indicated that population-overall viral non-suppression (D) and relative non-suppression among FSW and clients (d_i) each had strong and positive effects on 2030 CAI and AIR outcomes (p < 10⁻⁵). These effects corroborate the results of Objective 1. Regression residuals are shown in Figure B.21, while Figure 3 plots the estimated effects of subpopulation-specific non-suppression d_i, plus effect modification by epidemic conditions C_j. The effect of relative non-suppression on CAI and AIR among both FSW and clients generally increased with: larger FSW and client population sizes, and larger HIV incidence ratio among clients vs other men. Effects of relative non-suppression among clients further increased with duration buying sex (slower turnover), while effects of relative non-suppression among both groups on CAI specifically further increased with HIV incidence ratio among FSW vs other women.

• Download figure
• Open in new tab

Figure 3:

Estimated effects on relative additional infections of disproportionate viral non-suppression (d) among FSW and clients vs population overall, plus effect modification by epidemic conditions

(a) cumulative additional infections, (b) additional incidence rate by 2030 vs base case; FSW: female sex workers; Clients: of FSW; IR: incidence ratio in 2020; di: difference in subpopulation-i-specific viral non-suppression vs population overall within counterfactual scenario; points and error bars: mean and 95% CI for each effect estimated via Eq. (2).

4 Discussion

We sought to explore how inequalities in the ART cascade that intersect with HIV risk heterogeneity may influence the model-estimated prevention impacts of ART. In our applied analysis of Eswatini, we found that slower cascade scale-up that left behind female sex workers (FSW) and their clients could have resulted in 40–149% more infections by 2030 vs slower scale-up alone. We also found that the impact of leaving behind FSW and/or clients was largely determined by their population sizes and HIV incidence relative to the wider population.

Eswatini has recently surpassed 95-95-95 for the population overall [15]. These targets were achieved through numerous initiatives coordinated across sectors, including those led by the MaxART program [31,32]. Multiple stakeholders, including people living with HIV, healthcare providers, traditional and religious leaders, community groups, and researchers were engaged via multiple channels, such as Technical Working Groups, Community Advisory Boards, and specific meetings for prioritized groups (men and adolescents) [31,32]. Drawing on this engagement and social science research to understand barriers to care, cascade services were comprehensively strengthened via investments in training, infrastructure, anti-stigma communication, demand creation, and monitoring [31,32].

Among FSW living with HIV in Eswatini, data suggest that 88% were diagnosed and 86% were on ART in 2021 (i.e., ART coverage was 98% among those diagnosed) [27]. Data on viral suppression among FSW were lacking, and so assumed to be similar to other women in the simulated epidemic (see § A.3.8 for more details). Although lower than 95-95-95, this cascade among FSW living with HIV is higher than in many other regions [7,10]. Strong programs are required to attain high cascades among FSW, considering that women enter and exit sex work (turnover) and likely experience highest risk of HIV acquisition during sex work. That is, programs must ensure higher rates of HIV testing, ART initiation, and retention among FSW vs other women to achieve similar cascades. For example, we inferred that rates of HIV testing in 2016 was 80–227% higher among FSW vs other women to reproduce observed cascade data during model calibration.

In Eswatini, programs for key populations include safe access to tailored services via drop-in centers (locally known as TRUE), mobile outreach, venue-based, and one-on-one options [27]. Health and clinical services are also integrated with efforts to reduce structural vulnerabilities, including experiences of harassment, violence, and fear of seeking healthcare, through community empowerment, psycho-social and legal supports, and sensitization and training for police and healthcare workers [27]. These programs have been designed and refined with ongoing community leadership and engagement, allowing them to better meet the specific needs of key populations, for whom barriers to engagement in HIV care often intersect with drivers of HIV risk, including economic insecurity, mobility, stigma, discrimination, and criminalization [4,8–11,33]. Our data-informed modeling of cascade scale-up in Eswatini confirms that such an equity-focused approach to ART cascade scale-up can maximize prevention impacts, and accelerate overall reductions in HIV incidence.

Our study highlights the importance of reaching both FSW and their clients, echoing recent modelling studies of South Africa and Cameroon [34,35]. These studies found that gaps in HIV prevention and treatment for clients were among the largest contributors to onward transmission in recent years. Such findings reiterate the need for improved data on both FSW and clients, including estimates of population size, sexual behaviour, and ART cascade attainment. These estimates may be difficult to obtain because individuals are unlikely to report buying or selling sex in population-level surveys due to stigma and criminalization [36] and because many clients are highly mobile (including transient seasonal/occupational migration) [33]. Thus, innovative study designs, bias adjustments, and services may be needed to understand and meet clients’ needs.

While numerous modelling studies have examined the potential prevention impacts of ART cascade scale-up [13] (Appendix B.3), our study is the first to explore the impact of inequalities in ART cascade across subpopulations with consistent population overall cascade across scenarios. Similar work by Marukutira et al. [37] illustrated the limited impact of achieving 95-95-95 for only citizens and not immigrants in Botswana, while Maheu-Giroux et al. [38] illustrated the high cost-effectiveness of prioritizing key populations (including clients) for ART in Côte d’Ivoire. Indeed, our findings are likely generalizable to other epidemic contexts. HIV prevalence ratios between key populations and the population overall are relatively low in Eswatini vs elsewhere [39,40]; thus, the impact of cascade inequalities among key populations in other contexts would likely be even greater than we found for Eswatini. Moreover, as HIV incidence declines in many settings, transmissions may become concentrated among key populations [41,42], further magnifying the impact of cascade inequalities.

A primary strength of our analysis is the use of observed ART cascade scale-up to 95-95-95 in Eswatini as the base case, with plausible cascade inequalities explored in counterfactual scenarios. As noted above, the available data suggest that Eswatini has minimized cascade inequalities which persist elsewhere [7]. Thus, our counterfactual scenarios directly estimate the consequences of failing to address these inequalities. Second, drawing on our conceptual framework for risk heterogeneity [13, Table 1] and multiple sources of context-specific data [20–24,27], we captured several dimensions of risk heterogeneity, including: heterosexual anal sex, four types of sexual partnerships, sub-stratification of FSW and clients into higher and lower risk strata, and subpopulation turnover (Table 1, Appendix A). Accurate modelling of risk heterogeneity has been shown to mediate model-estimated ART prevention impacts [43], and is especially important when considering differential ART scale-up across subpopulations. Finally, our analytic approach to Objective 2, in which epidemic conditions are conceptualized as potential effect modifiers represents a unique methodological contribution to the HIV modelling literature.

Our study also has limitations. First, we did not model pre-exposure prophylaxis (PrEP). However, our analyses focus on the time period prior to widespread PrEP availability in Eswatini. Second, we did not consider transmitted drug resistance (TDR). However, drug resistance is more likely to emerge in the context of barriers to viral suppression [44]; thus, lower cascade among those at higher risk would likely accelerate emergence of transmitted drug resistance, and thereby magnify our findings. Finally, our model structure did not include age, and we only considered heterosexual HIV transmission in Eswatini. Future work can explore adaptation of the model to consider PrEP, TDR, age stratification, and additional modes of HIV transmission. While the magnitude of our results may change with such adaptations, we do not expect that the qualitative interpretation would change. In fact, our findings would likely generalize to other transmission networks and determinants of risk heterogeneity, including other key populations and subpopulations such as highly mobile populations and young women [33,45].

In conclusion, the HIV response must remain rooted in context-specific understandings of inequalities in HIV risk and in access to HIV services, which often stem from common upstream factors. Thus, differences in ART cascade within and between subpopulations at higher risk of HIV must be monitored, characterized, and addressed to fully realize the anticipated benefits of ART at both the individual and population levels.