SimpactCyan is an open-source simulator for individual-based models in HIV epidemiology. Its core algorithm is written in C++ for computational efficiency, while the R and Python interfaces aim to make the tool accessible to the fast-growing community of R and Python users. Transmission, treatment and prevention of HIV infections in dynamic sexual networks are simulated by discrete events. A generic “intervention” event allows model parameters to be changed over time, and can be used to model medical and behavioural HIV prevention programmes. First, we describe a more efficient variant of the modified Next Reaction Method that drives our continuous-time simulator. Next, we outline key built-in features and assumptions of individual-based models formulated in SimpactCyan, and provide code snippets for how to formulate, execute and analyse models in SimpactCyan through its R and Python interfaces. Lastly, we give two examples of applications in HIV epidemiology: the first demonstrates how the software can be used to estimate the impact of progressive changes to the eligibility criteria for HIV treatment on HIV incidence. The second example illustrates the use of SimpactCyan as a data-generating tool for assessing the performance of a phylodynamic inference framework.

In Burundi, malaria infection has been increasing in the last decade despite efforts to increase access to health services, and several intervention programs. The use of heterogeneous data can help to build predictive models of malaria cases. We built predictive frameworks: the generalized linear model (GLM), and artificial neural network (ANN), to predict malaria cases in four sub-groups and the overall general population. Descriptive results showed that more than half of malaria infections are observed in pregnant women and children under 5 years, with high burden to children between 12 and 59 months. Modelling results showed that, ANN model performed better in predicting total cases compared to GLM. Both model frameworks showed that education rates and Insecticide Treated Bed Nets (ITNs) had decreasing effects on malaria cases, some other variables had an increasing effect. Thus, malaria control and prevention interventions program are encouraged to understand those variables, and take appropriate measures such as providing ITNs, sensitization in schools and the communities, starting within high dense communities, among others. Early prediction of cases can provide timely information needed to be proactive for intervention strategies, and it can help to mitigate the epidemics and reduce its impact on populations and the economy.

Initial investigation into the emerging mpox outbreak of novel clade Ib in eastern Democratic Republic of the Congo has identified signs of sustained human-to-human transmission and epidemiological links to sexual contacts involving female sex workers (FSWs) ¹ , which have not been observed in previous clade Ia outbreaks. Using mathematical models incorporating age-dependent contact patterns, we quantified the role of frequent sexual interactions as opposed to community contacts in clade Ib’s dynamics and found that this additional mode of transmission could explain its increased outbreak potential compared with clade Ia. As with the globally-circulating clade IIb, transmitted predominantly among men who have sex with men ² , our findings reinforce the importance of protecting key population groups, specifically FSWs for clade Ib, in controlling ongoing mpox outbreaks.

(1) Background: Calibration of Simpact Cyan can help to improve estimates related to the transmission dynamics of the Human Immunodeficiency Virus (HIV). Age-mixing patterns in sexual partnerships, onward transmissions, and temporal trends of HIV incidence are determinants which can inform the design of efficient prevention, and linkage-to-care programs. Using an agent-based model (ABM) simulation tool, we investigated, through a simulation study, if estimates of these determinants can be obtained with high accuracy by combining summary features from different data sources. (2) Methods: With specific parameters, we generated the benchmark data, and calibrated the default model in three scenarios based on summary features for comparison. For calibration, we used Latin Hypercube Sampling approach to generate parameter values, and Approximation Bayesian Computation to choose the best fitting ones. In all calibration scenarios the mean square root error was used as a measure to depict the estimates accuracy. (3) Results: The accuracy measure showed relatively no difference between the three scenarios. Moreover, we found that in all scenarios, age and gender strata incidence trends were poorly estimated. (4) Conclusions: Using synthetic benchmarks, we showed that it is possible to infer HIV transmission dynamics using an ABM of HIV transmission. Our results suggest that any type of summary feature provides adequate information to estimate HIV transmission network determinants. However, it is advisable to check the level of accuracy of the estimates of interest using benchmark data.

Understanding age-mixing patterns in Human Immunodeficiency Virus (HIV) transmission networks can enhance the design and implementation of HIV prevention strategies in sub-Saharan Africa. Due to ethical consideration, it is less likely possible to conduct a benchmark study to assess which sampling strategy, and sub-optimal sampling coverage which can yield best estimates for these patterns. We conducted a simulation study, using phylogenetic trees to infer estimates of age-mixing patterns in HIV transmission, through the computation of proportions of pairings between men and women, who were phylogenetically linked across different age groups (15–24 years, 25–39 years, and 40–49 years); and the means, and standard deviations of their age difference. We investigated also the uncertainty around these estimates as a function of the sampling coverage in four sampling strategies: when missing sequence data were missing completely at random (MCAR), and missing at random (MAR) with at most 30%—50%—70% of women in different age groups being in the sample. The results suggested that age-mixing patterns in HIV transmission can be unveiled from proportions of phylogenetic pairings between men and women across age groups; and the mean, and standard deviation of their age difference. A 55% sampling coverage was sufficient to provide the best values of estimates of age-mixing patterns in HIV transmission with MCAR scenario. But we should be cautious in interpreting proportions of men phylogenetically linked to women because they may be overestimated or underestimated, even at higher sampling coverage. The findings showed that, MCAR was the best sampling strategy. This means, it is advisable not to use sequence data collected in settings where we can find a systematic imbalance of age and gender to investigate age-mixing in HIV transmission. If not possible, ensure to take into consideration the imbalance in interpreting the results.

In 2024, mpox cases surged in the Democratic Republic of the Congo (DRC) with cross-border spread to Burundi. We developed a transmission-dynamic model calibrated against surveillance data to understand drivers in enzootic (Clade Ia) and non-enzootic (Clade Ib) areas, and the potential impact of vaccination. In non-enzootic areas we estimated that 58-84% of transmission occurred within sexual networks. MVA-BN vaccination of sex workers could have averted 91% (95% CrI 81%–98%) of infections in Sud Kivu (DRC) but only 35% (95% CrI 26%–47%) in Bujumbura (Burundi), due to later outbreak detection. In historically enzootic Equateur (DRC), ongoing zoonotic spillover best explained sustained incidence. There, pledged Lc18m8 vaccines could have averted 42% (95% CrI 40%–46%) of infections; prioritising children improved impact. Across all settings, doubling vaccine coverage by using a single dose of MVA-BN outperformed two-dose strategies. Timely detection and tailored vaccination strategies are critical to reducing mpox burden.