Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2901884

Formats

Article sections

Authors

Related links

Clin Microbiol Infect. Author manuscript; available in PMC 2010 December 1.

Published in final edited form as:

PMCID: PMC2901884

NIHMSID: NIHMS104983

Corresponding author: H. Gaff, Virginia Modelling, Analysis and Simulation Center School of Community and Environmental Health, 3133A Technology Building, Old Dominion University, Norfolk, VA 23529, USA E-mail: ude.udo@ffagh

The publisher's final edited version of this article is available at Clin Microbiol Infect

Human monocytic ehrlichiosis (Ehrlichia chaffeensis), HME, is a tick-transmitted, rickettsial disease that has recently increased substantially in the USA from 142 reported cases in 2001 to 506 reported cases in 2005 [1,2]. There have been increasing surveys of tick populations over the past 10 years that have in turn supported the development of models for tick-borne disease transmission. Resulting HME models [3] suggest the importance of metapopulation structures, landscape environment parameters and periodic climatic effects in predicting the dynamics of HME transmission and the efficacy of control efforts, such as the reduction of the tick population through acaricide use. On this note, we describe a spatially-explicit model for HME

$$\begin{array}{cc}\hfill \frac{d{N}_{i}}{dt}& ={\beta}_{i}\left(\frac{{K}_{i}-{N}_{i}}{{K}_{i}}\right){N}_{i}-{b}_{i}{N}_{i}+\sum _{j}{m}_{ij}({N}_{j}-{N}_{i})\hfill \\ \hfill \frac{d{V}_{i}}{dt}& ={\widehat{\beta}}_{i}{V}_{i}\left(\frac{{M}_{i}{N}_{i}-{V}_{i}}{{M}_{i}{N}_{i}}\right)-\left({\widehat{b}}_{i}+{\delta}_{i}\right){V}_{i}+\sum _{j}{m}_{ij}({V}_{j}-{V}_{i})\hfill \\ \hfill \frac{d{Y}_{i}}{dt}& ={A}_{i}\left(\frac{{N}_{i}-{Y}_{i}}{{N}_{i}}\right){X}_{i}-{\beta}_{i}\frac{{N}_{i}{Y}_{i}}{{K}_{i}}-{b}_{i}{Y}_{i}+\sum _{{N}_{j}<{N}_{i}}{m}_{ij}\frac{{Y}_{i}}{{N}_{i}}({N}_{j}-{N}_{i})+\sum _{{N}_{j}>{N}_{i}}{m}_{ij}\frac{{Y}_{j}}{{N}_{j}}({N}_{j}-{N}_{i})\hfill \\ \hfill \frac{d{X}_{i}}{dt}& ={\widehat{A}}_{i}\left(\frac{{Y}_{i}}{{N}_{i}}\right)({V}_{i}-{X}_{i})-{\widehat{\beta}}_{i}\frac{{V}_{i}{X}_{i}}{{M}_{i}{N}_{i}}-\left({\widehat{b}}_{i}+{\delta}_{i}\right){X}_{i}+\sum _{{V}_{j}<{V}_{i}}{m}_{ij}\frac{{X}_{i}}{{V}_{i}}({V}_{j}-{V}_{i})+\sum _{{V}_{j}>{V}_{i}}{m}_{ij}\frac{{X}_{j}}{{V}_{j}}({V}_{j}-{V}_{i})\hfill \end{array}$$

transmission, and give a result illustrating the importance of migration in the dynamics of HME risk.

Our initial model is a simplified representation of the tick-host-HME system using a single host population, a single life stage and a single pathogen. The dynamics are modelled using differential equations on a spatial grid consisting of environmental patches. White-tailed deer (Odocoileus virginianus) serve as our generic host population and lone-star ticks (Amblyomma americanum) as the tick population. The population densities of the host and tick population in each patch are denoted by *N _{i}* and

The first equation above describes the population dynamics of the host population in patch *i*, for which there is logistic growth with carrying capacity *K _{i}* as well as an external death rate,

In initial simulations, the model was allowed to run for 1000 months on a single patch. We found that with the absence of control, HME prevalence reached between 5 and 15% infected for ticks with hosts at about 35% infected. With the application of acaricide from month 500 forward, the tick population is significantly reduced, and HME is eliminated from both populations in the model. If the acaracide is only applied for 3 years, we find that the tick population recovers to pretreatment levels within a year, while HME doesn't return until nearly 10 years later.

We then expanded the model with multiple patches to include environmental differences, including external death rates that vary due to exposure or predation. Our initial work has two habitat classifications: grass and woods. Each geographic area is linked to others using migration that depends on distance between the two areas and existence of any barriers that would prevent movement between them. The figure shows the effect of migration rates on the average percentage of ticks that are infected in the two patches, one where HME would be endemic in isolation (solid blue) and one where it would die out in isolation (dotted red) (Fig. 1). These results show that even a small amount of connection to an endemic area will quickly increase the prevalence of HME in an area that would be free of disease in isolation.

Average percentage of infected ticks in the endemic (solid blue) and disease-free (dotted red) patches after 1000 months. Migration rates must be very low in order for the individual patches to retain their natural disease dynamics.

Future work on this project will further explore the importance of metapopulation considerations in modelling of rickettsial disease. In particular, we consider the importance of environmental variability, host diversity, habitat fragmentation and targeted application of acaricide.

H. Gaff is supported by Grant Number K25AI067791 from the National Institute of Allergy and Infectious Diseases. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Allergy and Infectious Diseases or the National Institutes of Health.

No conflicts of interest declared.

1. Anonymous Centers for Disease Control and Prevention, Summary of Notifiable Diseases – United States, 2005. MMWR. 2007 March 30;54(53):2–92. [PubMed]

2. Anonymous Centers for Disease Control and Prevention, Summary of Notifiable Diseases – United States, 2001. MMWR. 2003 May 2;50(53):1–108. [PubMed]

3. Gaff H, Gross LJ. Analysis of a tick-borne disease model with varying population sizes in various habitats, Bulletin of Mathematical Biology. 2007;69:265–288. [PubMed]