PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2013; 8(6): e67267.
Published online Jun 28, 2013. doi:  10.1371/journal.pone.0067267
PMCID: PMC3696119
A Model for Chagas Disease with Oral and Congenital Transmission
Daniel J. Coffield, Jr.,1 Anna Maria Spagnuolo,2* Meir Shillor,2 Ensela Mema,3 Bruce Pell,4 Amanda Pruzinsky,5 and Alexandra Zetye2
1Mathematics Department, University of Michigan-Flint, Flint, Michigan, United States of America
2Department of Mathematics and Statistics, Oakland University, Rochester, Michigan, United States of America
3Department of Mathematical Sciences, New Jerseys Science & Technology University, University Heights, Newark, New Jersey, United States of America
4School of Mathematical & Statistical Sciences, Arizona State University, Tempe, Arizona, United States of America
5Chesapeake Research Consortium, U.S. EPA Chesapeake Bay Program Office, Annapolis, Maryland, United States of America
Herbert B. Tanowitz, Editor
Albert Einstein College of Medicine, United States of America
* E-mail: spagnuol/at/oakland.edu
Competing Interests: The authors have declared that no competing interests exist.
Conceived and designed the experiments: DJC AMS MS. Performed the experiments: DJC AMS EM BP AP AZ. Analyzed the data: DJC AMS MS EM BP AP AZ. Contributed reagents/materials/analysis tools: DJC AMS MS EM BP AP AZ. Wrote the paper: DJC AMS MS.
Received December 13, 2012; Accepted May 16, 2013.
This work presents a new mathematical model for the domestic transmission of Chagas disease, a parasitic disease affecting humans and other mammals throughout Central and South America. The model takes into account congenital transmission in both humans and domestic mammals as well as oral transmission in domestic mammals. The model has time-dependent coefficients to account for seasonality and consists of four nonlinear differential equations, one of which has a delay, for the populations of vectors, infected vectors, infected humans, and infected mammals in the domestic setting. Computer simulations show that congenital transmission has a modest effect on infection while oral transmission in domestic mammals substantially contributes to the spread of the disease. In particular, oral transmission provides an alternative to vector biting as an infection route for the domestic mammals, who are key to the infection cycle. This may lead to high infection rates in domestic mammals even when the vectors have a low preference for biting them, and ultimately results in high infection levels in humans.
Chagas disease is caused by infection with the parasite Trypanosoma cruzi and is a major source of suffering throughout Latin America. The disease leads to organ deformity and early death in about one third of the 8–10 million individuals infected, [1], [2]. Vector transmission by reduviids is largely responsible for the spread of T. cruzi, with some particular species specialized in domestic infection cycles. Other modes of transmission include blood tranfusions, organ transplants, oral transmission, and congenital transmission [3][6].
Although various drugs are under development and testing, current control of the transmission of Chagas disease remains largely based on vector control and blood-bank screening [7][10]. In particular, the Southern Cone Initiative was implemented in the 1990s with the goal of interrupting the transmission of Chagas disease in South American countries through insecticide spraying and blood screening [10][12]. This program has led to a dramatic decrease in transmission in several countries in South America, with some regions now reporting a considerable drop in infections from Triatoma infestans, the primary vector, and transmission virtually at zero [13], [14]. Additional control measures are treatment for acute Chagas disease and for congenital transmission cases [15].
While insecticide spraying for Chagas vectors has led to a significant decrease in new infections, improved housing, better drugs, and an effective vaccine are needed [16]. In particular, T. cruzi infection is likely to remain endemic in sylvatic hosts despite spraying efforts and neither insect control nor current drug treatment is optimal for this disease because of the long life span of infected human hosts, triatomine insecticide resistance, and the ease with which protozoans develop drug resistance [17][20].
Mathematical models for studying Chagas disease dynamics with seasonal insecticide spraying were presented in [21], [22]. In this work, we enhance the model in [21], adding the effects of congenital transmission in infected humans and infected dogs and excluding spraying. We also account for oral transmission by allowing the domestic mammals to consume the vectors, as observed experimentally in [23]. The predation term involves a density-dependent consumption rate in the form of a Holling Type II response, similar to that in [4]. There are other likely routes of oral transmission in domestic mammals such as ingesting feces-contaminated food and water or licking feces-contaminated fur. Though the vector consumption is derived on the basis of the animals preying on the vectors, the consumption term can still in some sense account for these other modes since the consumption is dependent on the vector density.
The primary aim of this work is to investigate the significance of the following modes of disease transmission relative to vector biting: 1.) oral transmission due to predation, and 2.) congenital transmission. In particular, we want to know if these transmission modes play a significant role in human infection and have implications for disease control. To this end, we study the additional transmission modes as enhancements to the model in [21] for comparison.
The model we present here consists of four nonlinear differential equations that describe the domestic transmission of the disease by predicting the total number of vectors, infected vectors, infected humans, and infected domestic mammals. In [24], a mathematical model for a small population (one household) was described, whereas the model in this work considers a large population (one village). We note that models were recently used in [25], [26] to study control issues of the non-domiciliated Chagas disease vectors in the Yucatan Peninsula, Mexico. Those models use difference equations while our model is a classical nonlinear dynamical system. The effectiveness of different control measures was recently studied in [27] using a dynamical system model. Additional modeling and field results can be found in [28][30].
The Methods section gives a detailed description of the model and its parameters. In Results, the model is used to produce simulations of the populations in a hypothetical rural village over thirty years using a baseline parameter set, provided in Table 1. Additional simulations are performed to investigate the model's sensitivity to various parameters. The Results section also compares simulations of the model presented in this work with the model in [21], which does not consider oral transmission in dogs nor congenital transmission. We conclude with the Discussion.
Table 1
Table 1
The model parameters and the baseline simulation values.
In this section, we present a new model for Chagas disease dynamics in a village. Information on pertinent biological processes can be found in [24]. The model represents the overall dynamics of the populations of vectors, infected vectors, infected humans, and infected domestic animals (mammals only) – referred to as 'dogs' in our model. In a rural village, there are also non-mammals that act as sources of blood meals but that are not hosts for T. cruzi (i.e., cannot become infected). They will be referred to as 'chickens' in our model.
We consider a relatively large representative rural village in South America so that a differential equations model is appropriate. Let the number of humans in the village be An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e065.jpg, the number of domestic mammals (dogs) be An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e066.jpg, and the number of domestic non-mammals (chickens) be An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e067.jpg. These are taken as constants for the sake of simplicity, since we consider a modest time period of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e068.jpg years. We denote by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e069.jpg the number of carrier insects living in the houses in the village at time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e070.jpg, the number of infected insects by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e071.jpg, the number of infected humans by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e072.jpg, and the number of infected dogs by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e073.jpg. Each non-infected population, excluding the chickens, is assumed to be susceptible.
We now describe the rate of change of each population in the village per day. The growth rate of the vectors depends on the successful hatching of eggs. As in [22], the egg hatching rate, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e074.jpg, is delayed due to the gestation time of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e075.jpg days. The hatching rate is a product of the following terms: the ratio of adult female vectors to total vectors; the number of eggs laid by an adult female per bite; the successful hatching rate of eggs after An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e076.jpg days; the total blood meal supply (in human factors), An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e077.jpg; and the delayed seasonal biting rate An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e078.jpg days prior to hatching, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e079.jpg, in units of bites per vector per day per human factor (Table 1). By extracting seasonal data from Castanera et al. [38], we obtain the following values: the ratio of adults to the entire population of triatomines is approximately 11/365 (so we take half that number to be adult females), the fraction of eggs that survive is 0.83, and the number of (eggs/bite)/(fed female) is 20 [37], [38]. So, the form of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e080.jpg is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e081.jpg
Note that we are assuming An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e082.jpg follows the seasonality of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e083.jpg.
Following [24], we use 'human factors' as the unit for the total blood supply in the following way: each human represents one human factor, each mammal An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e084.jpg represents An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e085.jpg human factors, and each non-mammal An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e086.jpg represents An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e087.jpg human factors; then, the total blood supply is given by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e088.jpg We use the standard notation An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e089.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e090.jpg to denote the positive and negative parts, respectively, of a function An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e091.jpg. Vector growth is then modeled by the following delay logistic term
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e092.jpg
as in [21], where An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e093.jpg is the carrying capacity of vectors in the village houses. The term An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e094.jpg is the rate at which vectors hatch at time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e095.jpg from eggs laid at time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e096.jpg. The expression An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e097.jpg represents the fraction of the food supply that was available to the female vectors at time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e098.jpg. This assumes that if the vector population at time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e099.jpg days prior to the current time An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e100.jpg was above the carrying capacity, then no eggs are laid. The death rate of the vector population depends on the following three factors: natural mortality, death due to overcrowding or growth beyond the carrying capacity, and death due to being eaten by the dogs. The natural death rate coefficient of triatomines is An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e101.jpg and the coefficient of the death rate of vectors above the carrying capacity is An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e102.jpg. These rates are assumed to be periodic functions with a period of one year and are included by adding the following term to the vector growth equation:
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e103.jpg
We now consider the death rate of vectors due to consumption by dogs. We use a Holling Type II functional response for the consumption term, defining the per dog consumption rate as
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e104.jpg
Here, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e105.jpg is the maximum number of vectors consumed by one dog per day and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e106.jpg is the vector number at which dogs consume at the rate of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e107.jpg vectors per day. This term is similar to one used in [4] where vector consumption by wild animals was considered. Thus, the rate of change in the total vector population within the village is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e108.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e109.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e110.jpg
Next, we consider the growth rates of the infected populations An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e111.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e112.jpg, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e113.jpg. We denote by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e114.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e115.jpg the probabilities of a vector becoming infected by biting an infected human or an infected dog, respectively. Also, the number of bites per vector per day is given by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e116.jpg, where An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e117.jpg is the same seasonal biting term used to define An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e118.jpg. Since the proportions of those bites that occur on infected humans and infected dogs are An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e119.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e120.jpg, respectively, the growth rate of infected vectors is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e121.jpg
We assume that the natural death rate of the infected vectors is also An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e122.jpg, i.e., carrying the parasites does not affect their life span. Similarly, the death rate due to growth beyond the carrying capacity is An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e123.jpg.
The death rate of the infected vectors due to predation by dogs is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e124.jpg
Collecting the terms above, we get the following rate of change for the infected vectors:
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e125.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e126.jpg
We turn now to the infected humans. When bitten by an infected vector, a susceptible human becomes infected with probability An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e127.jpg. As before, each vector is biting at a rate of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e128.jpg bites per day, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e129.jpg is the fraction of bites that are on susceptible humans. Thus, the growth rate of infected humans is given by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e130.jpg. The natural mortality rates for infected humans and susceptible human are denoted by An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e131.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e132.jpg, respectively. We assume that human reproduction is independent of infection status, since infected humans typically live beyond reproductive ages. Therefore, the assumption that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e133.jpg is constant implies that the birth rate for all humans is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e134.jpg
Here, the birth rate depends on An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e135.jpg, which is a consequence of the two preceding assumptions. However, we note that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e136.jpg is small so that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e137.jpg does not vary enough for this dependence to meaningfully affect the model. Furthermore, we investigated multiple other assumptions for the human birth rate, including scenarios for different birth rates of infected humans versus uninfected humans, and the simulation results were virtually indistinguishable in all cases.
We define An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e138.jpg to be the ratio, amongst babies born to infected mothers, of infected babies to the total number of babies. Therefore, the rate of change of infected humans is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e139.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e140.jpg
Under the assumption that both infected and uninfected dogs reproduce at the same rate, the rate of change of infected dogs is similar to that of infected humans, except that we must account for infection caused by the consumption of infected vectors. To this end, let An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e141.jpg be the probability that an uninfected dog becomes infected when bitten by an infected vector, let An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e142.jpg be the probability that an uninfected dog becomes infected after eating an infected vector, and note that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e143.jpg be the susceptible dog population. We also use An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e144.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e145.jpg to denote the natural mortality rates for infected dogs and susceptible dogs, respectively. Then, the rate of change of infected dogs is
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e146.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e147.jpg
where An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e148.jpg is the probability that an infected dog passes the infection to its offspring congenitally, and the congenital transmission term is similar to the term in the infected humans equation because An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e149.jpg is constant.
To complete the model we prescribe the initial values of the respective populations: An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e150.jpg, together with.
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e151.jpg
.
These equations and conditions form a mathematical model for the domestic dynamics of Chagas disease with oral and congenital transmission:
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e152.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e153.jpg
(1)
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e154.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e155.jpg
(2)
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e156.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e157.jpg
(3)
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e158.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e159.jpg
(4)
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e160.jpg
A mathematical equation, expression, or formula.
 Object name is pone.0067267.e161.jpg
(5)
The coefficient functions An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e162.jpg, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e163.jpg are one-year periodic since they are seasonally dependent. Also, it is natural to assume that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e164.jpg. Note that the delay differential equation (1) for the total vector population is not coupled to the other equations, so it can be solved independently.
In this section, we study the effects of the congenital and oral transmission terms, various model parameters, and compare the current model with the one in [21]. The Adams-Bashforth Fourth-Order Method was implemented in gfortran [31] and verified using Wolfram Mathematica [32]. The figures were generated with Wolfram Mathematica [32].
Baseline Case
We first compare the simulation results in a baseline case, similar to the one described in [21], with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e165.jpg houses, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e166.jpg humans, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e167.jpg dogs, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e168.jpg chickens in a respresentative village. We use a similar parameter set to define our baseline case with the only difference being a new value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e169.jpg for An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e170.jpg See Vector Biting Preference Studies for an explanation. We also add the parameters for congenital and oral transmission. See Table 1 for baseline parameters. The simulation timeframe is An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e171.jpg years. The values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e172.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e173.jpg significantly affect the model, so we investigate them further in the next section.
The probabilities for congenital transmission are taken from data in [5], [15], [40]. Since data for predation rates of domestic mammals is not currently available, we estimate the maximum value of the per dog vector consumption rate to be about An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e174.jpg vector per week. This estimate is within the range of values estimated for wild mammals in [40]. For the probability of dog infection via oral consumption of an infected vector, we use the value found in [40], which is a weighted average of values found from experiments with raccoons and opossums. An appropriate value of the parameter An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e175.jpg is not known. We chose it to be An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e176.jpg meaning that when a village house has An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e177.jpg vectors in it, then dogs consume vectors at a rate of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e178.jpg We note that although the model is sensitive to the total number of vectors, it is not sensitive to the value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e179.jpg Also, Figure 1 contains the graph of the seasonally dependent values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e180.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e181.jpg. See Table 1 and the references therein for the data used to estimate these values and note that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e182.jpg is proportional to An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e183.jpg.
Figure 1
Figure 1
Vector growth and mortality coefficients.
Figure 2 shows the seasonal oscillations of the total vectors, infected vectors, infected humans, and infected dogs for the models with and without congenital and oral transmission. Notice that the new model produces lower peaks for the total number of vectors whereas the infected human and infected dog populations take on larger values due to the new transmission modes. The infected vector population in the new model attains higher values early on, but within 10 years is nearly identical to the old model.
Figure 2
Figure 2
Simulation results comparing models.
Figure 3 shows model simulations for the number of vectors, infected vectors, infected humans, and infected dogs over 30 years using the the baseline parameter set with higher initial conditions. We note that all of the populations largely stabilize around a central value and then seasonally oscillate around that value. This is consistent with the endemic nature of the disease in rural villages that do not engage in control measures. Note that the oscillations in the vector populations are much larger than those of the humans and domestic mammals, because of the shorter life span of vectors.
Figure 3
Figure 3
Higher initial conditions.
Vector Biting Preference Studies
Recently, experimental work in [33] suggests that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e186.jpg is about seven times An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e187.jpg This strongly contradicts earlier estimates found in [34] and used in our previous work [21], [22]. We point out that this new study only compares vector preference between dogs and chickens, but does not consider vector preference for humans. Thus, vector preference between humans, dogs, and chickens remains unclear. Moreover, our work indicates that the preference factors substantially affect the dynamics of the infected populations. To address this, we perform simulations over various ranges of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e188.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e189.jpg See Figure 4.
Figure 4
Figure 4
Infected human and vector populations at year 30.
Note that for fixed An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e192.jpg and increasing An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e193.jpg infection in humans increases and then decreases, attaining its peak for An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e194.jpg between An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e195.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e196.jpg Although not shown, the number of infected dogs increases as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e197.jpg increases, as expected. This leads to a sharp increase in infected vectors for low An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e198.jpg values before leveling off later on. See Figure 4. The initial increase in human infections is due to the steep increase in infected vectors which more than compensates for the increased vector preference for dogs. Because the number of infected vectors levels off for large values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e199.jpg it follows that human infection decreases as the vectors primarily bite the dogs. Also, for each fixed An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e200.jpg human infections decrease as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e201.jpg increases. This is expected since chickens are not infective and do not contribute to the infection cycle. Hence a higher An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e202.jpg just diverts bites from dogs and humans.
In most of the following simulations, we allow An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e203.jpg to vary and consistently choose An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e204.jpg due to the relationship found in [33]. In simulations where An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e205.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e206.jpg are fixed, we choose An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e207.jpg to stay consistent with baseline studies in previous work and correspondingly set An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e208.jpg
Effects of Oral and Congenital Transmission
As expected, congenital transmission increases overall infection in humans and dogs. As a result, the number of infected vectors also increases, but not significantly. Estimates of congenital transmission probabilities are readily available and generally in the range of 2–10% [5], [15]. For values in this range, the effects of congenital transmission are modest and close to linear as a function of transmission probabilities, Figure 5. We note that a recent article reported dramatically higher congenital transmission probabilities in mice (33–66%) [35]. Simulations of our model with congenital transmission probabilities up to 50% reveal a continued, near-linear effect of vertical transmission on infected humans.
Figure 5
Figure 5
Infected humans at year 30 as a function of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e209.jpg, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e210.jpg, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e211.jpg.
The oral transmission is more complicated. In the case shown in Figure 5 with a baseline An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e215.jpg the number of infected humans at year 30 changes by at most 4 as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e216.jpg ranges from 0 to 0.15 with a peak attained at approximately 0.07. We note that initally a higher vector consumption rate results in more infections in both humans and dogs, as expected. However, as the consumption rate increases further, the number of human infections actually declines since there are fewer infected vectors feeding on humans, as depicted in Figure 6. Also, the total number of vectors steadily decreases, since more vectors are being consumed as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e217.jpg increases.
Figure 6
Figure 6
Populations at year 30 as functions of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e218.jpg.
Next, we consider oral transmission with higher and lower values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e220.jpg and, correspondingly, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e221.jpg Figures 7 and and88 show the four populations as functions of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e222.jpg with all other parameters except An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e223.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e224.jpg set to their baseline values. We point out that varying An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e225.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e226.jpg had the same effect on the populations as in Figure 5. In the case where An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e227.jpg the effects of oral transmission are more dramatic with the number of infected humans at year 30 increasing by approximately 30 as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e228.jpg ranges from 0 to 0.15. In comparison, when An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e229.jpg the number of infected humans starts at a significantly lower number and decreases. We observed similar trends with larger An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e230.jpg
Figure 7
Figure 7
Populations at year 30 as functions of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e231.jpg with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e232.jpg.
Figure 8
Figure 8
Populations at year 30 as functions of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e235.jpg with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e236.jpg.
Since the vector preference numbers do not affect the biting rate of the vectors, it follows that the total number of vectors is the same for a fixed value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e239.jpg So, in each of Figures 6, ,7,7, ,8,8, the graphs for the total number of vectors are identical. However, all infected populations have dramatically different outcomes. More specifically, in all cases, the number of infected dogs increases as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e240.jpg increases, because the dogs are eating more infected vectors. Also, dog infection increases with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e241.jpg since An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e242.jpg.
Now, we consider the number of infected humans for the different values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e243.jpg. It is interesting to note that the curves for the number of infected humans and infected vectors have the same shape in each instance of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e244.jpg We observe that for An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e245.jpg roughly less than An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e246.jpg (and therefore, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e247.jpg roughly less than An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e248.jpg), the curves for infected humans and infected vectors are increasing with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e249.jpg in An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e250.jpg. However, for An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e251.jpg roughly in the interval An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e252.jpg, the curves initially increase before decreasing with the peak moving to the left for higher values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e253.jpg. Also, it appears that the peak of each of the An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e254.jpg curves trails closely behind the peak of the corresponding An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e255.jpg curves since vector infection drives human infection.
We now explain the decline after the peak for infected vectors, and correspondingly, for infected humans. Recall that An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e256.jpg, so a higher An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e257.jpg value means a stronger vector preference for dogs and a resulting higher number of infected vectors. However, since the total number of vectors is independent of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e258.jpg, the ratio An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e259.jpg increases as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e260.jpg increases. Also, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e261.jpg is a decreasing function of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e262.jpg, making An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e263.jpg even higher for large An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e264.jpg. Thus, more of the vectors removed through dog predation are infected vectors for higher An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e265.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e266.jpg values. Initially as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e267.jpg increases, more dogs become infected which in turn leads to a higher number of infected vectors. But, as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e268.jpg increases further, more infected vectors are being eliminated through vector consumption than are being added through new infections.
Finally, for An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e269.jpg greater than about 1, the number of infected humans and infected vectors declines as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e270.jpg increases from 0 to 0.15. As in the latter part of the previous case, the high vector preference for dogs leads to a large infected vector population and a high An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e271.jpg ratio. The predation on vectors removes such a high proportion of infected vectors in this case that the number of infected vectors strictly decreases as a An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e272.jpg increases. Dog infection still increases with An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e273.jpg, leading to new vector infections, but not enough to overcome those being removed.
Past work has shown that dogs play an important role in the infection cycle. The new model, and particularly the inclusion of oral transmission, demonstrates an even more severe, negative impact of the dogs. In fact, Figures 7 and and88 reveal that high levels of human infection occur in both of the following cases: a high level of dog oral transmission coupled with a low vector biting preference for dogs and chickens; a low level of dog oral transmission coupled with a high vector preference for dogs and chickens. The key observation is that even if vector preference for dogs is low, dogs still become sufficiently infected through oral transmission to drive infection in vectors and humans. Furthermore, for low levels of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e274.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e275.jpg, dog infection is actually caused more by oral consumption of vectors than by vector biting. See Figure 9. The figure further reveals that human infection remains high at low An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e276.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e277.jpg values even when dogs cannot be infected through biting.
Figure 9
Figure 9
Infected humans and infected dog populations at year 30 in different scenarios.
Sensitivity to Consumption
As shown above and in Figure 10, the model is sensitive to the value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e279.jpg, but only at low values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e280.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e281.jpg. For example, when An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e282.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e283.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e284.jpg ranges from 0 to 0.15, the number of infected humans changes by about 30. However, the model is not very sensitive to An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e285.jpg for larger values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e286.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e287.jpg In fact, for any value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e288.jpg greater than the baseline value of 0.35, the number of infected humans at year 30 changes by at most 4 as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e289.jpg varies in the same range (0 to 0.15).
Figure 10
Figure 10
Infected humans at year 30 as a function of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e290.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e291.jpg.
Varying Dog and Chicken Levels
We now investigate the effects of changing the number of chickens and dogs in the village and their preference factors, while setting all other parameters to their baseline values. For fixed An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e295.jpg, the number of infected humans increases as the number of chickens increases, Figure 11. Although chickens are not infectious, a higher chicken population means a higher blood supply is available to the vectors, resulting in higher vector prevalence. This leads to more bites on dogs and humans and increased infected populations.
Figure 11
Figure 11
Infected humans at year 30 with different numbers of chickens and dogs.
Similarly, for a fixed value of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e297.jpg the number of infected humans increases as the number of dogs increases, Figure 11. The effect on An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e298.jpg of increasing the number of dogs is more dramatic than the effect of increasing the number of chickens because dogs are infectious. As can be seen in both cases, the number of infected humans initially increases with increasing An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e299.jpg before decreasing once An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e300.jpg is beyond about An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e301.jpg.
Larger Blood Supply Cases
We now consider a more realistic village, see [36], where each house in the village has An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e302.jpg humans, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e303.jpg dogs, and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e304.jpg chickens. All other parameters are taken to be the baseline values. The primary effect of the increased populations is a larger blood supply available to the vectors. Correspondingly, the total number of vectors is significantly higher than in the baseline case. See Figure 12. All of the infected populations are significantly higher and a greater percentage of all populations become infected.
Figure 12
Figure 12
Populations with different blood supplies.
This work presents a new model with seasonally dependent coefficients for the domestic transmission of Chagas disease, building upon the work in [21]. The model includes transmission through vector biting along with the new infection routes of congenital transmission in humans and domestic mammals as well as oral transmission in domestic mammals through consumption of infected vectors. Simulations indicate that oral transmission plays an important role in the infection cycle while the effect of congenital transmission is more limited.
The inclusion of congenital transmission in humans directly leads to more human infections. However, because both the birth rate and the probability of the infection passing from mother to child are relatively low, vertical transmission in humans leads to only a few new infections over 30 years in a village of 400. We note that the choice to have a constant human population limits the model's flexibility in choosing the birth rate and thus it is perhaps artificially low. Clearly, the effects of congenital transmission in humans would be more severe in villages with higher birth rates. On the other hand, since dogs have a relatively higher birth rate than humans, congenital transmission in dogs might be expected to substantially influence the number of infected dogs and indirectly increase the number of infected humans. However, the dogs become so easily infected by other transmission routes that the addition of congenital transmission to the model does not significantly affect dog infections. In particular, simulations show that the infected dog population quickly stabilizes at a high level of infection (usually around 60%–90% infected). The inclusion of congenital transmission does slightly increase the peak dog infection level, but not enough to substantially affect human infection.
The effects of oral transmission in dogs are more dramatic and complex than that of congenital transmission. Furthermore, the significance of the oral transmission is strongly tied to the vector biting preference numbers. For high values of An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e305.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e306.jpg, most of the dogs become infected through biting without oral transmission, though increasing the dog's consumption rate does moderately increase the number of infected dogs (Figure 8). However, in this case, the increased consumption causes a decline in the infected vector population, and correspondingly, in the infected human population. These declines are small though, which demonstrates that vector biting drives infection when the vectors strongly prefer to feed on dogs.
Alternatively, oral transmission is the driving force behind the infection cycle when An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e307.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e308.jpg are low (Figure 7). In this case, vector biting alone leads to only about An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e309.jpg of the dogs being infected after 30 years. However, adding oral transmission dramatically increases all of the infected populations and the level of infection is very sensitive to the dog's maximal consumption rate, An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e310.jpg. As noted in Results, the number of human infections in a representative village of 400 increases by about 30 infections over 30 years of simulation as An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e311.jpg is increased from 0 to the baseline value of one per week. This means that with significant oral transmission in dogs, human infection will remain high even if the vectors have a low preference for biting the dogs. Since the probability of transmission from dogs to vectors is significantly higher than the probability of transmission from humans to vectors, see Table 1, we know that the dogs are primarily responsible for infecting the vectors. In turn, the size of the infected vector population directly drives the number of human infections. So, when An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e312.jpg is low, oral transmission is the key route of dog infection, whereas biting is more important when An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e313.jpg is high. In either case, our simulations show that the level of dog infection remains high, resulting in a substantial number of human infections. This result is noteworthy because it suggests that the disease will persist at high levels even if measures are taken to deter the vectors from biting the dogs, e.g. using insecticide collars.
It is well-known and widely reported that domestic mammals are a major player in the infection cycle and the main reservoir of the parasite. Our work strengthens these conclusions and further demonstrates the need to remove mammals (dogs, cats, etc.) from the domestic settings. This is not a new control recommendation, but our simulations, and particularly the inclusion of oral transmission, show the fundamental, negative role the mammals play in causing human infections over a wide range of vector biting preferences. In fact, our simulations show that the infections persist endemically even with a small number of dogs (0.2 dogs/home). This is consistent with [34], where it was found that a 100% effective control method on at least 88% of the dogs would be needed to achieve a basic reproductive number smaller than 1.
In our model, the infection dies out with no dogs. And even if a total removal of domestic mammals is infeasible, reducing their numbers will likely lead to fewer human infections. As shown in Results, the number of infected humans only increases with the number of dogs in our model. This contrasts with the model in [24], where it was found that human infection declines when each household has more than two dogs, allowing the dogs to sufficiently divert vectors away from the humans. In our model, more dogs means a higher blood supply available to the vectors, and correspondingly, more vectors and higher infected populations. However, we note that we do see a similar decline in humans infections as the vector biting preference for dogs increases beyond about 2.5 human factors. See Figure 4.
We note that this work uses a predation term to model the oral transmission in dogs, and this term may not appropriately account for other likely routes of oral transmission such as licking of feces-contaminated fur and ingesting feces-contaminated food or water.
A weakness of this work is that the parameters are coming from different studies. However, data from the same studies do not currently exist. Although the simulation results are highly sensitive to the vector biting preferences for dogs and chickens, which are largely unknown, the primary control implication–eliminating domestic mammals–is independent of these parameters. That is, domestic mammals should be removed from the homes even if vector preference for them is low. Furthermore, due to recent work in [33], we have drastically changed the relationship between An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e314.jpg and An external file that holds a picture, illustration, etc.
Object name is pone.0067267.e315.jpg as compared to our previous work [21], [22]. Yet, the overall dynamics are very similar to our previous work and dogs remain the driving force of the infection cycle.
The dynamics of Chagas disease are indeed complex, so in addition to validation of the model, there are several open issues that are of interest for further study. First, we point out that human infection levels might be higher than in the simulations, because we did not include blood transfusions or oral transmission in humans, though recent outbreaks (e.g. feces in juice) indicate that the latter may be a significant source of human infection. Additionally, the effects of wild vectors and the disease in the wildlife were not investigated. We also did not directly consider vector mortality due to consumption by domestic non-mammals. Finally, the dynamics for the total human and domestic animal populations could be studied, which would allow for investigation of immigration and contact between neighboring villages.
Acknowledgments
The authors would like to thank the reviewers for their insightful comments and suggestions which improved the manuscript.
Funding Statement
National Science Foundation REU Grant DMS 0649099 (www.nsf.gov). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Organizacion Panamericana de la Salud (2006) Estimacion cuantitativa de la enfermedad de Chagas en las Americas. Washington DC: PAHO Publishing.
2. Schofield CJ, Jannin J, Salvatella R (2006) The future of chagas disease control. Trends Parasitol 22: 583–588. [PubMed]
3. Camandaroba ELP, Pinheiro Lima CM, Andrade SG (2002) Oral transmission of Chagas disease: Importance of Trypanosoma cruzi biodeme in the intragastric experimental infection. Rev Inst Med Trop Sao Paulo 44: 97–103. [PubMed]
4. Kribs-Zaleta CM (2006) Vector Consumption and Contact Process Saturation in Sylvatic Transmission of T. Cruzi. Math Popul Stud 13(2006): 135–152.
5. Kribs-Zaleta CM (2010) Alternative Transmission Modes for Trypanosoma Cruzi. Math Biosci Eng 7(3): 657–673. [PubMed]
6. Roellig DM, Ellis AE, Yabsley MJ (2009) Oral transmission of Trypanosoma cruzi with opposing evidence for the theory of carnivory. J Parasitol 95: 360–364. [PMC free article] [PubMed]
7. Garg N, Bhatia N (2005) Current status and future prospects for a vaccine against American trypanosomiasis. Expert Rev Vaccines 4(6): 867–880. [PubMed]
8. Rodriques CJ, de Castro SL (2002) A critical review on Chagas disease chemotherapy. Mem Inst Oswaldo Cruz 97(1): 3–24. [PubMed]
9. Schmunis GA (1999) Prevention of transfusional Trypanosoma cruzi infection in Latin America.Mem Inst Oswaldo Cruz (Suppl 1): 93–101. [PubMed]
10. Silveira A, Vinhaes M (1999) Elimination of vector-borne transmission of Chagas disease. Mem Inst Oswaldo Cruz (Suppl 1): 405–411. [PubMed]
11. Schmunis GA (1999) Prevention of transfusional Trypanosoma cruzi infection in Latin America. Mem Inst Oswaldo Cruz (Suppl 1): 93–101. [PubMed]
12. Schmunis GA, Cruz JR (2005) Safety of the blood supply in Latin America. Clin Microbiol Rev 18: 12–29. [PMC free article] [PubMed]
13. Dias JC (2007) Southern Cone Initiative for the elimination of domestic populations of Triatoma infestans and the interruption of transfusional Chagas disease. Historical aspects, present situation, and perspectives. Mem Inst Oswaldo Cruz. 102 (Suppl 1)11–18. [PubMed]
14. Massad E (2008) The elimination of Chagas disease from Brazil Epidemiol Infect. 136: 1153–1164. [PubMed]
15. Gürtler RE, Segura EL, Cohen JE (2003) Congenital transmission of Trypanosoma cruzi infection in Argentina, Emerg Infect Dis 9(1). Available: http://www.cdc.gov/ncidod/EID/vol9no1/02-0274.htm Accessed 29 May 2013. [PMC free article] [PubMed]
16. Moncayo A (2003) Chagas disease: current epidemiological trends after the interruption of vectorial and transfusional transmission in the Southern Cone countries. Mem Inst Oswaldo Cruz 98: 577–591. [PubMed]
17. Buckner FS, Wilson AJ, White TC, Van Voorhis WC (1998) Induction of resistance to azole drugs in Trypanosoma cruzi. Antimicrob Agents Chemother 42: 3245–3250. [PMC free article] [PubMed]
18. Diotaiuti L, Pereira AS, Loiola CF, Fernandes AJ, Schofield JC, et al. (1995) Inter-relation of sylvatic and domestic transmission of Trypanosoma cruzi in areas with and without domestic vectorial transmission in Minas Gerais. Mem Inst Oswaldo Cruz 90: 443–448. [PubMed]
19. Murta SM, Romanha AJ (1998) In vivo selection of a population of Trypanosoma cruzi and clones resistant to benznidazole. Parasitology 116(Pt 2): 165–171. [PubMed]
20. Vassena CV, Picollo MI, Zerba EN (2000) Insecticide resistance in Brazilian Triatoma infestans and Venezuelan Rhodnius prolixus. Med Vet Entomol 14: 51–55. [PubMed]
21. Spagnuolo AM, Shillor M, Kingsland L, Thatcher A, Toeniskoetter M, et al. (2012) A Logistic Delay Differential Equation Model for Chagas Disease with Interrupted Spraying Schedules. J Biol Dyn 6(2): 377–394doi:10.1080/17513758.2011.587896. [PubMed]
22. Spagnuolo AM, Shillor M, Stryker GA (2011) A Model for Chagas Disease with Controlled Spraying. J Biol Dyn 5(4): 299–317.
23. Reithinger R, Ceballos LA, Stariolo R, Davies CR, Grtler RE (2005) Chagas disease control: deltamethrin-treated collars reduce Triatoma infestans feeding success on dogs. Trans R Soc Trop Med Hyg 99: 502–508. [PubMed]
24. Cohen JE, Gürtler RE (2001) Modeling household transmission of American trypanosomiasis. Science 293: 694–698. [PubMed]
25. Barbu C, Dumonteil E, Gourbieŕe S (2009) Optimization of Control Strategies for Non-Domiciliated Triatoma dimidiata Chagas Disease Vector in the Yucata n Peninsula, Mexico. PLoS Negl Trop Dis 3(4): e416.doi:10.1371/journal.pntd.0000416. [PMC free article] [PubMed]
26. Gourbière S, Dumonteil E, Rabinovich JE, Minkoue R, Menu F (2008) Demographic and Dispersal Constraints for Domestic Infestation by Non-Domicilated Chagas Disease Vectors in the Yucatan Peninsula, Mexico. Am J Trop Med Hyg 78(1): 133–139. [PubMed]
27. Cruz-Pacheco G, Esteva L, Vargas C (2012) Control measures for Chagas disease. Math Biosci 237: 49–60. [PubMed]
28. Rabinovich JE, Himshoot P (1990) A population-dynamics simulation model of the main vectors of Chagas' Disease transmission, Rhodnius prolixus and Triatoma infestans. Ecol Modell 52: 249–266.
29. Rabinovich J, Schweigmann N, Yohai V, Wisnivesky-Colli C (2001) Probability of Trypanosoma cruzi transmission by Triatoma infestans (Hemiptera: Reduviidae) to the opossum Didelphis albiventris (Marsupialia: Didelphidae), Am J Trop Med Hyg. 65(2): 125–130. [PubMed]
30. Rabinovich JE, Wisnevesky-Colli C, Solarz ND, Gürtler RE (1990) Probability of Chagas disease by Triatomine infestans (Hemiptera: Reduviidae) in an endemic area of Santiago del Estero, Arentina. Bull World Health Organ 68(6): 737–746. [PubMed]
31. GNU Compiler Collection website. Available: http://gcc.gnu.org/fortran Accessed 2013 May 28.
32. Wolfram website. Available: http://www.wolfram.com Accessed 2013 May 28.
33. Gürtler RE, Ceballos LA, Ordóñez-Krasnowski P, Lanati LA, Stariolo R, et al. (2009) Strong Host-Feeding Preferences of the Vector Triatoma infestans Modified by Vector Density: Implications for the Epidemiology of Chagas Disease. PLoS Negl Trop Dis 3(5): e447.doi:10.1371/journal.pntd.0000447. [PMC free article] [PubMed]
34. Gürtler RE, Cecere MC, Lauricella MA, Cardinal MV, Kitron U, et al. (2007) Domestic dogs and cats as sources of Trypanosoma cruzi infection in rural northwestern Argentina. Parasitology 134: 69–82. [PMC free article] [PubMed]
35. Hall CA, Pierce EM, Wimsatt AN, Hobby-Dolbeer T, Meers JB (2010) Virulence and vertical transmission of two genotypically and geographically diverse isolates of Trypanosoma cruzi in mice. J Parasitol 96(2): 371–376. [PubMed]
36. Cardinal MV, Lauricella MA, Marcet PL, Orozco MM, Kitron U, et al. (2007) Impact of community-based vector control on house infestation and Trypanosoma cruzi infection in Triatoma infestans, dogs and cats in the Argentine Chaco. Acta Trop 103(3): 201–211. [PMC free article] [PubMed]
37. Gorla DE, Schofield CJ (1985) Analysis of egg mortality in experimental populations of Triatoma infestans under natural climatic conditions in Argentina. Bull Soc Vector Ecol 10: 107–117.
38. Castanera MB, Aparicio JP, Gürtler RE (2003) A stage-structured stochastic model of the population dynamics of Triatoma infestans the main vector of Chagas disease. Ecol Modell 162: 33–53.
39. Catalá S (1991) The biting rate of Triatoma infestans in Argentina. Med Vet Entomol 5(3): 325–333. [PubMed]
40. Kribs-Zaleta CM (2010) Estimating Contact Process Saturation in Sylvatic Transmission of Trypanosoma cruzi in the United States. PLoS Negl Trop Dis 4: e656. [PMC free article] [PubMed]
41. The World Factbook website. Available: https://www.cia.gov/library/publications/the-worldfactbook/geos/ar.html Accessed 29 May 2013.
42. Rassi Jr A, Rassi A, Marin-Neto JA (2009) Chagas heart disease: pathophysiologic mechanisms, prognostic factors and risk stratification. Mem Inst Oswaldo Cruz 104 (Suppl 1)152–158. [PubMed]
Articles from PLoS ONE are provided here courtesy of
Public Library of Science