Bayesian Estimation of MSM Population Size in Côte d'Ivoire

Population size, HIV prevalence, and antiretroviral therapy coverage among key populations in sub-Saharan Africa: collation and synthesis of survey data, 2010-23

BACKGROUND

Key population HIV programmes in sub-Saharan Africa require epidemiological information to ensure equitable and universal access to effective services. We aimed to consolidate and harmonise survey data among female sex workers, men who have sex with men, people who inject drugs, and transgender people to estimate key population size, HIV prevalence, and antiretroviral therapy (ART) coverage for countries in mainland sub-Saharan Africa.

METHODS

Key population size estimates, HIV prevalence, and ART coverage data from 39 sub-Saharan Africa countries between 2010 and 2023 were collated from existing databases and verified against source documents. We used Bayesian mixed-effects spatial regression to model urban key population size estimates as a proportion of the gender-matched, year-matched, and area-matched population aged 15-49 years. We modelled subnational key population HIV prevalence and ART coverage with age-matched, gender-matched, year-matched, and province-matched total population estimates as predictors.

FINDINGS

We extracted 2065 key population size data points, 1183 HIV prevalence data points, and 259 ART coverage data points. Across national urban populations, a median of 1·65% (IQR 1·35-1·91) of adult cisgender women were female sex workers, 0·89% (0·77-0·95) were men who have sex with men, 0·32% (0·31-0·34) were men who injected drugs, and 0·10% (0·06-0·12) were women who were transgender. HIV prevalence among key populations was, on average, four to six times higher than matched total population prevalence, and ART coverage was correlated with, but lower than, the total population ART coverage with wide heterogeneity in relative ART coverage across studies. Across sub-Saharan Africa, key populations were estimated as comprising 1·2% (95% credible interval 0·9-1·6) of the total population aged 15-49 years but 6·1% (4·5-8·2) of people living with HIV.

INTERPRETATION

Key populations in sub-Saharan Africa experience higher HIV prevalence and lower ART coverage, underscoring the need for focused prevention and treatment services. In 2024, limited data availability and heterogeneity constrain precise estimates for programming and monitoring trends. Strengthening key population surveys and routine data within national HIV strategic information systems would support more precise estimates.

FUNDING

UNAIDS, Bill & Melinda Gates Foundation, and US National Institutes of Health.

Challenges and Opportunities in Big Data Science to Address Health Inequities and Focus the HIV Response

PURPOSE OF REVIEW

Big Data Science can be used to pragmatically guide the allocation of resources within the context of national HIV programs and inform priorities for intervention. In this review, we discuss the importance of grounding Big Data Science in the principles of equity and social justice to optimize the efficiency and effectiveness of the global HIV response.

RECENT FINDINGS

Social, ethical, and legal considerations of Big Data Science have been identified in the context of HIV research. However, efforts to mitigate these challenges have been limited. Consequences include disciplinary silos within the field of HIV, a lack of meaningful engagement and ownership with and by communities, and potential misinterpretation or misappropriation of analyses that could further exacerbate health inequities. Big Data Science can support the HIV response by helping to identify gaps in previously undiscovered or understudied pathways to HIV acquisition and onward transmission, including the consequences for health outcomes and associated comorbidities. However, in the absence of a guiding framework for equity, alongside meaningful collaboration with communities through balanced partnerships, a reliance on big data could continue to reinforce inequities within and across marginalized populations.

An introduction to bayesian spatial smoothing methods for disease mapping: modeling county firearm suicide mortality rates

This article introduces bayesian spatial smoothing models for disease mapping-a specific application of small area estimation where the full universe of data is known-to a wider audience of public health professionals using firearm suicide as a motivating example. Besag, York, and Mollié (BYM) Poisson spatial and space-time smoothing models were fitted to firearm suicide counts for the years 2014-2018. County raw death rates in 2018 ranged from 0 to 24.81 deaths per 10 000 people. However, the highest mortality rate was highly unstable, based on only 2 deaths in a population of approximately 800, and 80.5% of contiguous US counties experienced fewer than 10 firearm suicide deaths and were thus suppressed. Spatially smoothed county firearm suicide mortality estimates ranged from 0.06 to 4.05 deaths per 10 000 people and could be reported for all counties. The space-time smoothing model produced similar estimates with narrower credible intervals as it allowed counties to gain precision from adjacent neighbors and their own counts in adjacent years. bayesian spatial smoothing methods are a useful tool for evaluating spatial health disparities in small geographies where small numbers can result in highly variable rate estimates, and new estimation techniques in R software have made fitting these models more accessible to researchers.

Leveraging Routinely Collected Program Data to Inform Extrapolated Size Estimates for Key Populations in Namibia: Small Area Estimation Study

BACKGROUND

Estimating the size of key populations, including female sex workers (FSW) and men who have sex with men (MSM), can inform planning and resource allocation for HIV programs at local and national levels. In geographic areas where direct population size estimates (PSEs) for key populations have not been collected, small area estimation (SAE) can help fill in gaps using supplemental data sources known as auxiliary data. However, routinely collected program data have not historically been used as auxiliary data to generate subnational estimates for key populations, including in Namibia.

OBJECTIVE

To systematically generate regional size estimates for FSW and MSM in Namibia, we used a consensus-informed estimation approach with local stakeholders that included the integration of routinely collected HIV program data provided by key populations' HIV service providers.

METHODS

We used quarterly program data reported by key population implementing partners, including counts of the number of individuals accessing HIV services over time, to weight existing PSEs collected through bio-behavioral surveys using a Bayesian triangulation approach. SAEs were generated through simple imputation, stratified imputation, and multivariable Poisson regression models. We selected final estimates using an iterative qualitative ranking process with local key population implementing partners.

RESULTS

Extrapolated national estimates for FSW ranged from 4777 to 13,148 across Namibia, comprising 1.5% to 3.6% of female individuals aged between 15 and 49 years. For MSM, estimates ranged from 4611 to 10,171, comprising 0.7% to 1.5% of male individuals aged between 15 and 49 years. After the inclusion of program data as priors, the estimated proportion of FSW derived from simple imputation increased from 1.9% to 2.8%, and the proportion of MSM decreased from 1.5% to 0.75%. When stratified imputation was implemented using HIV prevalence to inform strata, the inclusion of program data increased the proportion of FSW from 2.6% to 4.0% in regions with high prevalence and decreased the proportion from 1.4% to 1.2% in regions with low prevalence. When population density was used to inform strata, the inclusion of program data also increased the proportion of FSW in high-density regions (from 1.1% to 3.4%) and decreased the proportion of MSM in all regions.

CONCLUSIONS

Using SAE approaches, we combined epidemiologic and program data to generate subnational size estimates for key populations in Namibia. Overall, estimates were highly sensitive to the inclusion of program data. Program data represent a supplemental source of information that can be used to align PSEs with real-world HIV programs, particularly in regions where population-based data collection methods are challenging to implement. Future work is needed to determine how best to include and validate program data in target settings and in key population size estimation studies, ultimately bridging research with practice to support a more comprehensive HIV response.

Comparison of Empirically Derived and Model-Based Estimates of Key Population HIV Incidence and the Distribution of New Infections by Population Group in Sub-Saharan Africa

BACKGROUND

The distribution of new HIV infections among key populations, including female sex workers (FSWs), gay men and other men who have sex with men (MSM), and people who inject drugs (PWID) are essential information to guide an HIV response, but data are limited in sub-Saharan Africa (SSA). We analyzed empirically derived and mathematical model-based estimates of HIV incidence among key populations and compared with the Joint United Nations Programme on HIV/AIDS (UNAIDS) estimates.

METHODS

We estimated HIV incidence among FSW and MSM in SSA by combining meta-analyses of empirical key population HIV incidence relative to the total population incidence with key population size estimates (KPSE) and HIV prevalence. Dynamic HIV transmission model estimates of HIV incidence and percentage of new infections among key populations were extracted from 94 country applications of 9 mathematical models. We compared these with UNAIDS-reported distribution of new infections, implied key population HIV incidence and incidence-to-prevalence ratios.

RESULTS

Across SSA, empirical FSW HIV incidence was 8.6-fold (95% confidence interval: 5.7 to 12.9) higher than total population female 15-39 year incidence, and MSM HIV incidence was 41.8-fold (95% confidence interval: 21.9 to 79.6) male 15-29 year incidence. Combined with KPSE, these implied 12% of new HIV infections in 2021 were among FSW and MSM (5% and 7% respectively). In sensitivity analysis varying KPSE proportions within 95% uncertainty range, the proportion of new infections among FSW and MSM was between 9% and 19%. Insufficient data were available to estimate PWID incidence rate ratios. Across 94 models, median proportion of new infections among FSW, MSM, and PWID was 6.4% (interquartile range 3.2%-11.7%), both much lower than the 25% reported by UNAIDS.

CONCLUSION

Empirically derived and model-based estimates of HIV incidence confirm dramatically higher HIV risk among key populations in SSA. Estimated proportions of new infections among key populations in 2021 were sensitive to population size assumptions and were substantially lower than estimates reported by UNAIDS.

Summarizing methods for estimating population size for key populations: a global scoping review for human immunodeficiency virus research

BACKGROUND

Estimating the population sizes of key populations(people who inject drugs, men who have sex with men, transgender persons, and commercial sex workers) is critical for understanding the overall Human Immunodeficiency Virus burden. This scoping review aims to synthesize existing methods for population size estimation among key populations, and provide recommendations for future application of the existing methods.

METHODS

Relevant studies published from 1st January 2000 to 4th August 2020 and related to key population size estimation were retrieved and 120 of 688 studies were assessed. After reading the full texts, 81 studies were further excluded. Therefore, 39 studies were included in this scoping review. Estimation methods included five digital methods, one in-person method, and four hybrid methods.

FINDING

We summarized and organized the methods for population size estimateion into the following five categories: methods based on independent samples (including capture-recapture method and multiplier method), methods based on population counting (including Delphi method and mapping method), methods based on the official report (including workbook method), methods based on social network (including respondent-driven sampling method and network scale-up method) and methods based on data-driven technologies (Bayesian estimation method, Stochastic simulation method, and Laska, Meisner, and Siegel estimation method). Thirty-six (92%) articles were published after 2010 and 23 (59%) used multiple methods. Among the articles published after 2010, 11 in high-income countries and 28 in low-income countries. A total of 10 estimated the size of commercial sex workers, 14 focused on men who have sex with men, and 10 focused on people who inject drugs.

CONCLUSIONS

There was no gold standard for population size estimation. Among 120 studies that were related to population size estimation of key populations, the most commonly used population estimation method is the multiplier method (26/120 studies). Every method has its strengths and biases. In recent years, novel methods based on data-driven technologies such as Bayesian estimation have been developed and applied in many surveys.

The effect of men who have sex with men (MSM) on the spread of sexually transmitted infections

Sexually transmitted infections (STIs) have remained a worldwide public health threat. It is difficult to control the spread of STIs, not only because of heterogeneous sexual transmission between men and women but also because of the complicated effects of sexual transmission among men who have sex with men (MSM) and mother-to-child transmission. Many studies point to the existence of a ‘bisexual bridge’, where STIs spread from the MSM network via bisexual connections. However, it is unclear how the MSM network affects heterosexual networks as well as mother-to-child transmission. To analyse the effect of MSM on the spread of STIs, we divided the population into four subpopulations: (i) women, (ii) men who have sex with women only (MSW), (iii) men who have sex with both men and women (MSMW), (iv) men who have sex with men exclusively (MSME). We calculated the type-reproduction numbers of these four subpopulations, and our analysis determined what preventive measures may be effective. Our analysis shows the impact of bisexual bridge on the spread of STIs does not outweigh their population size. Since MSM and mother-to-child transmission rates do not have a strong synergistic effect when combined, complementary prevention measures are needed. The methodologies and findings we have provided here will contribute greatly to the future development of public health.

Estimating the Population Size of Female Sex Workers in Zimbabwe: Comparison of Estimates Obtained Using Different Methods in Twenty Sites and Development of a National-Level Estimate

National-level population size estimates (PSEs) for hidden populations are required for HIV programming and modelling. Various estimation methods are available at the site-level, but it remains unclear which are optimal and how best to obtain national-level estimates.

Using factor analyses to estimate the number of female sex workers across Malawi from multiple regional sources

PURPOSE

Human immunodeficiency virus (HIV) risks are heterogeneous in nature even in generalized epidemics. However, data are often missing for those at highest risk of HIV, including female sex workers. Statistical models may be used to address data gaps where direct, empiric estimates do not exist.

METHODS

We proposed a new size estimation method that combines multiple data sources (the Malawi Biological and Behavioral Surveillance Survey, the Priorities for Local AIDS Control Efforts study, and the Malawi Demographic Household Survey). We used factor analysis to extract information from auxiliary variables and constructed a linear mixed effects model for predicting population size for all districts of Malawi.

RESULTS

On average, the predicted proportion of female sex workers among women of reproductive age across all districts was about 0.58%. The estimated proportions seemed reasonable in comparing with a recent study Priorities for Local AIDS Control Efforts II (PLACE II). Compared with using a single data source, we observed increased precision and better geographic coverage.

CONCLUSIONS

We illustrate how size estimates from different data sources may be combined for prediction. Applying this approach to other subpopulations in Malawi and to countries where size estimate data are lacking can ultimately inform national modeling processes and estimate the distribution of risks and priorities for HIV prevention and treatment programs.

Evolving HIV epidemics: the urgent need to refocus on populations with risk

Purpose of review

To explore the comparative importance of HIV infections among key populations and their intimate partners as HIV epidemics evolve, and to review implications for guiding responses.

Recent findings

Even as concentrated epidemics evolve, new infections among current and former key population members and their intimate partners dominate new infections. Prevalent infections in the general population grow primarily because of key population turnover and infections among their intimate partners. In generalized epidemic settings, data and analysis on key populations are often inadequate to assess the impact of key population-focused responses, so they remain limited in coverage and under resourced. Models must incorporate downstream infections in comparing impacts of alternative responses.

Summary

Recognize that every epidemic is unique, moving beyond the overly simplistic concentrated/generalized epidemic paradigm that can misdirect resources. Guide HIV responses by gathering and using locally relevant data, understanding risk heterogeneity, and applying modeling at both national and sub-national levels to optimize resource allocations among different populations for greatest impact. Translate this improved understanding into clear, unequivocal advice for policymakers on where to focus for impact, breaking them free of the generalized/concentrated paradigm limiting their thinking and affecting their decisions.

The uptake of population size estimation studies for key populations in guiding HIV responses on the African continent

Background

There has been a heightened emphasis on prioritizing data to inform evidence-based HIV responses, including data focused on both defining the content and scale of HIV programs in response to evidence-based need. Consequently, population size estimation (PSE) studies for key populations have become increasingly common to define the necessary scale of specific programs for key populations. This study aims to assess the research utilization of these size estimates in informing HIV policy and program documents across the African continent.

Methods

This study included two phases; Phase 1 was a review of all PSE for key populations, including men who have sex with men (MSM), female sex workers (FSW), people who use drugs (PWUD), and transgender persons in the 54 countries across Africa published from January 2009—December 2017. Phase 2 was a review of 23 different types of documents released between January 2009 –January 2019, with a focus on the US President’s Emergency Plan for AIDS Relief (PEPFAR) and The Global Fund to Fight AIDS, Tuberculosis and Malaria investments, for evidence of stakeholder engagement in PSE studies, as well as key population PSE research utilization to inform HIV programming and international HIV investments.

Results

Of 118 size estimates identified in 39 studies, less than 15% were utilized in PEPFAR Country Operational Plans or national strategic health plan documents, and less than 2% in Global Fund Concept Notes. Of 39 PSE studies, over 50% engaged stakeholders in study implementation and identified target population stakeholders, a third of studies identified policy or program stakeholders, and 15% involved stakeholders in study design.

Conclusion

The past decade has seen an increase in PSE studies conducted for key populations in more generalized HIV epidemic settings which involve significant investments of finances and human resources. However, there remains limited evidence of sustained uptake of these data to guide the HIV responses. Increasing uptake necessitates effective stakeholder engagement and data-oriented capacity building to optimize research utilization and facilitate data-driven and human rights-affirming HIV responses.

Estimating Sizes of Key Populations at the National Level: Considerations for Study Design and Analysis

Background:

National estimates of the sizes of key populations, including female sex workers, men who have sex with men, and transgender women are critical to inform national and international responses to the HIV pandemic. However, epidemiologic studies typically provide size estimates for only limited high priority geographic areas. This article illustrates a two-stage approach to obtain a national key population size estimate in the Dominican Republic using available estimates and publicly available contextual information.

Methods:

Available estimates of key population size in priority areas were augmented with targeted additional data collection in other areas. To combine information from data collected at each stage, we used statistical methods for handling missing data, including inverse probability weights, multiple imputation, and augmented inverse probability weights.

Results:

Using the augmented inverse probability weighting approach, which provides some protection against parametric model misspecification, we estimated that 3.7% (95% CI = 2.9, 4.7) of the total population of women in the Dominican Republic between the ages of 15 and 49 years were engaged in sex work, 1.2% (95% CI = 1.1, 1.3) of men aged 15–49 had sex with other men, and 0.19% (95% CI = 0.17, 0.21) of people assigned the male sex at birth were transgender.

Conclusions:

Viewing the size estimation of key populations as a missing data problem provides a framework for articulating and evaluating the assumptions necessary to obtain a national size estimate. In addition, this paradigm allows use of methods for missing data familiar to epidemiologists.