SimpactCyan 1.0: An Open-source Simulator for Individual-Based Models in HIV Epidemiology with R and Python Interfaces

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.

[1]  Maurice S. Bartlett,et al.  Stochastic Processes or the Statistics of Change , 1953 .

[2]  Wim Delva,et al.  Connecting the dots: network data and models in HIV epidemiology , 2016, AIDS.

[3]  J J Goedert,et al.  A simple relationship between viral load and survival time in HIV-1 infection. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[4]  S. Tavaré.  Some probabilistic and statistical problems in the analysis of DNA sequences , 1986 .

[5]  David A. Rasmussen,et al.  Phylodynamic Inference for Structured Epidemiological Models , 2014, PLoS Comput. Biol..

[6]  Niel Hens,et al.  Lessons from a decade of individual-based models for infectious disease transmission: a systematic review (2006-2015) , 2017, BMC Infectious Diseases.

[7]  Samuel Alizon,et al.  Within-host and between-host evolutionary rates across the HIV-1 genome , 2013, Retrovirology.

[8]  Chaitra Gopalappa,et al.  Progression and Transmission of HIV/AIDS (PATH 2.0) , 2017, Medical decision making : an international journal of the Society for Medical Decision Making.

[9]  Till Bärnighausen,et al.  Impact of early initiation versus national standard of care of antiretroviral therapy in Swaziland’s public sector health system: study protocol for a stepped-wedge randomized trial , 2017, Trials.

[10]  David F Anderson,et al.  A modified next reaction method for simulating chemical systems with time dependent propensities and delays. , 2007, The Journal of chemical physics.

[11]  John Young,et al.  The ethics of ‘Trials within Cohorts’ (TwiCs): 2nd international symposium , 2017, Trials.

[12]  Sean Lucio Tolentino,et al.  Effective and efficient algorithms for simulating sexually transmitted diseases , 2014 .

[13]  Niema Moshiri,et al.  FAVITES: simultaneous simulation of transmission networks, phylogenetic trees and sequences , 2019, Bioinform..

[14]  Sikhulile Moyo,et al.  Impact of sampling density on the extent of HIV clustering. , 2014, AIDS research and human retroviruses.

[15]  Faryad Darabi Sahneh,et al.  GEMFsim: A stochastic simulator for the generalized epidemic modeling framework , 2016, J. Comput. Sci..

[16]  Fei Meng,et al.  A simulation-based method for efficient resource allocation of combination HIV prevention , 2013, SimuTools.

[17]  Erik M. Volz,et al.  Scalable relaxed clock phylogenetic dating , 2017 .

[18]  Klaus Peter Schliep,et al.  phangorn: phylogenetic analysis in R , 2010, Bioinform..

[19]  Guillaume Deffuant,et al.  Adaptive approximate Bayesian computation for complex models , 2011, Computational Statistics.

[20]  Benjamin Armbruster,et al.  Sexual Networks, Partnership Mixing, and the Female-to-Male Ratio of HIV Infections in Generalized Epidemics: An Agent-Based Simulation Study , 2015 .

[21]  Tanja Stadler,et al.  Phylogenetic Tools for Generalized HIV-1 Epidemics: Findings from the PANGEA-HIV Methods Comparison , 2016, Molecular biology and evolution.

[22]  Leigh F. Johnson,et al.  MicroCOSM: a model of social and structural drivers of HIV and interventions to reduce HIV incidence in high-risk populations in South Africa , 2018, bioRxiv.

[23]  Ye Chen,et al.  Implementation and applications of EMOD, an individual-based multi-disease modeling platform , 2018, Pathogens and disease.

[24]  Ramón Doallo,et al.  CircadiOmics: integrating circadian genomics, transcriptomics, proteomics and metabolomics , 2012, Nature Methods.

[25]  Andrew Rambaut,et al.  Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees , 1997, Comput. Appl. Biosci..

[26]  Marc Lipsitch,et al.  Development, Calibration and Performance of an HIV Transmission Model Incorporating Natural History and Behavioral Patterns: Application in South Africa , 2014, PloS one.

[27]  Deborah Donnell,et al.  Swaziland HIV Incidence Measurement Survey (SHIMS): a prospective national cohort study. , 2017, The lancet. HIV.

[28]  Chris Kenyon,et al.  Higher risk sexual behaviour is associated with unawareness of HIV-positivity and lack of viral suppression – implications for Treatment as Prevention , 2017, Scientific Reports.

[29]  Fernando González-Candelas,et al.  The substitution rate of HIV-1 subtypes: a genomic approach , 2017, Virus evolution.

[30]  Sebastian Bonhoeffer,et al.  Virulence and Pathogenesis of HIV-1 Infection: An Evolutionary Perspective , 2014, Science.

[31]  Larry Sawers,et al.  HIV and concurrent sexual partnerships: modelling the role of coital dilution , 2011, Journal of the International AIDS Society.

[32]  T. Déirdre Hollingsworth,et al.  Variation in HIV-1 set-point viral load: Epidemiological analysis and an evolutionary hypothesis , 2007, Proceedings of the National Academy of Sciences.

[33]  Jan C. Thiele.  R Marries NetLogo: Introduction to the RNetLogo Package , 2014 .

[34]  John Hargrove,et al.  How Should We Best Estimate the Mean Recency Duration for the BED Method? , 2012, PloS one.

[35]  Christophe Fraser,et al.  HIV-1 transmission, by stage of infection. , 2008, The Journal of infectious diseases.

[36]  Roel Bakker,et al.  STDSIM: A Microsimulation Model for Decision Support in the Control of HIV and Other STDs , 2000 .

[37]  Martina Morris,et al.  EpiModel: An R Package for Mathematical Modeling of Infectious Disease over Networks , 2017, bioRxiv.

[38]  Viktor Müller,et al.  Comparison of two individual-based models for simulating HIV epidemics in a population with HSV-2 using as case study Yaoundé-Cameroon, 1989-1998 , 2019, bioRxiv.

[39]  Alex Welte,et al.  Coital frequency and condom use in monogamous and concurrent sexual relationships in Cape Town, South Africa , 2013, Journal of the International AIDS Society.

[40]  Lyle McKinnon,et al.  Genital immunology and HIV susceptibility in young women , 2013, American journal of reproductive immunology.

[41]  Tulio de Oliveira,et al.  Transmission networks and risk of HIV infection in KwaZulu-Natal, South Africa: a community-wide phylogenetic study. , 2017, The lancet. HIV.

[42]  Huldrych F. Günthard,et al.  Using an Epidemiological Model for Phylogenetic Inference Reveals Density Dependence in HIV Transmission , 2013, Molecular biology and evolution.

[43]  Matthew H Todd,et al.  Open science is a research accelerator. , 2011, Nature chemistry.

[44]  Tanja Stadler,et al.  Assessment of Overlap of Phylogenetic Transmission Clusters and Communities in Simple Sexual Contact Networks: Applications to HIV-1 , 2016, PloS one.

[45]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[46]  Wim Delva,et al.  Calibration of individual-based models to epidemiological data: A systematic review , 2019, PLoS Comput. Biol..

[47]  J. Doob.  Topics in the theory of Markoff chains , 1942 .

[48]  Richard M Grimes,et al.  Re-engagement in HIV Care: A Clinical and Public Health Priority. , 2016, Journal of AIDS & clinical research.

[49]  Jonathan M. Mudge,et al.  Evidence for Transcript Networks Composed of Chimeric RNAs in Human Cells , 2012, PloS one.

[50]  Steven F. Railsback,et al.  Agent-Based and Individual-Based Modeling: A Practical Introduction , 2011 .

[51]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[52]  A. Rambaut,et al.  Episodic Sexual Transmission of HIV Revealed by Molecular Phylodynamics , 2008, PLoS medicine.

[53]  Kalliopi Kravari,et al.  A Survey of Agent Platforms , 2015, J. Artif. Soc. Soc. Simul..