The epidemiology of HIV population viral load in twelve sub-Saharan African countries

Abstract

Background

We examined the epidemiology and transmission potential of HIV population viral load (VL) in 12 sub-Saharan African countries.

Methods

We analyzed data from Population-based HIV Impact Assessments (PHIAs), large national household-based surveys conducted between 2015 and 2019 in Cameroon, Cote d’Ivoire, Eswatini, Kenya, Lesotho, Malawi, Namibia, Rwanda, Tanzania, Uganda, Zambia, and Zimbabwe. Blood-based biomarkers included HIV serology, recency of HIV infection, and VL. We estimated the number of people living with HIV (PLHIV) with suppressed viral load (<1,000 HIV-1 RNA copies/mL) and with unsuppressed viral load (viremic), the prevalence of unsuppressed HIV (population viremia), sex-specific HIV transmission ratios (number female incident HIV-1 infections/number unsuppressed male PLHIV per 100 persons-years [PY] and vice versa) and examined correlations between a variety of VL metrics and incident HIV. Country sample sizes ranged from 10,016 (Eswatini) to 30,637 (Rwanda); estimates were weighted and restricted to participants 15 years and older.

Results

The proportion of female PLHIV with viral suppression was higher than that among males in all countries, however, the number of unsuppressed females outnumbered that of unsuppressed males in all countries due to higher overall female HIV prevalence, with ratios ranging from 1.08 to 2.10 (median: 1.43). The spatial distribution of HIV seroprevalence, viremia prevalence, and number of unsuppressed adults often differed substantially within the same countries. The 1% and 5% of PLHIV with the highest VL on average accounted for 34% and 66%, respectively, of countries’ total VL. HIV transmission ratios varied widely across countries and were higher for male-to-female (range: 2.3–28.3/100 PY) than for female-to-male transmission (range: 1.5–10.6/100 PY). In all countries mean log₁₀ VL among unsuppressed males was higher than that among females. Correlations between VL measures and incident HIV varied, were weaker for VL metrics among females compared to males and were strongest for the number of unsuppressed PLHIV per 100 HIV-negative adults (R² = 0.92).

Conclusions

Despite higher proportions of viral suppression, female unsuppressed PLHIV outnumbered males in all countries examined. Unsuppressed male PLHIV have consistently higher VL and a higher risk of transmitting HIV than females. Just 5% of PLHIV account for almost two-thirds of countries’ total VL. Population-level VL metrics help monitor the epidemic and highlight key programmatic gaps in these African countries.

Introduction

Globally, an estimated 38 million people were living (PLHIV) with HIV in 2021 [1]. Sub-Saharan Africa (SSA), and in particular Eastern and Southern Africa, is the epicenter where most HIV epidemics are generalized and severe [1]. The massive scale-up of antiretroviral therapy (ART) beginning in the early-to-mid-2000s in SSA led to a dramatic reduction in HIV-related mortality and a slower decline in adult HIV incidence [1]. In the absence of a vaccine or cure, Treatment as Prevention, a concept informed by clinical trials demonstrated that ART-mediated viral load suppression (VLS) minimizes the risk of onward sexual transmission of HIV [2]. Using modeled findings, UNAIDS suggests that by meeting first the global 90-90-90 targets by 2020 and by 2025 the global 95-95-95 targets (95% of all PLHIV are aware of their HIV-positive status, among PLHIV aware of their HIV positive status, 95% are on ART; and, among those on ART, 95% have VLS), will help “end AIDS” by 2030 [3]. VLS reduces the HIV-related mortality and sexual transmission of HIV. These expected dual population-level effects of ART on mortality and HIV incidence make monitoring viral load (VL) essential.

However, the utility of VL metrics go beyond monitoring progress toward the third 95 target; notably VL also lends itself to inform a population’s (remaining) HIV transmission potential. Population viremia, the population prevalence of unsuppressed HIV often is strongly correlated with higher HIV incidence in the larger population [4–6]. Unfortunately, VL data that consider only PLHIV in the denominator as well as clinic-based VL metrics that exclude undiagnosed or untreated PLHIV often show statistically weak or insignificant relationships with HIV incidence in the background population [4, 7]. Further, while dichotomizing VL into suppressed and unsuppressed has utility for programming and policy, the categorization also leads to information loss as the level of VL is strongly correlated with HIV disease progression and transmission probability. Numerous VL metrics are described in the literature, and efforts to use such metrics are increasing. Examples include the mean and total “community viral load” as “markers of access to care and treatment” and indicators of a population’s viral burden, [6, 8, 9] or, borrowing a graphical concept from economics, the Lorenz curve (and associated Gini coefficient) help to display the often highly unequal distribution of VL in a given population as exemplified in a San Francisco patient cohort [10]. Above all, VL-related metrics allow describing the HIV transmission risk from PLHIV to uninfected individuals [11, 12].

As population-level VL data appear key to understanding progress towards effective HIV control, such data are becoming increasingly common in SSA countries, e.g., in Uganda and Kenya [4, 13, 14]. We examined data from 12 national household-based surveys in SSA with the goal to describe the epidemiology of VL and to examine the utility of various VL metrics for surveillance and epidemic monitoring at the national population level.

Materials and methods

Data sources, setting and survey design

PEPFAR, the President’s Emergency Plan for AIDS Relief, a U.S. Government program assisting HIV response efforts in selected low- and middle-income countries, funded a series of national household-based surveys in 2015–2019 to assess the impact of HIV programs in SSA countries. Implementing partners included country governments, ICAP at Columbia University in New York, and the Centers for Disease Control and Prevention (CDC). The field methods of these Population-based HIV Impact Assessments (PHIAs) are described in detail elsewhere [15, 16]. Briefly, PHIAs were powered for national adult HIV incidence estimation as well as sub-national VLS. A two-stage cluster sampling design was used to select clusters and households. The PHIA data sets are in the public domain and available upon request [17].

Study population

We analyzed PHIA data from Cameroon (2017–18, number of adult respondents with an unambiguous HIV result in the survey: 26,031), Cote d’Ivoire (2017–18, n = 17,813), Eswatini (2016–17, n = 10,016), Kenya (2018–19, n = 27,745), Lesotho (2016–17, n = 11,682), Malawi (2015–16, n = 17,187), Namibia (2017, n = 16,939), Rwanda (2018–19, n = 30,637), Tanzania (2016–17, n = 29,369), Uganda (2016–17, n = 29,024), Zambia (2016, n = 19,115), and Zimbabwe (2015–16, n = 20,577). The PHIAs sampled PLHIV irrespective of their reported serostatus knowledge, linkage to care, or retention in care. Among participants enrolled in the surveys the proportion who consented to a blood draw ranged from 87.5% to 99.7% (median: 92.5%). Our study considered male and female survey participants aged 15–64 years (except Lesotho and Zambia: 15–59 years) for whom an HIV test was conducted during the survey. Additional survey eligibility criteria included having spent the last night in the sampled household and speaking one of the offered interview languages. We excluded children from the analysis due to their much-diminished HIV transmission potential and because the number of sampled HIV-infected children was small in each country. We excluded the Ethiopia PHIA as sampling was restricted to urban areas and so did not allow for the examination of national VL metrics or the relationship between viremia and HIV incidence.

Data collection

Interview data were collected through structured face-to-face interviews at the households using mobile tablets. Original data variables included household characteristics, demographics (including their sex at birth, hence the terms male and female), HIV-related service uptake, and HIV-related risk behaviors. Participants were asked about their HIV testing history, past test results, and ART status.

Biologic measures

Detailed descriptions of our biomarker procedures have been published elsewhere [18]. Venous blood collected with Ethylenediaminetetraacetic acid (EDTA) as an anticoagulant was tested at the household for HIV using the countries’ national rapid test algorithm. All HIV-positive specimens underwent confirmatory testing using the Geenius HIV 1/2 Supplemental Assay (Bio-Rad, Hercules, CA, USA) or PCR testing. Additional HIV quality control testing varied by country. To determine VL, HIV-positive plasma specimens (or, if unavailable, dried blood spots (DBS), accounting for 3% of all participants) were further tested in-country using a Roche or Abbott platform [19]. The lower limit of detection was 20–40 copies/mL for plasma and 700–800 copies/mL for DBS. HIV-positive participants with a VL<1,000 HIV-1 RNA copies/mL were classified as virally suppressed; those with a VL≥1,000 copies/mL were classified as unsuppressed or viremic.

HIV-positive participants (except in Kenya and Rwanda) received a CD4⁺ T-cell count measurement in the field using the Pima™ CD4 Analyzer (Abbott Molecular Inc., Chicago, IL, USA, formerly Alere).

DBS from HIV-positive participants were qualitatively tested for the presence of three antiretroviral (ARV) drugs (two from the countries’ prevailing 1^st line ARV regimens, plus one from the prevailing 2^nd line regimens) using high performance liquid chromatography coupled with tandem mass spectrometry at the University of Cape Town, South Africa.

For HIV-1 recency testing, HIV-seropositive specimens were tested in-country with the limiting antigen (LAg) avidity enzyme immunoassay [20]. Specimens with a normalized optical density (ODn) value of ≤2.0 were re-tested for HIV recency in triplicate; median ODn values of ≤1.5 were classified as LAg-recent. For DBS specimens, the Maxim HIV-1 LAg DBS EIA was used; median values of ≤1.0 were classified as LAg-recent. Participants who tested LAg-recent, had a VL≥1,000 copies/mL, and tested ARV-negative were classified as having recent HIV-1-infection. Participants with a LAg ODn of >1.5 or a VL<1,000 copies/mL, or had detectable ARV were classified as having long-term HIV-1 infection. More detail is published elsewhere [21].

Data management and analysis

The data presented for each survey are based on participants with a valid biomarker result for HIV-1 status (HIV-2 status was not considered for the definition for HIV status and prevalence). The analytic weights account for the selection probabilities of primary sampling units, households, and individuals and incorporate adjustments for non-response by households, individuals, and biomarker participation. The weights were further adjusted to match each country’s age and sex distribution from a recent population projection. Estimates across all countries combined were also weighted by population size. We classified HIV-infected participants as diagnosed or aware of their HIV status if they reported a prior HIV positive test result or if they had detectable ARV. HIV-positive participants who tested ARV-positive were classified as being on ART.

For each country, we estimated HIV incidence based on the number of HIV-recent test results, mean duration of recency (130 days for most countries), time cutoff of one year, and a proportion false recent of 0.00 for the LAg avidity assay. Other computed indicators included the arithmetic mean VL, the number of suppressed and unsuppressed PLHIV by sex and age, countries’ total VL (weighted national sum-total of all PLHIV’s VL in a country, using 1 c/mL for undetectable VL). We computed VLS for all PLHIV (population VLS), diagnosed PLHIV, and PLHIV on ART. We computed the proportion virally suppressed among PLHIV who tested ARV-negative as well as the proportion unsuppressed among PLHIV who tested ARV-positive. We estimated population viremia by multiplying HIV prevalence by the proportion PLHIV who were unsuppressed (viremic).

We mapped HIV and viremia prevalence using geomasked enumeration area centroids [22] and spatial interpolation methods to explore the spatial dimensions of viremia. Specifically, we used kernel density smoothing within PrevR [23, 24]. A surface was computed using a Gaussian kernel which depends on each enumeration area’s number of HIV (or unsuppressed) cases and all tested cases. The “intensity” surfaces (number of persons per area) were calculated separately for each (weighted number of HIV-infected, unsuppressed [numerators] and tested individuals [denominator]) and a ratio was calculated from these to produce smooth surfaces of viremia and HIV prevalence. The results of this modeling were reworked in QGIS 3.4 [25] and ArcMap 10.5.1 (Environmental Systems Research Institute, Redlands, CA, USA) and presented with country-specific scales to each country’s minimum and maximum values using percent stretch with 0.5% clipping; this country-specific strategy highlights local differences in spatial pattern in HIV and viremia prevalence [26, 27]. To visualize the spatial pattern of PLHIV and unsuppressed PLHIV we applied raster math of each prevalence surface and the 2020 spatial population models (unconstrained, UN adjusted, 100m resolution) for each country from the WorldPop [28] program. Minor adjustments were made to rescale the resultant product totals to PHIA weighted estimates of adult population, PLHIV, and number of unsuppressed persons at the national level (15–59 or 15–64 years [see Table 1]). Adjustments are presented in S1 Table. The maps displaying the number of unsuppressed persons per pixel are produced using a “histogram equalize” stretch to increase contrast for visualization.

[Figure omitted. See PDF.]

Population-based HIV Impact Assessments, 2015–19.

We constructed Lorenz curves [29] and computed Gini coefficients [30] for VL by country using SAS v9.4 (SAS Institute, Cary, NC, USA) to construct population-weighted cumulative percentiles of total VL that correspond to percentiles of the population. In ART-naïve HIV-infected populations one would expect a scattered distribution of mostly unsuppressed VL values, however, in the presence of widespread treatment one can expect a growing disparity ranging from undetectable to high VL values. Undetectable VL values were assigned a VL value of 1 for the construction of the Lorenz curves. Gini coefficients are a measure of statistical dispersion, ranging from 0 (perfect parity) to 1 (perfect disparity), whereas Lorenz curves graphically depict disparity. For the computation of Gini coefficients and construction of Lorenz curves, we included all HIV-infected participants and ordered the observations by the weighted number of VL copies.

For each country we estimated the annualized HIV transmission ratios using the following formula: (estimated number of annual new HIV-1 infections/number PLHIV)* 100, thus deriving estimates of the number of HIV transmissions per 100 PLHIV per year. Assuming that all transmissions stem from unsuppressed PLHIV we further estimated effective HIV transmission ratios using the same formula but restricting the number of PLHIV to unsuppressed only. Assuming that all transmission was between opposite sexes we computed sex-specific effective HIV transmission ratios (male-to-female and female-to-male). Using the same data, we also estimated the number of unsuppressed adults per incident HIV infection for each country by dividing the number of unsuppressed adults by the number of incident HIV infections, for both sexes combined as well as for each sex (No. unsuppressed males / No. female incident HIV; No. unsuppressed females / No. male incident HIV), thus deriving a measure of how many unsuppressed adults/males/females on average lead to one incident case of adult/female/male HIV per year.

We examined the relationship between incident HIV and various metrics across the 12 countries by computing the coefficient of determination (R²). The examined metrics included total VL (computed as a country’s estimated sum of all PLHIV’s VL, absolute and log₁₀ transformed), population viremia, mean VL (per capita VL, i.e., HIV-infected, and uninfected adults combined), Gini coefficients, VLS (for all PLHIV, for diagnosed PLHIV, and for PLHIV on ART) and, for comparison purposes, HIV seroprevalence. For all metrics we computed R² values for both sexes combined; for select metrics we also computed R² values for opposite sex correlations to examine male-to-female and female-to-male transmission potentials. For this sub-analysis we re-computed Gini coefficients by including all participants, both HIV-uninfected and infected, assigning HIV-uninfected adults a log₁₀ value of 0 (i.e., an absolute VL value of 1 copy/mL).

Ethics statement

The survey protocols were approved by appropriate institutional review boards in each country, at Columbia University, Westat, UCSF (for Namibia), and at CDC. Institutional review board names, protocol numbers, and approval dates for each survey are shown as Supporting Information. Written informed consent was obtained from adult participants; assent and parental permission was obtained for participants under the age of consent (typically 15–17 years). HIV tests were returned to participants (In Kenya, participants could also opt to receive their HIV results at a health facility of their choice.); VL results were communicated to the health facility chosen by the participants. Efforts were made to link HIV-infected participants not on ART to a health care provider of their choice.

Results

Descriptive viral load metrics

Table 1 shows selected HIV and VL-related metrics; corresponding measures of uncertainty are presented in S2 Table. The estimated population viremia (population prevalence of unsuppressed HIV) ranged from 0.7% (Rwanda) to 8.3% (Lesotho). Tanzania was estimated to have the largest total VL (101.8 billion copies/mL), whereas Rwanda’s total VL offered the smallest estimate at 5.3 billion copies/mL. Log₁₀ transformed mean VL among unsuppressed PLHIV ranged from 4.36 (Malawi) to 4.88 (Cameroon), corresponding to an absolute VL range of 22,900 to 75,900 copies/mL. In each country examined, the mean log₁₀ VL among unsuppressed males (range across countries: 4.47–5.15) exceeded that among unsuppressed females (range across countries: 4.25–4.75); however, 95% CI for some countries overlapped (S2 Table). The highest mean log₁₀ VL for both males and females were observed in Cameroon.

Fig 1 displays the frequency distribution of categorical VL among ARV-negative PLHIV by sex, combined for all 12 countries. Females more frequently showed a VL category below 50,000 copies/mL (peak frequency 1,000<10,000 c/ml) whereas males more often showed a VL category above 50,000 copies/mL (peak frequency 100,000<250,000 c/ml). Fig 2 displays the log₁₀ VL by CD4⁺ T cell count category and sex; male ARV-negative PLHIV tended to exhibit higher VL than females. Fig 3 shows the country-specific population pyramids for PLHIV by age band and viral suppression status, illustrating substantial variation in the distribution of viremia by sex and age. In most countries the estimated number of both male and female unsuppressed adults increased from age 15–19 years onwards, often peaking around age 30–35 years and declining again thereafter. In all countries, female unsuppressed PLHIV outnumbered male unsuppressed PLHIV (Table 1, median: 59%). The ratio of female:male unsuppressed adults varied from 1.08 (Malawi) to 2.10 (Cameroon).

[Figure omitted. See PDF.]

12 Population-based HIV Impact Assessments, 2015–19. ARV: antiretroviral; PLHIV: people living with HIV; ND: Not detectable; K denotes values in thousands, M in millions.

[Figure omitted. See PDF.]

10 Population-based HIV Impact Assessments, 2015–19. Bars represent the 25th and 75th percentiles. ARV: antiretroviral; PLHIV: people living with HIV. Across all CD4 strata, the median viral load among males was 4.766; among females, 4.481.

[Figure omitted. See PDF.]

Population-based HIV Impact Assessments, 2015–2019. PLHIV: (adult) people living with HIV. Dark blue color: Males unsuppressed, Light blue color: males suppressed, Dark orange color: females unsuppressed, Light orange color: females suppressed.

Fig 4 displays the spatial distribution of estimated HIV seroprevalence, viremia prevalence, and the number of unsuppressed PLHIV. Differences across maps should be appreciated within the same country (rather than across countries) and are a function of the proportion of PLHIV suppressed (affecting viremia prevalence) as well as population density (also affecting the number unsuppressed PLHIV).

[Figure omitted. See PDF.]

Viremia prevalence indicates the estimated percentage of the adult population (HIV-infected and uninfected) with unsuppressed viral load. No. unsuppressed PLHIV indicates the estimated number of adult unsuppressed PLHIV per pixel (100m by 100m square). Note: The coloring ramp may not correlate in linear fashion with the count values of unsuppressed PLHIV, and maps should be compared within, rather than across, countries.

Lorenz curves and Gini coefficients

Fig 5 shows the Lorenz curves (cumulative No. of PLHIV vs. ordered cumulative total VL) for each country. The curves’ shapes across the 12 countries were similar, reflecting the large proportions of suppressed PLHIV. The Gini coefficients (Table 1) ranged from 0.789 (Cote d’Ivoire) to 0.939 (Namibia and Rwanda). Drawing from the data informing the Lorenz curves, the proportion of the total VL corresponding to the 1% of PLHIV with the highest VL ranged from 17.1% (Cote d’Ivoire) to 47.2% (Rwanda); whereas the highest 5% of PLHIV accounted for 45.9% (Cote d’Ivoire) to 77.6% (Namibia).

[Figure omitted. See PDF.]

12 Population-based HIV Impact Assessments, 2015–2019. The x-axis displays the cumulative percent of all adult PLHIV ranked by viral load. The y-axis displays the cumulative percent of total VL. PLHIV: people living with HIV.

Viremia, VL, and HIV incidence

Table 2 displays the HIV transmission ratios by country for both sexes combined as well as for opposite-sex transmission; corresponding measures of uncertainty are presented in S3 Table. Transmission ratios for both sexes combined ranged from 1.8 transmission events per 100 unsuppressed PY in Cote d’Ivoire to 15.0 in Uganda. In all countries, female-to-male transmission ratios were lower (ranging from 1.5 in Cote d’Ivoire to 10.6 in Uganda) than male-to-female transmission ratios (ranging from 2.3 in Cote d’Ivoire to 28.3 in Cameroon). Rearranging the data to estimate the number of unsuppressed adults per incident HIV per year suggests that across these 12 countries, in Uganda only 6.6 unsuppressed adults corresponded to each case of incident HIV whereas in Cote d’Ivoire the estimate was 56.1. We observe equally wide ranges when assuming that all transmission was between opposite sexes: For male-to-female transmission, country estimates range from 3.5 unsuppressed males per female incident HIV (Cameroon) to 44.1 in Cote d’Ivoire, whereas for female-to-male transmissions our estimates ranged from 9.4 unsuppressed females in Uganda to 67.4 in Cote d’Ivoire. Similar to transmission ratios, in all countries we estimate fewer unsuppressed males for each female incident HIV than vice versa.

[Figure omitted. See PDF.]

Population HIV Impact Assessments, sub-Saharan Africa, 2015–19.

Table 3 shows the correlation between select VL metrics and incident HIV. Correlations between population viremia and HIV incidence were strong (i.e., R²>0.70) for both sexes combined, for opposite sexes, and between 30+ year-old unsuppressed males and 15–29-year-old HIV-incident females; only the correlation between younger unsuppressed females and older HIV-incident males was moderate (0.51). Gini coefficients also showed a strong correlation with an inverse relationship to incidence. A strong correlation was also seen with expressing a country’s VL as total VL (absolute and log₁₀) with the number of incident HIV infections for both sexes combined as well as between unsuppressed males and HIV-incident females. Whereas a modest correlation (i.e., between 0.50–0.70) was observed between unsuppressed females and HV-incident males. The correlations for both sexes combined weakened when restricting PLHIV to those diagnosed or on treatment. The number of unsuppressed PLHIV per 100 HIV-uninfected adults (No. unsuppressed/100 HIV-neg) as well as HIV seroprevalence showed the strongest correlations with HIV incidence (R² = 0.92).

[Figure omitted. See PDF.]

12 Population-based HIV Impact Assessments, 2015–19.

Discussion

We examined population VL in 12 SSA countries in an attempt to add to the understanding of HIV epidemiology in this region. Key observations include that female unsuppressed PLHIV continue to outnumber males while ART-naïve male PLHIV have consistently higher VL than females. HIV transmission ratios were higher for males than females and vary significantly geographically while the geospatial distribution of population viremia–for some countries—appears markedly different than that for HIV seroprevalence in select countries.

Our study benefitted from several strengths related to the PHIAs: We utilized standardized data from 12 large high-quality [31] national population-based surveys that allowed for the concurrent estimation of HIV incidence, prevalence, and VL-related metrics. The surveys’ incidence estimates also inform UNAIDS-led national HIV incidence estimates and treatment status was measured through ARV testing. Nevertheless, several limitations need to be noted. Our surveys’ antibody-based HIV testing algorithms could not detect acute HIV, a stage characterized by seronegativity and rapidly rising VL. We assumed that the overall level of viremia is long-term. We did not consider behavioral effects (e.g., partner change rates, condom use, male-to-male sex, sex work, or injection drug use) nor could we link HIV transmitting with HIV acquiring persons at the individual level. Further, we could not consider the contribution of HIV subtypes or HIV disease stage (primary, latent, late stage) to HIV transmission, although we did observe a few individuals with very high VL, suggesting primary HIV infection. Our estimated transmission ratios were accompanied by wide standard errors, mostly reflecting the small numbers of recent HIV infections detected in our surveys. Our surveys’ cross-sectional design does not allow us to infer a causal relationship between viremia and incidence; as well, all our surveys’ recent HIV infections were viremic as defined per our testing algorithm for recent HIV-1 infection.

As suggested in earlier work by others (e.g., [4, 32]) population viremia can be a more suitable indicator for monitoring HIV epidemics as it better reflects the current burden of disease as well as the remaining potential for HIV-related morbidity, mortality and onward transmission compared to conventional HIV prevalence estimates. Effective ART programs removed the bulk of these countries’ total population VL, “concentrating” the remaining viremia in a shrinking subpopulation of unsuppressed PLHIV. Looking at the remaining prevalence of viremia rather than HIV seroprevalence, allowed us to paint a partially different picture of SSA’s HIV epidemics, one where e.g., Lesotho and Eswatini no longer face a hyperendemic epidemic (i.e., <15% viremia prevalence among adults [33]), Rwanda no longer has a generalized epidemic (<1%) and Kenya (1.4%) is close to achieving the same milestone.

Similarly, our viremia-based maps offer a more informative display of the spatial distribution of viremia compared to HIV seroprevalence. In some countries, the spatial distribution of viremia prevalence simply reflected the spatial distribution of HIV infection (Cote d’Ivoire, Malawi, Uganda). In other countries (Kenya, Cameroon, Namibia) clear differences emerged, apparent HIV “hotspots” subsided and “viremia hotspots” surfaced as distinct from HIV prevalence hotspots, illustrating areas with poor viral suppression. Other country comparisons were more mixed, nuanced, and intermediate, showing some areas with varying viral suppression. Our third type of maps, displaying spatial counts of unsuppressed PLHIV, often differed markedly from the two proportional displays of HIV and viremia. These count maps are a function of population density, HIV prevalence, and the suppressive effect of ART programs on VL. It is likely this latter type of map that may prove most useful for program planning and evaluation as it reflects the absolute burden of remaining unsuppressed HIV. In contrast, HIV seroprevalence and prevalence of viremia do not inform about the absolute burden (i.e., number of people with) of HIV or unsuppressed HIV, respectively. Despite the surveys’ large sample sizes, the number of detected participants with recent HIV infections was too small to meaningfully examine the spatial distribution of recent HIV infection alongside viremia.

A population’s mean and sum total VL are rarely used metrics reflecting the absolute size of its transmission and morbidity potential [7, 8] and was shown to correlate with the number of new HIV diagnoses [9]. In our study, the 12 countries’ total VL varied 19-fold, a function of population size, number of PLHIV, and ART uptake, and was strongly correlated with HIV incidence (both sexes), albeit much less so when restricted to female total VL.

The log₁₀ transformed mean VL among unsuppressed PLHIV varied substantially across countries. We did not examine possible factors such as HIV subtypes or the prevalence of co-infections (e.g., herpes simplex virus type 2) [34] for this variation. The literature suggests that the per act HIV transmission risk is closely correlated to log₁₀ plasma VL [14] which may help explain some of the variation for HIV incidence across countries. Of note, unsuppressed males in all countries showed higher mean log₁₀ VL values compared to unsuppressed females, an observation made previously, albeit in smaller settings [35]. This is consistent with our estimated transmission ratios which in all countries were larger for male-to-female than for female-to-male transmission although without doubt other behavioral and biologic factors [36] disproportionately contribute to female HIV acquisition as well. The highest mean log₁₀ VL (5.15) was estimated for unsuppressed males in Cameroon, the same country for which we also estimated the highest male-to-female transmission ratio (28.3).

Although females generally achieve higher proportions of VLS compared to males in all PHIA countries, examining the absolute numbers of unsuppressed PLHIV suggests that in all countries unsuppressed females still outnumber unsuppressed males, due to the higher HIV prevalence observed among females. This finding may serve as a reminder that while more progress towards VLS among males is warranted, most of the absolute unsuppressed HIV burden in SSA continues to be shouldered by females. These data also suggest that estimating the number of unsuppressed PLHIV across demographic and behavioral strata has its own merit, beyond VLS monitoring.

Lorenz curves, commonly used to examine income inequality, are an uncommon way to examine VL distribution [10]. While perhaps less useful or needed for routine epidemic monitoring, the extreme disparities depicted in our Lorenz curves illustrate the large-scale successes of most national ART programs that resulted in VLS in many or most PLHIV. As a result, at the time of these surveys, the 1% of adult PLHIV with the highest VL on average accounted for one third of countries’ remaining total VL, and the highest 5% accounted for two-thirds of total VL. Identifying, linking, and improving adherence of this small population segment to ART programs would suggest disproportional gains in reducing total VL and the corresponding HIV transmission potential, albeit our data did not allow us to infer that PLHIV with (very) high VL are also at higher transmission risk behaviorally. The corresponding Gini coefficients’ large values (close to 1) express these extreme but desired inequalities numerically and display a strong correlation with HIV incidence.

Examining VLS against ARV status showed that between 9% and 21% of PLHIV who tested ARV-negative were nevertheless suppressed. Given SSA’s public health approach to ART provision, we believe that few if any of our survey respondents were taking ARVs that our testing did not cover. These rather large proportions may be due to short-term lapses in ARV intake among some participants whereas a minority may have been elite controllers [37] or transiently controlled VL [38]. Conversely, 6% to 21% of PLHIV who tested ARV-positive were unsuppressed, indicating poor adherence, drug resistance, or having recently initiated ART. These are sizeable proportions that call for further improvement in adherence and drug resistance monitoring.

Our estimated HIV transmission ratios are similar to those reported for the US [39]. Transmission ratios and number of unsuppressed adults per incident case of HIV highlight a very different scale compared to HIV incidence rates. Whereas incidence rates contain the entire HIV-uninfected population in its denominator (at risk for HIV acquisition), effective transmission ratios indicate the force of HIV transmission from the perspective of the much smaller population of unsuppressed adults (at risk for HIV transmission). HIV transmission ratios may be seen as aggregate estimates of annualized effective reproductive numbers. The much higher values for effective transmission ratios (compared to incidence rates) underscore the importance of this relatively small sub-population of unsuppressed adults for ongoing HIV transmission risk. Most countries showed somewhat similar effective transmission ratios (combined for both sexes) between approximately 10 to 11 HIV transmission events per 100 unsuppressed adults per year; outliers included Uganda (15.0), Zimbabwe (6.3) and Cote d’Ivoire (1.8). Looking at the opposite sex transmission ratios we note two general observations: The male-to-female transmission ratios were uniformly higher than female-to-male transmission ratios and the sex-specific variations across countries were more marked. Cameroon and Cote d’Ivoire specifically warrant mention: Cote d’Ivoire offered the lowest estimated HIV transmission ratios, whereas Cameroon saw the largest difference between male-to-female and female-to-male transmission ratios. The reasons for the marked variations for sex-specific transmission ratios across countries and the extremely high male-to-female effective transmission ratio in Cameroon remain unclear and are beyond the scope of this analysis. Data from the same surveys suggest that the self-reported prevalence of medical male circumcision varied from 9.8% (Malawi) to 58.1% (Cameroon) [17].

Our estimated transmission ratios assume that all transmission occurred 1) within the respective countries, 2) within the sampled adult age ranges (no transmission between HIV-infected children and adults, or, where older age groups remained unsampled, from older age groups to sampled adult age groups), and 3) for sex-specific ratios, between opposite sexes (i.e., male-to-male transmission was ignored). We believe these assumptions are acceptable for entire countries (compared to smaller, subnational studies) although we recognize that our female-to-male transmission ratios are overestimates as some male-to-male transmission will have contributed to the estimated male HIV incidence (UNAIDS estimates 6% and 21% of incident HIV in 2020 in Eastern/Southern and Western/Central Africa was due to male-to-male sex, respectively [40]). To the extent that our estimates for HIV incidence and population viremia are national-level estimates, the resulting HIV transmission ratios incorporate most/all incident HIV occurring inside and outside of stable sexual relationships. We could not account for the potential HIV transmission between sex partners residing in different countries, although the net effect (offsetting ‘exported’ transmissions by ‘imported’ ones) can be expected to be somewhat smaller than the absolute number of such transmissions. Even so, our surveys sampled all ‘de-facto’ household members and visitors who slept in the sampled dwelling the night before, regardless of nationality. For the same reasons, estimates of transmission ratios quickly become unreliable when restricted to subgroups, such as age (e.g., young females and older men). Unlike HIV acquisition, no detectable signal can reliably identify transmitting individuals (except for the use of phylogenetic analysis).

Across the 12 PHIA countries, we were able to confirm a strong correlation between population VL and incident HIV, as confirmed in previous studies elsewhere [4, 32]. The substantially higher R² values for transmission from older unsuppressed males to younger HIV-incident females (compared to the opposite scenario) is supported by other studies [41, 42] and highlights again a potentially important component in SSA’s HIV epidemics. Expressing total VL in absolute or log₁₀ values showed no discernible difference. The lower R² values shown when restricting the population to PLHIV who are diagnosed or on treatment serve as a reminder about the advantage of population-based data as compared to routine service data that by definition are restricted to PLHIV utilizing such services [4].

With the large-scale roll-out of ART in SSA, the declining utility of HIV prevalence to monitor the epidemic may be compensated using VL. Given that virtually all transmission stems from unsuppressed individuals, viremia-related metrics indicate a population’s HIV transmission potential, can be estimated with more precision than HIV incidence in the survey or surveillance setting, and are more amenable to programmatic action. As a continuous or categorical metric, VL can inform various indicators for surveillance, evaluation, and monitoring purposes and are an excellent way to describe the spatial and demographic distribution of HIV-related transmission, morbidity, and mortality risk. The population-level prevalence of unsuppressed HIV and the number of unsuppressed PLHIV, disaggregated as needed, stand out in their added utility for epidemic monitoring and HIV programming and should be reported and evaluated whenever feasible.

Supporting information

S1 Table. Adult population, PLHIV, and number unsuppressed at time of survey and adjustments in creation of raster-based population counts.

https://doi.org/10.1371/journal.pone.0275560.s001

(DOCX)

S2 Table. Standard errors and 95% confidence intervals for Table 1.

https://doi.org/10.1371/journal.pone.0275560.s002

(XLSX)

S3 Table. Standard errors and 95% confidence intervals for transmission ratios shown in Table 2.

https://doi.org/10.1371/journal.pone.0275560.s003

(XLSX)

S1 File. IRB approval documentation.

https://doi.org/10.1371/journal.pone.0275560.s004

(DOCX)

S2 File. PLOS questionnaire on inclusivity.

https://doi.org/10.1371/journal.pone.0275560.s005

(DOCX)

Acknowledgments

We thank the PHIA survey respondents for their participation and PHIA study staff for their work. We are grateful for the many constructive comments received from the reviewers of this manuscript.

Disclaimer: The findings and conclusions in this report are those of the author(s) and do not necessarily represent the official position of the U.S. Department of Health and Human Services, the Centers for Disease Control and Prevention.

Citation: Hladik W, Stupp P, McCracken SD, Justman J, Ndongmo C, Shang J, et al. (2023) The epidemiology of HIV population viral load in twelve sub-Saharan African countries. PLoS ONE 18(6): e0275560. https://doi.org/10.1371/journal.pone.0275560

About the Authors:

Wolfgang Hladik

Roles: Conceptualization, Data curation, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-3869-1711

Paul Stupp

Roles: Conceptualization, Data curation, Formal analysis, Validation, Visualization

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Stephen D. McCracken

Roles: Investigation, Methodology, Visualization, Writing – original draft

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-9180-0366

Jessica Justman

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

ORICD: https://orcid.org/0000-0003-4704-3115

Clement Ndongmo

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0003-0177-8458

Judith Shang

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Emily K. Dokubo

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-1499-0151

Elizabeth Gummerson

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Isabelle Koui

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health, Abidjan, Cote d’Ivoire

Stephane Bodika

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-4294-5518

Roger Lobognon

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Hermann Brou

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Caroline Ryan

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Kristin Brown

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0001-9256-533X

Harriet Nuwagaba-Biribonwoha

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Leonard Kingwara

Roles: Data curation, Writing – review & editing

Affiliation: National AIDS and STI’s Control Programme, Ministry of Health, Nairobi, Kenya

Peter Young

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-7432-2876

Megan Bronson

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Duncan Chege

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Optatus Malewo

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Yohannes Mengistu

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Frederix Koen

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Andreas Jahn

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health, Lilongwe, Malawi

Andrew Auld

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0001-5089-9163

Sasi Jonnalagadda

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Elizabeth Radin

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Ndapewa Hamunime

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health, Windhoek, Namibia

Daniel B. Williams

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0003-1609-5748

Eugenie Kayirangwa

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Veronicah Mugisha

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Rennatus Mdodo

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Stephen Delgado

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Wilford Kirungi

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health, Kampala, Uganda

Lisa Nelson

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Christine West

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Samuel Biraro

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Kumbutso Dzekedzeke

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Danielle Barradas

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Owen Mugurungi

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health and Child Care, Harare, Zimbabwe

Shirish Balachandra

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Peter H. Kilmarx

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0001-6464-3345

Godfrey Musuka

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Hetal Patel

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Bharat Parekh

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-1088-4263

Katrina Sleeman

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

Robert A. Domaoal

Roles: Data curation, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-0629-9459

George Rutherford

Roles: Data curation, Writing – review & editing

Affiliation: University of California San Francisco, San Francisco, California, United States of America

ORICD: https://orcid.org/0000-0002-0707-8287

Tsietso Motsoane

Roles: Data curation, Writing – review & editing

Affiliation: Ministry of Health, Maseru, Lesotho

Anne-Cécile Zoung-Kanyi Bissek

Roles: Data curation, Writing – review & editing

Affiliations: Division of Operational Research for Health, Ministry of Health, Yaoundé, Cameroon, Faculty of Medicine and Biomedical Sciences, University of Yaoundé, Yaoundé, Cameroon

Mansoor Farahani

Roles: Data curation, Writing – review & editing

Affiliation: ICAP at Columbia University, New York, New York, United States of America

Andrew C. Voetsch

Roles: Data curation, Project administration, Writing – review & editing

Affiliation: Division of Global HIV and TB, US Centers for Disease Control and Prevention (CDC), Atlanta, GA, United States of America

ORICD: https://orcid.org/0000-0002-0783-1312

References

1. UNAIDS Fact Sheet 2022. Joint United Nations Programme on AIDS, Geneva, 2022. Accessed at https://www.unaids.org/en/resources/fact-sheet.

2. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Antiretroviral Therapy for the Prevention of HIV-1 Transmission. N Engl J Med. 2016;375(9):830–9. Epub 2016/07/19. pmid:27424812; PubMed Central PMCID: PMC5049503.

3. Prevailing Against Pandemics. WAD Report 2020. Joint United Nations Programme on HIV/AIDS. Geneva. 2020. Accessed 2-Jan-2023 at https://aidstargets2025.unaids.org/assets/images/prevailing-against-pandemics_en.pdf.

4. Tanser F, Vandormael A, Cuadros D, Phillips AN, de Oliveira T, Tomita A, et al. Effect of population viral load on prospective HIV incidence in a hyperendemic rural African community. Sci Transl Med. 2017;9(420). Epub 2017/12/15. pmid:29237762; PubMed Central PMCID: PMC6435884.

5. Havlir D, Lockman S, Ayles H, Larmarange J, Chamie G, Gaolathe T, et al. What do the Universal Test and Treat trials tell us about the path to HIV epidemic control? Journal of the International AIDS Society. 2020;23(2):e25455. pmid:32091179

6. Montaner JS, Lima VD, Barrios R, Yip B, Wood E, Kerr T, et al. Association of highly active antiretroviral therapy coverage, population viral load, and yearly new HIV diagnoses in British Columbia, Canada: a population-based study. The Lancet. 2010;376(9740):532–9. pmid:20638713

7. Miller WC, Powers KA, Smith MK, Cohen MS. Community viral load as a measure for assessment of HIV treatment as prevention. The Lancet infectious diseases. 2013;13(5):459–64. Epub 2013/03/30. pmid:23537801; PubMed Central PMCID: PMC4512165.

8. Castel AD, Befus M, Willis S, Griffin A, West T, Hader S, et al. Use of the community viral load as a population-based biomarker of HIV burden. AIDS. 2012;26(3):345–53. Epub 2011/10/20. pmid:22008660.

9. Das M, Chu PL, Santos GM, Scheer S, Vittinghoff E, McFarland W, et al. Decreases in community viral load are accompanied by reductions in new HIV infections in San Francisco. PLoS One. 2010;5(6):e11068. Epub 2010/06/16. pmid:20548786; PubMed Central PMCID: PMC2883572.

10. Christopoulos KA, Hartogensis W, Glidden DV, Pilcher CD, Gandhi M, Geng EH. The Lorenz curve: a novel method for understanding viral load distribution at the population level. J AIDS. 2017;31(2):309. pmid:27831945

11. Gopalappa C, Farnham PG, Chen YH, Sansom SL. Progression and Transmission of HIV/AIDS (PATH 2.0). Med Decis Making. 2017;37(2):224–33. Epub 2016/09/21. pmid:27646567.

12. Bonacci RA, Bradley H, Rosenberg ES, Holtgrave DR. Evaluating HIV Transmission Rates for the US National HIV/AIDS Strategy, 2010–2015. J Acquir Immune Defic Syndr. 2019;82(2):e37–e9. Epub 2019/06/11. pmid:31180996.

13. Conan N, Badawi M, Chihana ML, Wanjala S, Kingwara L, Mambula C, et al. Two‐fold increase in the HIV viral load suppression rate along with decreased incidence over six years in Ndhiwa sub‐county, Kenya. Tropical Medicine & International Health. 2021;26(12):1609–15.

14. Quinn TC, Wawer MJ, Sewankambo N, Serwadda D, Li C, Wabwire-Mangen F, et al. Viral load and heterosexual transmission of human immunodeficiency virus type 1. Rakai Project Study Group. N Engl J Med. 2000;342(13):921–9. Epub 2000/03/30. pmid:10738050.

15. Justman JE, Mugurungi O, El-Sadr WM. HIV Population Surveys—Bringing Precision to the Global Response. N Engl J Med. 2018;378(20):1859–61. Epub 2018/05/17. pmid:29768142.

16. PHIA Project—A Drop That Counts—Resources. ICAP Columbia University. New York. Accessed on 21-Jul-2020 at https://phia.icap.columbia.edu/resources/.

17. Population HIV Impact Assessment Data Sets. Accessed on 2-Jan-2023 at https://phia-data.icap.columbia.edu/datasets.

18. Patel HK, Duong YT, Birhanu S, Dobbs T, Lupoli K, Moore C, et al. A comprehensive approach to assuring quality of laboratory testing in HIV surveys: lessons learned from the Population-Based HIV Impact Assessment Project. JAIDS Journal of Acquired Immune Deficiency Syndromes. 2021;87:S17–S27. pmid:34166309

19. Saito S, Duong YT, Metz M, Lee K, Patel H, Sleeman K, et al. Returning HIV‐1 viral load results to participant‐selected health facilities in national Population‐based HIV Impact Assessment (PHIA) household surveys in three sub‐Saharan African Countries, 2015 to 2016. Journal of the International AIDS Society. 2017;20:e25004. pmid:29171193

20. Duong YT, Qiu M, De AK, Jackson K, Dobbs T, Kim AA, et al. Detection of recent HIV-1 infection using a new limiting-antigen avidity assay: potential for HIV-1 incidence estimates and avidity maturation studies. PLoS One. 2012;7(3):e33328. Epub 2012/04/06. pmid:22479384; PubMed Central PMCID: PMC3314002.

21. Voetsch AC, Duong YT, Stupp P, Saito S, McCracken S, Dobbs T, et al. HIV-1 recent infection testing algorithm with antiretroviral drug detection to improve accuracy of incidence estimates. JAIDS Journal of Acquired Immune Deficiency Syndromes. 2021;87:S73–S80. pmid:34166315

22. Population-based HIV Impact Assessment (PHIA) Geospatial Data Use Manual. New York, NY. October 2019.

23. Larmarange J, Bendaud V. HIV estimates at second subnational level from national population-based surveys. AIDS. 2014;28 Suppl 4:S469–76. Epub 2014/11/20. pmid:25406750; PubMed Central PMCID: PMC4247267.

24. Larmarange J, Vallo R, Yaro S, Msellati P, Méda N. Methods for mapping regional trends of HIV prevalence from Demographic and Health Surveys (DHS). Cybergeo: European Journal of Geography. Document 558. 2011.

25. QGIS Geographic Information System. Open Source Geospatial Foundation Project. http://qgis.org.

26. Spatial Data Repository, The Demographic and Health Surveys Program. ICF International. Available from http://spatialdata.dhsprogram.com/home/. Accessed 17 Feb 2021 (Uganda) and 10 Apr 2019 (Tanzania).

27. Database of Global Administrative Areas (GADM). Accessed 17 Mar 2021 at https://gadm.org/index.html.

28. WorldPop (School of Geography and Environmental Science UoSDoGaG, University of Louisville; Departement de Geographie, Universite de Namur) and Center for International Earth Science Information Network (CIESIN), Columbia University (2018). Global High Resolution Population Denominators Project—Funded by The Bill and Melinda Gates Foundation (OPP1134076). Accessed 17 March 2021 https://doi.org/10.5258/SOTON/WP00660.

29. Lorenz MO. Methods of measuring the concentration of wealth. J Publications of the American Statistical Association. 1905;9(70):209–19.

30. Gini C. Concentration and dependency ratios. J Rivista di Politica Economica. 1997;87:769–92.

31. Sachathep K, Radin E, Hladik W, Hakim A, Saito S, Burnett J, et al. Population-based HIV Impact Assessments Survey methods, response, and quality in Zimbabwe, Malawi, and Zambia. JAIDS Journal of Acquired Immune Deficiency Syndromes. 2021;87:S6–S16. pmid:34166308

32. Solomon SS, Mehta SH, McFall AM, Srikrishnan AK, Saravanan S, Laeyendecker O, et al. Community viral load, antiretroviral therapy coverage, and HIV incidence in India: a cross-sectional, comparative study. Lancet HIV. 2016;3(4):e183–90. Epub 2016/04/03. pmid:27036994; PubMed Central PMCID: PMC4863069.

33. UNAIDS Terminology Guidelines. Geneva, 2015. Accessed 21 May 2023 at https://www.unaids.org/sites/default/files/media_asset/2015_terminology_guidelines_en.pdf.

34. Baggaley RF, Hollingsworth TD. Brief report: HIV-1 transmissions during asymptomatic infection: exploring the impact of changes in HIV-1 viral load due to coinfections. J Acquir Immune Defic Syndr. 2015;68(5):594–8. Epub 2015/01/15. pmid:25585299; PubMed Central PMCID: PMC4387204.

35. Gray RH, Li X, Wawer MJ, Serwadda D, Sewankambo NK, Wabwire-Mangen F, et al. Determinants of HIV-1 load in subjects with early and later HIV infections, in a general-population cohort of Rakai, Uganda. J Infect Dis. 2004;189(7):1209–15. Epub 2004/03/20. pmid:15031789.

36. Ramjee G, Daniels B. Women and HIV in sub-Saharan Africa. J AIDS Research. 2013;10(30):1–9. pmid:24330537

37. Okulicz JF, Lambotte O. Epidemiology and clinical characteristics of elite controllers. Curr Opin HIV AIDS. 2011;6(3):163–8. Epub 2011/04/20. pmid:21502920.

38. Madec Y, Boufassa F, Porter K, Meyer L, Collaboration C. Spontaneous control of viral load and CD4 cell count progression among HIV-1 seroconverters. AIDS. 2005;19(17):2001–7. Epub 2005/11/02. pmid:16260907.

39. Li Z, Purcell DW, Sansom SL, Hayes D, Hall HI. Vital Signs: HIV Transmission Along the Continuum of Care—United States, 2016. MMWR Morb Mortal Wkly Rep. 2019;68(11):267–72. Epub 2019/03/22. pmid:30897075; PubMed Central PMCID: PMC6478059 potential conflicts of interest. No potential conflicts of interest were disclosed.

40. UNAIDS Data 2020. Joint United Nations Programme on HIV/AIDS. Geneva. 2020. Accessed 2-Sep-2021 at https://www.unaids.org/sites/default/files/media_asset/2020_aids-data-book_en.pdf.

41. Beauclair R, Helleringer S, Hens N, Delva W. Age differences between sexual partners, behavioural and demographic correlates, and HIV infection on Likoma Island, Malawi. J Scientific Reports. 2016;6:36121. pmid:27805053

42. Schaefer R, Gregson S, Eaton JW, Mugurungi O, Rhead R, Takaruza A, et al. Age-disparate relationships and HIV incidence in adolescent girls and young women: evidence from Zimbabwe. J AIDS. 2017;31(10):1461–70. pmid:28426534