An epidemiological network contains all the organisms involved (types) in the transmission of a parasite. The nodes of the network represent reservoirs, hosts, and vectors, while the links between the nodes represent the strength and direction of parasite movement. Networks that contain humans are of special interest because they are of concern to public health authorities. Under these circumstances, it is possible, in principle, to identify cycles (closed paths in the network) that include humans and select the ones that carry the maximum probability of human infection. The basic reproduction number R0 in such a network gives the average number of new infections of any type after the introduction of one individual infected by any type. To obtain R0 for complex networks, one can use the next-generation matrix (NGM) approach. Every entry in NGM will average the contribution of each link that connects two types. To tease the contribution of every cycle apart, we define the virulence as the geometric mean of the NGM entries corresponding to the links therein. This approach allows for the quantification of specific cycles of interest while it also makes the computation of the sensitivity and elasticity of the parameters easier. In this work, we compute the virulence for the transmission dynamics of Chagas disease for a typical rural area in Colombia incorporating the effect of environmental changes on the vector population size. We concluded that the highest contribution to human infection comes from humans themselves, which is a surprising and interesting result. In addition, sensitivity analysis revealed that increasing vector population size increases the risk of human infection.
What changes can be expected in human infection risk produced by shifts in land use and climate? The answer is obviously complicated as some diseases may actually shrink their distributions1,2 while others can potentially become more common.3–45 We are particularly interested in investigating infection risk of vector-borne tropical diseases, and among them Chagas disease. Chagas disease in Colombia is transmitted to humans primarily by vectors of the species Rhodnius prolixus6,7 that feed from reservoirs in what has been called the sylvatic cycle. The main reservoir in the wild is Didelphis marsupialis, while dogs, pigs, and cats may act as domiciliary reservoirs.8,9 Humans became part of the transmission network by constructing dwellings close to palm trees, the preferred habitat of the vector. Disease transmission from vectors to humans and reservoirs requires contact with insect feces infected with the parasite Trypanosoma cruzi.10 Transmission to insects from humans and reservoirs occurs when insects seek blood meals.11 In Colombia, there are 2 people million infected and 8 million at risk.7
Current changes in land use and possible shifts in temperature, humidity, and rain seasons can potentially alter the transmission of T. cruzi.9 Indeed, increased temperature has been linked to shorter insect breeding cycles,12 while palm tree mono-crops can provide huge areas for insect development. But are these changes in insect population size going to translate into an increased risk of human infection? After all, changes in land use while providing more land for insects could alter the composition of species in the environment (including ecological relations13,14) resulting in a different parasite transmission network, which could be either less or more efficient in spreading the disease.5,15–161718
Unfortunately, measuring the effect of climate or land use change by field experiments might be too slow or extremely complicated.4 However, epidemiological models can be used to understand the role of different species in disease transmission and the potential effect of environmental changes. Mathematical epidemiological models aim at describing or predicting the average dynamic of the transmission of a disease among members of a particular population.19–2021 An important quantity associated with these models is the basic reproductive number R0: a threshold value that indicates whether a disease is going to invade the population.22,23 If multiple populations are involved, R0 can be computed as the spectral radius of the next generation matrix (NGM). The NGM contains the number of secondary cases of each pair of types involved. Although, given a specific model, obtaining R0 is straightforward (cf.ref. 23), it does not tell us the contribution of each of the populations involved in the transmission of the parasite.21,22 In addition, its sensitivity to model parameters is often complicated to derive,24,25 making it less an attractive tool to investigate the effect of environmental changes.
To overcome these difficulties, we recently proposed an alternative quantity that we called the virulence.26 To compute the virulence, we start by considering the environment as a network where nodes correspond to species and the edges represent relationships between species (ie, who interacts with who). The virulence is defined for a path in the network that starts and finishes at the same node (a cycle). It is clear that changes in the environment can alter the network structure, for example, it can change node densities or it could change the overall connectivity. The virulence of a cycle is the geometric mean of the entries of the NGM for the edges involved in the cycle. The critical virulence, which is the greatest cycle virulence, is telling us that the reproductive capacity of a disease is concentrated on the critical circuit. The critical virulence bounds the basic reproductive number, and being easier to compute, it is more malleable for sensitivity analysis.
In this work, we present a simple mathematical model for Chagas disease transmission and find expressions for the virulence and the critical virulence for the system. In addition, using parameters that are biologically feasible, we identify key cycles in the network. Finally, we investigate the possible impact of an environmental change, modeled as a change in vector population size, by computing the sensitivity and elasticity of the critical virulence to this variable.
First we consider a general epidemiological framework for vector-borne diseases to derive the governing equations, but in the Results section, we will derive the parameters specific for Chagas disease in Colombia. The epidemiological system can be described by four different types of carriers of a pathogen: two different kinds of hosts, which we will call R1 and R2, a vector V, and humans, H. An epidemiological model for such a system is given by the equations:
• nj, SJ, Ij are the abundance of the total population, the population of susceptible and the infected individuals, respectively. j = V for vectors, j = 1 for hosts of type R1 j = 2 for hosts of type R2, and j = H for humans. We also have the relation: nj = Sj + Ij
• mj are death rates and recruitment rates in each population; they are assumed to be equal so the population size remains constant. We consider an increased mortality for humans that have acquired the pathogen represented by dH, so the natural death rate and recruitment rate for humans is given by: mH + dH;
• βj and bj are the infection rates from host or humans to vectors and vice versa, respectively.
Using the assumption that the population size remains constant (so birth and death rates are comparable), we put SJ = nj - Ij and the system 1-8 collapses to:
The transmission network is represented by the graph in Figure 1. Arrows contain the number of new infections of type j that are caused by the introduction of an infected individual of every type j. The network is characterized by nodes (populations) and arrows that connect them. A cycle is a path in the network that starts and finishes at the same node. A simple cycle is a cycle that does not enter and leave more than once through any node. Note that we do not include transmission within populations, and vectors move parasites between populations of reservoirs and hosts.
The NGM contains the number of new infections produced in each type in the system by an infected individual of every other type. Thus, in each entry of the NGM, we have the term associated with the arrow in the directed graph. The NGM for the system described in Figure 1 is given by:
The basic reproductive number R0 is the average number of infections produced by the contribution of the whole network. If R0 is greater than 1, the disease will establish itself in the population.27 Therefore, R0 is a quantitative approximation for long-term infection risk of the network. R0 can be computed from the NGM by finding the largest eigenvalue of the NGM.28 The contribution of every cycle to R0 is an important aspect of disease transmission: it will tell us the relative weight of a path for this network's infection risk.
Based on previous work26 and using ref. 29, we established that the virulence of a cycle µ(NGM), or the contribution to R0 from a given cycle, can be computed by taking the geometric mean of every entry of the NGM involved in the cycle. Furthermore, there is a cycle that contributes the most (µ (NGM)) to the reproductive capacity of the disease, and we call its virulence the critical virulence.
The simple circuits in our network are:
with respective virulence equal to:
The value µmax(NGM) is obtained by finding among all these simple circuits in the graph the greatest one.
Chagas disease in Colombia could be grossly characterized by the transmission of T. cruzi to vectors of the species R. prolixus (V) and sylvatic reservoirs such as D. marsupialis (R1) and domestic animals, like dogs or pigs (R2) and humans (H). The amount of parasite transmission between the nodes in this network is quantified by the virulences defined in the previous section. In order to compare these expressions, we use data from the literature (actual number choices detailed in the Discussion section) and some simple manipulations to reduce the amount of unknowns. To this end, we start by dividing NGM by nv to get the relative abundance to vectors, instead of the absolute abundance. Using an estimate of 1 domestic animal per 1000 vectors, 1 human per 1000 vectors, and 1 sylvatic reservoir every 2000 vectors we get:
Mortality rates can be estimated from life expectancy. We use and estimate a 2-year life expectancy for R1, 1 year for V, 3 years for R2, and 70 years for H. Thus,
Disease-induced death rate in humans can be estimated by assuming a decrease in life expectancy of an average of 20 years. Therefore:
In the model, disease transmission rates have units of new infections per susceptible per unit of time. Because they are difficult to estimate, it is more convenient to express them as relative odds of infection. Thus, we further divide NGM by βH (infection rate of vectors after a human encounter) and determine the relative difficulty of transmission compared to this route. We estimate that infection of vectors after a contact with an R2 is equally probable that an infection after a human contact, but can be twice as probable from R1 since they interact continuously in the sylvatic cycle. Thus,
Infection of H, R1, and R2 from V tends to be less efficient since it requires contact with insect feces. We also assumed that bH < b1 < b2 based on the number of contacts and reported incidences. Thus,
Virulence and critical virulence
Using the expressions for virulence and the set of parameters derived in the previous section and varying βH between 10-5 and 1, we found that the virulence and R0 vary linearly between:
for nv = 1 (ie, per vector). The critical virulence (see ref. 26, 29) is µmax(NGM) = µ(H→V→H)
Elasticity and sensitivity
E(µmax(NGM), nv) and S(µmax(NGM), nv) for βH between 10-5 and 1, and nV = 1 are:
In this work, we presented a simple mathematical model that contains some of the features that characterize the transmission of T. cruzi between sylvatic, domestic, and human populations. We consider three simple cycles: (i) R1 → V → R1, (ii) R2 → V → R2, and (iii) H → V → H. Other cycles are missing, for example, humans can have transmission via blood transfusion,30 insect nymphs could eat contaminated feces from adults,30 or domestic animals could eat insects.31
In addition, there are multiple other sylvatic populations involved, including reptiles and birds that do not host the parasite but feed the insects.32 However, the cycles considered are the best known and might be responsible for the greatest portion of parasite transmission.30,32
With the simple structure shown in Figure 1, we had to estimate a total of 10 parameters and the densities for each population involved. Population densities were normalized by insect population, so we estimated relative abundances. Reports suggest that D. marsupialis can have densities of two to three individuals per km2.31 For the same area, we found reports of 2000 insects (100 insects per palm times 20 palms).33 Similarly, we estimated that in rural areas, the density of humans and domestic animals is around six per km2. The values for the virulence for every cycle and R0 are reported per insect.
Mortality rates came from life expectancy reported in the literature. However, the experimental reports vary even within species; thus, we use approximate life spans in years. We found that R. prolixus lives for about a year12,34 and D. marsupialis for about two years.8 For humans, we used life expectancy of 70 years and assume the disease would shorten the life span by 20 years on average.7,10
Transmission rates were expressed as proportions to the transmission rate from humans to vectors. There are some reports that suggest that β's (parasite transmission via blood meals) are higher than b's (parasite transmission via contaminated feces).30 Because the rates of transmission are difficult to measure and they are critical for model results, we vary βH between 10-5 (1 out of 100.000 encounters results in an infected vector) and 1 (1 out of 1 encounter results in an infected vector) and kept the proportions constant to explore the behavior of the critical virulence. We believe that this range contains the biologically feasible values for the transmission rate.
Virulence and critical virulence
We found with the model that the cycle H → V → H is the critical cycle because it has the maximum virulence. The virulence for this cycle is 12% higher than the virulence of R1 → V → R1 and 30% higher compared to R2 → V → R2, for all values of βH. We also computed R0 and found that extremely low values of βH produce R0 = 3.4 X 10-7. Thus, a population in the order of millions of insects per km2 is needed to maintain the disease. On the other hand, when βH = 1, R0 = 3.4 X 10-2, and only 100 insects per km2 would be enough to maintain the disease in the population. This last number is greatly exceeded in many endemic areas.
This finding is interesting because very often health authorities target insect eradication, which is very important, but this study suggests that additional efforts should be made to identify people infected and provide them with treatment, not only for ethical reasons, but also because having humans infected is what contributes the most to disease establishment, at least in this simple model.
Elasticity and sensitivity
Climate and many other local environmental conditions are likely to have a large impact on vector-borne diseases, as survival, development, and physiological rates of vectors and hosts are often related to abiotic variables.35 In addition, climatic factors will, to a large extent, determine where a vector species can persist, and the same applies for many host species, therefore climate change could expand or diminish the areas where the disease can establish.25
Sensitivity analysis provides a way to measure how small changes in the parameters translate into variations in the critical virulence.24,36 We found that vector densities increase the critical virulence linearly for 10-5 ± βH ± 1. Sensitivity analysis showed that a 50% change in critical virulence can be expected after a small perturbation of insect population's size. The elasticity corroborated that this relation is maintained for the whole range of βH. This suggests that insect population densities play a prominent role in human infection. If they change because of land use shifts or climate variations, we could expect different patterns of disease transmission.
We conclude that insect population densities, as expected, play an important role in human Chagas disease infection risk. Sensitivity and elasticity analysis revealed that a small change in insect densities could translate into 50% increase in the number of human secondary cases. If changes in land use or climate produce changes in insect population sizes, then human infection risk is expected to change in a steep manner. In addition, we found that the simple cycle H → V → H contributes the most to Chagas disease establishment. These findings, taken together, suggested to us that it is important to actively screen human populations for infection while continuing the efforts to keep insect densities low. This is important for Colombia because only a small proportion of the population infected is diagnosed when donating blood or by health systems and historically strong effort has been directed to house improvement for insect eradication.
Conceived and designed the experiments: JMC, CS. Contributed to the writing of the manuscript: JMC, CS. Agree with manuscript results and conclusions: JMC, CS. Jointly developed the structure and arguments for the paper: JMC, CS. Made critical revisions and approved final version: JMC, CS. Both authors reviewed and approved of the final manuscript.