MGDrivE 2: A simulation framework for gene drive systems incorporating seasonality and epidemiological dynamics

Interest in gene drive technology has continued to grow as promising new drive systems have been developed in the lab and discussions are moving towards implementing field trials. The prospect of field trials requires models that incorporate a significant degree of ecological detail, including parameters that change over time in response to environmental data such as temperature and rainfall, leading to seasonal patterns in mosquito population density. Epidemiological outcomes are also of growing importance, as: i) the suitability of a gene drive construct for release will depend on its expected impact on disease transmission, and ii) initial field trials are expected to have a measured entomological outcome and a modeled epidemiological outcome. We present MGDrivE 2 (Mosquito Gene Drive Explorer 2): an extension of and development from the MGDrivE 1 simulation framework that investigates the population dynamics of a variety of gene drive architectures and their spread through spatially-explicit mosquito populations. Key strengths and improvements of the MGDrivE 2 framework are: i) the ability of parameters to vary with time and induce seasonal population dynamics, ii) an epidemiological module accommodating reciprocal pathogen transmission between humans and mosquitoes, and iii) an implementation framework based on stochastic Petri nets that enables efficient model formulation and flexible implementation. Example MGDrivE 2 simulations are presented to demonstrate the application of the framework to a CRISPR-based homing gene drive system intended to drive a disease-refractory gene into a population, incorporating time-varying temperature and rainfall data, and predict impact on human disease incidence and prevalence. Further documentation and use examples are provided in vignettes at the project’s CRAN repository. MGDrivE 2 is an open-source R package freely available on CRAN. We intend the package to provide a flexible tool capable of modeling gene drive constructs as they move closer to field application and to infer their expected impact on disease transmission.

[1]  P. Glynn.  A GSMP formalism for discrete event systems , 1989, Proc. IEEE.

[2]  Corrado Priami,et al.  Simulation Algorithms for Computational Systems Biology , 2017, Texts in Theoretical Computer Science. An EATCS Series.

[3]  Andrew G. Clark,et al.  A toxin-antidote CRISPR gene drive system for regional population modification , 2019, bioRxiv.

[4]  Peter G. Fennell,et al.  Limitations of discrete-time approaches to continuous-time contagion dynamics , 2016, Physical review. E.

[5]  Penny A Hancock,et al.  Malaria Journal BioMed Central Methodology , 2007 .

[6]  T. Kurtz.  Solutions of ordinary differential equations as limits of pure jump markov processes , 1970, Journal of Applied Probability.

[7]  Darren J. Wilkinson.  Stochastic Modelling for Systems Biology , 2006 .

[8]  Pejman Rohani,et al.  Ecological and immunological determinants of dengue epidemics. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Austin Burt,et al.  Requirements for effective malaria control with homing endonuclease genes , 2011, Proceedings of the National Academy of Sciences.

[10]  David L Smith,et al.  Quantitative, model-based estimates of variability in the generation and serial intervals of Plasmodium falciparum malaria , 2016, Malaria Journal.

[11]  Junru Liu,et al.  Transforming insect population control with precision guided sterile males with demonstration in flies , 2019, Nature Communications.

[12]  Martin A. Nowak,et al.  Daisy-chain gene drives for the alteration of local populations , 2016, Proceedings of the National Academy of Sciences.

[13]  John M. Marshall,et al.  Transforming insect population control with precision guided sterile males with demonstration in flies , 2018, Nature Communications.

[14]  E. Seneta,et al.  On Quasi-Stationary distributions in absorbing discrete-time finite Markov chains , 1965, Journal of Applied Probability.

[15]  Neil M Ferguson,et al.  Modelling the impact of vector control interventions on Anopheles gambiae population dynamics , 2011, Parasites & Vectors.

[16]  P. Gething,et al.  malariaAtlas: an R interface to global malariometric data hosted by the Malaria Atlas Project , 2018, Malaria Journal.

[17]  John M. Marshall,et al.  Toward the Definition of Efficacy and Safety Criteria for Advancing Gene Drive-Modified Mosquitoes to Field Testing , 2020, Vector borne and zoonotic diseases.

[18]  T. Scott,et al.  House-to-house human movement drives dengue virus transmission , 2012, Proceedings of the National Academy of Sciences.

[19]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[20]  A. James,et al.  Next-generation gene drive for population modification of the malaria vector mosquito, Anopheles gambiae , 2020, Proceedings of the National Academy of Sciences.

[21]  John M. Marshall,et al.  Overcoming evolved resistance to population-suppressing homing-based gene drives , 2016, Scientific Reports.

[22]  John C. Carlson,et al.  A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission , 2004, Malaria Journal.

[23]  Christian P. Robert,et al.  Stochastic Modelling for Systems Biology (second edition) , 2012 .

[24]  David L Smith,et al.  The Risk of a Mosquito-Borne Infectionin a Heterogeneous Environment , 2004, PLoS biology.

[25]  Darren J. Wilkinson,et al.  Stochastic Modelling for Systems Biology, Third Edition , 2018 .

[26]  J. Crawford,et al.  The distribution of hatching time in Anopheles gambiae , 2006, Malaria Journal.

[27]  Daniel T. Gillespie,et al.  The multivariate Langevin and Fokker–Planck equations , 1996 .

[28]  Héctor M. Sánchez C.,et al.  MGDrivE: A modular simulation framework for the spread of gene drives through spatially-explicit mosquito populations , 2018, bioRxiv.

[29]  K. Paaijmans,et al.  The importance of temperature fluctuations in understanding mosquito population dynamics and malaria risk , 2017, Royal Society Open Science.

[30]  E. Seneta,et al.  On Quasi-Stationary distributions in absorbing discrete-time finite Markov chains , 1965, Journal of Applied Probability.

[31]  Wes Hinsley,et al.  Reducing Plasmodium falciparum Malaria Transmission in Africa: A Model-Based Evaluation of Intervention Strategies , 2010, PLoS medicine.

[32]  Arno Solin,et al.  Applied Stochastic Differential Equations , 2019 .

[33]  Douglas M. Bates,et al.  Sparse and Dense Matrix Classes and Methods [R package Matrix version 1.2-18] , 2019 .

[34]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[35]  Diego Latella,et al.  Continuous approximation of collective system behaviour: A tutorial , 2013, Perform. Evaluation.

[36]  Corrado Priami,et al.  Incorporating extrinsic noise into the stochastic simulation of biochemical reactions: A comparison of approaches. , 2018, The Journal of chemical physics.

[37]  John Odentrantz,et al.  Markov Chains: Gibbs Fields, Monte Carlo Simulation, and Queues , 2000, Technometrics.

[38]  Jared B. Bennett,et al.  Efficient population modification gene-drive rescue system in the malaria mosquito Anopheles stephensi , 2020, Nature Communications.

[39]  Scott Meyers,et al.  Effective modern C++: 42 specific ways to improve your use of C++11 and C++14 , 2014 .

[40]  J. Kingman.  Markov population processes , 1969, Journal of Applied Probability.

[41]  Phillip Cassey,et al.  Dodging silver bullets: good CRISPR gene-drive design is critical for eradicating exotic vertebrates , 2017, Proceedings of the Royal Society B: Biological Sciences.

[42]  E. Seneta,et al.  On quasi-stationary distributions in absorbing continuous-time finite Markov chains , 1967, Journal of Applied Probability.

[43]  Raúl Toral,et al.  Stochastic Numerical Methods: An Introduction for Students and Scientists , 2014 .

[44]  Karline Soetaert,et al.  Solving Ordinary Differential Equations in R , 2012 .

[45]  Sebastian Funk,et al.  Transmission Dynamics of Zika Virus in Island Populations: A Modelling Analysis of the 2013–14 French Polynesia Outbreak , 2016, bioRxiv.

[46]  Peter Buchholz,et al.  Input Modeling with Phase-Type Distributions and Markov Models: Theory and Applications , 2014 .

[47]  P J Goss,et al.  Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[48]  D. Damiens,et al.  Field evaluation of seasonal trends in relative population sizes and dispersal pattern of Aedes albopictus males in support of the design of a sterile male release strategy , 2019, Parasites & Vectors.

[49]  J. Botella de Maglia,et al.  [Prevention of malaria]. , 1999, Revista clinica espanola.

[50]  Maia Martcheva,et al.  An Introduction to Mathematical Epidemiology , 2015 .

[51]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[52]  Bruce A. Hay,et al.  Cleave and Rescue, a novel selfish genetic element and general strategy for gene drive , 2019, Proceedings of the National Academy of Sciences.

[53]  Jared B. Bennett,et al.  Active Genetic Neutralizing Elements for Halting or Deleting Gene Drives. , 2020, Molecular cell.

[54]  Peter J. Haas,et al.  Stochastic Petri Nets: Modelling, Stability, Simulation , 2002 .

[55]  Anja Gronewold,et al.  Event-based modelling of ecological systems with asynchronous cellular automata , 1998 .

[56]  A. Burt,et al.  A CRISPR–Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes , 2018, Nature Biotechnology.

[57]  Paul Bratley,et al.  A guide to simulation , 1983 .

[58]  Corrado Priami,et al.  Hybrid Simulation Algorithms , 2017 .

[59]  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.

[60]  Andrew G. Clark,et al.  Reducing resistance allele formation in CRISPR gene drive , 2018, Proceedings of the National Academy of Sciences.

[61]  Fredros O. Okumu,et al.  Pathway to Deployment of Gene Drive Mosquitoes as a Potential Biocontrol Tool for Elimination of Malaria in Sub-Saharan Africa: Recommendations of a Scientific Working Group , 2018, The American journal of tropical medicine and hygiene.

[62]  Paul J. Hurtado,et al.  Generalizations of the ‘Linear Chain Trick’: incorporating more flexible dwell time distributions into mean field ODE models , 2018, Journal of Mathematical Biology.

[63]  David L Smith,et al.  Statics and dynamics of malaria infection in Anopheles mosquitoes , 2004, Malaria Journal.

[64]  Jared B. Bennett,et al.  Development of a confinable gene drive system in the human disease vector Aedes aegypti , 2020, eLife.

[65]  E. Bier,et al.  The dawn of active genetics , 2016, BioEssays : news and reviews in molecular, cellular and developmental biology.

[66]  S W Lindsay,et al.  Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae) , 2003, Bulletin of Entomological Research.

[67]  L. Allen.  Some discrete-time SI, SIR, and SIS epidemic models. , 1994, Mathematical biosciences.

[68]  Norbert Perrimon,et al.  Temperature-sensitive control of protein activity by conditionally splicing inteins , 2004, Nature Biotechnology.

[69]  C. Mbogo,et al.  Larval density dependence in Anopheles gambiae s.s., the major African vector of malaria. , 2013, The Journal of animal ecology.

[70]  Karline Soetaert,et al.  Solving Differential Equations in R: Package deSolve , 2010 .

[71]  S. Lal,et al.  Epidemiology and control of malaria , 1999, Indian journal of pediatrics.

[72]  L. Johnson,et al.  Thermal biology of mosquito‐borne disease , 2019, Ecology letters.

[73]  X. Rodó,et al.  A Model for a Chikungunya Outbreak in a Rural Cambodian Setting: Implications for Disease Control in Uninfected Areas , 2014, PLoS neglected tropical diseases.