- Open Access
A reduction in ecological niche for Trypanosoma cruzi-infected triatomine bugs
Parasites & Vectorsvolume 12, Article number: 240 (2019)
Theory predicts that parasites can affect and thus drive their hosts’ niche. Testing this prediction is key, especially for vector-borne diseases including Chagas disease. Here, we examined the niche use of seven triatomine species that occur in Mexico, based on whether they are infected or not with Trypanosoma cruzi, the vectors and causative parasites of Chagas disease, respectively. Presence data for seven species of triatomines (Triatoma barberi, T. dimidiata, T. longipennis, T. mazzottii, T. pallidipennis, T. phyllosoma and T. picturata) were used and divided into populations infected and not infected by T. cruzi. Species distribution models were generated with Maxent 3.3.3k. Using distribution models, niche analysis tests of amplitude and distance to centroids were carried out for infected vs non-infected populations within species.
Infected populations of bugs of six out of the seven triatomine species showed a reduced ecological space compared to non-infected populations. In all but one case (T. pallidipennis), the niche used by infected populations was close to the niche centroid of its insect host.
Trypanosoma cruzi may have selected for a restricted niche amplitude in triatomines, although we are unaware of the underlying reasons. Possibly the fact that T. cruzi infection bears a fitness cost for triatomines is what narrows the niche breadth of the insects. Our results imply that Chagas control programmes should consider whether bugs are infected in models of triatomine distribution.
Pathogens and hosts are co-evolving with each other. This interaction implies that the life history traits of either actor is evolutionarily and ecologically driven by the life history traits of the opposite actor . At the phylogenetic level, for example, this has been widely evidenced by the close evolutionary radiations between pathogens and hosts in many different taxa [2,3,4].
Parasites can modulate their host’s niche in two ways . First, parasites can alter the behavior and/or morphology of the host, thereby driving its niche specialization . For example, for cases where parasites manipulate their hosts, infected and non-infected host individuals will show trophic divergence as they are likely to differ in spatial and/or resource use . An illustrative study case is that of acanthocephalans that parasitize freshwater amphipods and whose effects induce behavioral changes by amphipods to facilitate parasite transmission to the definitive bird host [7, 8]. In this example, infected amphipod hosts make use of different food resources compared to non-infected amphipods leading to different infection-dependent trophic roles in freshwater ecosystems . Secondly, by modulating a host’s micro-ecological changes, parasites can lead to a host’s niche restriction so that the latter cannot maximize their niche breadth . A related example to this restriction applies to explain distribution of parasites with complex life histories (for example, a parasite with an intermediate and definitive host), and whose distribution is expected to be narrower than that of their hosts . In this case, the requirement for multiple hosts renders parasites more susceptible to local extinction (e.g. if hosts in some areas are more resistant) compared to parasites that require only one host to complete their life-cycle. However, a niche restriction for parasite–host pairs [5, 10] has not been tested (a close related example is that of ). According to a niche restriction test, one would predict that infected hosts will have a more restricted ecological niche than non-infected hosts.
Members of the subfamily Triatominae (Hemiptera: Reduviidae) are vectors of Trypanosoma cruzi, the causal agent of Chagas disease. Chagas disease is endemic and highly important in Latin American countries as it causes around 12,000 deaths annually  and is difficult to control. Trypanosoma cruzi is a typical parasite with a complex life-cycle; it requires bugs to carry and transfer it to humans and other mammalian reservoirs where it is transmitted through blood-feeding and defecation on the reservoir’s skin . Triatomines make use of different blood sources that include not only mammals, but birds and even reptiles . However, T. cruzi can only replicate in mammals so these animals play a key role in the survival of the parasite [15, 16]. Different local and international programmes and resources have been devoted to controlling this disease, one of which is to look at the ecological niche of triatomines to predict ecological conditions where infections are more likely to occur [17,18,19,20]. However, estimates based on ecological niche constructions usually assume that all bug individuals are equally likely to bear the parasite. That is, studies of triatomine niche do not consider that infected and non-infected individuals may differ in niche characteristics and so researchers have carried out niche predictions using both types of individuals within the same group (e.g. [21, 22]). Despite this, several sources of information indicate this assumption is unlikely as bugs defend themselves against the parasite, which may result in individual variation in parasite infection: (i) even when high numbers of parasites can enter the bug, there is a massive mortality of the former  and immune responses by the latter , suggesting an active reaction of the bug against the parasite; (ii) there is variation in the number of bugs that become infected in natural conditions (e.g. 9–39% infection rates in Triatoma phyllosoma [25,26,27]; in T. dimidiata, differential T. cruzi prevalence (0–100%) and discrete typing unit (DTU) specificity for individual haplogroups has been reported  implying that some bugs are more susceptible to infection; and (iii) T. cruzi infection has a strong impact on life history traits [22, 29, 30] and leads to a reduction triatomine survival in a temperature-dependent manner [31,32,33,34], implying that such fitness costs explain resistance by bugs as indicated above. Given the above, modelling of triatomine ecological niche should consider whether bugs are infected or not as, according to the theoretical background, infected bugs may have a different niche use.
In this study, we tested if there are niche use differences between T. cruzi-infected and non-infected bug populations of seven triatomine species. In particular, we explored the environmental associations of infected vs non-infected bugs’ niche uses to evaluate whether niche centroids (i.e. the environmental space where populations are expected to be better adapted to their environments) were closer than expected given their niche amplitude. Additionally, we compared the niche size (i.e. range of the environmental conditions) of infected vs non-infected (per vector species). The analysis of the environmental space allowed us to understand the environmental determinants of the geographical distribution, i.e. whether infected vectors occurred preferentially in a restricted (and common) region of the ecological space that could restrict their distribution in the geographical space. We discuss our results in terms of both our understanding of factors driving distribution and niche usage in this triatomine–trypanosomatid system, as well as in the prevention of Chagas disease.
Models of current potential distribution showed that the suitability areas for all triatomine species were variable (Fig. 1). Based on the environmental PCA analysis (Fig. 2), it could be determined that infected populations use a reduced ecological space when compared with the ecological niche generated from non-infected populations (see Student’s t-test results in Fig. 3), except for the case of T. pallidipennis in which the reverse was observed (Fig. 3). By examining the average distance between the centroids of infected populations (21 combinations; mean distance = 1.6) compared to the average distance between infected and non-infected pairs (mean distance = 1.3), we could observe that no particular environmental envelope for infected populations exists (i.e. there is not a reduced parasite niche); furthermore, the average inter-centroid distance between all combinations of infected populations is not significantly different that the inter-centroid distance of all triatomines (one-tailed Student’s t-test, t = 1.12, P = 0.12; Table 1). In most cases (5 out of 7), the distances between niche centroids for infected and non-infected pairs were lower than the average distance of the non-infected niche centroids to their corresponding random points (Fig. 3). The two species whose inter-centroid distances were higher were T. barberi and T. mazzottii, the latter exceeding a standard deviation from the mean of centroid-random point distances.
As predicted, there were intra-specific differences in niche characteristics in T. cruzi-non-infected triatomines compared to T. cruzi-infected triatomines, with the latter showing a more restricted niche breadth. This relationship could not be explained by a restricted niche for T. cruzi as a whole, suggesting that parasite occurrence is more likely where triatomines are closer to their niche centroid. This difference may be a consequence of the cost exerted by parasites on triatomines, which can limit triatomines’ niche usage. As indicated before, a fitness cost for the host has been already shown in triatomines in terms of reduced survival in the face of infection . This would select for defense mechanisms by the host which would explain, for example, the observed differences in triatomine’s immune ability against the parasite along altitude  and a reduction in T. cruzi population soon after infection . Furthermore, since not all bug species show the same trend in niche use for infected and non-infected animals, we hypothesize that this pattern may be seen as the different aspects of a co-evolutionary arms race between T. cruzi and triatomines. In this regard, either some populations or even species may have a reduced fitness costs (e.g. the case of species whose infected and non-infected populations have similar niche “traits”). Towards testing these ideas, it would be interesting to compare fitness differences derived from T. cruzi infection in species where infected and non-infected populations do not differ in niche (e.g. T. longipennis) vs those population that do differ (e.g. T. picturata). We predict that fitness costs will be higher in species where niches of infected and non-infected populations differ in comparison to species where such niches for the same groups are the same.
A second important issue to explain our results is that there are factors other than fitness costs driving a niche use difference, such as, for example, preference for a specific food source [35,36,37,38]. Given the relatively large area extensions occupied by triatomines (both infected and non-infected populations, especially for T. dimidiata, T. longipennis and T. phyllosoma), other abiotic factors may also play a role. For example, Langford & Janovi  found that soil conditions [loam soils (comprising mainly sand, silt and clay), temperature (5–35 °C) and soil moisture (< 50%)] were key to determining infection by lungworms in frogs that occupy large distribution ranges. Although T. cruzi is not a free-living parasite like the infective stages of lungworms, it certainly depends on several species of mammals to complete its life-cycle . Such plasticity in host use may imply that the parasite is unable to cope with all possible niche variations in terms of finding a host or a vector and reproducing within it. In general, host specificity determines how successful a parasite can be in invading or widening its habitat . In fact, Chagas disease risk in humans is not associated with mammal host diversity  which provides support to the suggestion that T. cruzi is not fully able to reach as many hosts and vector species as possible. A third topic to understand niche differences and similarities in infected and non-infected populations is related to temperature . In general, triatomines are heavily affected by temperature . For example, hatching rates in infected Rhodnius prolixus triatomines are temperature-dependent . A fourth topic that may determine parasite–triatomine distribution is that of anthropogenic factors. In this regard, studies of predictors of domestic areas have rarely included both actors simultaneously. For example, we know that predictors of triatomine occurrence include exterior (e.g. house age, upkeep, spatial location in the town ); public street lights ) and interior housing characteristics (number of inhabitants (over seven), overhead storage space, grain shed, cats, pigs and dogs ), altitude and mean annual precipitation . As for the parasite, a likely predictor includes unplastered mud-stick houses . When both actors have been found, the predictors for both were the accumulation of woodpiles , the periphery of rural villages  and presence of domestic mammals . In summary, food preferences, temperature and anthropogenic factors are three candidate factors that may determine triatomine distribution and niche. In this regard, laboratory experiments should be combined with distribution to assess the niche properties to then locate which conditions better predict the geographical location of infected triatomines.
It should be noted that we did not control the type of T. cruzi discrete typing unit (DTU) in our analysis as we considered a broad category of infected animals. In Mexico, it is predicted that all six described DTUs (TCI-TCVI) exist . The predicted distribution ranges overlap considerably for all six types within the country . This implies that is likely that we have more than one DTU in our database. In fact, mixed infections by different T. cruzi genotypes have been found in T. dimidiata populations from the Yucatan Peninsula . Unfortunately, the way infection was determined in our database (via direct microscopic observation), makes parasite genotyping impossible. Finally, a key conclusion is that our results should make people refine Chagas prevention programmes. Currently, for example, many of these programmes assume that all triatomines found in the wild are equally likely to harbor T. cruzi parasites which our results refute. These studies predict triatomine or mammal distributions with the aim of elucidating risk indices [16, 22, 54]. Thus, future studies should consider whether bugs bear T. cruzi parasites and, if possible, T. cruzi lineage to outline Chagas risks. In this regard, the present study can be used to guide the design strategies of triatomine control based on the geographical regions we depict according to our present models, and with this, optimize the allocation of resources.
The present study provides evidence that Trypanosoma cruzi infection correlates with niche characteristics used by six triatomine bug species that occur in Mexico. Two implications of this are: (i) T. cruzi may restrict the triatomine ecological niche so that triatomines cannot maximize their ability to colonize their full ecological spectrum; and (ii) predictions of triatomine distribution and triatomine control strategies should take into account that T. cruzi-infected populations may be geographically located in more particular environments.
We used data of presence of seven Triatominae species occurring within the Mexican territory: T. barberi (Guanajuato, Guerrero, Hidalgo, Jalisco, Mexico state, Michoacan, Morelos and Oaxaca states; Additional file 1: Table S1), T. dimidiata (Campeche, Chiapas, Colima, Guerrero, Hidalgo, Jalisco, Mexico state, Michoacan, Oaxaca, Puebla, Queretaro, Quintana Roo, Tabasco, Tamaulipas, Veracruz and Yucatan states; Additional file 1: Table S2), T. longipennis (Aguascalientes, Colima, Chihuahua, Guanajuato, Jalisco, Michoacan, Nayarit, Sinaloa and Zacatecas states; Additional file 1: Table S3), T. mazzottii (Guerrero and Oaxaca states; Additional file 1: Table S4), T. pallidipennis (Guanajuato, Guerrero, Jalisco, Mexico state, Michoacan, Morelos, Nayarit, Oaxaca and Puebla states; Additional file 1: Table S5), T. phyllosoma (Oaxaca and Guerrero states; Additional file 1: Table S6) and Triatoma picturata (Nayarit and Jalisco states; Additional file 1: Table S7). Data were gathered from diverse sources: (i) Entomology Laboratory, Institute of Epidemiological Diagnosis and Reference (InDRE/ Health Secretary; Mexico City); (ii) technicians and/or entomologists from the health services of each Mexican state and locality where there is a problem of domiciliary infestation. Note that the effort and periodicity of collection depended on each particular state and sanitary jurisdiction; (iii) Entomology Laboratory, National School of Biological Sciences, Instituto Politécnico Nacional (Mexico City) whose members collected specimens from intradomiciliary areas; and (iv) literature references that include collections from urbanized and non-urbanized areas (see Additional file 1: Table S8). Records used were those falling between 1999 and 2013 for all species as collection was more intense in this time window. When collected, samples were checked for the presence of T. cruzi via direct microscopic observation of fecal content at 40× magnification. Feces were obtained by abdominal pressure and diluted in phosphate buffered saline prior to microscope observation . Although this method implies a detection error, it has been safely validated so the error is minor . This validation is based on the fact that the high motility of the parasite makes it very conspicuous during microscopic surveillance .
For bioclimatic variables, we used WorldClim v.1.4 (www.worldclim.org) data  at a cell size of 0.041666669. To establish a set of uncorrelated climatic variables, we intersected the variables with presence points, but removed some variables [BIO1, annual mean temperature; BIO2, mean diurnal range [mean of monthly (maximum temp − minimum temp)]; BIO3, isothermality [(BIO2/BIO7) × 100)]; BIO5, maximum temperature of warmest month; BIO6, minimum temperature of coldest month; BIO8, mean temperature of wettest quarter; BIO9, mean temperature of driest quarter; BIO10, mean temperature of warmest quarter; BIO12, annual precipitation; BIO13, precipitation of wettest month; BIO14, precipitation of driest month; BIO17, precipitation of driest quarter; BIO18, precipitation of warmest quarter; and BIO19, precipitation of coldest quarter] with an exploratory data analysis and a Pearson correlation analysis (i.e. any value > 0.7), a jack-knife contribution and biological knowledge of species. After this, the final data set included temperature seasonality (BIO4), temperature annual range (BIO7), mean temperature of warmest quarter (BIO11), precipitation seasonality (BIO15) and precipitation of coldest quarter (BIO16).
Potential species distribution modeling
Species distribution models (SDM) were generated with Maxent 3.3.3k . Among software that makes use of presence-only data, Maxent has proven to be reliable under different datum conditions and applications . In fact, despite having faced a large theoretical and empirical discussion, it has become one of the most widespread programs for species distribution models . Briefly, Maxent uses georeferenced data points where the species has been found together with environmental variables to estimate the distribution (geographical range) of a species by finding the distribution which has maximum entropy (i.e. is closest to geographically uniform) subject to constraints derived from their predicting variables . A crucial component of model calibration is the delimitation of an area accessible for species dispersal , over a relevant time period, the M (mobility), in the BAM model . As such, we assumed that the biogeographical provinces of Mexico  occupied by each triatomine species are the boundaries of the areas accessible to colonize. Therefore, all biogeographical regions containing at least one presence record of the presence database were used in combination to mask the environmental predictors for each species. To corroborate that the calibration area was appropriated (i.e. covering both suitable and unsuitable environmental combinations ), we visually inspected the distribution of species data points within the background area (i.e. the “M” region) in a biplot comprised by annual mean temperature and annual precipitation. Final models were constructed on a species-specific basis after testing predictive ability under different parameter settings that modify model complexity. For this task, we used the same data partition (random seed of 70% for training and 30% for testing), the regularization number (i.e. an inverse penalization parameter) by using 0.5, 1, 2, 3 and 4 and the model features (i.e. quadratic, linear, etc.), in order to have all possible combinations. Other parameters were kept to default (convergence = 10–5, maximum number of iterations = 500, background = 10,000). Duplicate records were removed. To generate an estimation of predictive accuracy, we ran bootstrap replicated type using 10 replicates. Finally, extrapolation and clamping were not allowed. All of this was done to find which combination of settings and variables generated the best outcomes, while minimizing the number of model parameters, as well as producing ‘closed’, bell-shaped response curves which guaranteed model calibration. The models were evaluated based on the highest area under the curve value (AUC), and lower AUC difference between training and testing outputs .
The values of each of the bioclimatic variables used in niche models were extracted to every pixel in the Mexican territory to produce a matrix where a principal component analysis was applied on the correlation matrix. Next, the eigenvalues for the two first principal components were aggregated, as they accounted for the 73.378% of the environmental variance. Using the potential distribution models of triatomine species (the final presence map) and the correspondent infected data points, a final matrix was produced that included eigenvalues for every pixel of the presence maps (distribution model of species i and infected data points of species i). These matrices were converted to shapefiles in ARCMap using the principal components environment as the space. Then, a minimum convex polygon was created to generate: (i) its geometric centroid and (ii) 500 random points; with that, a multi-distance matrix between the centroid and each of their correspondent 500 random points was created for each vector and infected populations. We observed whether the inter-centroid distance differed significantly (∞ = 0.05) of the average centroid-random points distances. Additionally, we tested whether these two groups (infected and non-infected populations), had the same range and position using Student’s t-test.
species distribution models
area under the curve
discrete typing units
Schmid-Hempel P. Evolutionary parasitology: the integrated study of infections, immunology, ecology, and genetics. Oxford: Oxford University Press; 2011.
Brooks DR, McLennan DA. Phylogeny, ecology, and behavior. Chicago: University of Chicago Press; 1991.
Clayton DH, Bush SE, Johnson KP. Ecology of congruence: past meets present. Syst Biol. 2004;53:165–73.
Martinsen ES, Perkins SL, Schall JJ. A three-genome phylogeny of malaria parasites (Plasmodium and closely related genera): evolution of life-history traits and host switches. Mol Phylogenet Evol. 2008;47:261–73.
Poulin R. Evolutionary ecology of parasites. Princeton: Princeton University Press; 2011.
Britton JR, Andreou D. Parasitism as a driver of trophic niche specialisation. Trends Parasitol. 2016;32:437–45.
Bauer A, Haine E, Perrot-Minnot M, Rigaud T. The acanthocephalan parasite Polymorphus minutus alters the geotactic and clinging behaviours of two sympatric amphipod hosts: the native Gammarus pulex and the invasive Gammarus roeseli. J Zool. 2005;267:39–43.
Médoc V, Bollache L, Beisel JN. Host manipulation of a freshwater crustacean (Gammarus roeseli) by an acanthocephalan parasite (Polymorphus minutus) in a biological invasion context. Int J Parasitol. 2006;36:1351–8.
Médoc V, Piscart C, Maazouzi C, Simon L, Beisel JN. Parasite-induced changes in the diet of a freshwater amphipod: field and laboratory evidence. Parasitology. 2011;138:537–46.
Rohde K. Niche restriction in parasites: proximate and ultimate causes. Parasitology. 1994;109:S69–84.
Morand S, Bordes F, Pisanu B, de Bellocq JG, Krasnov B. The geography of defence. In: Morand S, Krasnov BR, editors. The biogeography of host–parasite interactions. Oxford: Oxford University Press; 2010. p. 159–72.
Maher SP, Ellis C, Gage KL, Enscore RE, Peterson AT. Range-wide determinants of plague distribution in North America. Am J Trop Med Hyg. 2010;83:736–42.
Bonney KM. Chagas disease in the 21st century: a public health success or an emerging threat? Parasite. 2014;21:11.
Tyler KM, Olson CL, Engman DM. The life cycle of Trypanosoma cruzi. In: Tyler KM, Miles MA, editors. American trypanosomiasis World Class Parasites. Boston: Kluwer Academic Publishers; 2003.
Noireau F, Diosque P, Jansen AM. Trypanosoma cruzi: adaptation to its vectors and its hosts. Vet Res. 2009;40:26.
Jansen AM, Xavier SC, Roque AL. Ecological aspects of Trypanosoma cruzi: wild hosts and reservoirs. In: Telleria J, Tibayrenc A, editors. American Trypanosomiasis Chagas disease. London: Elsevier; 2017. p. 243–64.
Guhl F. Geographical distribution of Chagas disease. In: Tellería J, Tibnayrenc M, editors. American trypanosomiasis Chagas disease (second edition): one hundred years of research. London: Elsevier; 2017. p. 89–112.
Rengifo-Correa L, Stephens C, Morrone JJ, Téllez-Rendón JL, González-Salazar C. Understanding transmissibility patterns of Chagas disease through complex vector-host networks. Parasitology. 2017;144:760–72.
Ramsey JM, Peterson AT, Carmona-Castro O, Moo-Llanes DA, Nakazawa Y, Butrick M, et al. Atlas of Mexican Triatominae (Reduviidae: Hemiptera) and vector transmission of Chagas disease. Mem Inst Oswaldo Cruz. 2015;110:339–52.
Gurgel-Gonçalves R, Galvao C, Costa J, Peterson AT. Geographic distribution of Chagas disease vectors in Brazil based on ecological niche modeling. J Trop Med. 2012;2012:705326.
Ibarra-Cerdeña CN, Valiente-Banuet L, Sánchez-Cordero V, Stephens CR, Ramsey JW. Trypanosoma cruzi reservoir-triatomine vector co-occurrence networks reveal meta-community effects by synanthropic mammals on geographic dispersal. PeerJ. 2017;5:e3152.
Peterson AT, Sánchez-Cordero V, Beard CB, Ramsey JW. Ecological niche modeling and potential reservoirs for Chagas disease, Mexico. Emerg Infect Dis. 2002;8:662–7.
Ferreira RC, Kessler RL, Lorenzo MG, Paim RM, Ferreira Lde L, Probst CM, et al. Colonization of Rhodnius prolixus gut by Trypanosoma cruzi involves an extensive parasite killing. Parasitology. 2016;143:434–43.
De Fuentes-Vicente JA, Cabrera-Bravo M, Enríquez-Vara JN, Flores-Villegas AL, Vidal-López DG, Gutiérrez-Cabrera AE, et al. Relationships between altitude, triatomine (Triatoma dimidiata) immune response and virulence of Trypanosoma cruzi, the causal agent of Chagasʼ disease. Med Vet Entomol. 2016;31:63–71.
Ramsey JM, Ordoñez R, Cruz-Celis A, Alvear AL, Chávez V, López R, et al. Distribution of domestic Triatominae and stratification of Chagas disease transmission in Oaxaca, Mexico. Med Vet Entomol. 2000;14:19–30.
Vidal-Acosta V, Ibáñez S, Bernal R, Martínez-Campos C. Infección natural de chinches Triatominae con Trypanosoma cruzi asociada a la vivienda humana en México. Salud Publ Mex. 2000;42:496–503.
Villalobos G, Martínez-Hernández F, de la Torre P, Laclette JP, Espinoza B. Entomological indices, feeding sources, and molecular identification of Triatoma phyllosoma (Hemiptera: Reduviidae) one of the main vectors of Chagas disease in the Istmo de Tehuantepec, Oaxaca, Mexico. Am J Trop Med Hyg. 2011;85:449–97.
Pech-May A, Mazariegos-Hidalgo CJ, Izeta-Alberdi A, López-Cancino SA, Tun-Ku E, De la Cruz-Félix K, et al. Genetic variation and phylogeography of the Triatoma dimidiata complex evidence a potencial center of origin and recent divergence of haplogroups having differential Trypanosoma cruzi and DTU infections. PLoS Negl Trop Dis. 2019;13:e0007044.
Alejandre-Aguilar R, Nogueda-Torres B, Calvo-Méndez ML, Cortés-Jiménez M. Estudio comparativo de la susceptibilidad de cinco especies de triatominos (Insecta: Reduviidae) a la infección con Trypanosoma cruzi. Rev Latinoam Microbiol. 1993;35:201–6.
Botto-Mahan C. Trypanosoma cruzi induces life-history trait changes in the wild kissing bug Mepraia spinolai: implications for parasite transmission. Vector Borne Zoonotic Dis. 2009;9:505–10.
Fellet MR, Lorenzo MG, Elliot SL, Carrasco D, Guarneri AA. Effects of infection by Trypanosoma cruzi and Trypanosoma rangeli on the reproductive performance of the vector Rhodnius prolixus. PLoS ONE. 2014;9:e105255.
Elliot SL, Rodrigues JO, Lorenzo MG, Martins-Filho OA, Guarneri AA. Trypanosoma cruzi, etiological agent of Chagas disease, is virulent to its triatomine vector Rhodnius prolixus in a temperature-dependent manner. PLoS Negl Trop Dis. 2015;9:e0003646.
Luz C, Fargues J, Grunewald J. The effect of fluctuating temperature and humidity on the longevity of starved Rhodnius prolixus (Hem., Triatominae). J Appl Ent. 1998;122:219–22.
Garcia AR, de Paula-Rocha A, Moreira CC, Rocha SL, Guarneri AA, Elliot SL. Screening of fungi for biological control of a triatomine vector of chagas disease: temperature and trypanosome infection as factors. PLoS Negl Trop Dis. 2016;10:e0005128.
Gurtler RE, Ceballos LA, Ordóñez-Krasnowski P, Lanati L, Kitron U. Strong host-feeding preferences of the vector Triatoma infestans modified by vector density: implications for the epidemiology of Chagas disease. PLoS Negl Trop Dis. 2009;3:e447.
Rabinovich JE, Kitron UD, Obed Y, Yoshioka M, Gottdenker N, Chaves LF. Ecological patterns of blood-feeding by kissing-bugs (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2011;106:479–94.
Otálora-Luna F, Pérez-Sánchez AJ, Sandoval C, Aldana E. Evolution of hematophagous habit in Triatominae (Heteroptera: Reduviidae). Rev Chil His Nat. 2015;88:4.
De Fuentes-Vicente JA, Gutiérrez-Cabrera AE, Flores-Villegas AL, Lowenberger C, Benelli G, Salazar-Schettino PM, et al. What makes an effective Chagas disease vector? Factors underlying Trypanosoma cruzi-triatomine interactions. Acta Trop. 2018;183:23–31.
Langford GJ, Janovy J Jr. Ecological factors responsible for the geographic distribution of Rhabdias joaquinensis: where do lungworms infect anurans in nature? Parasitol Res. 2016;115:1305–13.
WHO. Chagas disease (American trypanosomiasis). Geneva: World Health Organization; 2017. http://www.who.int/chagas/en/. Accessed 6 Apr 2017.
Poulin R, Mouillot D. Parasite specialization from a phylogenetic perspective: a new index of host specificity. Parasitology. 2003;126:473–80.
Oda E, Solari A, Botto-Mahan C. Effects of mammal host diversity and density on the infection level of Trypanosoma cruzi in sylvatic kissing bugs. Med Vet Entomol. 2014;28:384–90.
Carcavallo RU. Climatic factors related to Chagas disease transmission. Mem Inst Oswaldo Cruz. 1999;94:367–9.
Canals M, Solis R, Valderas J, Ehrenfeld M, Cattan PE. Preliminary studies on temperature selection and activity cycles of Triatoma infestans and T. spinolai (Heteroptera: Reduviidae), Chilean vectors of Chagasʼ disease. J Med Entomol. 1997;34:11–7.
Cohen JM, Wilson ML, Cruz-Celis A, Ordoñez R, Ramsey JM. Infestation by Triatoma pallidipennis (Hemiptera: Reduviidae: Triatominae) is associated with housing characteristics in rural Mexico. J Med Ent. 2006;43:1252–60.
Pacheco-Tucuch FS, Ramirez-Sierra MJ, Gourbière S, Dumonteil E. Public street lights increase house infestation by the Chagas disease vector Triatoma dimidiata. PLoS ONE. 2012;7:e36207.
Campbell-Lendrum DH, Angulo VM, Esteban L, Tarazona Z, Parra GJ, Restrepo M, et al. House-level risk factors for triatomine infestation in Colombia. Int J Epidemiol. 2007;36:866–72.
Mott KE, Muniz TM, Lehman JS Jr, Hoff R, Morrow RH Jr, de Oliveira T, et al. House construction, triatomine distribution, and household distribution of seroreactivity to Trypanosoma cruzi in a rural community in northeast Brazil. Am J Trop Med Hyg. 1978;27:1116–22.
Coutinho CF, Souza-Santos R, Lima MM. Combining geospatial analysis and exploratory study of triatomine ecology to evaluate the risk of Chagas disease in a rural locality. Acta Trop. 2012;121:30–3.
Ramirez-Sierra MJ, Herrera-Aguilar M, Gourbière S, Dumonteil E. Patterns of house infestation dynamics by non-domiciliated Triatoma dimidiata reveal a spatial gradient of infestation in rural villages and potential insect manipulation by Trypanosoma cruzi. Trop Med Int Health. 2010;15:77–86.
Ramsey JM, Gutiérrez-Cabrera AE, Salgado-Ramírez L, Peterson AT, Sánchez-Cordero V, Ibarra-Cerdeña CN. Ecological connectivity of Trypanosoma cruzi reservoirs and Triatoma pallidipennis hosts in an anthropogenic landscape with endemic Chagas disease. PLoS ONE. 2012;7:e46013.
Izeta-Alberdi A, Ibarra-Cerdeña CN, Moo-Llanes DA, Ramsey JW. Geographical, landscape and host associations of Trypanosoma cruzi DTUs and lineages. Parasit Vectors. 2016;9:631.
López-Cancino SA, Tun-Ku E, De la Cruz-Felix HK, Ibarra-Cerdeña CN, Izeta-Alberdi A, Pech-May A, et al. Landscape ecology of Trypanosoma cruzi in the southern Yucatan Peninsula. Acta Trop. 2015;151:58–72.
Guhl F, Aguilera G, Pinto N, Vergara D. Updated geographical distribution and ecoepidemiology of the triatomine fauna (Reduviidae: Triatominae) in Colombia. Biomedica. 2007;27:143–62.
Schenone H. Xenodiagnosis. Mem Inst Oswaldo Cruz. 1999;94:289–94.
Velasco CO, Guzmán BC. Manual de técnicas de laboratorio. Instituto Nacional de Salud Pública (Mexico City), Mexico: INDRE; 1994.
Villagrán ME, Marín C, Hurtado A, Sánchez-Moreno M, de Diego J. Natural infection and distribution of triatomines (Hemiptera: Reduviidae) in the state of Queretaro, Mexico. Trans R Soc Trop Med. 2008;102:833–8.
Hijmans R, Cameron S, Parra J, Jones P, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25:1965–78.
Phillips S, Anderson R, Schapire R. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.
Elith J, Michael K, Steven P. The art of modelling range shifting species. Methods Ecol Evol. 2010;1:330–42.
Phillips S, Anderson R, Dudík M, Schapire R, Blair M. Opening the black box: an open-source release of Maxent. Ecography. 2017;40:887–93.
Barve N, Barve V, Jiménez-Valverde A, Lira-Noriega A, Maher P, Peterson A, et al. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol Model. 2011;222:1810–9.
Soberón J, Peterson A. Interpretation of models of fundamental ecological niches and species distributional areas. Biodiver Inform. 2005. https://doi.org/10.17161/bi.v2i0.4.
Portal de Geoinformación Sistema Nacional de Información Sobre Biodiversidad. CONABIO. http://www.conabio.gob.mx/informacion/gis/. Accessed 10 Mar 2018.
Owens H, Campbell L, Dornak L, Saupe E, Barve N, Soberon J, et al. Constraints on interpretation of ecological niches models by limited environment ranges on calibration areas. Ecol Model. 2013;263:10–8.
This work was supported by a UNAM PAPIIT grant (IN206618) to ACA. GV and MRO were supported by postdoctoral fellowships from Consejo Nacional de Ciencia y Tecnología. Funding bodies had no role in the design of the study and collection, analysis or interpretation of data, or in writing the manuscript.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional file. The datasets used and/or analyzed during the present study are available from the corresponding author upon request.
CIC, GV, ANB, MRO, AEGC and ACA conceived the study. GV, ANB, JADFV, JLTR, HH, FMH, AEGC and ACA collected data. CIC, GV, ANB, JADFV, MRO, AEGC and ACA analyzed and interpreted data. CIC, GV, ANB, JADFV, MRO, AEGC and ACA contributed to writing for intellectual content. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.