Developing a temperature-driven map of the basic reproductive number of the emerging tick vector of Lyme disease Ixodes scapularis in Canada.

A mechanistic model of the tick vector of Lyme disease, Ixodes scapularis, was adapted to a deterministic structure. Using temperature normals smoothed by Fourier analysis to generate seasonal temperature-driven development rates and host biting rates, and a next generation matrix approach, the model was used to obtain values for the basic reproduction number (R(0)) for I. scapularis at locations in southern Canada where the tick is established and emerging. The R(0) at Long Point, Point Pelee and Chatham sites where I. scapularis are established, was estimated at 1.5, 3.19 and 3.65, respectively. The threshold temperature conditions for tick population survival (R(0)=1) were shown to be the same as those identified using the mechanistic model (2800-3100 cumulative annual degree days >0°C), and a map of R(0) for I. scapularis, the first such map for an arthropod vector, was drawn for Canada east of the Rocky Mountains. This map supports current risk assessments for Lyme disease risk emergence in Canada. Sensitivity analysis identified host abundance, tick development rates and summer temperatures as highly influential variables in the model, which is consistent with our current knowledge of the biology of this tick. The development of a deterministic model for I. scapularis that is capable of providing values for R(0) is a key step in our evolving ability to develop tools for assessment of Lyme disease risk emergence and for development of public health policies on surveillance, prevention and control.

[1]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[2]  A. Pugliese,et al.  Effects of Tick Population Dynamics and Host Densities on the Persistence of Tick-borne Infections Article in Press , 2006 .

[3]  Lisette de Pillis,et al.  Mathematical Biology at an Undergraduate Liberal Arts College , 2010, CBE life sciences education.

[4]  M. Begon,et al.  Persistence of tick-borne virus in the presence of multiple host species: tick reservoirs and parasite mediated competition. , 1999, Journal of theoretical biology.

[5]  H. Artsob,et al.  The rising challenge of Lyme borreliosis in Canada. , 2008, Canada communicable disease report = Releve des maladies transmissibles au Canada.

[6]  D. Rogers,et al.  Incidence from coincidence: patterns of tick infestations on rodents facilitate transmission of tick-borne encephalitis virus , 1999, Parasitology.

[7]  S W Lindsay,et al.  Climate change and vector-borne diseases: a regional analysis. , 2000, Bulletin of the World Health Organization.

[8]  L. R. Lindsay,et al.  Predicting the speed of tick invasion: an empirical model of range expansion for the Lyme disease vector Ixodes scapularis in Canada , 2012 .

[9]  R. May,et al.  Infectious Diseases of Humans: Dynamics and Control , 1991, Annals of Internal Medicine.

[10]  J. Řeháček,et al.  Persistence of Tick-Borne Encephalitis Virus in Hibernating Bats. , 1961 .

[11]  P. Hudson,et al.  Thresholds for disease persistence in models for tick-borne infections including non-viraemic transmission, extended feeding and tick aggregation. , 2003, Journal of theoretical biology.

[12]  D. Waltner-Toews,et al.  Investigation of Relationships Between Temperature and Developmental Rates of Tick Ixodes scapularis (Acari: Ixodidae) in the Laboratory and Field , 2004, Journal of medical entomology.

[13]  M. Lankester,et al.  Spring migratory birds (Aves) extend the northern occurrence of blacklegged tick (Acari:Ixodidae). , 1996, Journal of medical entomology.

[14]  O. Diekmann,et al.  Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation , 2000 .

[15]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[16]  G. Beauchamp,et al.  Investigation of Ground Level and Remote-Sensed Data for Habitat Classification and Prediction of Survival of Ixodes scapularis in Habitats of Southeastern Canada , 2006, Journal of medical entomology.

[17]  J A P Heesterbeek,et al.  The Basic Reproduction Number for Complex Disease Systems: Defining R0 for Tick‐Borne Infections , 2008, The American Naturalist.

[18]  D J Rogers,et al.  Climate change and vector-borne diseases. , 2006, Advances in parasitology.

[19]  M Bigras-Poulin,et al.  A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. , 2005, International journal for parasitology.

[20]  T. Gillespie,et al.  Survival and development of the different life stages of Ixodes scapularis (Acari: Ixodidae) held within four habitats on Long Point, Ontario, Canada. , 1998, Journal of medical entomology.

[21]  T. Porco A mathematical model of the ecology of Lyme disease. , 1999, IMA journal of mathematics applied in medicine and biology.

[22]  M. E. Alexander,et al.  Modelling strategies for controlling SARS outbreaks , 2004, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[23]  Thomas Caraco,et al.  Stage‐Structured Infection Transmission and a Spatial Epidemic: A Model for Lyme Disease , 2002, The American Naturalist.

[24]  R. Ostfeld,et al.  BIODIVERSITY AND THE DILUTION EFFECT IN DISEASE ECOLOGY , 2001 .

[25]  Nicholas H. Ogden,et al.  Active and Passive Surveillance and Phylogenetic Analysis of Borrelia burgdorferi Elucidate the Process of Lyme Disease Risk Emergence in Canada , 2010, Environmental health perspectives.

[26]  A. Pugliese,et al.  Seasonal population dynamics of ticks, and its influence on infection transmission: A semi-discrete approach , 2004, Bulletin of mathematical biology.

[27]  E. Lyons,et al.  Pandemic Potential of a Strain of Influenza A (H1N1): Early Findings , 2009, Science.

[28]  S. Randolph,et al.  A modified matrix model to describe the seasonal population ecology of the European tick Ixodes ricinus , 2011 .

[29]  Xiao-Qiang Zhao,et al.  Threshold Dynamics for Compartmental Epidemic Models in Periodic Environments , 2008 .

[30]  Martijn Gough Climate change , 2009, Canadian Medical Association Journal.

[31]  D. Cummings,et al.  Strategies for containing an emerging influenza pandemic in Southeast Asia , 2005, Nature.

[32]  E. Daniels,et al.  Simulation of blacklegged tick (Acari:Ixodidae) population dynamics and transmission of Borrelia burgdorferi. , 1997, Journal of medical entomology.

[33]  B. Grenfell Revisiting the basic reproductive number for malaria and its implications for malaria control , 2007 .

[34]  M. Blaser,et al.  Critical role of nosocomial transmission in the toronto sars outbreak. , 2004, Mathematical biosciences and engineering : MBE.

[35]  D. Haile,et al.  Computer simulation of population dynamics of the American dog tick (Acari: Ixodidae). , 1989, Journal of medical entomology.

[36]  M. Bigras-Poulin,et al.  Climate change and the potential for range expansion of the Lyme disease vector Ixodes scapularis in Canada. , 2006, International journal for parasitology.

[37]  D. Dennis,et al.  Reported distribution of Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae) in the United States. , 1998, Journal of medical entomology.

[38]  G. Yohe,et al.  A globally coherent fingerprint of climate change impacts across natural systems , 2003, Nature.

[39]  S. Mahamdallie,et al.  Integrated Mapping of Establishment Risk for Emerging Vector-Borne Infections: A Case Study of Canine Leishmaniasis in Southwest France , 2011, PloS one.

[40]  M. Bigras-Poulin,et al.  Role of Migratory Birds in Introduction and Range Expansion of Ixodes scapularis Ticks and of Borrelia burgdorferi and Anaplasma phagocytophilum in Canada , 2008, Applied and Environmental Microbiology.

[41]  R Meiswinkel,et al.  Mapping the basic reproduction number (R₀) for vector-borne diseases: a case study on bluetongue virus. , 2009, Epidemics.

[42]  David L Smith,et al.  Revisiting the Basic Reproductive Number for Malaria and Its Implications for Malaria Control , 2007, PLoS biology.

[43]  A. Steere,et al.  The emergence of Lyme disease. , 2004, The Journal of clinical investigation.

[44]  A. Dobson Population Dynamics of Pathogens with Multiple Host Species , 2004, The American Naturalist.

[45]  M Bigras-Poulin,et al.  Vector seasonality, host infection dynamics and fitness of pathogens transmitted by the tick Ixodes scapularis , 2006, Parasitology.

[46]  P. Mead,et al.  Surveillance for Lyme disease--United States, 1992-2006. , 2008, Morbidity and mortality weekly report. Surveillance summaries.

[47]  M. Bigras-Poulin,et al.  Risk maps for range expansion of the Lyme disease vector, Ixodes scapularis, in Canada now and with climate change , 2008, International journal of health geographics.

[48]  Nicolas Bacaër,et al.  The epidemic threshold of vector-borne diseases with seasonality , 2006, Journal of mathematical biology.

[49]  Thomas Smith,et al.  A Research Agenda for Malaria Eradication: Modeling , 2011, PLoS medicine.

[50]  D. Kirschner,et al.  A methodology for performing global uncertainty and sensitivity analysis in systems biology. , 2008, Journal of theoretical biology.

[51]  I. Foppa The basic reproductive number of tick–borne encephalitis virus , 2005, Journal of mathematical biology.

[52]  Nicholas H. Ogden,et al.  The emergence of Lyme disease in Canada , 2009, Canadian Medical Association Journal.