Distributional potential of the Triatoma brasiliensis species complex at present and under scenarios of future climate conditions
Parasites & Vectors volume 7, Article number: 238 (2014)
The Triatoma brasiliensis complex is a monophyletic group, comprising three species, one of which includes two subspecific taxa, distributed across 12 Brazilian states, in the caatinga and cerrado biomes. Members of the complex are diverse in terms of epidemiological importance, morphology, biology, ecology, and genetics. Triatoma b. brasiliensis is the most disease-relevant member of the complex in terms of epidemiology, extensive distribution, broad feeding preferences, broad ecological distribution, and high rates of infection with Trypanosoma cruzi; consequently, it is considered the principal vector of Chagas disease in northeastern Brazil.
We used ecological niche models to estimate potential distributions of all members of the complex, and evaluated the potential for suitable adjacent areas to be colonized; we also present first evaluations of potential for climate change-mediated distributional shifts. Models were developed using the GARP and Maxent algorithms.
Models for three members of the complex (T. b. brasiliensis, N = 332; T. b. macromelasoma, N = 35; and T. juazeirensis, N = 78) had significant distributional predictivity; however, models for T. sherlocki and T. melanica, both with very small sample sizes (N = 7), did not yield predictions that performed better than random. Model projections onto future-climate scenarios indicated little broad-scale potential for change in the potential distribution of the complex through 2050.
This study suggests that T. b. brasiliensis is the member of the complex with the greatest distributional potential to colonize new areas: overall; however, the distribution of the complex appears relatively stable. These analyses offer key information to guide proactive monitoring and remediation activities to reduce risk of Chagas disease transmission.
Chagas disease affects millions of people across Latin America, yet no vaccine has been developed [1, 2]. This disease is caused by a protozoan parasite, Trypanosoma cruzi, transmitted to humans mainly by blood-sucking bugs of the subfamily Triatominae . Most triatomine species are distributed in Central and South America. Brazil holds the most diverse triatomine fauna, several of which are recognized as constituting complexes of species [4–6].
Multidisciplinary studies have indicated that the Triatoma brasiliensis complex is a monophyletic group [7–9] comprising four species, one of which includes two distinct subspecies, such that the relevant taxa in the group are T. b. brasiliensis, T. b. macromelasoma, T. juazeirensis, T. melanica, and T. sherlocki. Each member of the group can be identified by external morphological characteristics and a taxonomic key was recently published . This complex is widespread in Brazil, occurring in 12 states, primarily within the caatinga and cerrado biomes . Members of the complex present distinct epidemiological importance, morphological characteristics, natural history traits, ecological requirements, genetic characteristics, and dispersal abilities [11–19].
Triatoma b. brasiliensis is the best-studied member of the complex in terms of epidemiology and spatial distribution: it occurs in six states, uses diverse food sources, and is distributed across several ecotopes. Because this species presents one of T. cruzi the highest infection rates, as well as high intradomiciliary infestation rates, it is considered the most important Chagas disease vector in the region [6, 10]. Despite significant reductions of Chagas disease incidence in Brazil in recent decades, T. b. brasiliensis and the parasites it transmits remain priority public health concerns [20, 21]. However, the transitory nature of epidemiological importance is illustrated by the case of T. sherlocki, a species described as new to science only in 2002 , and subsequently recognized as a fifth member of the T. brasiliensis complex based on genetic and ecological characteristics [8, 9]. Recently, however, this species was found colonizing human domiciles in a mining area in Bahia state, which highlights the quick changes that may occur in triatomine ecology, particularly in the face of drastic anthropogenic environmental changes . Hence, monitoring vector capacity and epidemiological importance is crucial for understanding and controlling Chagas disease transmission [23, 24].
One of the most important concerns regarding Chagas disease transmission is the degree to which vector species have the potential to colonize new areas, expanding risk areas for transmission. Some authors have suggested that degradation of sylvatic environments promotes invasion of new ecotopes by triatomines [25–27]. Another study, however, did not support the hypothesis that environmental degradation per se is the major factor in domiciliary invasions; instead, landscape patterns and characteristics of human dwellings may drive the probability that triatomines invade and colonize homes . As several studies have shown that climatic characteristics shape geographic distributions of triatomines [29–32], tools that improve mapping of present-day triatomine distributions, as well as how they respond to climate change, are crucial to proactive measures in controlling and monitoring this disease. Costa and Peterson  recently reviewed ecological niche models (ENMs) as a tool permitting better understanding along these lines; and Diniz-Filho and colleagues  analyzed geographic patterns and environmental correlates of species richness in triatomine distributions.
ENMs offer a powerful approach to characterizing and differentiating ecological niches of vectors, broadening knowledge of relevant vector and parasite-reservoir distributions, defining areas at risk for transmission, and assessing overall potential geographic distributions of vector species. ENMs have been used to understand evolutionary processes involved in differentiating members of the T. brasiliensis complex , with patterns of ecological similarity matching patterns of molecular differentiation within the complex [7, 11]. ENMs have also been useful in identifying vector-host associations  and mapping geographic distributions of triatomine species in detail [36, 37]. In spite all this research attention, however, important distributional questions remain unaddressed.
Global climate change is anticipated to have critical impacts on biodiversity, with synergistic effects producing significant changes in distributions, natural history, and behavior of many organisms [38, 39]. Climatic change has already influenced natural environments and human activities, and is likely already molding disease transmission [40–42]. ENMs have been used to make ‘forecasts’ of climate change effects on biodiversity . Hence, here, we use ENM approaches (1) to develop estimates of potential geographic distributions for each member of the T. brasiliensis complex, and (2) to evaluate potential for colonizing new areas in response to climate change based on modeled future-climate scenarios for 2020 and 2050.
Data were compiled from multiple sources: the entomological collections of Fundação Oswaldo Cruz, secretarias de saúde (health departaments) from states where the T. brasiliensis species complex occurs, and field collections by the research team of the Laboratório de Biodiversidade Entomológica - IOC - FIOCRUZ (2004–2008 [7, 8, 10, 11, 13, 14, 16]). In total, we assembled 459 occurrence points across 12 states on which to base ENMs for the complex; all specimens collected at these localities are deposited in the Entomological Collection of Instituto Oswaldo Cruz (CEIOC), spanning 1911–2008, including type specimens deposited by the authors of the species. All modern field work employed the same capture methods (i.e., capture by exhaustion).
Of the 459 occurrence points, most were readily identifiable based on external morphology and geographic range [9, 10]. However, to confirm distributional limits for members of the complex, we selected 41 localities falling between ranges of different taxa, and made additional field captures of ~2160 specimens, all of which were identified carefully via reference to taxonomic studies [4, 9, 44–46], and by comparison with specimens in the collections of the Instituto Oswaldo Cruz, U.S. National Museum of Natural History, and American Museum of Natural History. As such, we can state confidently that identifications along these range borders were made correctly.
All records were plotted on maps, and checked further for consistency with species’ known ranges. Triatoma b. brasiliensis (N = 332), T. b. macromelasoma (N = 35), and T. juazeirensis (N = 78) had sufficient points for standard niche modeling techniques (described below); however, we had only limited occurrence data for T. melanica and T. sherlocki (each N = 7), which reflects their rarity and the scarcity of their populations. Occurrence data were divided randomly into calibration and evaluation data sets using 80% and 20% of points, respectively; because of small sample sizes for T. melanica and T. sherlocki, an alternative evaluation method was used (see Model evaluation, below). For models developed to evaluate implications of scenarios of future climates for the species complex as a whole, all occurrence data for the complex were combined into a single occurrence data set; this step maximizes sample sizes and avoids dispersal-limited situations in which distributional limits are set by environmental characteristics that are difficult to model .
We used two sets of environmental data to construct models for each species. First, to estimate fine details of the current geographic distribution of each species, we used multitemporal Normalized Difference Vegetation Index (NDVI) from the Advanced Very High Resolution Radiometer (AVHRR; ~1 km native spatial resolution) sensor. These monthly NDVI images covered April 1992–March 1993. We chose imagery from this period because it best matched the bulk of our occurrence data (1992–2008). We assumed that NDVI patterns across the landscape in 1992–1993 was representative of the broader period.
Second, to develop coarse-resolution projections of the potential distributions of species under scenarios of global climate change for 2020 and 2050, we used the WorldClim data set at 1 km spatial resolution . We included seven “bioclimatic” variables: annual mean temperature, mean diurnal temperature range, maximum temperature of the warmest month, minimum temperature of the coldest month, annual precipitation, and precipitation of the wettest and driest months . Among the 19 available climatic variables within the Worldclim data set , these seven were relatively uncorrelated across the study region.
Following recent recommendations , we restricted model calibration to areas hypothesized to fall within the likely long-term movement and colonization potential (area M) of each species. To this end, for present-day distributional predictions of suitable habitat, we buffered occurrence points for each species by one half of the longest diameter of that species’ range. For climate change projections, we buffered around all points by a broader area of 500 km (i.e., combining species). After model calibration, major geographic features (e.g., São Francisco and Jequitinhonha rivers) were used to reduce M areas for each species, as members of this complex likely are not easily able to cross such barriers; we further considered neighboring distributions of species as barriers to dispersal (e.g., future zone between distributions of T. b. brasiliensis and T. b. macromelasoma; see summary in Figure 1).
We used two algorithms and applied a consensus method to model ecological niches: a genetic algorithm (GARP [51–53]) and a maximum-entropy approach (Maxent ). Both algorithms relate environmental conditions to known occurrence records, and then compare those relationships to the landscape over which observations were derived (see detailed treatment in ). We combined the resulting models by mapping areas where predictions agreed from both algorithms. This final model can be used to identify areas presenting favorable environmental conditions across the landscape. Details of each approach follow.
GARP divides occurrence data into subsets for model calibration (rule training and intrinsic testing) and model evaluation (extrinsic testing). Through an iterative process, GARP uses calibration data to develop a set of rules from a variety of methods (BIOCLIM, logistic regression, etc.). GARP then modifies rules to maximize model predictivity, testing rules for accuracy, using the intrinsic testing data to determine whether rules should be incorporated into or rejected from the model. This process continues for 1000 iterations or until convergence, and models are evaluated using the extrinsic testing data . We developed 100 model runs, with consensus via soft extrinsic omission thresholds ; all other parameters were left to default. The 10 “best” models were retained and summed to produce a final grid.
Maxent fits a probability distribution to species’ occurrences, such that entropy of the probability distribution is maximized (i.e., the probability distribution is maximally spread out), but subject to constraints imposed by the environmental characteristics of the occurrence data. Response curves (linear, quadratic, etc.) are constructed and then smoothed with regularization parameters to reduce model overfitting . For Maxent models for estimating current distributions, we used default settings for all parameters, except that we selected a random seed. For climate change scenarios, we also used default settings, except that we specified a random seed with 50% of points set aside for model evaluation.
We used a minimum training threshold  to convert raw model outputs (i.e., continuous outputs) into actual distributional estimates (i.e., a binary representation). Such thresholds , also allowed us to scale Maxent and GARP models similarly. However, we modified this threshold to consider potential for error in the occurrence data. That is, instead of the threshold that includes 100% of the training data, we sought the highest threshold that includes (100 - E)% of the training data, where E is an estimate of the proportion of occurrence data that may include locational error ; we assumed E = 5%, based on our experience with these occurrence data.
In view of recent concerns  regarding robustness of the popular receiver operating characteristic (ROC) approach , we used a partial ROC area under the curve (AUC ) for T. b. brasiliensis, T. b. macromelasoma, and T. juazeirensis. This approach is based on traditional ROC approaches , but takes into account the degree of coverage of the commission error spectrum (x-axis) by model predictions, and prioritizes omission error over commission error in calculating model robustness. All tests were developed based on the random subsets of occurrence data described above. We used a Visual Basic program  to compare model AUCs to null expectations via bootstrapping 50% of the test points 1000 times, allowing for E = 5% error among the occurrence data.
For T. melanica and T. sherlocki, the partial ROC method was not appropriate because of small sample sizes (N = 7). Instead, we used a jackknife approach, designed for models based on limited occurrence data . In this method, a set of models is built, each based on the occurrence but data omitting a different, single data point; these excluded points are used to test the predictive ability of the model, given a particular proportion of area predicted present . For each algorithm and species, then, we constructed 7 models, each based on 6 training points, setting aside one point for model testing. All resulting models were thresholded as described above. These values and proportions of suitable areas for each model were used within the program provided by Pearson and colleagues .
Captures were carried out with technical support of technicians of the Secretaria de Vigilância em Saúde and Fundação Nacional de Saúde, which are the official governmental institutions charged with monitoring and controlling triatomines in domiciles in Brazil. All field collections followed established and approved procedures of those institutions, including obtaining permission from owners of domiciles.
In all, 459 occurrence sites for the T. brasiliensis species complex were included in the analysis. Taxonomic confirmations based in detailed comparison of color patterns matched with known geographic distributions of members of the complex [9, 10]. In this sense, our data were in line with earlier results [8–10] indicating parapatric geographic distributions of members of the T. brasiliensis species complex, with the exception of marginal sympatry between T. juazeirensis and T. melanica, and between T. juazeirensis and T. sherlocki. Figure 1 shows distributional data available for the five species of the complex, as well as the M hypotheses for each species.
Distributional estimates for the three species for which adequate sample sizes were available were statistically robust (P < 0.001). Predictive power of models for T. melanica and T. sherlocki, for which only minimal occurrence data were available, however, was not significant (all P > 0.05), using the Pearson small-sample test.
Triatoma b. brasiliensis
Known occurrence data comes from the states of Maranhão, Piauí, Ceará, Rio Grande do Norte, and Paraíba, encompassing an extent of occurrence of about 1.6 × 106 km2 (Figure 1). Model predictions included parts of those states, but also areas of Tocantins, Pará, Goiás, Minas Gerais, Distrito Federal, and Bahia, as falling within the species’ distributional potential. Occurrence of T. b. brasiliensis in Pernambuco is rare , although suitable conditions extend into this state. The southeastern edge of the predicted area for this species was restricted by the distribution of T. b. macromelasoma, which likely confines or obscures any possible dispersal potential in that direction.
Triatoma b. macromelasoma
The distribution of this subspecies is limited to the state of Pernambuco [9, 10], with an extent of occurrence of about 120,000 km2 (Figure 1). This state (and the species’ distributional area) is bordered by Piauí, Ceará, Rio Grande do Norte, and Paraíba (occupied by T. b. brasiliensis), and by Sergipe and Bahia (occupied by T. juazeirensis). These separations may be mediated by the São Francisco River and/or the presence of T. juazeirensis to the south; separation from T. b. brasiliensis depends on phenomena yet to be identified [14, 19]. Triatoma b. macromelasoma is therefore constrained from population expansion on most sides by presence of other species, or via broadly unsuitable areas (to the north and southwest); however, within this defined region, a suitable area is restricted predominantly to the west of the Atlantic forest but east of the drier, more inland region of the caatinga (Figure 1).
This species has been recorded almost exclusively in Bahia, in a broad area of about 500,000 km2 (Figure 1), with only a few records from southern Pernambuco [9, 10]. The northern and western limits of its distributional area are set by the São Francisco River, and the southern limit may be set by the Jequitinhonha River. Suitable areas are limited almost completely to the caatinga region, with only minor extension into the cerrado in Minas Gerais or into the wetter Atlantic forest to the east. We note that two occurrences have been recorded north of the São Francisco River. These points may represent established resident populations; however, these bugs may have been transported there passively by humans, and may not be the result of natural colonization across the river (Figure 1).
Occurrence of this species is limited primarily to Minas Gerais, in an area of about 1.1 × 106 km2 (Figure 1). However, it is known from one location in southern Bahia, where it is sympatric with T. juazeirensis. Because this species is known from only seven sites, model accuracy was poor, and predictions regarding its distribution should be considered with caution. Our models show sparse areas of suitability across seven states: Bahia, Distrito Federal, Goiás, Minas Gerais, Mato Grosso do Sul, São Paulo, and Rio de Janeiro. In Bahia, areas possibly suitable for T. melanica overlap T. juazeirensis east of the São Francisco River (Figure 1).
Records for this species are limited to north-central Bahia, in a circumscribed area of about 40,000 km2 (Figure 1). Triatoma sherlocki is the most restricted member in the group as regards potential distributional area, and appears to be confined to the northern part of the Chapada Diamantina. It remains unclear whether T. sherlocki and T. juazeirensis are parapatric or sympatric, because both occur under the same general conditions (cracks, crevices, and holes in exposed rocks) within the caatinga. However, the area of transition between the two remains unknown, awaiting in-depth, on-the-ground surveys. Model predictions were poor for this species, so interpretation of results should be made with caution (Figure 1).
Climate change scenarios
Within the distribution of the T. brasiliensis complex, the climate change scenario explored here suggests an increase in mean temperature of 1.72°C and decrease in mean precipitation of 55.6 mm by 2050 across the region. Projections of niche model rule sets onto future climate conditions for the complex suggested little likely change in species’ distributions with changing climates (Figure 2). Near-term projections (2020) anticipate areas of retraction in the northwest, along the western edge of the distribution of T. b. brasiliensis, as well as along the eastern margin where the cerrado meets Atlantic forest. Some possible expansion may be expected in central Bahia and northern Minas Gerais. Overall, however, the distribution of the complex remains largely stable under these projections. The longer-term projections (2050) yielded more or less similar patterns, with expansion anticipated in the same areas, but without retraction (Figure 2).
Additional file 1 shows biomes and political boundaries of Brazilian states across the distribution at areas of members of the T. brasiliensis species complex.
This study presents a first comprehensive round of niche modeling and mapping of potential distributions for members of the T. brasiliensis complex, although several of our previous studies have touched on similar topics [14, 33, 35, 37]. It worth mentioning that our assumption that neighboring distributions may serve as barriers to dispersal was based on background information posed of periodic field work since 1996 [7, 8, 10, 11, 13, 14, 16]. Until now, we have never collected two species in the same location. Species distribution seems to have remained stable until now. Given our knowledge of these species and their dispersal capacity, we show all the potential areas that are suitable given the restrictions to the dispersal capacity of each. The general approach and reasoning are also presented in greater detail by authors . Model predictions for T. b. brasiliensis, T. b. macromelasoma, and T. juazeirensis were statistically significantly robust in their distributional predictions; however, models for T. sherlocki and T. melanica were unable to make predictions that were better than random. Because the five members of the complex show distinct epidemiologic importance [10, 13, 16] and dispersal abilities , these distributional patterns have considerable importance for human exposure patterns.
For example, T. melanica is a highly differentiated member of the complex [7, 9, 11, 12], and occupies sylvatic ecotopes, invading homes only sporadically; indeed, we know of no records of broad domiciliary colonization by this species [10, 45]. On the other hand, T. sherlocki, the most distinct species on all grounds [8, 10] described recently as a sylvatic species new to science , was recently found colonizing domiciles constructed relatively near natural rock piles in a mining area [16, 17]. Models for both of these species identified only restricted suitable areas, yet the epidemiological potential in real life appears to differ markedly.
Despite the fact that T. sherlocki occurs in areas close to the distribution of T. juazeirensis, it may occupy distinct microhabitats among rocky outcrops [16, 46]. Triatoma sherlocki is the only species in the complex that is flightless, although its dispersal capacity may be compensated by longer legs and possibly better walking ability . A recent study found support for high reproductive compatibility between T. sherlocki and the remaining members of T. brasiliensis species complex . What is more, laboratory progeny between the parapatric T. juazeirensis and T. sherlocki were more vigorous in terms of longevity and dispersal capacity [17, 65]. However, despite the physical proximity and genetic compatibility of these two species, no hybrids or intermediate forms have yet been collected in the field, such that factors involved in their isolation remain unknown [17–19]. Recent experiments under laboratory conditions showed that T. juazeirensis presents higher reproductive potential as regards fecundity and fertility when compared to T. sherlocki.
Triatoma juazeirensis is restricted almost entirely to Bahia, although a few specimens with this external phenotype have been collected recently in Pernambuco. These occurrences may represent genetic reassortments in a hybrid zone between T. juazeirensis and T. b. macromelasoma, leading to elevated phenotypic variability in Pernambuco. The alternative hypothesis is that T. juazeirensis phenotypes reflect recent dispersal of the species into new areas.
Our models suggest that the member of the complex with greatest opportunity to colonize new areas is T. b. brasiliensis. Model projections suggest that it could find suitable conditions in the central part of the cerrado, in Tocantins and Goiás, northwestern Minas Gerais, and western Bahia. This species is known to have high genetic variability  and broad feeding habits [13, 66], such that it may be able to explore new environmental and ecological possibilities . Therefore, the areas mentioned above should be monitored for colonization by T. b. brasiliensis. In fact, a few records for this species state have already been accumulated in Tocantins by technicians of the Fundação Nacional de Saúde/Secretaria de Saúde .
Climate and environmental changes
Models that estimate climate change effects on species’ distributions can be important tools in anticipating distributional shifts . Triatomines may be able to respond rapidly to environmental changes, which may in turn influence their biological cycles. The effects of climate on triatomines and T. cruzi transmission have been studied almost since the original description of Chagas disease: in T. infestans, warmer temperatures both accelerate embryonic development  and allow for additional generations per year . Therefore, anticipation of distributional shifts with respect to future climate scenarios is essential for effective control of Chagas disease transmission.
The T. brasiliensis complex is distributed chiefly within two biomes [9, 10, 45, 46], in habitats that are anticipated to see rather strong climate change effects, at least on vegetation . Indeed, beyond simple distributional effects that are the object of the niche models, increasing temperatures may have more direct physiological effects on triatomines, such as elevated metabolic rates, elevated egg production, etc. [72–77]. Such shifts could have important impacts on vectorial transmission of Chagas disease by members of the complex. Nonetheless, the distributional consequences anticipated by our models appear to be subtle or nearly negligible, suggesting that this complex in this region is unlikely to show large-scale distributional shifts in reaction to changing climates.
Several factors potentially may have affected the models that we have explored. A first comment concerns T. melanica and T. sherlocki. First, five of the seven occurrence points for T. melanica were collected within small towns and cities; however, the satellite products (NDVI) we used measured greenness of the landscape. In this case, then, the algorithms were likely relating the data points to reflectance patterns characteristic of urban areas. Although it is important to identify all potential habitats (including urban environments), we do not feel that our environmental data permitted our models to characterize the sylvatic environment for this species adequately. In larger data sets, as for T. b. brasiliensis or T. juazeirensis, urban-collected points would be joined by others that give a clearer view of all environmental conditions used by the species.
Second, four data points for T. sherlocki were georeferenced precisely, but perhaps incorrectly, as falling in areas uncharacteristic of their natural environments (e.g., in the center of agricultural fields away from large rock outcrops). It is possible that these samples represent coarse-spatial-resolution representations of otherwise correct locations, or passive human transport; regardless, these points probably do not correspond to established colonies, consequently affecting model results. Clearly, further documentation of the distributional ecology of this species is needed.
Third, both T. melanica[9, 10, 45] and T. sherlocki[16, 17] are essentially sylvatic (the domiciliation process for T. sherlocki may be a new event ). For the latter species, the proximity of its natural ecotope (rocky outcrops) and human domiciles (also sometimes built with rocks), combined with sylvatic hosts might have pressured this change in behavior, as has been suggested previously for other triatomines . Our data for these species may have been biased geographically because most records came from long-term (~20 years) surveys by the Secretaria de Saúde, which only surveys human domiciles. Broader searches in sylvatic environments would likely increase numbers of occurrence points for these species, and would therefore improve model predictions in the future.
One final concern relates to the two T. juazeirensis points from north of the São Francisco River. Passive transportation of bugs has been implicated in dispersal of T. infestans across the Southern Cone of Latin America. After domiciliation, insects can easily be transported by humans (e.g., in clothes or luggage [78, 79]). We suggest that T. juazeirensis occurrences north of the São Francisco River may be attributed to such human-mediated passive transportation. Other possible explanations for these occurrence points include misidentification, poor georeferencing, or unexpected phenotypes generated by hybridization and reassortment . Regardless, this river likely functions at least as a partial natural barrier to northward expansion of this species, since records north of the river are very uncommon.
Chagas disease affects millions of people across Latin America, yet this disease is still largely neglected. In eastern Brazil, the T. brasiliensis species complex represents the group of vectors with the greatest potential to transmit this disease. Until now, the extent of the potential distribution of habitat for each member of this complex was not fully characterized. Our results indicated that, of all members of this complex, T. b. brasiliensis exhibits the greatest potential to colonize new areas; however, our model projections do not anticipate much change in the distribution of the complex as a consequence of climate change. Hence, here, we present information useful in guiding proactive monitoring and remediation activities to reduce risk of further Chagas disease transmission in northeastern Brazil.
Schofield CJ, Kabayo JP: Trypanosomiasis vector control in Africa and Latin America. Parasit Vectors. 2008, 1: 24-10.1186/1756-3305-1-24. doi:10.1186/1756-3305-1-24
Hashimoto K, Schofield CJ: Elimination of Rhodnius prolixus in Central America. Parasit Vectors. 2012, 22: 45-doi:10.1186/1756-3305-5-45
Organization WH: World Health Report: Changing History. 2004, Geneva: World Health Organization, Available: http://www.who.int/tdr/research/ntd/chagas/en/. Accessed 28 August 2012
Lent H, Wygodzinsky P: Revision of the Triatominae (Hemiptera, Reduviidae), and their significance as vectors of Chagas disease. Bull Am Mus Nat Hist. 1979, 163: 123-520.
Galvão C: A checklist of the current valid species of the subfamily Triatominae Jeannel, 1919 (Hemiptera, Reduviidae) and their geographical distribution, with nomenclatural and taxonomic notes. Zootaxa. 2003, 202: 1-36.
Costa J, Lorenzo M: Biology, diversity and strategies for the monitoring and control of triatomines–Chagas disease vectors. Mem Inst Oswaldo Cruz. 2009, 1: 46-51.
Monteiro FA, Donnelly MJ, Beard CB, Costa J: Nested clade and phylogeographic analyses of the Chagas disease vector Triatoma brasiliensis in northeast Brazil. Mol Phylogenet Evol. 2004, 32: 46-56. 10.1016/j.ympev.2003.12.011.
Mendonça VJ, da Silva MT, de Araújo RF, Júnior JM, Júnior MB, Almeida CE, Costa J, Graminha MA, Cicarelli RM, da Rosa JA: Phylogeny of Triatoma sherlocki (Hemiptera: Reduviidae: Triatominae) inferred from two mitochondrial genes suggests its location within the Triatoma brasiliensis complex. Am J Trop Med Hyg. 2009, 81: 858-864. 10.4269/ajtmh.2009.08-0664.
Costa J, Correia NC, Neiva VL, Gonçalves TCM, Felix M: Revalidation and redescription of Triatoma brasiliensis macromelasoma Galvão, 1956 and an identification key for the Triatoma brasiliensis complex (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2013, 108: 785-789. 10.1590/0074-0276108062013016.
Costa J, Almeida CE, Dotson EM, Lins A, Vinhaes M, Silveira AC, Beard CB: The epidemiologic importance of Triatoma brasiliensis as a Chagas disease vector in Brazil: a revision of domiciliary captures during 1993–1999. Mem Inst Oswaldo Cruz. 2003, 98: 443-449. 10.1590/S0074-02762003000400002.
Costa J, Freitas-Sibajev MG, Marchon-Silva V, Pires MQ, Pacheco RS: Isoenzymes detect variation in populations of Triatoma brasiliensis (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 1997, 92: 459-464. 10.1590/S0074-02761997000400002.
Costa J, Orthrud MA, Marchon-Silva V, Almeida CE, Freitas-Sibajev MG, Panzera F: Morphological studies on the Triatoma brasiliensis Neiva, 1911 (Hemiptera, Reduviidae, Triatominae) genital structures and eggs of different chromatic forms. Mem Inst Oswaldo Cruz. 1997, 92: 493-498. 10.1590/S0074-02761997000400009.
Costa J, de Almeida JR, Britto C, Duarte R, Marchon-Silva V, Pacheco RS: Ecotopes, natural infection and trophic resources of Triatoma brasiliensis (Hemiptera, Reduviidae, Triatominae). Mem Inst Oswaldo Cruz. 1998, 93: 7-13.
Costa J, Peterson AT, Beard CB: Ecologic niche modeling and differentiation of populations of Triatoma brasiliensis Neiva, 1911, the most important Chagas’ disease vector in northeastern Brazil (Hemiptera, Reduviidae, Triatominae). Am J Trop Med Hyg. 2002, 67: 516-520.
Costa J, Almeida CE, Dujardin JP, Beard CB: Crossing experiments detect genetic incompatibility among populations of Triatoma brasiliensis Neiva, 1911 (Heteroptera, Reduviidae, Triatominae). Mem Inst Oswaldo Cruz. 2003, 98: 637-639. 10.1590/S0074-02762003000500009.
Almeida CE, Folly-Ramos E, Peterson AT, Lima-Neiva V, Gumiel M, Duarte R, Lima MM, Locks M, Beltrão M, Costa J: Could the bug Triatoma sherlocki be vectoring Chagas disease in small mining communities in Bahia, Brazil?. Med Vet Entomol. 2009, 23: 410-417. 10.1111/j.1365-2915.2009.00822.x.
Almeida CE, Oliveira HL, Correia N, Dornak LL, Gumiel M, Neiva VL, Harry M, Mendonça VJ, Costa J, Galvão C: Dispersion capacity of Triatoma sherlocki, Triatoma juazeirensis and laboratory-bred hybrids. Acta Trop. 2012, 122: 71-79. 10.1016/j.actatropica.2011.12.001.
Correia N, Almeida CE, Lima-Neiva V, Gumiel M, Dornak LL, Lima MM, Medeiros LM, Mendonça VJ, da Rosa JA, Costa J: Cross-mating experiments detect reproductive compatibility between Triatoma sherlocki and other members of the Triatoma brasiliensis species complex. Acta Tropica. 2013, 128: 162-167. 10.1016/j.actatropica.2013.06.019.
Costa J, Peterson AT, Dujardin JP: Morphological evidence suggests homoploid hybridization as a possible mode of speciation in the Triatominae (Hemiptera, Heteroptera, Reduviidae). Infect Genet Evol. 2009, 9: 263-70. 10.1016/j.meegid.2008.12.005.
Silveira A, Vinhaes M: Elimination of vector-borne transmission of Chagas disease. Mem Inst Oswaldo Cruz. 1999, 1: 405-411.
Borges-Pereira J, Sarquis O, Zauza PL, Britto C, Lima MM: Epidemiologia da doença de Chagas em quatro localidades rurais de Jaguaruana, estado do Ceará: soroprevalência da infecção, parasitemia e aspectos clínicos. Rev Soc Bras Med Trop. 2008, 41: 345-351. 10.1590/S0037-86822008000400005.
Papa AR, Jurberg J, Carcavallo RU, Cerqueira RL, Barata JMS: Triatoma sherlocki sp.n. coletada na Bahia, Brasil (Hemiptera, Reduviidae, Triatominae). Entomol Vect. 2002, 9: 133-146.
Luitgards-Moura JF, Vargas AB, Almeida CE, Magno-Esperança G, Agapito-Souza R, Folly-Ramos E, Costa J, Tsouris P, Rosa-Freitas MG: A Triatoma maculata (Hemiptera, Reduviidae, Triatominae) population from Roraima, Amazon region, Brazil, has some bionomic characteristics of a potential Chagas disease vector. Rev Inst Med Trop Sao Paulo. 2005, 47: 131-137.
Luitgards-Moura JF, Borges-Pereira J, Costa J, Zauza PL, Rosa-Freitas MG: On the possibility of autochthonous Chagas disease in Roraima, Amazon region, Brazil, 2000-2001. Rev Inst Med Trop Sao Paulo. 2005, 47: 45-54.
Forattini OP, Ferreira OA, Silva E, Rabello EX: Aspectos ecológicos da Tripanossomíase americana: XII - Variação regional da tendência de Panstrongylus megistus à domiciliação. Rev Saude Publica. 1978, 12: 209-233.
Romana CA, Brunstein D, Collin-Delavaud A, Sousa O, Ortega-Barria E: Public policies of development in Latin America and Chagas’ disease. Lancet. 2003, 362 (9383): 579-
Dias FBS, Quartier M, Diotaiuti L, Mejía G, Harry M, Lima ACL, Davidson R, Mertens F, Lucotte M, Romaña CA: Ecology of Rhodnius robustus Larrousse, 1927 (Hemiptera, Reduviidae, Triatominae) in Attalea palm trees of the Tapajós River Region (Pará State, Brazilian Amazon). Parasites & Vectors. 2014, 7: 154-10.1186/1756-3305-7-154.
Leite GR, dos Santos CB, Falqueto A: Influence of the landscape on dispersal of sylvatic triatomines to anthropic habitats in the Atlantic Forest. J Biogeogr. 2011, 38: 651-663. 10.1111/j.1365-2699.2010.02442.x.
Ramsey JM, Ordoñez R, Cruz-Celis A, Alvear AL, Chavez V, Lopez R, Pintor JR, Gama F, Carrillo S: Distribution of domestic triatominae and stratification of Chagas disease transmission in Oaxaca, Mexico. Med Vet Entomol. 2000, 14: 19-30. 10.1046/j.1365-2915.2000.00214.x.
Zeledón R, Ugalde JA, Paniagua LA: Entomological and ecological aspects of six sylvatic species of triatomines (Hemiptera, Reduviidae) from the collection of the National Biodiversity Institute of Costa Rica, Central America. Mem Inst Oswaldo Cruz. 2001, 96: 757-764. 10.1590/S0074-02762001000600002.
Galvão C, Jurberg J, Carcavallo RU, Segura CAM, Girón IG, Casas SIC: Distribuição geográfica e dispersão alti-latitudinal de alguns gêneros e espécies da tribo Triatomini Jeannel, 1919 (Hemiptera, Reduviidae, Triatominae). Mem Inst Oswaldo Cruz. 1998, 93: 33-38.
de Casas SI C, Carcavallo RU, Mena Segura CA, Galíndez Girón I: Bioclimatic factors of Triatominae distribution: useful techniques for studies on climate change. Entomol Vect. 1994, 1: 51-67.
Costa J, Peterson AT: Ecological niche modeling as a tool for understanding distributions and interactions of vectors, hosts, and etiologic agents of Chagas disease. Adv Exp Med Biol. 2012, 710: 59-70. 10.1007/978-1-4419-5638-5_7.
Diniz-Filho JAF, Ceccarelli S, Hasperué W, Rabinovich J: Geographical patterns of Triatominae (Heteroptera: Reduviidae) richness and distribution in the Western Hemisphere. Insect Conservation and Diversity. 2013, doi:10.1111/icad.12025
Peterson AT, Sanchez-Cordero V, Beard CB, Ramsey JM: Ecologic niche modeling and potential reservoirs for Chagas disease, Mexico. Emerg Infect Dis. 2002, 8: 662-667. 10.3201/eid0807.010454.
Sandoval-Ruiz CA, Zumaquero-Rios JL, Rojas-Soto OR: Predicting geographic and ecological distributions of triatomine species in the southern Mexican state of Puebla using ecological niche modeling. J Med Entomol. 2008, 45: 540-546. 10.1603/0022-2585(2008)45[540:PGAEDO]2.0.CO;2.
Gurgel-Gonçalves R, Galvão C, Costa J, Peterson AT: Geographic distribution of Chagas disease vectors in Brazil based on ecological niche modeling. J Trop Med. 2012, doi:10.1155/2012/705326
Hannah L, Midgley GF, Hughes G, Bomhard B: The view from the Cape: extinction risk, protected areas and climate change. BioScience. 2005, 55: 231-42. 10.1641/0006-3568(2005)055[0231:TVFTCE]2.0.CO;2.
Parmesan C, Yohe G: A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003, 421: 37-42. 10.1038/nature01286.
Epstein PR, Ferber D: Changing planet, changing health: how the climate crisis threatens our health and what we can do about it. 2011, Berkeley: University of California Press, 355-
Reiter P: Climate change and mosquito-borne disease. Environ Health Perspect. 2001, 1: 141-161.
Chaves LF, Koenraadt CJ: Climate change and highland malaria: fresh air for a hot debate. Q Rev Biol. 2010, 85: 27-55. 10.1086/650284.
Botkin D, Saxe H, Sobel M, Stockwell D: Forecasting the effects of global warming on biodiversity. BioScience. 2007, 57: 227-236. 10.1641/B570306.
Galvão AB: Triatoma brasiliensis macromelasoma n. subsp. (Reduviidade, Hemíptera). Rev Bras Malariol Doenças Trop. 1956, 7: 455-457.
Costa J, Argolo AM, Felix M: Redescription of Triatoma melanica Neiva & Lent, 1941, new status (Hemiptera: Reduviidae: Triatominae). Zootaxa. 2006, 1385: 47-58.
Costa J, Felix M: Triatoma juazeirensis sp. nov. from the state of Bahia, northeastern Brazil (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2007, 102: 87-90.
Saupe EE, Barve V, Myers CE, Soberón J, Barve N, Hensz CM, Peterson AT, Owens HL, Lira-Noriega A: Variation in niche and distribution model performance: the need for a priori assessment of key causal factors. Ecol Modell. 2012, 237–238: 11-22.
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005, 25: 1965-1978. 10.1002/joc.1276.
Jiménez-Valverde A, Nakazawa Y, Lira-Noriega A, Peterson AT: Environmental correlation structure and ecological niche model projections. Biodiversity Informatics. 2009, 6: 28-35.
Barve N, Barve V, Jiménez-Valverde A, Lira-Noriega A, Maher SP, Peterson AT, Soberón J, Villalobos F: The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol Mod. 2011, 222: 1810-1819. 10.1016/j.ecolmodel.2011.02.011.
Stockwell DRB, Noble IR: Induction of sets of rules from animal distribution data: a robust and informative method of analysis. Math Comp Simul. 1992, 33: 385-90. 10.1016/0378-4754(92)90126-2.
Pearson RG, Raxworthy C, Nakamura M, Peterson AT: Predicting species’ distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J Biogeography. 2007, 34: 102-117.
Pearson RG, Thuiller W, Araújo MB, Martinez-Meyer E, Brotons L, McClean C, Miles L, Segurado P, Dawson TP, Lees DC: Model-based uncertainty in species range prediction. J Biogeography. 2006, 33: 1704-1711. 10.1111/j.1365-2699.2006.01460.x.
Phillips S, Anderson R, Schapire R: Maximum entropy modeling of species geographic distributions. Ecol Mod. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.
Peterson AT, Soberón J, Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, Araújo MB: Ecological Niches and Geographic Distributions (MPB-49). 2011, Princeton: Princeton University Press, 336-
Stockwell SDRB, Peters D: The GARP modeling system: problems and solutions to automated spatial prediction. Int J Geogr Inf Sci. 1999, 13: 143-158. 10.1080/136588199241391.
Anderson RP, Lew D, Peterson AT: Evaluating predictive models of species’ distributions: criteria for selecting optimal models. Ecol Mod. 2003, 162: 211-232. 10.1016/S0304-3800(02)00349-6.
Phillips SJ, Dudík M: Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008, 31: 161-175. 10.1111/j.0906-7590.2008.5203.x.
Peterson AT, Papes M, Eaton M: Transferability and model evaluation in ecological niche modeling: a comparison of GARP and Maxent. Ecography. 2007, 30: 550-560. 10.1111/j.0906-7590.2007.05102.x.
Lobo JM, Jiménez-Valverde A, Real R: AUC: a misleading measure of the performance of predictive distribution models. Global Ecol Biogeogr. 2008, 17: 145-151. 10.1111/j.1466-8238.2007.00358.x.
Elith J, Graham CH, Anderson RP, Dudík M, Ferrier S, Guisan A, Hijmans RJ, Huettmann F, Leathwick JR, Lehmann A, Li J, Lohmann LG, Loiselle BA, Manion G, Moritz C, Nakamura M, Nakazawa Y, Overton JMcC M, Peterson AT, J Phillips S, Richardson K, Scachetti-Pereira R, Schapire RE, Soberón J, Williams S, Wisz MS, Zimmermann NE: Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006, 29: 129-151. 10.1111/j.2006.0906-7590.04596.x.
Peterson AT, Papeş M, Soberón J: Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol Mod. 2008, 213: 63-72. 10.1016/j.ecolmodel.2007.11.008.
Fielding AH, Bell JF: A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997, 24: 38-49. 10.1017/S0376892997000088.
Pearson RG, Raxworthy CJ, Nakamura M, Peterson AT: Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J Biogeogr. 2007, 34: 102-117.
Neiva VL, Gumiel M, Lima MM, Gonçalves TCM, Provance DW, Almeida CE, Costa J: Deposition, incubation period and hatching of eggs from Triatoma juazeirensis Costa & Felix and Triatoma sherlocki Papa, Jurberg, Carcavallo, Cerqueira & Barata (Hemiptera: Reduviidae) under laboratory conditions. EntomoBrasilis. 2012, 5: 130-136.
Alencar JE: História natural da doença de Chagas no estado do Ceará. Fortaleza. Imprensa Universitária da Universidade Federal do Ceará. 1987, 341-
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-494. 10.1590/S0074-02762011000400016.
Siqueira MF, Peterson AT: Consequences of global climate change for geographic distributions of cerrado tree species. Biota Neotrop. 2003, 3: 2-
Rabinovich JE: Vital statistics of Triatominae (Hemiptera: Reduviidae) under laboratory conditions. I. Triatoma infestans Klug. J Med Entomol. 1972, 9: 351-370.
Hack W: Estudios sobre biología del Triatoma infestans (Klug, 1834) (Hemiptera, Reduviidae). An Inst Med Regional. 1955, 4: 125-147.
Dobson A, Jolly A, Rubenstein D: The greenhouse effect and biological diversity. Trends Ecol Evol. 1989, 4: 64-68. 10.1016/0169-5347(89)90150-X.
Lucena DT: Evolução do Triatoma brasiliensis Neiva, 1911, no laboratório. Rev Bras Biol. 1960, 20: 295-302.
Perondini ALP, Costa MC, Brasileiro VLF:Biologia do Triatoma brasiliensis. II observações sobre a autogenia. Rev Saúde Públ. 1975, 9: 363-70.
Juarez E, Silva EP: Comportamento do Triatoma sordida em condições de laboratório. Rev Saúde Pública. 1982, 3: 1-36.
Silva IG: Influência da temperatura na biologia de triatomíneos. I. Triatoma rubrovaria (Blanchad, 1843) (Hemiptera, Reduviidae). Rev Goiana Med. 1985, 31: 1-37.
Almeida CE, Francischetti CN, Pacheco RS, Costa J: Triatoma rubrovaria (Blanchard, 1843) (Hemiptera-Reduviidae-Triatominae) III: patterns of feeding, defecation and resistance to starvation. Mem Inst Oswaldo Cruz. 2003, 98: 367-372. 10.1590/S0074-02762003000300012.
Almeida CE, Folly-Ramos E, Agapito-Souza R, Magno-Esperança G, Pacheco RS, Costa J: Triatoma rubrovaria (Blanchard, 1843) (Hemiptera-Reduviidae-Triatominae) IV: bionomic aspects on the vector capacity of nymphs. Mem Inst Oswaldo Cruz. 2005, 100: 231-235. 10.1590/S0074-02762005000300002.
Aragão MB: Sobre a dispersão do Triatoma infestans. Rev Soc Bras Med Trop. 1971, 5: 183-91. 10.1590/S0037-86821971000400002.
Forattini OP: Biogeografia, origem e distribuição da domiciliação de triatomíneos no Brasil. Rev Saúde Pública. 1980, 14: 265-299.
We thank the curators of the entomological collections of the U.S. National Museum of Natural History (Dr. Tom Henry) and American Museum of Natural History (Dr. Randall Schuh), for allowing us to work with their collections. FUNASA and Secretaria de Vigilância em Saúde kindly provided part of the dataset of triatomine occurrence points and support for field work. Vanessa Lima-Neiva, Elaine Folly Ramos, Nathalia Correia, Martha Locks, and Maria Beltrão (coordinator of Projeto Central, Bahia) provided essential help in the field.
Financial support was provided by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, process numbers 2010/17027-0 and 2011/22378-0). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
Conceived the study: JC and ATP. Collected data: JC and CEA. Identified samples: JC; Model and figure production: LLD and ATP. Wrote the manuscript: JC, LLD, CEA, and ATP. All authors read and approved the final version of the manuscript.
Electronic supplementary material
Additional file 1: Map showing biomes and political boundaries of Brazilian states across the distribution at areas of members of the Triatoma brasiliensis species complex. AL = Alagoas, AP = Amapá, BA = Bahia, CE = Ceará, DF = Distrito Federal, ES = Espírito Santo, GO = Goiás, MA = Maranhão, MT = Mato Grosso, MS = Mato Grosso do Sul, MG = Minas Gerais, PR = Paraná, PB = Paraíba, PA = Pará, PE = Pernambuco, PI = Piauí, RJ = Rio de Janeiro, RN = Rio Grande do Norte, SE = Sergipe, SP = São Paulo, TO = Tocantins. (PDF 609 KB)
About this article
Cite this article
Costa, J., Dornak, L.L., Almeida, C.E. et al. Distributional potential of the Triatoma brasiliensis species complex at present and under scenarios of future climate conditions. Parasites Vectors 7, 238 (2014). https://doi.org/10.1186/1756-3305-7-238