Dynamics of non-cohabiting sex partnering in sub-Saharan Africa: a modelling study with implications for HIV transmission

Dynamics of Neisseria gonorrhoeae transmission among female sex workers and clients: A mathematical modeling study

BACKGROUND

This study aimed to examine the transmission dynamics of Neisseria gonorrhoeae (NG) in heterosexual sex work networks (HSWNs) and the impact of variation in sexual behavior and interventions on NG epidemiology.

METHODS

The study employed an individual-based mathematical model to simulate NG transmission dynamics in sexual networks involving female sex workers (FSWs) and their clients, primarily focusing on the Middle East and North Africa region. A deterministic model was also used to describe NG transmission from clients to their spouses.

RESULTS

NG epidemiology in HSWNs displays two distinct patterns. In the common low-partner-number HSWNs, a significant proportion of NG incidence occurs among FSWs, with NG prevalence 13 times higher among FSWs than clients, and three times higher among clients than their spouses. Interventions substantially reduce incidence. Increasing condom use from 10 % to 50 % lowers NG prevalence among FSWs, clients, and their spouses from 12.2 % to 6.4 %, 1.2 % to 0.5 %, and 0.4 % to 0.2 %, respectively. Increasing symptomatic treatment coverage among FSWs from 0 % to 100 % decreases prevalence from 10.6 % to 4.5 %, 0.8 % to 0.4 %, and 0.3 % to 0.1 %, respectively. Increasing asymptomatic treatment coverage among FSWs from 0 % to 50 % decreases prevalence from 8.2 % to 0.4 %, 0.6 % to 0.1 %, and 0.2 % to 0.0 %, respectively, with very low prevalence when coverage exceeds 50 %. In high-partner-number HSWNs, prevalence among FSWs saturates at a high level, and the vast majority of incidence occurs among clients and their spouses, with a limited impact of incremental increases in interventions.

CONCLUSION

NG epidemiology in HSWNs is typically a "fragile epidemiology" that is responsive to a range of interventions even if the interventions are incremental, partially efficacious, and only applied to FSWs.

Understanding dynamics and overlapping epidemiologies of HIV, HSV-2, chlamydia, gonorrhea, and syphilis in sexual networks of men who have sex with men

Introduction

We aimed to investigate the overlapping epidemiologies of human immunodeficiency virus (HIV), herpes simplex virus type 2 (HSV-2), chlamydia, gonorrhea, and syphilis in sexual networks of men who have sex with men (MSM), and to explore to what extent the epidemiology of one sexually transmitted infection (STI) relates to or differs from that of another STI.

Methods

An individual-based Monte Carlo simulation model was employed to simulate the concurrent transmission of STIs within diverse sexual networks of MSM. The model simulated sexual partnering, birth, death, and STI transmission within each specific sexual network. The model parameters were chosen based on the current knowledge and understanding of the natural history, transmission, and epidemiology of each considered STI. Associations were measured using the Spearman's rank correlation coefficient (SRCC) and maximal information coefficient (MIC).

Results

A total of 500 sexual networks were simulated by varying the mean and variance of the number of partners for both short-term and all partnerships, degree correlation, and clustering coefficient. HSV-2 had the highest current infection prevalence across the simulations, followed by HIV, chlamydia, syphilis, and gonorrhea. Threshold and saturation effects emerged in the relationship between STIs across the simulated networks, and all STIs demonstrated moderate to strong associations. The strongest current infection prevalence association was between HIV and gonorrhea, with an SRCC of 0.84 (95% CI: 0.80-0.87) and an MIC of 0.81 (95% CI: 0.74-0.88). The weakest association was between HSV-2 and syphilis, with an SRCC of 0.54 (95% CI: 0.48-0.59) and an MIC of 0.57 (95% CI, 0.49-0.65). Gonorrhea exhibited the strongest associations with the other STIs while syphilis had the weakest associations. Across the simulated networks, proportions of the population with zero, one, two, three, four, and five concurrent STI infections were 48.6, 37.7, 11.1, 2.4, 0.3, and < 0.1%, respectively. For lifetime exposure to these infections, these proportions were 13.6, 21.0, 22.9, 24.3, 13.4, and 4.8%, respectively.

Conclusion

STI epidemiologies demonstrate substantial overlap and associations, alongside nuanced differences that shape a unique pattern for each STI. Gonorrhea exhibits an "intermediate STI epidemiology," reflected by the highest average correlation coefficient with other STIs.

Countering the rise of syphilis: A role for doxycycline post-exposure prophylaxis?

Doxycycline post-exposure prophylaxis (PEP) holds the potential to mitigate increasing rates of syphilis among sexual minority men (SMM) in the US yet has received limited attention. Since evaluation of this intervention in actual populations is not currently feasible, we used agent-based models (ABM) to assess the population-level impact of this strategy. We adapted ABM of HIV and HPV transmission, representing a population of 10,230 SMM in Philadelphia, Pennsylvania, US. Parameter inputs were derived from the literature, and ABM outputs during the pre-intervention period were calibrated to local surveillance data. Intervention scenarios varied doxycycline uptake by 20, 40, 60, 80 and 100%, while assuming continued condom use and syphilis screening and treatment. Under each intervention scenario, we incorporated treatment adherence at the following levels: 0, 20, 40, 60, 80 and 100%. Long-term population impact of prophylactic doxycycline was measured using the cumulative incidence over the 10-year period and the percentage of infections prevented attributable to doxycycline at year 10. An uptake scenario of 20% with an adherence level of 80% would reduce the cumulative incidence of infections by 10% over the next decade, translating to 57 fewer cases per 1000 SMM. At year 10, under the same uptake and adherence level, 22% of infections would be prevented due to doxycycline PEP in the instances where condoms were not used or failed. Findings suggest that doxycycline PEP will have a modest impact on syphilis incidence when assuming a reasonable level of uptake and adherence. Doxycycline PEP may be most appropriate as a secondary prevention measure to condoms and enhanced syphilis screening for reducing infections among SMM.

Sexual Behavior Surveys Should Ask More: Covering the Diversity of Sexual Behaviors That May Contribute to the Transmission of Pathogens

Supplemental digital content is available in the text.

Modeling HPV vaccination scale-up among urban young men who have sex with men in the context of HIV

INTRODUCTION

Young men who have sex with men (YMSM) are at high risk to contract human papillomavirus (HPV). While an effective vaccination exists, its use among YMSM is markedly lower compared to non-MSM and women. This study compares scaling up HPV vaccination in conjunction with other prevention strategies.

METHODS

An agent-based model of urban YMSM (≤26 years of age) reflective of the demography of Philadelphia, PA, simulated for up to ten years of follow-up to examine anal and oral transmission of the HPV genotypes covered in the nonavalent (9v) vaccine: 6, 11, 16, 18, 31, 33, 45, 52, 58. Starting HPV prevalences ranged from a high of 18% (type 6) to a low of 6% (type 31); overall 65% of individuals carried any HPV genotype. Simulated levels of vaccination were ranged from 0% to 13% (present-day level), 25%, 50%, 80% (Healthy People 2020 target), and 100% in conjunction with condom use and HIV seroadaptive practices. The primary outcome was the relative reduction in HPV infection.

RESULTS

Compared to present-day vaccination levels (13%), scaling-up vaccination led to expected declines in 10-year post-simulation HPV prevalence. Anal HPV (any 9v types) declined by 9%, 27%, 46%, and 58% at vaccination levels of 25%, 50%, 80%, and 100%, respectively. Similarly, oral HPV (any 9v types) declined by 11%, 33%, 57%, and 71% across the same levels of vaccine uptake. Comparing the prevention strategies, condoms blocked the greatest number of anal transmissions when vaccination was at or below present-day levels. For oral transmission, vaccination was superior to condom use at all levels of coverage.

CONCLUSIONS

Public health HPV preventions strategies should continue to emphasize the complementary roles of condoms and vaccination, especially for preventing oral infection. Improving vaccination coverage will ultimately have the greatest impact on reducing HPV infection among YMSM.

Does infection with Chlamydia trachomatis induce long-lasting partial immunity? Insights from mathematical modelling

OBJECTIVES

To explore whether existence of long-lasting partial immunity against reinfection with Chlamydia trachomatis is necessary to explain C. trachomatis prevalence patterns by age and sexual risk, and to provide a plausible estimate for the effect size, defined here as a reduction in susceptibility to reinfection.

METHODS

A population-based mathematical model was constructed to describe C. trachomatis natural history and transmission dynamics by age and sexual risk. The model was parameterised using natural history, and epidemiological and sexual behaviour data, and applied for UK and US data. Sensitivity analyses were conducted to assess the robustness of predictions to variations in model structure and to examine the impact of alternative assumptions for the mechanism underlying partial immunity.

RESULTS

Partial immunity against reinfection was found necessary to explain observed C. trachomatis prevalence patterns by age and sexual risk. The reduction in susceptibility to reinfection was estimated at 93% using UK data (95% uncertainty interval (UI)=88%-97%) and at 67% using US data (95% UI=24%-88%). The model-structure sensitivity analyses affirmed model predictions. The immunity-mechanism sensitivity analyses suggested a mechanism of susceptibility reduction against reinfection or a mechanism of infectious-period duration reduction upon reinfection.

CONCLUSIONS

A strong long-lasting partial immunity against C. trachomatis reinfection should be present to explain observed prevalence patterns. The mechanism of immunity could be either a reduction in susceptibility to reinfection or a reduction in duration of infection on reinfection. C. trachomatis infection appears to naturally elicit a strong long-lasting immune response, supporting the concept of vaccine development.

Hepatitis C virus infection spontaneous clearance: Has it been underestimated?

OBJECTIVES

Hepatitis C virus (HCV) clearance rate (f_clearance) is defined as the proportion of infected persons who will spontaneously clear their infection after acute infection. We aimed to estimate f_clearance using a novel approach that avoids limitations in existing estimates, and to clarify the link between f_clearance and HCV viremic rate-the latter being the proportion of RNA positivity among those antibody positive.

METHODS

A mathematical model was developed to describe HCV transmission. f_clearance was estimated by fitting the model to probability-based and nationally representative population-based data for Egypt (Egypt 2008 and Egypt 2015) and USA (NHANES A and NHANES B). Uncertainty and sensitivity analyses were conducted.

RESULTS

f_clearance was estimated at 39.9% (95% uncertainty interval (UI): 34.3%-46.4%) and 33.5% (95% UI: 29.2%-38.3%) for Egypt 2008 and Egypt 2015 data, respectively; and at 29.6% (23.0%-37.1%) and 39.9% (31.2%-51.0%) for NHANES A and NHANES B data, respectively. f_clearance was found related to HCV viremic rate through (approximately) the formula f_clearance=1.16 (1-HCV viremic rate). HCV viremic rate was higher with higher risk of HCV exposure. Robustness of results was demonstrated in uncertainty and sensitivity analyses.

CONCLUSION

One-third of HCV-infected persons clear their infection spontaneously, higher than earlier estimates-the immune-system capacity to clear HCV infection may have been underestimated.

Characterizing the temporal evolution of the hepatitis C virus epidemic in Pakistan

Pakistan has the second largest number of HCV infections in the world. We assessed past, present and future levels and trends of the HCV epidemic in Pakistan. An age-structured mathematical model was developed and analysed to describe transmission dynamics over 1980-2050. The model was fitted to a nationally representative survey and a comprehensive database of systematically gathered HCV Ab prevalence data. HCV Ab and chronic infection prevalences peaked at 5.3% and 3.9% in 2000 but were projected to decline to 4.3% and 3.2% by 2017, 3.4% and 2.6% by 2030 and 2.6% and 1.9% by 2050, respectively. The number of chronically infected individuals was estimated at 6 663 906 in 2017 and was projected to peak at 6 665 900 in 2018 and decline to 6 372 100 in 2030 and 5 131 500 in 2050. Annual number of new infections peaked at 346 740 in 1992 but was projected to decline to 198 320 in 2017, 151 090 in 2030 and 98 120 in 2050. Incidence rate per 100 000 person-year peaked at 343 in 1988 but was projected to decline to 99 in 2017, 62 in 2030 and 36 in 2050. Prevalence and incidence varied by age, and the majority of new infections occurred in the 20-39 age group. Prevalence and incidence of HCV in Pakistan have been slowly declining for two decades-Pakistan is enduring a large epidemic that will persist for decades if not controlled. Nearly, 10% of global infections are in Pakistan, with about 200 000 additional infections every year. Rapid and mass scale-up of prevention and treatment programmes are critically needed.

Temporal evolution of HIV sero-discordancy patterns among stable couples in sub-Saharan Africa

INTRODUCTION

Objective was to examine the temporal variation of HIV sero-discordancy in select representative countries (Kenya, Lesotho, Mali, Niger, Tanzania, and Zimbabwe) in sub-Saharan Africa at different HIV epidemic scales. A sero-discordant couple is defined as a stable couple (SC) in which one partner is HIV-positive while the other is HIV-negative.

METHODS

A deterministic compartmental mathematical model was constructed to describe HIV transmission dynamics. The model was pair-based, that is explicitly modeling formation of SCs and infection dynamics in both SCs and in single individuals. The model accommodated for different forms of infection statuses in SCs. Using population-based nationally-representative epidemiologic and demographic input data, historical (1980-2014) and future (2015-2030) trends of sero-discordancy and other demographic and epidemiologic indicators were projected throughout HIV epidemic phases.

RESULTS

As the epidemics emerged, about 90% of SCs affected by HIV were sero-discordant. This proportion declined to 45%-88% at epidemic peak and stabilized as the epidemics started their natural decline. The largest reductions in sero-discordancy were in high HIV-prevalence countries. As the epidemics further declined with antiretroviral therapy (ART) scale-up, the proportion of sero-discordant couples among HIV-affected couples was projected to increase to 70%-92% by 2030. The proportion of sero-discordant couples among all SCs increased as the epidemics emerged and evolved, then peaked at 2%-20% as the epidemics peaked, and then declined as the epidemics declined to reach 0.3%-16% by 2030.

CONCLUSIONS

Sero-discordancy patterns varied with the evolution of the epidemics, and were affected by both epidemic phase and scale. The largest variations were found in high HIV-prevalence countries. The fraction of stable couples that are sero-discordant, as opposed to being sero-concordant positive, was projected to increase with ART scale-up and further HIV incidence decline over the coming two decades. These findings inform strategic planning and resource allocation for interventions among sero-discordant couples.

HIV incidence among people who inject drugs in the Middle East and North Africa: mathematical modelling analysis

INTRODUCTION

Emerging HIV epidemics have been documented among people who inject drugs (PWID) in the Middle East and North Africa (MENA). This study estimates the HIV incidence among PWID due to sharing needles/syringes in MENA. It also delineates injecting drug use role as a driver of the epidemic in the population, and estimates impact of interventions.

METHODS

A mathematical model of HIV transmission among PWID was applied in seven MENA countries with sufficient and recent epidemiological data and HIV prevalence ≥1% among PWID. Estimations of incident and/or prevalent infections among PWID, ex-PWID and sexual partners of infected current and ex-PWID were conducted.

RESULTS

The estimated HIV incidence rate for 2017 among PWID ranged between 0.7% per person-year (ppy) in Tunisia and 7.8% ppy in Pakistan, with Libya being an outlier (24.8% ppy). The estimated number of annual new infections was lowest in Tunisia (n = 79) and Morocco (n = 99), and highest in Iran and Pakistan (approximately n = 6700 each). In addition, 20 to 2208 and 5 to 837 new annual infections were estimated across the different countries among sexual partners of PWID and ex-PWID respectively. Since epidemic emergence, the number of total ever acquired incident infections across countries was 706 to 90,015 among PWID, 99 to 18,244 among sexual partners of PWID, and 16 to 4360 among sexual partners of ex-PWID. The estimated number of prevalent infections across countries was 341 to 23,279 among PWID, 119 to 16,540 among ex-PWID, 67 to 10,752 among sexual partners of PWID, and 12 to 2863 among sexual partners of ex-PWID. Increasing antiretroviral therapy (ART) coverage to the global target of 81% - factoring in ART adherence and current coverage - would avert about half of new infections among PWID and their sexual partners. Combining ART with harm reduction could avert over 90% and 70% of new infections among PWID and their sexual partners respectively.

CONCLUSIONS

There is considerable HIV incidence among PWID in MENA. Of all new infections ultimately due to injecting drug use, about 75% are among PWID and the rest among sexual partners. Of all prevalent infections ultimately attributed to injecting drug use as epidemic driver, about half are among PWID, 30% among ex-PWID and 20% among sexual partners of PWID and ex-PWID. These findings call for scale-up of services for PWID, including harm reduction as well as testing and treatment services.

The Effect of PrEP on HIV Incidence Among Men Who Have Sex With Men in the Context of Condom Use, Treatment as Prevention, and Seroadaptive Practices

BACKGROUND

HIV preexposure prophylaxis (PrEP) is an effective tool in preventing HIV infection among high-risk men who have sex with men (MSM). It is unknown how effective PrEP is in the context of other implemented HIV prevention strategies, including condom use, seroadaption, and treatment as prevention (TasP). We evaluate the impact of increasing uptake of PrEP in conjunction with established prevention strategies on HIV incidence in a high-risk population of MSM through simulation.

METHODS

Agent-based simulation models representing the sexual behavior of high-risk, urban MSM in the United States over the period of 1 year were used to evaluate the effect of PrEP on HIV infection rates. Simulations included data for 10,000 MSM and compared increasing rates of PrEP uptake under 8 prevention paradigms: no additional strategies, TasP, condom use, seroadaptive behavior, and combinations thereof.

RESULTS

We observed a mean of 103.2 infections per 10,000 MSM in the absence of any prevention method. PrEP uptake at 25% without any additional prevention strategies prevented 30.7% of infections. In the absence of PrEP, TasP, condom use, and seroadaptive behavior independently prevented 27.1%, 48.8%, and 37.7% of infections, respectively, and together prevented 72.2%. The addition of PrEP to the 3 aforementioned prevention methods, at 25% uptake, prevented an additional 5.0% of infections.

CONCLUSIONS

To achieve a 25% reduction in HIV infections by 2020, HIV prevention efforts should focus on significantly scaling up access to PrEP in addition to HIV testing, access to antiretroviral therapy, and promoting condom use.

HIV and herpes simplex virus type 2 epidemiological synergy: misguided observational evidence? A modelling study

Objectives

To investigate whether observational studies of HIV and herpes simplex virus type 2 (HSV-2) infections have the capacity to assess the HIV/HSV-2 epidemiological synergy.

Methods

An individual-based Monte Carlo model was used to simulate HIV/HSV-2 epidemics in two scenarios: no HIV/HSV-2 biological interaction and HSV-2 seropositivity enhancing HIV acquisition. Cross-sectional observational studies were simulated by sampling individuals from the population to assess resulting crude and adjusted ORs of the HIV/HSV-2 association. Meta-analyses were conducted to estimate the pooled mean ORs. Impact of under-reporting of sexual behaviour and miscapture of high-risk individuals was assessed through sensitivity analyses.

Results

Assuming no HIV/HSV-2 biological interaction, the crude HIV/HSV-2 OR ranged between 1.38 and 9.93, with a pooled mean of 6.45 (95% CI 5.81 to 7.17). Adjustment for the number of sexual partners over last year, over lifetime and for both partner numbers simultaneously reduced the mean OR to 5.45 (95% CI 4.90 to 6.06), 3.70 (95% CI 3.32 to 4.12) and 3.54 (95% CI 3.17 to 3.94), respectively. Assuming HIV/HSV-2 biological interaction, the crude OR ranged between 3.44 and 9.95, with a pooled mean of 8.05 (95% CI 7.14 to 9.07). The adjustments reduced the mean OR to 7.00 (95% CI 6.21 to 7.90), 3.76 (95% CI 3.32 to 4.25) and 3.68 (95% CI 3.25 to 4.17), respectively. Under-reporting of partners reduced the confounder-adjustment effects. Miscapture of high-risk individuals considerably lowered the estimated ORs.

Conclusions

It is difficult to control for sexual-behaviour confounding in observational studies. The observed HIV/HSV-2 association appears more consistent with two infections sharing the same mode of transmission, rather than with HSV-2 enhancing HIV acquisition.

Nonpaternity and Half-Siblingships as Objective Measures of Extramarital Sex: Mathematical Modeling and Simulations

Background

Understanding the epidemiology of HIV and other sexually transmitted infections (STIs) requires knowledge of sexual behavior, but self-reported behavior has limitations. We explored the reliability and validity of nonpaternity and half-siblings ratios as biomarkers of current and past extramarital sex.

Methods

An individual-based Monte Carlo simulation model was constructed to describe partnering and conception in human populations with a focus on Sub-Saharan Africa (SSA). The model was parameterized with representative biological, behavioral, and demographic data.

Results

Nonpaternity and half-siblings ratios were strongly correlated with extramarital sex, with Pearson correlation coefficients (PCC) of 0.79 (95% CI: 0.71–0.86) and 0.77 (0.68–0.84), respectively. Age-specific nonpaternity ratios correlated with past extramarital sex at time of conception for different scenarios: for example, PCC, after smoothing by moving averages, was 0.75 (0.52–0.89) in a scenario of steadily decreasing nonmarital sex and 0.39 (0.01–0.73) in a scenario of transient drops in nonmarital sex. Simulations assuming self-reported levels of extramarital sex from Kenya yielded nonpaternity levels lower than global nonpaternity data, suggesting sizable underreporting of extramarital sex.

Conclusions

Nonpaternity and half-siblings ratios are useful objective measures of extramarital sex that avoid limitations in self-reported sexual behavior.

Sexual network drivers of HIV and herpes simplex virus type 2 transmission

OBJECTIVES

HIV and herpes simplex virus type 2 (HSV-2) infections are sexually transmitted and propagate in sexual networks. Using mathematical modeling, we aimed to quantify effects of key network statistics on infection transmission, and extent to which HSV-2 prevalence can be a proxy of HIV prevalence.

DESIGN/METHODS

An individual-based simulation model was constructed to describe sex partnering and infection transmission, and was parameterized with representative natural history, transmission, and sexual behavior data. Correlations were assessed on model outcomes (HIV/HSV-2 prevalences) and multiple linear regressions were conducted to estimate adjusted associations and effect sizes.

RESULTS

HIV prevalence was one-third or less of HSV-2 prevalence. HIV and HSV-2 prevalences were associated with a Spearman's rank correlation coefficient of 0.64 (95% confidence interval: 0.58-0.69). Collinearities among network statistics were detected, most notably between concurrency versus mean and variance of number of partners. Controlling for confounding, unmarried mean/variance of number of partners (or alternatively concurrency) were the strongest predictors of HIV prevalence. Meanwhile, unmarried/married mean/variance of number of partners (or alternatively concurrency), and clustering coefficient were the strongest predictors of HSV-2 prevalence. HSV-2 prevalence was a strong predictor of HIV prevalence by proxying effects of network statistics.

CONCLUSION

Network statistics produced similar and differential effects on HIV/HSV-2 transmission, and explained most of the variation in HIV and HSV-2 prevalences. HIV prevalence reflected primarily mean and variance of number of partners, but HSV-2 prevalence was affected by a range of network statistics. HSV-2 prevalence (as a proxy) can forecast a population's HIV epidemic potential, thereby informing interventions.

Impact of treatment on hepatitis C virus transmission and incidence in Egypt: A case for treatment as prevention

Egypt has launched a hepatitis C virus (HCV) treatment programme using direct-acting antivirals (DAAs). Our aim was to assess the impact of five plausible programme scale-up and sustainability scenarios for HCV treatment as prevention in Egypt. We developed and analysed a mathematical model to assess programme impact using epidemiologic, programming and health economics measures. The model was parametrized with current and representative natural history, HCV prevalence and programme data. HCV incidence in Egypt is declining, but will persist at a considerable level for decades unless controlled by interventions. Across the five programme scenarios, 1.75-5.60 million treatments were administered by 2030. Reduction in incidence (annual number of new infections) by 2030 ranged between 29% and 99%, programme-attributed reduction in incidence rate (new infections per susceptible person per year) ranged between 18% and 99%, number of infections averted ranged between 42 393 and 469 599, and chronic infection prevalence reached as low as 2.8%-0.1%. Reduction in incidence rate year by year hovered around 7%-15% in the first decade of the programme in most scenarios. Treatment coverage in 2030 ranged between 24.9% and 98.8%, and number of treatments required to avert one new infection ranged between 9.5 and 12.1. Stipulated targets for HCV by 2030 could not be achieved without scaling-up treatment to 365 000 per year and sustaining it for a decade. In conclusion, DAA scale-up will have an immense and immediate impact on HCV incidence in Egypt. Elimination by 2030 is feasible if sufficient resources are committed to programme scale-up and sustainability. HCV treatment as prevention is a potent and effective prevention approach.

Effects of contact structure on the transient evolution of HIV virulence

Early in an epidemic, high densities of susceptible hosts select for relatively high parasite virulence; later in the epidemic, lower susceptible densities select for lower virulence. Thus over the course of a typical epidemic the average virulence of parasite strains increases initially, peaks partway through the epidemic, then declines again. However, precise quantitative outcomes, such as the peak virulence reached and its timing, may depend sensitively on epidemiological details. Fraser et al. proposed a model for the eco-evolutionary dynamics of HIV that incorporates the tradeoffs between transmission and virulence (mediated by set-point viral load, SPVL) and their heritability between hosts. Their model used implicit equations to capture the effects of partnership dynamics that are at the core of epidemics of sexually transmitted diseases. Our models combine HIV virulence tradeoffs with a range of contact models, explicitly modeling partnership formation and dissolution and allowing for individuals to transmit disease outside of partnerships. We assess summary statistics such as the peak virulence (corresponding to the maximum value of population mean log₁₀ SPVL achieved throughout the epidemic) across models for a range of parameters applicable to the HIV epidemic in sub-Saharan Africa. Although virulence trajectories are broadly similar across models, the timing and magnitude of the virulence peak vary considerably. Previously developed implicit models predicted lower virulence and slower progression at the peak (a maximum of 3.5 log₁₀ SPVL) compared both to more realistic models and to simple random-mixing models with no partnership structure at all (both with a maximum of ≈ 4.7 log₁₀ SPVL). In this range of models, the simplest random-mixing structure best approximates the most realistic model; this surprising outcome occurs because the dominance of extra-pair contact in the realistic model swamps the effects of partnership structure.

Pathogens such as HIV can evolve rapidly when the environment changes. One important aspect of a pathogen’s environment is the probability that an infectious contact (a sneeze for a respiratory disease, or an unprotected sex act for a sexually transmitted disease) encounters an uninfected person and thus has a chance to transmit the pathogen. As an epidemic grows the number of uninfected people shrinks, changing evolutionary pressures on the pathogen. While researchers have used models to explore pathogen evolution during epidemics, their models usually neglect important processes such as people entering and leaving sexual partnerships. We compared several evolutionary models for HIV that include partnership dynamics as well as sexual contact outside of stable partnerships. Models of intermediate complexity predicted lower virulence midway through the epidemic (a minimum of 15 years to progress to AIDS) than either more realistic models or simple models with no partnership structure (both with a minimum of 7.25 years to progress to AIDS), because random sexual contacts tended to wash out the effects of stable partnerships. Researchers trying to predict the evolution of pathogens must try to understand the implications of their modeling choices; models of intermediate complexity may not produce intermediate conclusions.

Could Circumcision of HIV-Positive Males Benefit Voluntary Medical Male Circumcision Programs in Africa? Mathematical Modeling Analysis

Background

The epidemiological and programmatic implications of inclusivity of HIV-positive males in voluntary medical male circumcision (VMMC) programs are uncertain. We modeled these implications using Zambia as an illustrative example.

Methods and Findings

We used the Age-Structured Mathematical (ASM) model to evaluate, over an intermediate horizon (2010–2025), the effectiveness (number of VMMCs needed to avert one HIV infection) of VMMC scale-up scenarios with varying proportions of HIV-positive males. The model was calibrated by fitting to HIV prevalence time trend data from 1990 to 2014. We assumed that inclusivity of HIV positive males may benefit VMMC programs by increasing VMMC uptake among higher risk males, or by circumcision reducing HIV male-to-female transmission risk. All analyses were generated assuming no further antiretroviral therapy (ART) scale-up.

The number of VMMCs needed to avert one HIV infection was projected to increase from 12.2 VMMCs per HIV infection averted, in a program that circumcises only HIV-negative males, to 14.0, in a program that includes HIV-positive males. The proportion of HIV-positive males was based on their representation in the population (e.g. 12.6% of those circumcised in 2010 would be HIV-positive based on HIV prevalence among males of 12.6% in 2010). However, if a program that only reaches out to HIV-negative males is associated with 20% lower uptake among higher-risk males, the effectiveness would be 13.2 VMMCs per infection averted. If improved inclusivity of HIV-positive males is associated with 20% higher uptake among higher-risk males, the effectiveness would be 12.4. As the assumed VMMC efficacy against male-to-female HIV transmission was increased from 0% to 20% and 46%, the effectiveness of circumcising regardless of HIV status improved from 14.0 to 11.5 and 9.1, respectively. The reduction in the HIV incidence rate among females increased accordingly, from 24.7% to 34.8% and 50.4%, respectively.

Conclusion

Improving inclusivity of males in VMMC programs regardless of HIV status increases VMMC effectiveness, if there is moderate increase in VMMC uptake among higher-risk males and/or if there is moderate efficacy for VMMC against male-to-female transmission. In these circumstances, VMMC programs can reduce the HIV incidence rate in males by nearly as much as expected by some ART programs, and additionally, females can benefit from the intervention nearly as much as males.

HCV prevalence can predict HIV epidemic potential among people who inject drugs: mathematical modeling analysis

BACKGROUND

Hepatitis C virus (HCV) and HIV are both transmitted through percutaneous exposures among people who inject drugs (PWID). Ecological analyses on global epidemiological data have identified a positive association between HCV and HIV prevalence among PWID. Our objective was to demonstrate how HCV prevalence can be used to predict HIV epidemic potential among PWID.

METHODS

Two population-level models were constructed to simulate the evolution of HCV and HIV epidemics among PWID. The models described HCV and HIV parenteral transmission, and were solved both deterministically and stochastically.

RESULTS

The modeling results provided a good fit to the epidemiological data describing the ecological HCV and HIV association among PWID. HCV was estimated to be eight times more transmissible per shared injection than HIV. A threshold HCV prevalence of 29.0% (95% uncertainty interval (UI): 20.7-39.8) and 46.5% (95% UI: 37.6-56.6) were identified for a sustainable HIV epidemic (HIV prevalence >1%) and concentrated HIV epidemic (HIV prevalence >5%), respectively. The association between HCV and HIV was further described with six dynamical regimes depicting the overlapping epidemiology of the two infections, and was quantified using defined and estimated measures of association. Modeling predictions across a wide range of HCV prevalence indicated overall acceptable precision in predicting HIV prevalence at endemic equilibrium. Modeling predictions were found to be robust with respect to stochasticity and behavioral and biological parameter uncertainty. In an illustrative application of the methodology, the modeling predictions of endemic HIV prevalence in Iran agreed with the scale and time course of the HIV epidemic in this country.

CONCLUSIONS

Our results show that HCV prevalence can be used as a proxy biomarker of HIV epidemic potential among PWID, and that the scale and evolution of HIV epidemic expansion can be predicted with sufficient precision to inform HIV policy, programming, and resource allocation.

Population sexual behavior and HIV prevalence in Sub-Saharan Africa: missing links?

OBJECTIVES

Patterns of sexual partnering should shape HIV transmission in human populations. The objective of this study was to assess empirical associations between population casual sex behavior and HIV prevalence, and between different measures of casual sex behavior.

METHODS

An ecological study design was applied to nationally representative data, those of the Demographic and Health Surveys, in 25 countries of Sub-Saharan Africa. Spearman rank correlation was used to assess different correlations for males and females and their statistical significance.

RESULTS

Correlations between HIV prevalence and means and variances of the number of casual sex partners were positive, but small and statistically insignificant. The majority of correlations across means and variances of the number of casual sex partners were positive, large, and statistically significant. However, all correlations between the means, as well as variances, and the variance of unmarried females were weak and statistically insignificant.

CONCLUSIONS

Population sexual behavior was not predictive of HIV prevalence across these countries. Nevertheless, the strong correlations across means and variances of sexual behavior suggest that self-reported sexual data are self-consistent and convey valid information content. Unmarried female behavior seemed puzzling, but could be playing an influential role in HIV transmission patterns.

Characterizing HIV epidemiology in stable couples in Cambodia, the Dominican Republic, Haiti, and India

Using a set of statistical methods and HIV mathematical models applied on nationally representative Demographic and Health Survey data, we characterized HIV serodiscordancy patterns and HIV transmission dynamics in stable couples (SCs) in four countries: Cambodia, the Dominican Republic, Haiti, and India. The majority of SCs affected by HIV were serodiscordant, and about a third of HIV-infected persons had uninfected partners. Overall, nearly two-thirds of HIV infections occurred in individuals in SCs, but only about half of these infections were due to transmissions within serodiscordant couples. The majority of HIV incidence in the population occurred through extra-partner encounters in SCs. There is similarity in HIV epidemiology in SCs between these countries and countries in sub-Saharan Africa, despite the difference in scale of epidemics. It appears that HIV epidemiology in SCs may share similar patterns globally, possibly because it is a natural ‘spillover’ effect of HIV dynamics in high-risk populations.