Lyme Disease Models of Tick-Mouse Dynamics with Seasonal Variation in Births, Deaths, and Tick Feeding

Lyme disease is the most common vector-borne disease in the United States impacting the Northeast and Midwest at the highest rates. Recently, it has become established in southeastern and south-central regions of Canada. In these regions, Lyme disease is caused by Borrelia burgdorferi, which is transmitted to humans by an infected Ixodes scapularis tick. Understanding the parasite-host interaction is critical as the white-footed mouse is one of the most competent reservoir for B. burgdorferi. The cycle of infection is driven by tick larvae feeding on infected mice that molt into infected nymphs and then transmit the disease to another susceptible host such as mice or humans. Lyme disease in humans is generally caused by the bite of an infected nymph. The main aim of this investigation is to study how diapause delays and demographic and seasonal variability in tick births, deaths, and feedings impact the infection dynamics of the tick-mouse cycle. We model tick-mouse dynamics with fixed diapause delays and more realistic Erlang distributed delays through delay and ordinary differential equations (ODEs). To account for demographic and seasonal variability, the ODEs are generalized to a continuous-time Markov chain (CTMC). The basic reproduction number and parameter sensitivity analysis are computed for the ODEs. The CTMC is used to investigate the probability of Lyme disease emergence when ticks and mice are introduced, a few of which are infected. The probability of disease emergence is highly dependent on the time and the infected species introduced. Infected mice introduced during the summer season result in the highest probability of disease emergence.