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

**|**HHS Author Manuscripts**|**PMC2746865

Formats

Article sections

- Abstract
- 1 Introduction
- 2 Empirical Data and Motivation for the Model
- 3 Model Derivation
- 4 Model Analysis
- 5 Numerical Examples
- 6 Discussion
- References

Authors

Related links

J Theor Biol. Author manuscript; available in PMC 2010 October 21.

Published in final edited form as:

Published online 2009 July 16. doi: 10.1016/j.jtbi.2009.07.009

PMCID: PMC2746865

NIHMSID: NIHMS133046

Linda J. S. Allen,^{a,}^{*} Curtis L. Wesley,^{b} Robert D. Owen,^{c,}^{d} Douglas G. Goodin,^{e} David Koch,^{e} Colleen B. Jonsson,^{f} Yong-Kyu Chu,^{f} J. M. Shawn Hutchinson,^{e} and Robert L. Paige^{a}

The publisher's final edited version of this article is available at J Theor Biol

See other articles in PMC that cite the published article.

New habitat-based models for spread of hantavirus are developed which account for interspecies interaction. Existing habitat-based models do not consider interspecies pathogen transmission, a primary route for emergence of new infectious diseases and reservoirs in wildlife and man. The modeling of interspecies transmission has the potential to provide more accurate predictions of disease persistence and emergence dynamics. The new models are motivated by our recent work on hantavirus in rodent communities in Paraguay. Our Paraguayan data illustrate the spatial and temporal overlap among rodent species, one of which is the reservoir species for Jabora virus and others which are spillover species. Disease transmission occurs when their habitats overlap. Two mathematical models, a system of ordinary differential equations (ODE) and a continuous-time Markov chain (CTMC) model, are developed for spread of hantavirus between a reservoir and a spillover species. Analysis of a special case of the ODE model provides an explicit expression for the basic reproduction number, _{0}, such that if _{0} < 1, then the pathogen does not persist in either population but if _{0} > 1, pathogen outbreaks or persistence may occur. Numerical simulations of the CTMC model display sporadic disease incidence, a new behavior of our habitat-based model, not present in other models, but which is a prominent feature of the seroprevalence data from Paraguay. Environmental changes that result in greater habitat overlap result in more encounters among various species that may lead to pathogen outbreaks and pathogen establishment in a new host.

Successful transmission of a directly transmitted pathogen requires opportunities for contact between species. These opportunities often occur when the preferred habitat of a species overlaps or is invaded by a second species. Interspecies interactions, especially among species competitively utilizing the same resources, often result in aggressive encounters. If a pathogen is present in a reservoir host, the encounter may result in pathogen transmission to a naive host or adaptation of the pathogen to create a new reservoir. The reservoir population, the carrier of the pathogen and the long-term host, often does not exhibit disease symptoms or experience any additional mortality.

In this investigation, we develop, analyze, and numerically simulate solutions to two new habitat-based models for the spread of a directly transmitted pathogen between two species. Our goal is to model the process of interspecies pathogen transmission based on species habitat preferences. The motivation for the models comes from our recent study of hantavirus in Paraguay. Hantavirus (Family *Bunyaviridae*) is a genus of viruses, each generally associated with a specific rodent species (i.e., mice and rats). Approximately 30 different hantaviruses exist throughout the world, some of which cause human infection (Mills et al., 1997). Human infection is incidental, generally due to indirect transmission from contact with infectious rodent excreta, but may result in hantavirus pulmonary syndrome with a mortality rate as high as 37% (CDC, 2002). One of the reservoir species for hantavirus in Paraguay is *Akodon montensis* (Montane Akodont found in Eastern Paraguay, Northeastern Argentina and Southeastern Brazil) carrier of Jabora virus (JABV, GenBank # EF492471). Our empirical data show that although these species exhibit different habitat preferences, the combination of partial habitat flexibility and temporally variable climatic and soil and vegetation conditions, results in periodic microgeographic sympatry of *Akodon* with one or both of the spillover species.

Mathematical models for the spread of hantavirus in rodents have concentrated primarily on the dynamics of the reservoir population (Abramson and Kenkre, 2003; Abramson et al., 2003; Allen et al, 2006a; Allen et al., 2003; Allen et al., 2006b; Sauvage et. al, 2007; Sauvage et al., 2003; Wesley, 2008; Wesley et al., 2009; Wolf et al., 2006). A multi-species epizootic model for susceptible and infected hosts was formulated and analyzed by McCormack and Allen (2007) but this model was not spatially explicit and did not account for differences in epizootiology of reservoir and spillover species. Our new models take into account habitat partitioning and important differences in the epizootiology of the reservoir and spillover populations. The role of the spillover species in pathogen and disease emergence is not well understood. It has been speculated that the spillover species may contribute to maintenance of the pathogen in the wild, provided there is spillback infection (McCormack and Allen, 2007) or the spillover species may be instrumental in the evolution of new hantaviruses (Chu et al., 2006). Spillover infections occur in hantavirus (Delfraro et al., 2008; Palma et al., 2009; Klingstrom et al., 2002; Torrez-Martinez et al., 2005; Wiedmann et al., 2005) but are not unique to hantaviruses (Daszak et al, 2000); they have been documented in other zoonotic diseases including rabies (Nadin-Davis and Loza-Ruio, 2006; Nel et al., 1997), Nipah virus (Chua, 2003), canine distemper, parvovirus (Fiorello et. al, 2006), and the SARS coronavirus (Holmes, 2003).

The habitat-based models consist of three regions: a preferred habitat for each of the reservoir and the spillover populations, and a third region of overlap (or boundary region) where interspecies encounters and pathogen transmission may occur. We formulate two models, the first model is a deterministic model, a system of ordinary differential equations (ODE), whereas the second model is a stochastic model, a continuous-time Markov chain (CTMC) model. The ODE system is analytically tractable in the case that encounters in the boundary region are brief. In this case an explicit expression for the basic reproduction number, _{0}, can be calculated, the threshold for disease outbreaks. The basic reproduction number is one of the most important parameters in the study of disease ecology. Specifically, it is the number of secondary infections caused by introduction of one infectious individual into an entirely susceptible population (see Anderson and May, 1991; Hethcote, 2000). The ODE system does not capture the few cases that occur due to interspecies interactions in the region of overlap. Therefore, we formulate a CTMC model for this purpose. Numerical solutions of the stochastic model illustrate sporadic infection in the spillover species when habitats overlap, a prominent feature of the seroprevalence data from Paraguay. Analysis and simulation of our new models show that as the number of encounters in the overlap region and the time spent in the overlap region increase (which may be triggered by habitat change), there is greater likelihood of pathogen outbreaks and disease persistence in the reservoir and spillover populations (through increase in _{0}). Interspecies encounters and pathogen transmission in the region of overlap may be the first step in the evolution of a new hantavirus strain.

Recent data collected in Paraguay (2005–2007) have shown cases of hantaviral infection in *A. montensis,* the reservoir species for JABV. Spillover infection, presumably of JABV, has been found in several other mouse species, including *Necromys lasiurus* (Hairy-tailed Akodont in Central Brazil, Southeastern Peru, Eastern Paraguay and Northeastern Argentina) and *Oxymycterus delator* (Paraguayan Hocicudo in Eastern Paraguay and South Central Brazil).

Our field work was conducted in the Mbaracayú Biosphere Reserve in eastern Paraguay, which lies in the western-most portion of the Interior Atlantic Forest. Vegetation composition in this area is typical of the mixture of intact, disturbed, and deforested areas found in eastern Paraguay (Fernández Soto and Mata Olmo, 2001).

Within this landscape of mixed habitat types, rodents were sampled on two mark-recapture grids, R3A and R3B, representing contrasting potential habitat for *A. montensis* (see Figure 1). Site R3A is largely deforested, with its natural cover replaced by native and exotic graminoids and forbs. This site is highly disturbed by human activities and is intensively managed for pasturing and grazing. Reforestation is suppressed and graminoid cover maintained by frequent prescribed burning. Large ungulate grazers, primarily domestic cattle (*Bos taurus*) are present on this site year-round. Vegetation in the site is dominated by species of the genus *Andropogon,* warm temperature/tropical grasses used as grazing forage. Other common vegetation genera include *Merostachys* (bamboo) along the fringes of pastures and *Xyris* (a forb) in lower, wetter areas. Although dominated by herbaceous species, R3A retains a few islands of woody vegetation and trees, especially along the edges of the deforested areas. These forest remnants are better microhabitats for *Akodon,* and most captures of this species were along the northwest and southeast corners of the grid, the edges closest to the woodlands.

Natural color satellite imagery showing location and general land cover characteristics of the two data collection grids. Larger image (acquired Feb. 2003 from Landsat ETM+), shows the general setting of the collection grids within and near the Mbaracayú **...**

Site R3B contrasts with R3A in that its dominant cover consists of native forest and its associated vegetation community. Although native cover remains, the site shows evidence of recent human disturbance, especially selective logging and nearby road construction. These disturbances have resulted in fragmentation of the native vegetation cover, producing numerous internal edges and gaps in the forest canopy. These edges and gaps are associated with a dense understory at the forest floor, favorable habitat for *A. montensis* (Pardiñas and D’Elía, 2003). Dominant vegetation genera in areas of intact or mostly intact forest include *Cedrela* and *Balfourodendron.* In disturbed areas, a lower canopy dominated by *Sorocea bonplandii* (a lower, woody shrub-like tree) is common, with herbaceous species and *Bromelia* common in the understory, and *Merostachys* frequently in forest gaps.

Each sampling grid consisted of an 11×11 square array of trap stations, set 10 m apart. One standard Sherman Live Trap was set on the ground at each station. Where vegetation structure permitted (R3B only), another trap was placed 2–3 m above ground, in branches or vines, to sample other more arboreal species. However, although *A. montensis is* preferentially a forest-dweller, it only rarely climbs in the vegetation (651 of 661 [98.5 %] of our captures of this species on R3B were on the ground). Both grids were sampled in nine sessions from February 2005 through May 2007, using standard mark-recapture techniques for small mammals (Wilson et al., 1996). Each session included eight consecutive nights (seven in February 2005). Animals were individually marked with a subdermally-implanted Passive Integrated Transponder (PIT) tag, which can be read by passing the electronic reader near the animal’s body. Date, grid, trap station, and the PIT tag number were recorded, animals were identified to species, and sex and weight were recorded. Blood, saliva, urine, and feces were collected the first time an animal was captured during each sampling session, for assay of hantavirus antibody and viral RNA presence. Animals were then released at the site of capture. All field protocols were approved by the Texas Tech University Animal Care and Use Committee.

A total of 582 captures was recorded on R3B, and 332 on R3A. Species captured on R3B included *A. montensis, Calomys callosus, Oligoryzomys fornesi, O. nigripes, and Oryzomys megacephalus;* on R3A, *A. montensis, C. callosus, C. tener, N. lasiurus, O. nigripes,* and *O. delator* were encountered. In this report we consider only the populations of *A. montensis, N. lasiurus,* and *O. delator,* as exemplifying the scenario being modeled. *A. montensis* has been described as the primary reservoir of JABV. RNA-positive individuals of *Akodon* were encountered on R3B and R3A, and seropositive individuals of N. *lasiurus* and *O. delator* were trapped and identified on R3A. Each of these species is widely distributed in the central Southern Cone of South America (Carleton and Musser, 2005), and may be locally abundant in their preferred habitat. Habitat preferences differ somewhat among the three, which is critical to this field situation and to the model which we present herein. *N. lasiurus* and *O. delator* prefer grasslands, with *Necromys* preferring dry soil and *Oxymycterus* preferring wet or even saturated soils. *A. montensis* preferentially inhabits disturbed woodlands, also venturing into old-fields and grasslands which include forbs and brushy growth (Redford and Eisenberg, 1992; Goodin et al., manuscript). Our data from nine sampling sessions through 27 months support these descriptions of habitat preferences, and further indicate that the microgeographic separation among these three species is partial and temporally variable and often incomplete, with overlap (microsympatry) occurring sporadically between *A. montensis* and one or both of the spillover species. Figure 2 illustrates this situation on site R3A, and also indicates the presence of seropositive individuals for each of these species in close proximity to the others.

Data for nine trapping sessions in 2005–2007 for three species *A. montensis* (AKMO), *N. lasiurus* (NLAS) and *O. delator* (ODEL) at site R3A. Symbols indicate numbers of animals captured at each station during eight nights of trapping, and the number **...**

At site R3B, *A. montensis* was captured but neither *N. lasiurus* nor *O. delator.* Based on 2005–2007 data, 21 out of 84 *A. montensis* males (25%) tested for hantavirus showed positive titers for antibodies or RNA, including 12 (14.3%) that were only antibody-positive, eight (9.5%) that were antibody-positive and RNA-positive, and one (1.2%) that was only RNA-positive. Only four out of 68 *A. montensis* females (5.9%) tested for hantavirus showed positive titers for antibodies with no detectable viral RNA. In other studies of hantavirus ecology, male seroprevalence was higher than female seroprevalence (Bernshtein et al., 1999; Childs, et al., 1994; Glass et al., 1998; Klein et al., 2001; McIntyre et al., 2005; Mills et al., 1997; Yahnke et al., 2001).

Based on the empirical data, we model only male rodents in the spread of hantavirus and use two infectious stages for the reservoir host, a highly infectious stage and a persistent stage, *I* and *P.* The highly infectious stage represents animals that are RNA-positive and may or may not have antibodies, and the persistent stage represents animals that are only antibody-positive. This latter persistent stage is less infectious than the highly infectious stage. Two infectious stages were assumed in models for Puumala hantavirus in bank voles (Sauvage et al., 2007; Sauvage et al., 2003; Wolf et al., 2006).

First, two basic models for the reservoir and the spillover species, each within their own habitat, are formulated. Then these two models are merged into a habitat-based model, where interspecies pathogen transmission may occur in a region of overlap of the two habitats (such as R3A).

The disease stages for the reservoir species include susceptible, *S _{r},* exposed or latently infected,

$$\begin{array}{l}\frac{d{S}_{r}}{dt}={b}_{r}{N}_{r}-{S}_{r}({\beta}_{I}{I}_{r}+{\beta}_{P}{P}_{r})-{S}_{r}{d}_{r}({N}_{r})\\ \frac{d{E}_{r}}{dt}={S}_{r}({\beta}_{I}{I}_{r}+{\beta}_{P}{P}_{r})-{\delta}_{r}{E}_{r}-{E}_{r}{d}_{r}({N}_{r})\\ \frac{d{I}_{r}}{dt}={\delta}_{r}{E}_{r}-{\gamma}_{r}{I}_{r}-{I}_{r}{d}_{r}({N}_{r})\\ \frac{d{P}_{r}}{dt}={\gamma}_{r}{I}_{r}-{P}_{r}{d}_{r}({N}_{r}).\end{array}$$

(3.1)

The total population density satisfies the differential equation

$$\frac{d{N}_{r}}{dt}={N}_{r}\phantom{\rule{0.16667em}{0ex}}[{b}_{r}-{d}_{r}({N}_{r})].$$

(3.2)

From the assumptions on *d _{r},* it follows that there exists a unique positive constant

The spillover species responds differently to hantavirus infection. Presumably, the infection is only short-term, an acute stage, *A*, and therefore, we assume no disease-related deaths occur. But this assumption can be modified. The model is an SEAR model, where animals pass through the stages of being susceptible, latent, infectious, and finally recovered. A subscript *s* is used to identify the spillover species and distinguish it from the reservoir species. The differential equations for hantaviral infection in the spillover species are similar to the reservoir species but transmission of hantavirus occurs only from the infectious stage *A _{s}* and

$$\begin{array}{l}\frac{d{S}_{s}}{dt}={b}_{s}{N}_{s}-{\beta}_{A}{S}_{s}{A}_{s}-{S}_{s}{d}_{s}({N}_{s})\\ \frac{d{E}_{s}}{dt}={\beta}_{A}{S}_{s}{A}_{s}-{\delta}_{s}{E}_{s}-{E}_{s}{d}_{s}({N}_{s})\\ \frac{d{A}_{s}}{dt}={\delta}_{s}{E}_{s}-{\gamma}_{s}{A}_{s}-{A}_{s}{d}_{s}({N}_{s})\\ \frac{d{R}_{s}}{dt}={\gamma}_{s}{A}_{s}-{R}_{s}{d}_{s}({N}_{s}).\end{array}$$

(3.3)

We assume *d _{s}*(

The reservoir and spillover species generally have preferred habitats as shown by the data. The spillover species is rarely found in the habitat where the reservoir species is dominant and vice versa (see Figure 2). However, contact between these two species occurs in a boundary or overlap region adjacent to their habitats, where densities of the two species may be relatively low. Encounters between infectious and susceptible animals in this boundary region may result in interspecies transmission of hantavirus. Figure 3 is a schematic of the three regions, the preferred habitats for the reservoir and the spillover species and the boundary region.

A schematic of the three regions representing preferred habitats for the reservoir and spillover species and the boundary or overlap region adjacent to the two habitats.

These habitats are connected via movement to and from the preferred habitat and the boundary region. Time spent in this boundary region is short for both species. The majority of the population is susceptible, especially in the case of the spillover species.

Suppose the per capita rate of movement *p _{i}* into the boundary region is low i.e.,

Based on the preceding assumptions, the differential equations for the reservoir species in its preferred habitat take the following form:

$$\begin{array}{l}\frac{d{S}_{r}}{dt}={b}_{r}{N}_{r}-{S}_{r}({\beta}_{I}{I}_{r}+{\beta}_{P}{P}_{r})-{S}_{r}{d}_{r}({N}_{r})-{p}_{i}{S}_{r}+{p}_{o}{S}_{a}\\ \frac{d{E}_{r}}{dt}={S}_{r}({\beta}_{I}{I}_{r}+{\beta}_{P}{P}_{r})-{\delta}_{r}{E}_{r}-{E}_{r}{d}_{r}({N}_{r})-{p}_{i}{E}_{r}+{p}_{o}{E}_{a}\\ \frac{d{I}_{r}}{dt}={\delta}_{r}{E}_{r}-{\gamma}_{r}{I}_{r}-{I}_{r}{d}_{r}({N}_{r})-{p}_{i}{I}_{r}+{p}_{o}{I}_{a}\\ \frac{d{P}_{r}}{dt}={\gamma}_{r}{I}_{r}-{P}_{r}{d}_{r}({N}_{r})-{p}_{i}{P}_{r}+{p}_{o}{P}_{a}\end{array}$$

(3.4)

Similar differential equations apply to the spillover species, where terms for movement into and out of the preferred habitat are added to the differential equations in (3.3).

Because rodents are in the boundary region for a short period of time, on the order of days, no births nor deaths occur in this region. We assume that the carrying capacity in the preferred habitat remains constant. That is, the stable preferred habitat density for the reservoir species is *K _{r}* which can be thought of as a population source for the boundary region. The carrying capacity in the boundary region may increase or decrease relative to

$$\begin{array}{l}\frac{d{S}_{a}}{dt}=-{S}_{a}({\beta}_{a1}{I}_{a}+{\beta}_{a2}{P}_{a}+{\beta}_{a3}{A}_{b})+{p}_{i}{S}_{r}-{p}_{o}{S}_{a}\\ \frac{d{E}_{a}}{dt}={S}_{a}({\beta}_{a1}{I}_{a}+{\beta}_{a2}{P}_{a}+{\beta}_{a3}{A}_{b})-{\delta}_{a}{E}_{a}+{p}_{i}{E}_{r}-{p}_{o}{E}_{a}\\ \frac{d{I}_{a}}{dt}={\delta}_{a}{E}_{a}-{\gamma}_{a}{I}_{a}+{p}_{i}{I}_{r}-{p}_{o}{I}_{a}\\ \frac{d{P}_{a}}{dt}={\gamma}_{a}{I}_{a}+{p}_{i}{P}_{r}-{p}_{o}{P}_{a}\end{array}$$

(3.5)

and for the spillover species they are

$$\begin{array}{l}\frac{d{S}_{b}}{dt}=-{S}_{b}({\beta}_{b1}{I}_{a}+{\beta}_{b2}{P}_{a}+{\beta}_{b3}{A}_{b})+{p}_{i}{S}_{s}-{p}_{o}{S}_{b}\\ \frac{d{E}_{b}}{dt}={S}_{b}({\beta}_{b1}{I}_{a}+{\beta}_{b2}{P}_{a}+{\beta}_{b3}{A}_{b})-{\delta}_{b}{E}_{b}+{p}_{i}{E}_{s}-{p}_{o}{E}_{b}\\ \frac{d{A}_{b}}{dt}={\delta}_{b}{E}_{b}-{\gamma}_{b}{A}_{b}+{p}_{i}{A}_{s}-{p}_{o}{A}_{b}\\ \frac{d{R}_{b}}{dt}={\gamma}_{b}{A}_{b}+{p}_{i}{R}_{s}-{p}_{o}{R}_{b}.\end{array}$$

(3.6)

Rodents may change their disease status while in the boundary, e.g., *E _{a}* →

$${\int}_{0}^{\infty}{p}_{o}exp(-{p}_{o}t){\int}_{0}^{t}{\delta}_{a}exp(-{\delta}_{a}s)\phantom{\rule{0.16667em}{0ex}}ds\phantom{\rule{0.16667em}{0ex}}dt=\frac{{\delta}_{a}}{{p}_{o}+{\delta}_{a}}.$$

(3.7)

Other models based on more general probability distributions such as the gamma distribution provide alternative formulations (Feng et al., 2007; Lloyd, 2001a, 2001b). More data are required to determine the form of the distributions. In this investigation, we consider the simplest form, an exponential distribution. If *p _{o}* max

The total male population densities in the preferred habitat and in the boundary region are *N _{r}* and

$$\begin{array}{l}\frac{d{N}_{r}}{dt}={p}_{o}{N}_{a}+{N}_{r}[{b}_{r}-{p}_{i}-{d}_{r}({N}_{r})]\\ \frac{d{N}_{a}}{dt}={p}_{i}{N}_{r}-{p}_{o}{N}_{a}\\ \frac{d{N}_{s}}{dt}={p}_{o}{N}_{b}+{N}_{s}[{b}_{s}-{p}_{i}-{d}_{s}({N}_{s})]\\ \frac{d{N}_{b}}{dt}={p}_{i}{N}_{s}-{p}_{o}{N}_{b}.\end{array}$$

(3.8)

Initial conditions are nonnegative and strictly positive in the preferred habitat; *N _{r}*(0) > 0 and

The ODE model can be easily extended to a CTMC model which includes variability in the birth, death, transmission, and movement processes. In the CTMC model, the 16 random variables are integer-valued taking on values in the set {0, 1, 2,…}. There are 38 different events, that include births, deaths, transmission, and movement. Let (*t*) be a vector of 16 discrete random variables associated with the CTMC process:

$$\mathcal{X}=({S}_{r},{E}_{r},{I}_{r},{P}_{r},{S}_{a},{E}_{a},{I}_{a},{P}_{a},{S}_{b},{E}_{b},{A}_{b},{R}_{b},{S}_{s},{E}_{s},{A}_{s},{R}_{s}).$$

Based on the 38 events, the infinitesimal transition probabilities can be defined, Prob{Δ(*t*)*|*(*t*)}, where Δ (*t*) = (*t* + Δ*t*) − (*t*) for Δ*t* sufficiently small (Allen, 2003; Karlin and Taylor, 1975). For example, the probability of a birth in the reservoir population is

$$\text{Prob}\{\mathrm{\Delta}\mathcal{X}(t)=(1,0,\dots ,0)\mathcal{X}(t)\}={b}_{r}{N}_{r}(t)\mathrm{\Delta}t+o(\mathrm{\Delta}t)={b}_{r}\sum _{i=1}^{4}{X}_{i}(t)\mathrm{\Delta}t+o(\mathrm{\Delta}t),$$

where the time step Δ*t* is chosen so that the possibility of more than one transition or change in Δ*t* units of time is negligible. The 38 events and their corresponding transition probabilities are described in the Appendix.

Due to climatic variations within the year-dry, wet, and transitional (D, W, T) periods-there may be greater overlap of the habitats during certain periods of the year. One way of modeling this variability in the overlap region is to modify the rate of movement into or out of the boundary region for each species, depending on their habitat preferences during each of these periods. For example, *A. montensis* is seen in the overlap region more frequently in period T than in periods D or W, when densities of *O. delator* are high. In the model, we do not specifically include this seasonal variability but we do consider the effects of changes in *p _{i}* and

There are three types of equilibria for the ODE habitat-based model: an extinction equilibrium, where the population density is zero, a unique disease-free equilibrium (DFE), where all the infectious and recovered states are zero but the susceptible states are positive, and enzootic equilibria (EE), where some infectious states have positive values. It can be easily shown that the extinction equilibrium is unstable; the population persists. For example, it follows from equations (3.8) that in the preferred habitats, the population densities approach a constant value, their respective carrying capacities,

$$\underset{t\to \infty}{lim}{N}_{r}(t)={K}_{r}\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.38889em}{0ex}}\underset{t\to \infty}{lim}{N}_{s}(t)={K}_{s}.$$

The preferred habitats serve as a population source for the boundary region. The densities in the boundary region depend on the densities of these source populations and the movement rates into and out of this region. In particular, in the boundary region, the reservoir and spillover population densities are

$$\underset{t\to \infty}{lim}{N}_{a}(t)=\frac{{p}_{i}}{{p}_{o}}{K}_{r}={K}_{a}\text{and}\phantom{\rule{0.38889em}{0ex}}\underset{t\to \infty}{lim}{N}_{b}(t)=\frac{{p}_{i}}{{p}_{o}}{K}_{s}={K}_{b},$$

respectively. As the ratio *p _{i}/p_{o}* increases, so do the population densities in the boundary region, whereas the densities in the preferred habitats will approach their respective carrying-capacities. Hence, the total reservoir population density (preferred + boundary) is

We calculate reproduction numbers for each of the preferred habitats and an approximation to the overall basic reproduction number for the habitat-based model (3.4)–(3.6), _{0}. If _{0} > 1, then it is likely that the disease persists in the reservoir and spillover populations. The reproduction numbers can be calculated using the next generation matrix approach (van den Driessche and Watmough, 2002). For the general system (3.4)–(3.6), it is possible to show that the basic reproduction number is a positive root of a fourth degree polynomial but it is difficult to obtain a simple analytical expression for _{0}. The simplifying assumption

$${\delta}_{i}=0={\gamma}_{i},\phantom{\rule{0.38889em}{0ex}}i=a,b,$$

(4.9)

leads to an explicit expression for _{0} (shown below). This explicit expression is a close approximation to the overall basic reproduction number, if the time spent in the boundary is short relative to the time spent in each of the disease states. This expression is very useful in interpreting the contributions to disease outbreaks by the reservoir and the spillover species and in making comparisons to other reproduction numbers.

Assume the condition (4.9) holds. First, the reproduction number for the reservoir species (assuming the spillover species is not present) is

$${\mathcal{R}}_{0}^{r}=\frac{{p}_{o}{K}_{r}({\beta}_{I}{\delta}_{r}{b}_{r}+{\beta}_{P}{\delta}_{r}{\gamma}_{r})+{p}_{i}{K}_{a}({\beta}_{a1}{\delta}_{r}{b}_{r}+{\beta}_{a2}{\delta}_{r}{\gamma}_{r})}{{p}_{o}({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r}){b}_{r}}.$$

Second, the reproduction number for the spillover species (assuming the reservoir species is not present) is

$${\mathcal{R}}_{0}^{s}=\frac{{p}_{o}{K}_{s}{\beta}_{A}{\delta}_{s}+{p}_{i}{K}_{b}{\beta}_{b3}{\delta}_{s}}{{p}_{o}({\delta}_{s}+{b}_{s})({\gamma}_{s}+{b}_{s})}.$$

If there are no interspecies interactions so that *p _{i}* = 0 =

$$\begin{array}{l}{\overline{S}}_{r}=\frac{{K}_{r}}{{\mathcal{R}}^{r}}\\ {\overline{E}}_{r}=\frac{{b}_{r}{K}_{r}}{{\delta}_{r}+{b}_{r}}\left(1-\frac{1}{{\mathcal{R}}^{r}}\right)\\ {\overline{I}}_{r}=\frac{{\delta}_{r}}{{\gamma}_{r}+{b}_{r}}{\overline{E}}_{r}\\ {\overline{P}}_{r}=\frac{{\gamma}_{r}}{{b}_{r}}{\overline{I}}_{r},\end{array}$$

(4.10)

whenever the reproduction number for the SEIP model (3.1) is

$${\mathcal{R}}^{r}=\frac{{K}_{r}({\beta}_{I}{\delta}_{r}{b}_{r}+{\beta}_{P}{\delta}_{r}{\gamma}_{r})}{({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r}){b}_{r}}>1.$$

At equilibrium, the proportion of animals that are RNA- or antibody-positive, (*Ī _{r}* +

$$\frac{{\delta}_{r}}{{\delta}_{r}+{b}_{r}}\left(1-\frac{1}{{\mathcal{R}}^{r}}\right).$$

(4.11)

To derive an expression for the basic reproduction number for the habitat-based model when (4.9) holds, we first define an expression
${\mathcal{R}}_{0}^{c}$ which depends on interspecies or crossover transmission. That is, let the intraspecies transmission parameters be zero, *β _{I}* =

$${\mathcal{R}}_{0}^{c}=\left[\frac{{p}_{i}{\delta}_{r}{K}_{b}({\beta}_{b1}{b}_{r}+{\beta}_{b2}{\gamma}_{r})}{{p}_{o}({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r}){b}_{r}}\right]\phantom{\rule{0.38889em}{0ex}}\left[\frac{{p}_{i}{\delta}_{s}{\beta}_{a3}{K}_{a}}{{p}_{o}({\gamma}_{s}+{b}_{s})({\delta}_{s}+{b}_{s})}\right].$$

Note that
${\mathcal{R}}_{0}^{c}$ depends on the ratio *p _{i}/p_{o}* directly and indirectly through the carrying capacities in the overlap region,

$${\mathcal{R}}_{0}=\frac{1}{2}\left({\mathcal{R}}_{0}^{r}+{\mathcal{R}}_{0}^{s}+\sqrt{{({\mathcal{R}}_{0}^{r}-{\mathcal{R}}_{0}^{s})}^{2}+4{\mathcal{R}}_{0}^{c}}\right).$$

(4.12)

(Derivation of this formula is given in the Appendix.) It follows that ${\mathcal{R}}_{0}\ge max\{{\mathcal{R}}_{0}^{r},{\mathcal{R}}_{0}^{s}\}$; interspecies pathogen transmission increases the basic reproduction number. A similar relationship was shown in a multi-species SI model of McCormack and Allen (2007).

The local stability of the DFE follows directly from the results of van den Driessche and Watmough (2002). Global stability of the DFE when _{0} < 1 and condition (4.9) holds can be verified by construction of a Liapunov function. The full system (3.4)–(3.6) consists of 16 differential equations which makes it difficult to find an explicit closed form solution for an enzootic equilibrium (EE). However, existence and uniqueness of a positive EE can be verified when _{0} > 1 and condition (4.9) holds. It is shown, in the Appendix, that the EE for the full system (3.4)–(3.6) is a fixed point of *E _{r}* =

A basic reproduction number _{0} exists for system (3.4)–(3.6) such that if _{0} < 1, the DFE is locally asymptotically stable. If condition (4–9) holds for system (3.4)–(3.6), then

_{0}has the form given in (4.12),- if
_{0}< 1, then the DFE is globally asymptotically stable, and - if
_{0}> 1; then the DFE is unstable and there exists a unique positive enzootic equilibrium.

Selection of parameters values is based on estimates from the literature for hantavirus (*δ _{i}* and

$${d}_{r}({N}_{r})={b}_{r}{N}_{r}/{K}_{r}\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.38889em}{0ex}}{d}_{s}({N}_{s})={b}_{s}{N}_{s}/{K}_{s},$$

so that *d _{r}*(

The transmission parameters *β _{j}* cannot be estimated directly. Instead, we make some reasonable assumptions about their relationship to disease transmission in the infectious stages,

$${\beta}_{{a}_{1}}={\beta}_{I},\phantom{\rule{0.38889em}{0ex}}{\beta}_{{a}_{2}}={\beta}_{P},\phantom{\rule{0.38889em}{0ex}}{\beta}_{{b}_{3}}={\beta}_{A},$$

(5.13)

$${\beta}_{{b}_{1}}=2{\beta}_{I},\phantom{\rule{0.38889em}{0ex}}{\beta}_{{b}_{2}}=2{\beta}_{P},\phantom{\rule{0.38889em}{0ex}}{\beta}_{{a}_{3}}=2{\beta}_{A}.$$

(5.14)

The basic parameter values are given in Table 1.

If (4.9) holds and the remaining parameter values are as in Table 1 with *p _{i}* = 8 and

$${\mathcal{R}}_{0}^{r}=1.42,\phantom{\rule{0.38889em}{0ex}}{\mathcal{R}}_{0}^{s}=0.04,\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.16667em}{0ex}}{\mathcal{R}}_{0}=\mathrm{1.42.}$$

The overall basic reproduction number for the parameters in Table 1 is _{0}= 1.38 which is close to the approximation 1.42. The disease persists in the habitat-based model. For the parameter values in Table 1, one sample path of the CTMC model and the solution to the ODE model are graphed for the two infectious stages of the reservoir species (see Figure 4).

Although the pathogen persists, the infection in the spillover population is very low (straight line is the ODE equilibrium value); only sporadic infection occurs in the sample path for the spillover species in the preferred habitat and in the boundary region (Figure 5).

Solution to the ODE model and one sample path of the CTMC model for the spillover species in its preferred habitat (*A*_{s}) and in the boundary region (*A*_{b}). In the ODE model, *Ā*_{s} = 0.07 and *Ā*_{b} = 0.02, whereas, in the CTMC model, the values **...**

For the ODE model with parameter values given in Table 1, *p _{i}* = 8, and

$$\begin{array}{c}({\overline{S}}_{r},{\overline{E}}_{r},{\overline{I}}_{r},{\overline{P}}_{r};\phantom{\rule{0.38889em}{0ex}}{\overline{S}}_{a},{\overline{E}}_{a},{\overline{I}}_{a},{\overline{P}}_{a})=(72.4,2.6,10.0,15.2;\phantom{\rule{0.38889em}{0ex}}11.1,0.3,1.5,2.5)\\ ({\overline{S}}_{s},{\overline{E}}_{s},{\overline{A}}_{s},{\overline{R}}_{s};\phantom{\rule{0.38889em}{0ex}}{\overline{S}}_{b},{\overline{E}}_{b},{\overline{A}}_{b},{\overline{R}}_{b})=(49.1,0.06,0.07,0.8;\phantom{\rule{0.38889em}{0ex}}7.5,0.04,0.02,0.1).\end{array}$$

With no interspecies transmission and no overlap region (*p _{i}* = 0 =

The CTMC simulation with interspecies transmission illustrates the sporadic infection in the spillover population (as in site R3A) and provides information about the variability in number of cases. A quasistationary probability distribution is reached in the CTMC model (conditional on nonextinction). Approximations (estimated from 10,000 sample paths) to the quasistationary probability distributions for the two infectious stages in the reservoir species are graphed in Figure 6. The mean values for *I _{r}* and

Encounters that lead to interspecies pathogen transmission can be measured by the magnitude of
${\mathcal{R}}_{0}^{c}$. The greater the habitat overlap, the greater the number of interspecies and intraspecies encounters which in turn increase the likelihood of pathogen outbreaks and disease persistence. Changes that affect the overlap region will have the greatest impact on the parameters *p _{i}* and

Biologically-motivated models for pathogen spread between two species were formulated, an ODE model (3.4)–(3.6) and a CTMC model. The models are based on the fact that spatial overlap of habitats leads to greater numbers of interspecies encounters. From the ODE model, an explicit expression for the basic reproduction number _{0} was calculated based on the assumption (4.9), as well as reproduction numbers for the preferred habitats,
${\mathcal{R}}_{0}^{r}$ and
${\mathcal{R}}_{0}^{s}$, and for crossover or interspecies transmission,
${\mathcal{R}}_{0}^{c}$. In this case, we showed global stability of the disease-free equilibrium when _{0} < 1 and existence of an enzootic equilibrium when _{0} > 1. Greater numbers of interactions among species allow the pathogen to be transmitted more frequently from an infectious host to a susceptible host. This, in turn, increases _{0} so that it exceeds the reproduction number in the preferred habitat, _{0} > * ^{r}* (Figure 7), which ultimately results in greater likelihood of outbreaks and disease persistence. As illustrated in Figure 2, the overlap region is spatially- and temporally-dependent. We did not consider temporal variability of this overlap region which may depend on seasonal variations. But we did include demographic variability due to births, deaths, transmission, and movement in the CTMC model. Seasonal variations, in general, will cause additional variability in the solution behavior (e.g., Allen et al., 2005). As more data are collected, the effects of seasonal and climatic variations on the reservoir and spillover species will be studied. In addition, controlled studies are needed to obtain data on the duration and shape of the rodents’ disease stage distributions.

Interspecies pathogen transmission, where a known virus “jumps” into a new host, is one of the primary reasons for the large increase in emerging diseases in wildlife in recent years (Daszak et al., 2000; Parrish et al., 2008; Richomme et. al., 2006). Our mathematical models illustrate the first step in this emergence and the role that the spillover species may play in emerging diseases. Our models were developed for spread of hantavirus in rodents but can be modified and applied to other species, where spatial spread results in spillover infection.

This research was supported by a grant from the Fogarty International Center #R01TW006986-02 under the NIH NSF Ecology of Infectious Diseases initiative. We thank R. K. McCormack for preliminary discussions on this work, the Fundación Moises Bertoni for facilitating access to the field sites, and the Vendramini family for allowing us to work in Estancia Rama III. The Secretaría de Ambiente provided necessary permits for working with wildlife. In addition, we thank the referees for their helpful suggestions.

The basic reproduction number _{0} for system (3.4)–(3.6) is obtained via the next generation matrix approach (van den Driessche and Watmough, 2002). Reorder the vector of 16 state variables as follows:

$$\mathcal{Y}=({E}_{r},{E}_{a},{E}_{b},{E}_{s},{I}_{r},{I}_{a},{A}_{b},{A}_{s},{P}_{r},{P}_{a},{S}_{r},{S}_{a},{S}_{b},{S}_{s},{R}_{b},{R}_{s}).$$

Then ˙ = − , where is the changes due to new infections and is the other transitions. The Jacobian matrix of and , evaluated at the DFE, is

$${D}_{\mathcal{F}}=\left(\begin{array}{cc}F& \mathbf{O}\\ \mathbf{O}& \mathbf{O}\end{array}\right)\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.38889em}{0ex}}{D}_{\mathcal{V}}=\left(\begin{array}{cc}V& \mathbf{O}\\ {J}_{1}& {J}_{2}\end{array}\right),$$

where *F* and *V* are 10 × 10 matrices, corresponding to the Jacobian matrices for the first 10 variables in . The symbol **O** represents a zero matrix. Matrix −*J*_{2}, a 6 × 6 matrix, has eigenvalues with negative real part. The basic reproduction number, _{0}, is the spectral radius of

$$F{V}^{-1}=\left(\begin{array}{cc}A& B\\ \mathbf{O}& \mathbf{O}\end{array}\right),$$

where matrix *A* is a 4 × 4 nonnegative matrix. Thus,

$${\mathcal{R}}_{0}=\rho (A),$$

the largest positive root of the characteristic polynomial of matrix *A*, a fourth degree polynomial. Under condition (4.9),

$$A=\left(\begin{array}{cccc}{c}_{1}& {c}_{1}& 0& 0\\ {c}_{2}& {c}_{2}& {c}_{3}& {c}_{3}\\ {c}_{4}& {c}_{4}& {c}_{5}& {c}_{5}\\ 0& 0& {c}_{6}& {c}_{6}\end{array}\right).$$

The elements of matrix *A* are

$$\begin{array}{l}{c}_{1}=\frac{{K}_{r}{\delta}_{r}({\beta}_{I}{b}_{r}+{\beta}_{P}{\gamma}_{r})}{({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r}){b}_{r}},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{c}_{2}=\frac{{p}_{i}{K}_{a}{\delta}_{r}({\beta}_{a1}{b}_{r}+{\beta}_{{a}_{2}}{\gamma}_{r})}{{p}_{o}({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r}){b}_{r}}\\ {c}_{3}=\frac{{p}_{i}{K}_{a}{\beta}_{a3}{\delta}_{s}}{{p}_{o}({\delta}_{s}+{b}_{s})({\gamma}_{s}+{b}_{s})},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{c}_{4}=\frac{{p}_{i}{K}_{b}{\delta}_{r}({\beta}_{b1}{b}_{r}+{\beta}_{b2}{\gamma}_{r})}{{p}_{o}({\delta}_{r}+{b}_{r})({\gamma}_{r}+{b}_{r})}\\ {c}_{5}=\frac{{p}_{i}{K}_{b}{\beta}_{b3}{\delta}_{s}}{{p}_{o}({\delta}_{s}+{b}_{s})({\gamma}_{s}+{b}_{s})},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{c}_{6}=\frac{{K}_{s}{\beta}_{A}{\delta}_{s}}{({\delta}_{s}+{b}_{s})({\gamma}_{s}+{b}_{s})},\end{array}$$

(6.15)

where *b _{r}* =

The characteristic polynomial of matrix *A* is

$${\lambda}^{2}[{\lambda}^{2}-({c}_{1}+{c}_{2}+{c}_{5}+{c}_{6})\lambda +({c}_{1}+{c}_{2})({c}_{5}+{c}_{6})-{c}_{3}{c}_{4}].$$

Note that
${c}_{1}+{c}_{2}={\mathcal{R}}_{0}^{r},{c}_{5}+{c}_{6}={\mathcal{R}}_{0}^{s}$ and
${c}_{3}{c}_{4}={\mathcal{R}}_{0}^{c}$. It follows that the largest positive root of the quadratic equation is _{0}, given in equation (4.12). Local stability of the DFE if _{0} < 1 and instability if _{0} > 1 follows directly from van den Driessche and Watmough (2002). This verifies the first part of Theorem 4.1.

Before we give the proof of Theorem 4.1 (ii)–(iii), some preliminary facts are stated in the two cases, _{0} < 1 and _{0} > 1. When _{0} < 1, to verify global stability, we use an expression equivalent to _{0},

$${\widehat{\mathcal{R}}}_{0}={\mathcal{R}}_{0}^{s}+{\mathcal{R}}_{0}^{r}+{\mathcal{R}}_{0}^{c}-{\mathcal{R}}_{0}^{r}{\mathcal{R}}_{0}^{s}.$$

It can be easily shown that _{0} < 1 (> 1, = 1) iff ○_{0} < 1 (> 1, = 1). Let

$${D}_{0}=\{({S}_{r},{I}_{r},\mathrm{..},{A}_{b},{R}_{b}){\mathbb{R}}_{+}^{16}:0\le {N}_{r}\le {K}_{r},0\le {N}_{a}\le {K}_{a},0\le {N}_{s}\le {K}_{s},0\le {N}_{b}\le {K}_{b}\}.$$

It is clear that *D*_{0} is forward invariant. As *t* → ∞, the total population density of each patch approaches their respective carrying capacities *K _{r}, K_{a}, K_{s},* and

$${D}_{1}=\{({S}_{r},{I}_{r},\mathrm{..},{A}_{b},{R}_{b}){\mathbb{R}}_{+}^{16}:{N}_{r}={K}_{r},{N}_{a}={K}_{a},{N}_{s}={K}_{s},{N}_{b}={K}_{b}\}$$

which is contained in *D*_{0}. We will show for _{0} < 1 that solutions to the limiting system (beginning in *D*_{1}) approach the DFE. Then global asymptotic stability of the DFE will follow from the theory of asymptotically autonomous systems.

The EE can be can be written as a fixed point problem. It can be easily verified that at an EE, the values of the variables *P _{r}, I_{r}, S_{r}, P_{a},* and

$${E}_{r}=\frac{{\mathrm{\Gamma}}_{r}}{({b}_{r}+{\delta}_{r})[{\mathrm{\Gamma}}_{r}+{b}_{r}{K}_{r}{J}_{r}]/({b}_{r}{K}_{r})}=f({E}_{r},{E}_{s})$$

(6.16)

and

$${E}_{s}=\frac{{\mathrm{\Gamma}}_{s}}{({b}_{s}+{\delta}_{s})[{\mathrm{\Gamma}}_{s}+{b}_{s}{K}_{s}{J}_{s}]/({b}_{s}{K}_{s})}=g({E}_{r},{E}_{s}),$$

(6.17)

where

$$\begin{array}{l}{\mathrm{\Gamma}}_{r}={N}_{1,r}{J}_{r}+{p}_{o}{N}_{2,r},\phantom{\rule{0.38889em}{0ex}}{\mathrm{\Gamma}}_{s}={N}_{1,s}{J}_{s}+{p}_{o}{N}_{2,s},\\ {J}_{r}={N}_{2,r}/{K}_{a}+{p}_{o},\phantom{\rule{0.38889em}{0ex}}{J}_{s}={N}_{2,s}/{K}_{b}+{p}_{o},\\ {N}_{1,r}={c}_{1}({b}_{r}+{\delta}_{r}){E}_{r},\phantom{\rule{0.38889em}{0ex}}{N}_{1,s}={c}_{6}({b}_{s}+{\delta}_{s}){E}_{s},\\ {N}_{2,r}={c}_{2}({b}_{r}+{\delta}_{r}){E}_{r}+{c}_{3}({b}_{s}+{\delta}_{s}){E}_{s},\phantom{\rule{0.38889em}{0ex}}{N}_{2,s}={c}_{s}({b}_{s}+{\delta}_{s}){E}_{s}+{c}_{4}({b}_{r}+{\delta}_{r}){E}_{r},\end{array}$$

and *c*_{1},…, *c*_{6} are defined in (6.15).

Assume _{0} < 1. Consider the limiting system of (3.4)–(3.6), where the total population densities are at their respective carrying capacities. Define the Liapunov function for the limiting system as follows:

$$V={z}_{1}{E}_{r}+{z}_{2}{E}_{a}+{z}_{3}{I}_{r}+{z}_{4}{I}_{a}+{z}_{5}{P}_{r}+{z}_{6}{P}_{a}+{z}_{7}{E}_{s}+{z}_{8}{E}_{b}+{z}_{9}{A}_{s}+{z}_{10}{A}_{b},$$

where

$$\begin{array}{l}{z}_{1}=({\gamma}_{s}+{b}_{s})({\delta}_{s}+{b}_{s}){b}_{r}{p}_{o}^{2}{\delta}_{r}(1-{\mathcal{R}}_{0}^{s})={z}_{2}\\ {z}_{3}=({\delta}_{r}+{b}_{r}){b}_{r}{p}_{o}^{2}({\gamma}_{s}+{b}_{s})({\delta}_{s}+{b}_{s})(1-{\mathcal{R}}_{0}^{s})\\ {z}_{4}=[{b}_{r}{\delta}_{r}{K}_{A}{\beta}_{a1}+{b}_{r}{p}_{o}^{2}({b}_{r}+{\delta}_{r})]({\gamma}_{s}+{b}_{s})({b}_{s}+{\delta}_{s})(1-{\mathcal{R}}_{0}^{s})+{b}_{r}{\delta}_{r}{\delta}_{s}{p}_{i}{K}_{A}{K}_{b}{\beta}_{a3}{\beta}_{b1}\\ {z}_{5}=({\delta}_{r}{p}_{i}{K}_{A}{\beta}_{a2}+{\delta}_{r}{K}_{r}{\beta}_{P}{p}_{o})({\gamma}_{s}+{b}_{s})({b}_{s}+{\delta}_{s}){p}_{o}(1-{\mathcal{R}}_{0}^{s})+{\delta}_{r}{\delta}_{s}{p}_{i}^{2}{K}_{A}{K}_{b}{\beta}_{a3}{\beta}_{b2}\\ {z}_{6}={\delta}_{r}[({p}_{i}{K}_{A}{\beta}_{a2}+{p}_{o}{K}_{r}{\beta}_{P}+{b}_{r}{K}_{A}{\beta}_{a2}){p}_{o}({b}_{s}+{\gamma}_{s})({b}_{s}+{\delta}_{s}(1-{\mathcal{R}}_{0}^{s})]+{\delta}_{r}{p}_{i}({b}_{r}+{p}_{i}){\delta}_{s}{K}_{A}{\beta}_{a3}{K}_{b}{\beta}_{b2}\\ {z}_{7}={p}_{i}{\delta}_{r}{\delta}_{s}{K}_{A}{\beta}_{a3}{p}_{o}{b}_{r}={z}_{8}\\ {z}_{9}=({b}_{s}+{\delta}_{s}){p}_{i}{\delta}_{r}{K}_{A}{\beta}_{a3}{p}_{o}{b}_{r}\\ {z}_{10}={K}_{A}{\beta}_{a3}{b}_{r}{p}_{o}{\delta}_{r}[({\delta}_{s}+{b}_{s})({\gamma}_{s}+{b}_{s})(1-{\mathcal{R}}_{0,1}^{s})+{p}_{i}({b}_{s}+{\delta}_{s})].\end{array}$$

The coefficients *z _{i}, i* = 1,…, 10, are positive if
${\mathcal{R}}_{0}^{s}<1$ and
${\mathcal{R}}_{0,1}^{s}={c}_{6}<1$. But

$$Y=({E}_{r},{E}_{a},{I}_{r},{I}_{a},{P}_{r},{P}_{a},{E}_{s},{E}_{b},{A}_{s},{A}_{b})$$

equals the zero vector. Differentiating *V* with respect to *t* along solution trajectories for the limiting system leads to

$$\frac{dV}{dt}=({\widehat{\mathcal{R}}}_{0}-1){p}_{o}({\gamma}_{r}+{b}_{r})({\delta}_{r}+{b}_{r})({\gamma}_{s}+{b}_{s})({\delta}_{s}+{b}_{s}){I}_{r}-\mathrm{\Omega}(Y),$$

where Ω (**0**) = **0**. A computer algebra system can be used to verify that Ω(*Y*) > 0 for *Y* ≥ **0**, *Y* ≠ **0**. It follows that *dV*/*dt* ≤ 0 for _{0} < 1 (or equivalently _{0} < 1). The Liapunov-Lasalle extension theorem (LaSalle, 1976) implies solutions to the limiting system in *D*_{1} approach the largest positively invariant subset of the set where *dV*/*dt* = 0. Hence, solutions to the limiting system converge to the DFE.

As noted previously, the *ω*-limit set of solutions to (3.4)–(3.6) is contained in *D*_{1}. The preceding argument shows that all solutions beginning in *D*_{1} converge to the DFE. It follows from the theory of asymptotically autonomous systems (Thieme, 1992, Theorem 4.1 or Castillo-Chavez and Thieme, 1995, p. 39) that if _{0} < 1, then the *ω*-limit set of (3.4)–(3.6) is the DFE.

Next, assume _{0} > 1. To prove existence and uniqueness of a positive EE, denote the vector-valued function *F*(*E _{r},E_{s}*) = (

$${F}^{\prime}(0,0)=\left(\begin{array}{cc}{c}_{1}+{c}_{2}& {c}_{3}\\ {c}_{4}& {c}_{5}+{c}_{6}\end{array}\right),$$

where *c*_{1},…,*c*_{6} are defined in (6.15). It is clear that *F*′(0, 0) is irreducible with spectral radius equal to _{0} > 1. With the aid of a computer algebra system, it can be shown that *F* is a monotone increasing function in each of its variables. These properties of *F* are sufficient to show the existence of a positive EE (Hethcote and Thieme, 1985, p. 209). To ensure a unique EE exists, a sufficient condition is that the function *F* is strictly sublinear: A vector-valued function *G*(*x*) = (*G _{i}*(

We define positive constants *ε _{r}* and

$${\epsilon}_{r}=\frac{(1-{h}^{2}){M}_{1,r}+(1-h){M}_{2,r}+h(1-h){N}_{2,r}^{2}{N}_{1,r}^{2}}{{M}_{3,r}({N}_{2,r}{N}_{1,r}{h}^{2}+{M}_{4,r}h+{p}_{i}{K}_{r}^{2}{b}_{r})}$$

and

$${\epsilon}_{s}=\frac{(1-{h}^{2}){M}_{1,s}+(1-h){M}_{2,s}+h(1-h){N}_{2,s}^{2}{N}_{1,s}^{2}}{{M}_{3,s}({N}_{2,s}{N}_{1,s}{h}^{2}+{M}_{4,s}h+{p}_{i}{K}_{s}^{2}{b}_{s})},$$

where

$$\begin{array}{l}{M}_{1,r}={p}_{i}{K}_{r}{N}_{1,r}{N}_{2,r}{N}_{3,r},\phantom{\rule{0.38889em}{0ex}}{M}_{1,s}={p}_{i}{K}_{s}{N}_{1,s}{N}_{2,s}{N}_{3,s},\\ {M}_{2,r}=[{P}_{i}{N}_{3,r}^{2}+{b}_{r}{N}_{2,r}^{2}]{p}_{i}{K}_{r}^{2},\phantom{\rule{0.38889em}{0ex}}{M}_{2,s}=[{p}_{i}{N}_{3,s}^{2}+{b}_{s}{N}_{2,s}^{2}]{p}_{i}{K}_{s}^{2},\\ {M}_{3,r}={p}_{i}{K}_{r}{N}_{3,r}+{N}_{2,r}{N}_{1,r},\phantom{\rule{0.38889em}{0ex}}{M}_{3,s}={p}_{i}{K}_{s}{N}_{3,s}+{N}_{2,s}{N}_{1,s},\\ {M}_{4,r}={b}_{r}{K}_{r}{N}_{2,r}+{p}_{i}{K}_{r}{N}_{3,r}\phantom{\rule{0.38889em}{0ex}}{M}_{4,s}={b}_{s}{K}_{s}{N}_{2,s}+{p}_{i}{K}_{s}{N}_{3,s},\end{array}$$

$${N}_{3,r}=({c}_{1}+{c}_{2})({b}_{r}+{\delta}_{r}){E}_{r}+{c}_{3}({b}_{s}+{\delta}_{s}){E}_{s},$$

and

$${N}_{3,s}=({c}_{5}+{c}_{6})({b}_{s}+{\delta}_{s}){E}_{s}+{c}_{4}({b}_{r}+{\delta}_{r}){E}_{r}.$$

Hence, the function *F* defined in (6.16)–(6.17) is strictly sublinear. All of the conditions of Theorem 2.1 inHethcote and Thieme (1985, p. 209) are satisfied. The function *F* has a unique fixed point which implies that a unique positive EE exists for the full system (3.4)–(3.6).

A CTMC model can be derived based on the ODE model (3.4)–(3.6), where variability is included in the birth, death, transmission, and movement processes (Allen, 2003). There are 38 events in the CTMC model, where a change occurs. Let *X*(*t*) = (*X*_{1}(*t*),…, *X*_{16}(*t*)) be a vector of 16 discrete random variables

$$X=({S}_{r},{E}_{r},{A}_{r},{P}_{r},{S}_{a},{E}_{a},{A}_{a},{P}_{a},{S}_{b},{E}_{b},{I}_{b},{R}_{b},{S}_{s},{E}_{s},{I}_{s},{R}_{s}),$$

where *X*_{1}(*t*) = *S _{r}*(

$$P(\mathrm{\Delta}XX)=\{\begin{array}{ll}{b}_{r}{\sum}_{i=1}^{4}{X}_{i}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{1}=1\hfill \\ {X}_{1}({\beta}_{I}{X}_{3}+{\beta}_{P}{X}_{4})\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{1}=-1,{a}_{2}=1\hfill \\ {d}_{r}{X}_{1}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{1}=-1\hfill \\ {\delta}_{r}{X}_{2}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{2}=-1,{a}_{3}=1\hfill \\ {d}_{r}{X}_{2}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{2}=-1\hfill \\ {\gamma}_{r}{X}_{3}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{3}=-1,{a}_{4}=1\hfill \\ {d}_{r}{X}_{3}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{3}=-1\hfill \\ {d}_{r}{X}_{4}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{r}=-1\hfill \\ {X}_{5}({\beta}_{a1}{X}_{7}+{\beta}_{a}{X}_{8}+{\beta}_{a3}{X}_{11})\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{5}=-1,{a}_{6}=1\hfill \\ {X}_{9}({\beta}_{b1}{X}_{7}+{\beta}_{b2}{X}_{8}+{\beta}_{b3}{X}_{11})\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{9}=-1,{a}_{10}=1\hfill \\ {b}_{s}{\sum}_{i=13}^{16}{X}_{i}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{13}=1\hfill \\ {\beta}_{A}{X}_{13}{X}_{15}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{13}=-1,{a}_{14}=1\hfill \\ {d}_{s}{X}_{13}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{13}=-1\hfill \\ {\delta}_{s}{X}_{14}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{14}=-1,{a}_{15}=1\hfill \\ {d}_{s}{X}_{14}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{14}=-1\hfill \\ {\gamma}_{s}{X}_{15}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{15}=-1,{a}_{16}=1\hfill \\ {d}_{s}{X}_{15}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{15}=-1\hfill \\ {d}_{s}{X}_{16}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{16}=-1\hfill \\ {\delta}_{a}{X}_{6}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{6}=-1,{a}_{7}=1\hfill \\ {\gamma}_{a}{X}_{7}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{7}=-1,{a}_{8}=1\hfill \\ {\delta}_{b}{X}_{10}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{10}=-1,{a}_{11}=1\hfill \\ {\gamma}_{b}{X}_{11}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{11}=-1,{a}_{12}=1\hfill \end{array},$$

where *d _{r}* =

$$P(\mathrm{\Delta}XX)=\{\begin{array}{ll}{p}_{i}{X}_{1}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{1}=-1,{a}_{5}=1\hfill \\ {p}_{o}{X}_{5}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{1}=1,{a}_{5}=-1\hfill \\ {p}_{i}{X}_{2}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{2}=-1,{a}_{6}=1\hfill \\ {p}_{o}{X}_{6}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{2}=1,{a}_{6}=-1\hfill \\ {p}_{i}{X}_{3}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{3}=-1,{a}_{7}=1\hfill \\ {p}_{o}{X}_{7}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{3}=1,{a}_{7}=-1\hfill \\ {p}_{i}{X}_{4}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{4}=-1,{a}_{8}=1\hfill \\ {p}_{o}{X}_{8}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{4}=1,{a}_{8}=-1\hfill \\ {p}_{i}{X}_{13}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{13}=-1,{a}_{9}=1\hfill \\ {p}_{0}{X}_{9}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{13}=1,{a}_{9}=-1\hfill \\ {p}_{i}{X}_{14}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{14}=-1,{a}_{10}=1\hfill \\ {p}_{o}{X}_{10}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{14}=1,{a}_{10}=-1\hfill \\ {p}_{i}{X}_{15}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{15}=-1,{a}_{11}=1\hfill \\ {p}_{o}{X}_{11}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{15}=1,{a}_{11}=-1\hfill \\ {p}_{i}{X}_{16}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{16}=-1,{a}_{12}=1\hfill \\ {p}_{o}{X}_{12}\mathrm{\Delta}t+o(\mathrm{\Delta}t),\hfill & {a}_{16}=1,{a}_{12}=-1\hfill \end{array}.$$

The probability of no change is one minus the sum of the 38 probabilities defined above and the probability of any other change is *o*(Δ*t*).

Differential equations for the joint probability function of *X* follow from these infinitesimal transition probabilities. The joint probability function is a solution of the forward Kolmogorov differential equations, an infinite system of differential equations. For multivariate processes, the form of these differential equations is often too complicated for analytical purposes. This is the case for our CTMC model with 16 random variables. The general form of the forward Kolmogorov differential equations depends on the set *S* = *U*×*U*×…×*U* = *U*^{16}, where *U* = {−1,0,1}. Let *h _{a}*(

$$\frac{d{\mathcal{P}}_{X}}{dt}=-{\mathcal{P}}_{X}\sum _{aS}$$

where
$X{\mathbb{Z}}_{+}^{16}$, _{+} = {0,1,2,…} (see e.g., Bailey, 1990).

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

- Abramson G, Kenkre VM. Spatiotemporal patterns in hantavirus infection. Phys Rev E. 2002;66:0011912-1–5. [PubMed]
- Abramson G, Kenkre VM, Yates TL, Parmenter RR. Traveling waves of infection in the hantavirus epidemics. Bull Math Biol. 2003;65:519–534. [PubMed]
- Allen LJS. An Introduction to Stochastic Processes with Applications to Biology. Prentice-Hall, Upper Saddle River; NJ: 2003.
- Allen LJS, Allen EJ, Jonsson CB. The impact of environmental variation on hantavirus infection in rodents. In: Gumel AB, Castillo-Chavez C, Mickens RE, Clemence DP, editors. Contemporary Mathematics Series, 410, Proceedings of the Joint Summer Research Conference on Modeling the Dynamics of Human Diseases: Emerging Paradigms and Challenges. AMS, Providence; RI: 2006a. pp. 1–15.
- Allen LJS, Langlais M, Phillips CJ. The dynamics of two viral infections in a single host population with applications to hantavirus. Math Biosci. 2003;186:191–217. [PubMed]
- Allen LJS, McCormack RK, Jonsson CB. Mathematical models for hantavirus infection in rodents. Bull Math Biol. 2006b;68:511–524. [PubMed]
- Anderson RM, May RM. Infectious Diseases of Humans, Dynamics and Control. Oxford Univ, Press; Oxford: 1991.
- Bailey NTJ. The Elements of Stochastic Processes with Applications to the Natural Sciences. John Wiley &; Sons; New York: 1990.
- Begon M, Hazel SM, Baxby D, Bown K, Cavanagh R, Chantrey J, Jones T, Bennett M. Transmission dynamics of a zoonotic pathogen within and between wildlife host species. Proc Roy Soc Lond B. 1999;266:1939–1945. [PMC free article] [PubMed]
- Bernshtein AD, Apekina NS, Mikhailova TV, Myasnikov YA, Khlyap LA, Korotkov YS, Gavrilovskaya IN. Dynamics of Puumala hantavirus infection in naturally infected bank voles (
*Clethrinomys glareolus*) Arch Virol. 1999;144:2415–2428. [PubMed] - Carleton MD, Musser GG. Order Rodentia. In: Wilson DE, Reeder DM, editors. Mammal Species of the World. 3. Vol. 2. Johns Hopkins University Press; Baltimore, MD: 2005. pp. 745–752.
- Castillo-Chavez C, Thieme HR. Asymptotically autonomous epidemic models. In: Arino O, Axelrod D, Kimmel M, Langlais M, editors. Mathematical Population Dynamics: Analysis of Heterogeneity, Theory of Epidemics. Vol. 1. 1995. pp. 33–49.
- CDC MMWR. Hantavirus pulmonary syndrome – United States: updated recommendations for risk reduction. 2002 July 26;51(RR09):1–12. [PubMed]
- Childs JE, Ksiazek TG, Spiropoulou CF, Krebs JW, Morzunov S, Maupin GO, Gage KL, Rollin PE, Sarisky J, Enscore RE, Frey JK, Peters CJ, Nichol ST. Serologic and genetic identification of
*Peromyscus maniculatus*as the primary rodent reservoir for a new hantavirus in the Southwestern United States. J Infect Dis. 1994;169:1271–280. [PubMed] - Chu YK, Milligan B, Owen RD, Goodin DG, Jonsson CB. Phylogenetic and geographical relationships of hantavirus strains in eastern and western Paraguay. Am J Trop Med Hyg. 2006;75:1127–1134. [PMC free article] [PubMed]
- Chua KB. Nipah virus outbreak in Malaysia. J Clin Virol. 2003;26:265–275. [PubMed]
- Daszak P, Cunningham AA, Hyatt AD. Emerging infectious disease of wildlife: threats to biodiversity and human health. Science. 2000;287:443–449. [PubMed]
- Delfraro A, Tome L, D’Elía G, Clara M, Achával F, Russi JC, Rodonz JR. Juquitiba-like hantavirus from two nonrelated rodent species, Uruguay. Emerg Infect Dis. 2008;14:1447–51. [PMC free article] [PubMed]
- Feng Z, Xu D, Zhao H. Epidemiological models for non-exponentially distributed disease states and applications to disease control. Bull Math Biol. 2007;69:1511–1536. [PubMed]
- Fernández Soto A, Mata Olmo R. Deforestatión y dinámica vegetal en un area de frontera agrícola del la Región Oriental del Paraguay. Revista Geonotas (Universidade Estadual de Maringá) 2001;5:1.
- Fiorello CV, Noss AJ, Deem SL. Demography, ecology, and pathogen exposure of domestic dogs in the Izozog of Bolivia. Conserv Biol. 2006;20:762–771. [PubMed]
- Glass GE, Livingston W, Mills JN, Hlady WG, Fine JB, Higgler W, Coke T, Frazier D, Atherley S, Rollin PE, Ksiazek TG, Peters CJ, Childs JE. Black Creek Canal Virus infection in
*Sigmodon hispidus*in southern Florida. Am J Trop Med Hyg. 1998;59:699–703. [PubMed] - Goodin DG, Paige R, Owen RD, Ghimire K, Koch DE, Chu Y-K, Jonsson CB. Microhabitat characteristics of
*Akodon montensis,*a reservoir for hantavirus, and hantaviral seroprevalence in an Atlantic forest site in eastern Paraguay. Journal of Vector Ecology Submitted to. [PubMed] - Hethcote HW. The mathematics of infectious diseases. SIAM Review. 2000;42:599–653.
- Hethcote HW, Thieme HR. Stability of the endemic equilibrium in epidemic models with subpopulations. Math Biosci. 1985;75:205–227.
- Holmes K. SARS-associated coronavirus. New Engl J Med. 2003;348:1948–1951. [PubMed]
- Karlin S, Taylor H. A First Course in Stochastic Processes. 2. Academic Press; NY: 1975.
- Klein SL, Bird BH, Glass GE. Sex differences in immune responses and viral shedding following Seoul virus infection in Norway rats. Am J Trop Med Hyg. 2001;65:57–63. [PubMed]
- Klingstrom J, Heyman P, Escutenaire S, Sjölander KB, De Jaegere F, Henttonen H, Lundkvist A. Rodent host specificity of European hantaviruses: evidence of Puumala virus interspecific spillover. J Med Virol. 2002;68:581–588. [PubMed]
- LaSalle JP. The Stability of Dynamical Systems. SIAM; Philadelphia: 1976.
- Lee HW, Lee PW, Baek LJ, Song CK, Seong IW. Intraspecific transmission of Hantaan virus, etiological agent of Korean hemorrhagic fever, in the rodent
*Apodemus agrarius*. Am J Trop Med Hyg. 1981;30:1106–1112. [PubMed] - Lloyd A. Realistic distributions of infectious periods in epidemic models. Theor Pop Biol. 2001a;60:59–71. [PubMed]
- Lloyd A. Destabilization of epidemic models with the inclusion of realistic distributions of infectious periods. Proc R Soc Lond B. 2001b;268:985–993. [PMC free article] [PubMed]
- McCormack RK, Allen LJS. Disease emergence in multi-host epidemic models. Math Med Biol. 2007;24:17–34. [PubMed]
- McIntyre NE, Chu Y-K, Owen RD, Abuzeineh A, De La Sancha N, Dick CW, Holsomback T, Nisbet RA, Jonsson C. A longitudinal study of Bayou virus, hosts, and habitat. Amer J Trop Med Hyg. 2005;73:1043–1049. [PubMed]
- Mills JN, Ksiazek TG, Ellis BA, Rollin PE, Nichol ST, Yates TL, Gannon WL, Levy CE, Engelthaler DM, Davis T, Tanda DT, Frampton JW, Nichols CR, Peters CJ, Childs JE. Patterns of association with mammals in the major biotic communities of the southwestern United States. Am J Trop Med Hyg. 1997;56:273–284. [PubMed]
- Mills JN, Yates TL, Ksiazek TG, Peters CJ, Childs JE. Long-term studies of hantavirus reservoir populations in the southwestern United States: Rationale, potential and methods. Emerg Infect Dis. 1999;5:95–101. [PMC free article] [PubMed]
- Nadin-Davis SA, Loza-Rubio E. The molecular epidemiology of rabies associated with chiropteran hosts in Mexico. Virus Res. 2006;117:215–226. [PubMed]
- Nel L, Jacobs J, Jaftha J, Meredith C. Natural spillover of a distinctly Canidae-associated biotype of rabies into an expanded wildlife host in southern Africa. Virus Genes. 1997;15:79–82. [PubMed]
- Padula P, Fogueroa R, Navarrette M, Pizarro E, Cadiz R, Bellomo C, Jofre C, Zaror L, Rodriguez E, Murua R. Transmission study of Andes hantavirus infection in wild sigmodontine rodents. J Virol. 2004;78:11972–11979. [PMC free article] [PubMed]
- Palma RE, Polop JJ, Owen RD. Hantavirus-host ecology in the Southern cone of South America: Argentina, Chile, Paraguay and Uruguay. Vector Borne and Zoonotic Disease 2009 Submitted to.
- Pardiňas UF, D’Elía G. The genus
*Akodon*(Muroidea: Sigmodontinae) in Misiones, Argentina. Mammalian Biology. 2003;68:129–143. - Parrish CR, Holmes EC, Morens DM, Park EC, Burke DS, Calisher CH, Laughlin CA, Saif LJ, Daszak P. Cross-species virus transmission and the emergence of new epidemic diseases. Microbiol Mol Biol Rev. 2008;72:457–470. [PMC free article] [PubMed]
- Redford KH, Eisenberg JF. Mammals of the Neotropics: The southern cone. University of Chicago Press; Chicago: 1992.
- Richomme C, Gauthier D, Fromont E. Contact rates and exposure to inter-species disease transmission in mountain ungulates. Epidemiol Inf. 2006;134:21–30. [PubMed]
- Sauvage F, Langlais M, Pontier D. Predicting the emergence of human han-tavirus disease using a combination of viral dynamics and rodent demographic patterns. Epidemiol Infect. 2007;135:46–56. [PubMed]
- Sauvage F, Langlais M, Yoccoz NG, Pontier D. Modelling hantavirus in fluctuating populations of bank voles: the role of indirect transmission on virus persistence. J Anim Ecol. 2003;72:1–13.
- Thieme HR. Convergence results and a Poincaré-Bendixson trichotomy for asymptotically autonomous differential equations. J Math Biol. 1992;30:755763.
- Torrez-Martinez N, Bharadwaj M, Goade D, Delury J, Moran P, Hicks B, Nix B, Davis JL, Hjelle B. Bayou virus-associated hantavirus pulmonary syndrome in Eastern Texas: identification of the rice rat,
*Oryzomys palustris*, as reservoir host. Emerg Infect Dis. 1998;4:105–111. [PMC free article] [PubMed] - van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002;180:29–48. [PubMed]
- Weidmann M, Schmidt P, Vackova M, Krivanec K, Munclinger P, Hufert FT. Identification of genetic evidence for dobrava virus spillover in rodents by nested reverse transcription (RT)-PCR and TaqMan RT-PCR. J din Microbiol. 2005;43:808–812. [PMC free article] [PubMed]
- Wesley CL. PhD Dissertation. Texas Tech University; Lubbock, Texas, U. S. A: 2008. Discrete-time and continuous-time models with applications to the spread of hantavirus in wild rodents and human populations.
- Wesley CL, Allen LJS, Jonsson CB, Chu YK, Owen RD. A discrete-time rodent-hantavirus model structured by infection and developmental stages. International Conference on Difference Equations and Applications July, 2006, Kyoto, Japan. Advanced Studies in Pure Mathematics. 2009;53:1–12.
- Wilson DE, Cole FR, Nichols JD, Rudran R, Foster MS. Standard Methods for Mammals. Smithsonian Institution Press; Washington: 1996. Measuring and Monitoring Biological Diversity.
- Wolf C, Langlais M, Sauvage F, Pontier D. A multi-patch epidemic model with periodic demography, direct and indirect transmission and variable maturation rate. Math Pop Studies. 2006;13:153–177.
- Yahnke CJ, Meserve PL, Ksiazek TG, Mills JN. Patterns of infection with Laguna Negra virus in wild populations of
*Calomys laucha*in the central Paraguayan chaco. Am J Trop Med Hyg. 2001;65:768–776. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |