- Open Access
Potential impact of climate change on the geographical distribution of two wild vectors of Chagas disease in Chile: Mepraia spinolai and Mepraia gajardoi
Parasites & Vectors volume 12, Article number: 478 (2019)
Mepraia gajardoi and Mepraia spinolai are endemic triatomine vector species of Trypanosoma cruzi, a parasite that causes Chagas disease. These vectors inhabit arid, semiarid and Mediterranean areas of Chile. Mepraia gajardoi occurs from 18° to 25°S, and M. spinolai from 26° to 34°S. Even though both species are involved in T. cruzi transmission in the Pacific side of the Southern Cone of South America, no study has modelled their distributions at a regional scale. Therefore, the aim of this study is to estimate the potential geographical distribution of M. spinolai and M. gajardoi under current and future climate scenarios.
We used the Maxent algorithm to model the ecological niche of M. spinolai and M. gajardoi, estimating their potential distributions from current climate information and projecting their distributions to future climatic conditions under representative concentration pathways (RCP) 2.6, 4.5, 6.0 and 8.5 scenarios. Future predictions of suitability were constructed considering both higher and lower public health risk situations.
The current potential distributions of both species were broader than their known ranges. For both species, climate change projections for 2070 in RCP 2.6, 4.5, 6.0 and 8.5 scenarios showed different results depending on the methodology used. The higher risk situation showed new suitable areas, but the lower risk situation modelled a net reduction in the future potential distribution areas of M. spinolai and M. gajardoi.
The suitable areas for both species may be greater than currently known, generating new challenges in terms of vector control and prevention. Under future climate conditions, these species could modify their potential geographical range. Preventive measures to avoid accidental human vectorial transmission by wild vectors of T. cruzi become critical considering the uncertainty of future suitable areas projected in this study.
Chagas disease, or American trypanosomiasis, is an endemic vector-borne disease affecting between six and eight million people, with an attributed number of deaths of approximately 12,000 per year worldwide . Its causative agent is the protozoan Trypanosoma cruzi, transmitted by hematophagous insects of the subfamily Triatominae to mammals . In Chile, the current human prevalence of Chagas disease is 0.7%, with 0.6% and 1.5% in urban and rural zones, respectively . The vector species present in Chile are Triatoma infestans, Mepraia gajardoi, M. parapatrica and M. spinolai . These triatomines occur in rural and suburban zones from 18°30′S to 34°36′S .
Mepraia gajardoi is currently detected on the northern coastal zones where the arid climate is common, whilst M. spinolai can be found in valleys in the Mediterranean-semiarid climate zones . The domestic vectorial transmission by the domestic vector T. infestans was declared interrupted in 1999; however, wild vectors are still a problem in rural areas of Chile . The sylvatic M. spinolai is very abundant in stone quarries of periurban zones  where it feeds on wild rodents, goats, dogs, cats, rabbits and humans [8, 9], and several home invasion complaints are notified to the authorities every year (data requested from http://www.portaltransparencia.cl). Mepraia gajardoi is abundant near seaweed collector settlements, where it preferably feeds on sea birds, marine mammals, lizards, dogs, cats and humans [10, 11]. These situations are epidemiologically relevant, especially considering that the prevalence of T. cruzi in M. spinolai populations can reach up to 76.1%  and 27.0% for M. gajardoi . Under this scenario, it is necessary to explore potential areas where these species can be detected, increasing information on habitat preferences, niche requirements and geographical distributions.
In general, a disease transmission system is composed of a set of species that interact in hosting and vectoring a pathogen in space , and all the component species have a unique biogeography related to its ecological niche (i.e. those conditions appropriate for its survival and reproduction), interactions with other species and accessible areas over time [14, 15]. Therefore, a large part of disease transmission risk corresponds to the intersection of the distributional areas of the species involved .
Species distribution models (SDMs) are proper tools to examine the potential geographical distribution of species . These models correlate the georeferenced occurrences of a species with environmental information (for instance, climatic information), identifying suitable areas for the survival of its populations given its niche requirements . Climate is key to understanding the geographical distribution of species at a large spatial scale [18,19,20]. Niche conservation, i.e. the tendency of species to maintain their niche requirements over space and time [21, 22], is frequently assumed in SDMs. If the niche is conserved, it is possible to project the niche requirements to other geographical regions or to other times, past or future [21, 22]. Under this assumption, SDMs have been extended to the study of infectious diseases and host, parasite, reservoir and/or vector modelling [20, 23]. For example, the geographical distributions of West Nile virus vectors have been modelled , the spatial dynamics of dengue vectors and human dengue cases , among many others. Chagas disease risk has also been modelled [26, 27], as well as the reservoirs of T. cruzi  and some of its vectors [29,30,31,32,33,34,35,36]. SDMs are also useful to project the potential distribution under climate change conditions, which could help to identify spatial changes of infectious diseases [37, 38]. Modelling under future climate conditions can be important for making decisions about control and disease surveillance, anticipating appropriate measures . The aim of this study is to estimate the potential geographical distribution of two wild vectors of Chagas disease in Chile, M. spinolai and M. gajardoi, under current and future climate scenarios.
This study was carried out considering the current distribution of both Mepraia species that includes desert, matorral and steppe ecoregions in the Pacific side of the Southern Cone of South America  (Additional file 1: Figure S1, modified from ). We chose this criterion because the wild vectors of Chagas disease included in this study are considered endemic for those ecoregions .
Triatomine species and occurrence data
We only used georeferenced occurrences data to allow their association with environmental data (Additional file 2: Datasets S1 and S2). Our sources were: field databases collected by different research groups between 2008–2016 (unpublished data); 19 published scientific articles with georeferenced locations (Additional file 3: Table S1); collections of the Museo Nacional de Historia Natural (Santiago, Chile) and the Museo Entomológico de la Universidad Metropolitana de Ciencias de la Educación (Santiago, Chile); and reported sightings of triatomines corresponding to notifications of house or peridomicile intrusion informed by rural communities to public health centres, which included the insect specimen, obtained through the governmental website (http://www.portaltransparencia.cl). All specimens were identified to the species level by researchers with entomological training using taxonomical keys [4, 6]. The compiled database included 790 occurrences for M. spinolai and 19 for M. gajardoi. Geographical duplicated occurrences, considering 1 km2 for each occurrence point, were removed from the database using NicheToolBox (http://shiny.conabio.gob.mx:3838/nichetoolb2/), so the available database for modelling included 151 occurrences for M. spinolai and 13 for M. gajardoi (Additional file 2: Datasets S1 and S2).
We used climate data from WorldClim v.1.4 because it includes both current and future climate conditions (http://worldclim.org/version1) . We selected a subset of five of the 19 available climate variables within the WorldClim dataset, based on expert knowledge of the biology of these vectors [27, 42,43,44]. In the models, we included annual trends (Bio12: annual precipitation), seasonality (Bio7: temperature annual range) and limiting or extreme environmental factors (Bio2: mean temperature diurnal range; Bio10: mean temperature of warmest quarter; Bio11: mean temperature of coldest quarter) . We used the data at 30 arc-second spatial resolution (approximately 1 km2 at the equator).
Based on the hypothesis of accessible areas by dispersal over relevant time periods (M, hereafter) for both Mepraia species [15, 45], we set 2° (~222 km) buffer areas around each occurrence. The genus Mepraia is a dispersal restricted group, given that most developmental stages disperse by walking, unlike winged males that can also fly [4, 43, 46]; therefore, the accessible areas used were conservative compared to previous reports [e.g. 35, 37]. This buffer was constrained by the Andes and the Pacific Ocean, which constitute geomorphological features that Mepraia species would not naturally overcome . These areas were used to extract the background data for modelling each species niche (Figs. 1, 2). We analysed the correlation among the selected variables within these calibration areas using 10,000 random points plus occurrences in R project v.3.4.1 (Stats and corrplot packages) .
The areas considered for projecting included the whole geographical area delimited by the administrative boundaries of Chile and Peru, using the same geomorphological criterion explained before . For climate change projections, general circulation models (GCMs) were used, which resulted from the fifth phase of the Coupled Model Intercomparison Project CMIP5 , part of the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC) . Each GCM exhibits different climate scenarios or representative concentration pathways (RCPs), which represent different ways of greenhouse gas concentrations resulting in radiative forcing, depending on human activity mitigation . Scenarios are: (i) RCP 2.6, considered an optimistic scenario for maximum mitigation; (ii) RCP 4.5 and (iii) RCP 6.0, considered stabilization scenarios, in which there is increased radiative forcing but it stabilizes; and (iv) RCP 8.5, the most pessimistic scenario, in which higher levels of radiative forcing are recorded and the increase in temperature does not stop . For future scenarios, we used five GCMs with the climatic projections for 2070: CCSM4 (CC, hereafter), GISS-E2-R (GS, hereafter), HadGEM-ES (HE, hereafter), IPSL-CM5A-LR (IP, hereafter) and MIROC-ESM (MR, hereafter). Each of the selected GCM has the four RCPs available . These RCPs have 19 bioclimatic variables in raster format. The RCPs were obtained from climate model simulations planned as part of the World Climate Research Programme’s Fifth Coupled Model Intercomparison Project (CMIP5) [48, 50].
Species distribution modelling
SDMs were constructed using Maxent v.3.3.3k [52, 53]. This software predicts species distribution, correlating sets of environmental predictors across a user-defined landscape that is divided into grid cells, with species occurrences . Maxent is a machine-learning method that assesses the probability of a species distribution, by estimating a probability of distribution function of maximum entropy . The method uses randomly selected pseudo-absences within an a priori-defined background area . Maxent generally performs better than other software commonly used for SDMs; it has been widely used [53, 55, 56], and its ability to make predictions is supported when occurrence records are scarce [57, 58], as in this case. Maxent’s average performance was among the first algorithms tested in the Classic BAM (Biotic Abiotic Mobility) scenario and also for species with narrow niche breadth, but its results were variable depending on the measuring test . No algorithm will be optimal under all circumstances because there is not a single best approach .
To run Maxent we used 10,000 background points, added samples to background, set auto features, and the regularization multiplier was set to 1. We allowed the creation of response curves, ran jackknife to measure variable importance (i.e. percent contribution and permutation importance ), activated Logscale raw/cumulative pictures, background predictions and verbose. Taking into account the number of occurrences of both species, cross-validation was replicated 15 times for M. spinolai and only three times for M. gajardoi; these numbers of replicates (k) corresponded to the number of independent subsets, used for training (k minus 1 subset), and evaluation of the model (using the remaining subset), with a maximum of 500 iterations. We used the average distribution map as output with logistic format for data values. Assuming niche conservation, the projections were performed using strict model transference between calibrated models and their current and projected future climate conditions. To this end, we deactivated clamping and extrapolation in Maxent, as recommended . The performance of the models was evaluated using the area under the receiver operating characteristic curve (AUC) . The AUC ranges from 0.5 for a model that performs no better than chance to 1.0 if the model fits perfectly. The AUC quantifies the degree to which the model identifies presences more accurately than random . The constructed model was projected to current and future conditions.
We performed multivariate environmental similarity surface (MESS), which represents how similar a point is to a reference set of points with respect to a set of predictor variables. MESS provides an index of environmental similarity between each pixel and the median of the most dissimilar variable in M [62, 63]. We also performed mobility-oriented parity (MOP) analysis, which identifies areas of strict extrapolation and calculates environmental similarity between the calibration (M) and projection regions. MOP calculates multivariate distances from environmental variables associated with points across the projection region to a user-specified proportion of the environmental variables associated with points in the calibration region ; in our case, we used 50%. Both MESS and MOP analyses were performed in NicheToolBox .
For binary prediction of suitability (presence/absence) we used the minimum training presence threshold ; that is, the lowest predicted suitability value corresponding to an occurrence of each species in their respective M . Once the binary predictions were obtained in QGIS Desktop v.2.18.19, we used two different approaches to calculate the suitability for each RCP. The first approach used the Boolean logic operator “OR”, allowing that if at least one of the GCM results for the same RCP resulted in a suitable cell, it was considered suitable for that RCP; this was the higher public health risk situation. The second approach involved the Boolean logic operator “AND”, forcing suitability in all five GCM results for the same RCP to retain the cell as suitable, conforming the lower public health risk situation. After this, we compared the resulting raster of the RCPs from each approach with the current potential distribution, obtaining two maps per RCP showing stable areas (present in the current potential distribution and in the resulting future climatic model), retraction areas (present in the current potential distribution but absent in the future projection) and expansion areas (absent in the current distribution but present in the modelled future scenario). The areas corresponding to each category were calculated using r.report tool in QGIS Desktop v.2.18.19. Maps were elaborated using publicly available shapefiles , our compiled occurrences, and the processed outputs of Maxent and NicheToolBox, in QGIS Desktop v.2.18.19. A complete summary of all used methods is depicted in a flow diagram included as Additional file 4: Figure S2.
Variables in the models
The final model for M. spinolai included the five environmental variables mentioned above. The correlation coefficients are shown in Additional file 5: Table S2. Bio12 had a percent contribution of 60.6 and a permutation importance of 41.7, followed by Bio10 (30.8; 39.7), Bio7 (6.0; 9.9), Bio11 (4.0; 7.6) and Bio2 (1.2; 1.2). The final model for M. gajardoi included all the above variables except for Bio10, since this variable showed 0% of contribution in the first model (i.e. first run including the same variables used for M. spinolai), so it was not included in the following analyses. Bio2 had a percent contribution of 94.4 and a permutation importance of 37.6, followed by Bio11 (3.0; 21.3), Bio12 (2.2; 40.7) and Bio7 (0.4; 0.4). The correlation coefficients are shown in Additional file 5: Table S3 and the response curves’ figures for these variables are available in the repository.
Validation of SDMs
The average AUC value for the model of M. spinolai was 0.878 with standard deviation 0.055. In the case of M. gajardoi, the average value of AUC was 0.984 and standard deviation 0.015.
Extrapolation risk assessment
The resulting figures of MESS and MOP analyses are available in the repository, which show that there are no areas with strict extrapolation (i.e. with climate values outside the range of those in the calibration region). Both MESS and MOP for M. gajardoi showed a thin coastal band in the projection areas of Peru similar to its calibration area. In Chile, the similarity of the projection is concentrated in the north-central part of the country. For M. spinolai, both analyses showed more extended similar areas encompassing the Andes not including the Peruvian Amazonia. In Chile, the similarity areas showed a restriction in the northern part related to the desert, extending its distribution further south.
Current and future SDMs for M. spinolai and M. gajardoi
The minimum suitability value of M. spinolai occurrences used as threshold for its binary prediction was 0.0126. All probability values below that number were considered non-suitable areas. The suitability area predicted under current climatic conditions for M. spinolai (317,580 km2) showed that the potential areas have similarities to its known distribution, including the inland valleys and, to a lesser extent, coastal areas and the Andes . However, the current potential distribution extended to the north. This species showed an increased suitability inland in the Atacama Region but a reduction of suitability in the Atacama Desert, with some potential presence north in Peru (Fig. 1). It slightly extended south, encompassing approximately from 7° to 35°S.
The threshold value for binary prediction of M. gajardoi was 0.0407. All probability values above or equal to that number were considered suitable areas. The current potential distribution projected for M. gajardoi (42,727 km2) included its known distribution but extended through the coast of Peru. Figure 2 shows that it is a coastal species, with some inland suitable areas, but always near the coast, ranging from 8° to 27°S.
The suitability areas in all future climate scenarios for M. spinolai and M. gajardoi are shown in Figs. 3, 4, 5, 6, 7, 8, 9 and 10, respectively. Tables 1 and 2 show in detail the stable, retraction and expansion areas in future scenarios, compared to the current potential distribution of M. spinolai and M. gajardoi, respectively.
For M. spinolai in the future higher risk public health situation, the different RCPs showed similar potential distributions, maintaining most of its known distribution, but with new areas expanding this species range to the north to the arid diagonal of the Atacama Desert (between 24° and 27°S) and to the south (c.36°S), with some retraction in the coast of Peru. On the other hand, the lower risk situation showed a retraction of potential suitable areas, with higher losses in the 6.0 scenario (Table 1). Retraction occurred to the south in all scenarios modelled, losing more areas while increasing the scenarios’ severity. The arid diagonal area appears to be a barrier to dispersion, showing retraction areas in all scenarios except for RCP 2.6. The current potential distribution in Peru was reduced in all scenarios but in different magnitudes (Figs. 3, 4, 5, 6).
The potential higher risk situation in the distribution of M. gajardoi in the future showed some new suitable inland areas near the coast that seem to be consistent among RCP, and other potential future distribution areas further south of its known distribution by the coast. In general, the potential future distributions in the four scenarios were similar to the current potential suitable area but showed retraction in the potential areas by the coast of Peru, especially in the most pessimistic scenario (RCP 8.5). In the lower risk situation, this retraction tendency in Peru was evident, even in the most optimistic scenario (RCP 2.6). An overall reduction trend of potential suitable Chilean areas was detected in all scenarios in the lower risk situation (with different magnitudes regarding RCP), showing a substantial loss of suitable areas in both countries in the most pessimistic scenario (Figs. 7, 8, 9, 10).
In this study, we collected data from different sources on the occurrence of M. spinolai and M. gajardoi, endemic triatomine species from Chile, to assess their current climate suitability and predict their potential distribution under future climate scenarios. Both models performed significantly different than chance.
The SDM built for M. spinolai was mainly associated with annual precipitation (Bio12) and mean temperature of the warmest quarter (Bio10), both contributing with more than 90% in the model. Precipitation above 400 mm decreased its occurrence probability, with the maximum suitability values about 100 mm. High rainfall might be creating an excessively humid microenvironment promoting entomopathogenic fungus growth ; therefore, this would probably decrease vector fitness. On the other hand, moderate precipitation may be indirectly related to the presence of M. spinolai, because rainfall promotes primary productivity (i.e. plant communities) and this, in turn, positively affects the abundance of small rodent species such as Octodon degus and Phyllotis darwini, hosts of M. spinolai [69, 70]. Then, the increase in host populations would have a positive effect on vector population demographic parameters.
Mean temperature of the warmest quarter is probably relevant because M. spinolai inhabits areas with marked seasonality that affects this ectothermic species, forcing reproduction during the warm season [7, 71]. Temperature is related to triatomine dispersal to locate hosts and mates, and diapause induction [44, 47]. The maximum suitability for this variable was at 24 °C; however, under 10 °C and above 24 °C, suitability drops abruptly. The suitability maps for T. infestans in Chile and Peru  are quite similar to ours for M. spinolai. Temperature annual range, mean diurnal range and mean temperature of the coldest quarter also contributed to M. spinolai model, showing similarities to the models reported for several triatomine species, including T. infestans .
Regarding the suitability areas shown in the model, the current potential distribution of M. spinolai showed two main areas: one limited to semiarid areas from 27° to 34°S, and another from 8° to 27°S including the Atacama Desert and some areas of Peru. The first area includes most of the known distribution of this species . In the second main area, no M. spinolai occurrences have been notified, but it is a suitable climate zone to invade according to our model. The arid diagonal of the Atacama Desert has probably prevented its colonization. Given that T. infestans colonized houses in this northern area  and considering that these two species have been detected coexisting under the same sylvatic environmental conditions , it is plausible that M. spinolai could be currently present but undetected or could eventually invade this area [73, 74]. In fact, our current potential distribution model for M. spinolai in Chile is very similar to that recently modelled for T. infestans in this country .
The model of M. gajardoi was mainly related to mean diurnal range (Bio2), which is the mean of the monthly minimum temperature subtracted from the maximum temperature, showing a decrease in the suitability under wide thermal oscillation. Then, the model is explained mainly by temperature. The current potential distribution of M. gajardoi included mainly coastal areas from 8° to 27°S, characterized by the Pacific anticyclone effect that prevents most rainfall and high temperature differences . Therefore, this vector species would be adapted to areas with scarce thermal and humidity fluctuations, preventing its spread to the Atacama Desert, characterized by wide daily thermal variations and excessive dryness.
When analysing the overlap between the current potential distribution from both species, we detected an overlap area encompassing around 24°40′ to 26°35′S by the coast, an area presenting another Mepraia species, M. parapatrica. In this zone, it has been reported that climatic and ecologic characteristics are intermediate from those where M. spinolai and M. gajardoi inhabit . In our study, we confirm that the climate of this area would be suitable for both M. gajardoi and M. spinolai, without discarding that hybridization could be occurring . There is another area where both species’ current potential distributions overlap, mainly near the coast of Peru, but neither of these species have been reported. We suggest that this area should be sampled, to assess if any member of the Mepraia genus are present, especially considering that this area is currently colonized by triatomine species from the Triatoma and Panstrongylus genus . Therefore, the ecological conditions, including biotic and abiotic factors, might be conducive for triatomine colonization.
Regarding RCPs 2.6, 4.5, 6.0 and 8.5, it is predicted that the temperature will increase an average of 1.0, 1.8, 2.2 and 3.7 °C, respectively . For Chile, the increase in temperature under a severe scenario will be modest in the coast (0.5 to 1.0 °C) but will increase towards the Andes, where it may rise up to 5 °C in north and central Chile . It is expected that these arid and semi-arid zones (that encompass most of the potential distribution areas of these species) will show modest precipitation changes, but the south will present a significant reduction in rainfall and an increase in temperature . These climate modifications would promote the future expansion of the potential distribution of M. spinolai. For the coastal areas used by M. gajardoi, low thermal oscillation and scarce precipitation change are expected in the future . A previous study has shown that insect species living at higher latitudes, as our study subjects, have broader thermal tolerance and live in climate cooler than their physiological optimal; therefore, climate warming may enhance their fitness performance . Our results for both M. spinolai and M. gajardoi show an increase in all future scenarios when modelling the higher public health risk situation, and a reduction in the suitable areas when modelling the lower risk situation, so there is consistency in the approaches utilized. Transmission risk models are useful and predictive if they can anticipate the effects of climate changes on spatial risk patterns, so this detected uncertainty must be taken into account in future predictions . In our study, we considered health risk using vector distribution modelling according to climate. To be more informative for public health risk assessments, future studies should include information on vector competence and/or infection prevalence, in case these parameters are affected by climate.
In this study, we used interpolated climate data. These data were obtained from available climate stations; however, the density of these is relatively low in the study area and therefore these interpolated data may not necessarily reflect its exact climate, particularly in mountainous areas . Remotely sensed data may be more informative than ground-based data sets ; in our case, in spite of its limitations, we used the Worldclim dataset because it includes future projections needed for future modelling, but for new studies of the potential current distribution of these vectors, data from remote sensors may be preferable. Our projections consider climate at a regional scale, but it is well known that microrefugia are present (i.e. sites that support locally favourable climate within unfavourable regional climate), which may allow populations to persist outside their main distributions . Then, it is possible our modelled distributions could be constrained and, therefore, underestimated.
A limitation of our models relates to their validation. To construct the SDMs, we used all the occurrences reported in the literature and all those provided by colleagues and the Ministry of Health; therefore, we were unable to validate them with additional independent occurrences (locations not used for modelling). Instead, we relied on the species identification accuracy, ensuring no false positive records were included as occurrences. Even though the performance of our models was evaluated by the AUC, a controversial metric for climatic niche modelling, we consider our results are relevant in terms of current and future public health risk. In addition, MESS and MOP analyses were used to minimize extrapolation errors to assess predictions . In our study, the suitability areas in the models were mostly restricted to environmental zones similar to the calibration areas of each species; therefore, supporting our models.
The areas indicated as suitable in the current potential distribution could be inspected to evaluate if these species are currently present. However, the lack of species presences could be the result of constraints in their connectivity to present or past populations, parameter directly dependent on dispersal . Both Mepraia species present adult alary polymorphism , so their ability to surpass unsuitable areas to reach further suitable ones would be impaired, compared to other winged triatomine species. This would prevent colonization when suitable patches are not contiguous unless they are passively transported, as occurred with T. infestans . An additional constraint to occupancy is microhabitat quality. These sylvatic species require rocky outcrops, bromeliads or rock crevices to establish, with the appropriate vertebrate feeding sources to maintain triatomine colonies [4, 8, 10, 13, 43, 70, 82, 83]. One last consideration when sampling to detect vectors in the field is the sensitivity of the used methods. Assuming detection is imperfect, areas could be determined as negative when in fact they may contain triatomines, so this should also be considered for model validation . This was not an issue in the construction of our models, because Maxent only requires presences as input. However, occurrences’ sampling bias may have partially influenced our results, as we cannot assure that the whole accessible areas were sampled to obtain the input occurrences for this study, and we did not explicitly model the probability of sampling a location in the calibration areas . Future studies should consider all these biases when projecting species distributions.
We assume that the construction of our models considered all the relevant climatic variables for Mepraia species; however, at this point we cannot discard that relevant variables for these SDMs were overlooked. Future studies on these vector species should also consider physiological data to predict species ranges. In this study, we used two approaches to assess potential suitability areas in the future. The lower risk situation is a conservative approach, because it would not overestimate but could underestimate the future potential distribution of both species. Conversely, the higher risk situation could show the opposite limitation. Furthermore, current weather observations (of 30 years as minimum) have not been compared with the baseline of the GCMs, at least for Chile, to assess which model fits better for this region. To account for this caveat, we included the results of five different GCMs for each RCP. Currently, species of the genus Mepraia seem to be constrained by the Pacific Ocean to the west and the Andes mountain range to the east. In these areas, several different climates are present in short stretches of land. For example, the area from 18° to 25°S is subdivided into three types of climates: coastal desert, interior desert and high desert , implying that climate change projections in face of this geographical complexity may present strong representation problems.
Finally, if triatomines conserve their climatic niche into the future, i.e. niche conservatism, species ought to migrate or become extinct when facing adverse conditions . However, other alternatives must be considered, whether species can modify its climatic niche due to phenotypic plasticity or adaptive microevolution . In our study that topic was not covered, but future studies should consider the fact that phenotypic plasticity can be an attenuating factor of climate change effects, and assist in short- and long-term adaptation [86, 87].
The climate requirements of M. spinolai are mostly related to annual precipitation and mean temperature of the warmest quarter. Its potential distribution was projected to semi-arid and Mediterranean climates, fitting the known distribution of this species, located in inland valleys of arid and semiarid climate and, to a lesser extent, along the coast. The SDM of M. gajardoi was built mainly with mean temperature diurnal range. The potential distribution of M. gajardoi was limited to areas with coastal desert climate of Chile and Peru. Additional sampling and surveillance are required to corroborate if these vectors are inhabiting areas detected as currently suitable by our models but still have no reports of their presence. Under climate change scenarios, M. spinolai and M. gajardoi would mainly conserve their current distributions, and expand to zones with similar climates. However, the biotic features (e.g. host and refuge availability) in these areas could be restricting factors that will need to be considered when generating vector control programs. Prevention campaigns should be anticipated in areas where historically no wild triatomines have been reported, and control measures must be established and/or reinforced in areas where these triatomines are expected to persist.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The data generated and/or analysed during the present study are available in the Figshare repository, https://figshare.com/s/79896af2b019b2eac0a7.
representative concentration pathways
species distribution model
Intergovernmental Panel on Climate Change
general circulation models
Fifth phase of the Coupled Model Intercomparison Project
Fifth Assessment Report
area under the curve
Cucunubá ZM, Okuwoga O, Basáñez MG, Nouvellet P. Increased mortality attributed to Chagas disease: a systematic review and meta-analysis. Parasit Vectors. 2016;9:42.
Panzera F, Dujardin JP, Nicolini P, Caraccio MN, Rose V, Tellez T, et al. Genomic changes of Chagas disease vector, South America. Emerg Infect Dis. 2004;10:438–46.
Canals M, González C, Canals L, Canals A, Caceres D, Alvarado S, et al. What do the numbers tell us about the temporal evolution of Chagasʼ disease? Rev Chil Infectol. 2017;34:120–7.
Frías-Lasserre D. A new species and karyotype variation in the bordering distribution of Mepraia spinolai (Porter) and Mepraia gajardoi Frías et al. (Hemiptera:Reduviidae: Triatominae) in Chile and its parapatric model of speciation. Neotrop Entomol. 2010;39:572–83.
Bacigalupo A, Torres-Pérez F, Segovia V, García A, Correa JP, Moreno L, et al. Sylvatic foci of the Chagas disease vector Triatoma infestans in Chile: description of a new focus and challenges for control programs. Mem Inst Oswaldo Cruz. 2010;105:633–41.
Frías DA, Henry AA, González CR. Mepraia gajardoi: a new species of Triatominae (Hemiptera: Reduviidae) from Chile and its comparison with Mepraia spinolai. Rev Chil Hist Nat. 1998;71:177–88.
Cattan PE, Pinochet A, Botto-Mahan C, Acuña MI, Canals M. Abundance of Mepraia spinolai in a periurban zone of Chile. Mem Inst Oswaldo Cruz. 2002;97:285–7.
Canals M, Cruzat L, Molina MC, Ferreira A, Cattan PE. Blood host sources of Mepraia spinolai (Heteroptera: Reduviidae), wild vector of Chagas disease in Chile. J Med Entomol. 2001;38:303–7.
Molina MC, Cattan PE, Canals M, Cruzat L, Aguillón JC, Ferreira A. A simple immunometric assay to assess the feeding habits of Mepraia spinolai, a Trypanosoma cruzi vector. Parasitol Res. 2004;92:375–9.
Sagua F, Araya J, González J, Neira I. Mepraia spinolai in the Southeast Pacific Ocean Coast (Chile)—first insular record and feeding pattern on the Pan de Azúcar Island. Mem Inst Oswaldo Cruz. 2000;95:167–70.
González CR, Reyes C, Canals A, Parra A, Muñoz X, Rodríguez K. An entomological and seroepidemiological study of the vectorial-transmission risk of Chagas disease in the coast of Northern Chile. Med Vet Entomol. 2015;29:387–92.
Coronado X, Rozas M, Botto-Mahan C, Ortíz S, Cattan PE, Solari A. Molecular epidemiology of Chagas disease in the wild transmission cycle: the evaluation in the sylvatic vector Mepraia spinolai from an endemic area of Chile. Am J Trop Med Hyg. 2009;81:656–9.
Toledo A, Vergara F, Campos R, Botto-Mahan C, Ortiz S, Coronado X, et al. Trypanosoma cruzi genotypes in Mepraia spinolai from wild ecotopes in Northern Chile. Am J Trop Med Hyg. 2012;88:285–8.
Peterson AT. Mapping disease transmission risk: enriching models using biogeography and ecology. Baltimore: Johns Hopkins University Press; 2014.
Barve N, Barve V, Jiménez-Valverde A, Lira-Noriega A, Maher SP, Peterson AT, et al. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol Model. 2011;222:1810–9.
Elith J, Leathwick JR. Species distribution models: ecological explanation and prediction across space and time. Annu Rev Ecol Evol Syst. 2009;40:677–977.
Soberón J, Nakamura M. Niches and distributional areas: concepts, methods, and assumptions. Proc Natl Acad Sci USA. 2009;106:19644–50.
Pearson RG, Dawson TP. Predicting the impacts of climate change on the distribution of species: are bioclimatic envelopes useful? Glob Ecol Biogeogr. 2003;12:361–71.
Soberón J. Grinnellian and Eltonian niches and geographic distributions of species. Ecol Lett. 2007;10:1115–23.
Peterson AT. Ecological niche conservatism: a time-structured review of evidence. J Biogeogr. 2011;38:817–27.
Wiens JJ, Graham CH. Niche conservatism: integrating evolution, ecology, and conservation biology. Annu Rev Ecol Evol Syst. 2005;36:519–39.
Wiens JJ, Ackerly DD, Allen AP, Anacker BL, Buckley LB, Cornell HV, et al. Niche conservatism as an emerging principle in ecology and conservation biology. Ecol Lett. 2010;13:1310–24.
Peterson AT. Ecologic niche modeling and spatial patterns of disease transmission. Emerg Inf Dis. 2006;12:1822–6.
Peterson AT, Komar N, Komar O, Navarro-Sigüenza AG, Robbins MB, Martínez-Meyer E. West Nile virus in the New World: potential impacts on bird species. Bird Conserv Int. 2004;14:215–32.
Peterson AT, Martínez-Campos C, Nakazawa Y, Martínez-Meyer E. Time-specific ecological niche modeling predicts spatial dynamics of vector insects and human dengue cases. Trans R Soc Trop Med Hyg. 2005;99:647–55.
Sarkar S, Strutz SE, Frank DM, Rivaldi C-L, Sissel B, Sánchez-Cordero V. Chagas disease risk in Texas. PLoS Negl Trop Dis. 2010;4:e836.
Tapia-Garay V, Figueroa DP, Maldonado A, Frías-Lasserre D, González CR, Parra A, et al. Assessing the risk zones of Chagas’ disease in Chile, in a world marked by global climatic change. Mem Inst Oswaldo Cruz. 2018;113:24–9.
Peterson AT, Sánchez-Cordero V, Beard CB, Ramsey JM. Ecologic niche modeling and potential reservoirs for Chagas disease, Mexico. Emerg Inf Dis. 2002;8:662–7.
Beard CB, Pye G, Steurer FJ, Rodríguez R, Campman R, Peterson AT, et al. Chagas disease in a domestic transmission cycle in southern Texas, USA. Emerg Inf Dis. 2003;9:103–5.
Bustamante DM, Monroy MC, Rodas AG, Juarez JA, Malone JB. Environmental determinants of the distribution of Chagas disease vectors in south-eastern Guatemala. Geospat Health. 2007;2:199–211.
Carbajal AL, Porcasi X, Noireau F, Diotaiuti L, Gorla DE. The association between the geographic distribution of Triatoma pseudomaculata and Triatoma wygodzinskyi (Hemiptera: Reduviidae) with environmental variables recorded by remote sensors. Infect Genet Evol. 2009;9:54–61.
López-Cárdenas J, Gonzalez FE, Salazar MP, Gallaga JC, Ramírez E, Martinez J, et al. Fine-scale predictions of distributions of Chagas disease vectors in the State of Guanajuato. Mexico. J Med Entomol. 2010;42:1068–81.
Ibarra-Cerdeña CN, Zaldívar-Riverón A, Peterson AT, Sánchez-Cordero V, Ramsey JM. Phylogeny and niche conservatism in North and Central American triatomine bugs (Hemiptera: Reduviidae: Triatominae), vectors of Chagas’ disease. PLoS Negl Trop Dis. 2014;8:e3266.
Ceccarelli S, Balsalobre A, Susevich ML, Echeverria MG, Gorla DE, Marti GA. Modelling the potential geographic distribution of triatomines infected by Triatoma virus in the southern cone of South America. Parasit Vectors. 2015;8:153.
Gómez-Palacio A, Arboleda S, Dumonteil E, Peterson AT. Ecological niche and geographic distribution of the Chagas disease vector, Triatoma dimidiata (Reduviidae: Triatominae): evidence for niche differentiation among cryptic species. Infect Genet Evol. 2015;36:15–22.
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.
Costa J, Dornak LL, Almeida CE, Peterson AT. Distributional potential of the Triatoma brasiliensis species complex at present and under scenarios of future climate conditions. Parasit Vectors. 2014;7:238.
Escobar LE, Craft ME. Advances and limitations of disease biogeography using ecological niche modeling. Front Microbiol. 2016;7:1174.
Myers MF, Rogers DJ, Cox J, Flahault A, Hay SI. Forecasting disease risk for increased epidemic preparedness in public health. Adv Parasitol. 2000;47:309–30.
Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Poewell GVN, Underwood EC, et al. Terrestrial ecoregions of the world: a new map of life on earth. BioScience. 2001;51:933–8.
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–78.
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.
Botto-Mahan C, Cattan PE, Canals M, Acuña M. Seasonal variation in the home range and host availability of the blood-sucking insect Mepraia spinolai in the wild environment. Acta Trop. 2005;95:160–3.
Ehrenfeld MJ, Canals M, Cattan PE. Population parameters of Triatoma spinolai under different environmental conditions and densities. J Med Entomol. 1998;35:740–4.
Campos C, Torres-Pérez F, Botto-Mahan C, Coronado X, Solari A. High phylogeographic structure in the sylvatic vectors of Chagas disease of the genus Mepraia (Hemiptera: Reduviidae). Infect Genet Evol. 2013;19:280–6.
Brémond P, Salas R, Waleckx E, Buitrago R, Aliaga C, Barnabé C, et al. Variations in time and space of an Andean wild population of T. infestans at a microgeographic scale. Parasit Vectors. 2014;7:164.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2017. https://www.R-project.org/. Accessed 24 Jul 2017.
Taylor KE, Stouffer RJ, Meehl GA. A summary of the CMIP5 experiment design; 2009. http://pcmdi.llnl.gov/mips/cmip5/docs/Taylor_CMIP5_design.pdf. Accessed 10 Sept 2019.
IPCC. Climate change 2013: the physical science basis. In: Stocker TF, Qin D, Plattner GK, Tignor M, Allen SK, Boschung J, et al., editors. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. Cambridge: Cambridge University Press; 2013. p. 708–20.
Meinshausen M, Smith SJ, Calvin K, Daniel JS, Kainuma MLT, Lamarque J-K, et al. The RCP greenhouse gas concentrations and their extensions from 1765 to 2300. Clim Change. 2011;109:213–41.
Martin GM, Bellouin N, Collins WJ, Culverwell ID, Halloran PR, Hardiman SC, et al. The HadGEM2 family of Met Office Unified Model climate configurations. Geosci Model Dev. 2011;4:723–57.
Phillips SJ, Dudik M, Schapire RE. A maximum entropy approach to species distribution modeling. In: Proceedings of the twenty-first international conference on machine learning. New York: ACM Press; 2004. p. 472–86.
Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.
Merow C, Smith MJ, Silander JA. A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography. 2013;36:1058–69.
Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51.
Ortega-Huerta MA, Peterson AT. Modeling ecological niches and predicting geographic distributions: a test of six presence-only methods. Rev Mex Biodivers. 2008;79:205–16.
Hernandez PA, Graham CH, Master LL, Albert DL. The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography. 2006;29:773–85.
Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A, et al. Effects of sample size on the performance of species distribution models. Divers Distrib. 2008;14:763–73.
Qiao H, Soberón J, Peterson AT. No silver bullets in correlative ecological niche modelling: insights from testing among many potential algorithms for niche estimation. Meth Ecol Evol. 2015;6:1126–36.
Searcy CA, Shaffer HB. Do ecological niche models accurately identify climatic determinant of species ranges? Am Nat. 2016;187:423–35.
Anderson RP. A framework for using niche models to estimate impacts of climate change on species distributions. Ann N Y Acad Sci. 2013;1297:8–28.
Elith J, Kearney M, Phillips S. The art of modelling range-shifting species. Methods Ecol Evol. 2010;1:330–42.
Owens HL, Campbell LP, Dornak L, Saupe EE, Barve N, Soberón J, et al. Constraints on interpretation of ecological niche models by limited environmental ranges on calibration areas. Ecol Modell. 2013;263:10–8.
Osorio-Olvera L, Barve V, Barve N, Soberón J, Falconi M. ntbox: From getting biodiversity data to evaluating species distribution models in a friendly GUI environment. R package version 0.2.5.4; 2018. https://github.com/luismurao/ntbox. Accessed 18 Dec 2018.
Garza M, Feria TP, Casillas EA, Sanchez-Cordero V, Rivaldi C-L, Sarkar S. Projected future distributions of vectors of Trypanosoma cruzi in North America under climate change scenarios. PLoS Negl Trop Dis. 2014;8:e2818.
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–17.
Porto CE. Orogénesis soluciones geográficas. Porlamar, Venezuela; 2015. http://tapiquen-sig.jimdo.com. Accessed 2 Nov 2018.
Luz C, Fargues J. Factors affecting conidial production of Beauveria bassiana from fungus-killed cadavers of Rhodnius prolixus. J Invertebr Pathol. 1998;72:97–103.
Lima M, Previtali MA, Meserve PL. Climate and small rodent dynamics in semi-arid Chile: the role of lateral and vertical perturbations and intra-specific processes. Clim Res. 2006;30:125–32.
Chacón F, Bacigalupo A, Quiroga JF, Ferreira A, Cattan PE, Ramírez-Toloza G. Feeding profile of Mepraia spinolai, a sylvatic vector of Chagas disease in Chile. Acta Trop. 2016;162:171–3.
Ordenes H, Ehrenfeld M, Cattan PE, Canals M. Tripano-triatomine infection index of Triatoma spinolai in a zone with epidemiological risk for Chagas disease. Rev Med Chil. 1996;124:1053–7.
Frías-Lasserre D, González CR, Reyes C, Blanco de Carvalho D, Oliveira J, Canals M, et al. Wing polymorphism and Trypanosoma cruzi infection in wild, peridomestic and domestic collections of Mepraia spinolai (Hemiptera: Reduviidae) from Chile. J Med Entomol. 2017;54:1061–6.
Abad-Franch F, Ferraz G, Campos C, Palomeque FS, Grijalva MJ, Aguilar HM, et al. Modeling disease vector occurrence when detection is imperfect: infestation of Amazonian palm trees by triatomine bugs at three spatial scales. PLoS Negl Trop Dis. 2010;4:e620.
Crawford BA, Kribs-Zaleta CM. Vector migration and dispersal rate for sylvatic Trypanosoma cruzi transmission. Ecol Complex. 2013;14:145–56.
Ancapichún S, Garcés-Vargas J. Variability of the Southeast Pacific Subtropical Anticyclone and its impact on sea surface temperature off north-central Chile. Cienc Mar. 2015;41:1–20.
Campos-Soto R, Panzera F, Pita S, Lages C, Solari A, Botto-Mahan C. Experimental crosses between Mepraia gajardoi and M. spinolai and hybrid chromosome analyses reveal the occurrence of several isolation mechanisms. Infect Genet Evol. 2016;45:205–12.
Chávez J. Contribution to the study of Triatominae in Peru: geographic distribution, nomenclature and taxonomic notes. An Fac Med. 2006;67:65–76.
Garreaud RD. Cambio Climático: Bases físicas e impactos en Chile. Revista Tierra Adentro – INIA. 2011;93:1–14.
Deutsch CA, Tewksbury JJ, Huey RB, Sheldon KS, Ghalambor CK, Haak DC, et al. Impact of climate warming in terrestrial ectotherms across latitude. Proc Natl Acad Sci USA. 2008;105:6668–72.
Dobrowski SZ. A climatic basis for microrefugia: the influence of terrain on climate. Global Change Biol. 2011;17:1022–35.
Torres-Pérez F, Acuna-Retamar M, Cook JA, Bacigalupo A, García A, Cattan PE. Statistical phylogeographic of Chagas disease vector Triatoma infestans: testing biogeographic hypotheses of dispersal. Infect Genet Evol. 2011;11:167–74.
Bacigalupo A, Segura JA, García A, Hidalgo J, Galuppo S, Cattan PE. First finding of Chagas disease vectors associated with wild bushes in the Metropolitan Region of Chile. Rev Med Chil. 2006;134:1230–6.
Botto-Mahan C, Sepúlveda M, Vidal M, Acuña-Retamar M, Ortiz S, Solari A. Trypanosoma cruzi infection in the wild kissing bug Mepraia gajardoi from the Chilean Southern Pacific Ocean coast. Acta Trop. 2008;105:166–9.
Di Castri F, Hajek ER. Bioclimatología de Chile. Santiago: Editorial de la Universidad Católica de Chile; 1976.
Valladares F, Matesanz S, Guilhaumon F, Araújo MB, Balaguer L, Benito-Garzón M, et al. The effects of phenotypic plasticity and local adaptation on forecast of species range shifts under climate change. Ecol Lett. 2014;17:1351–64.
Nicotra AB, Atkin OK, Bonser SP, Davidson AM, Finnegan EJ, Mathesius U, et al. Plant phenotypic plasticity in a changing climate. Trends Plant Sci. 2010;15:684–92.
Chevin L-M, Lande R, Mace GM. Adaptation, plasticity, and extinction in a changing environmental: towards a predictive theory. PLoS Biol. 2010;8:e1000357.
We thank Mario Elgueta and Daniel Frías for allowing us to examine the insect collections from the Museo Nacional de Historia Natural (Santiago, Chile) and the Museo Entomológico de la Universidad Metropolitana de Ciencias de la Educación (UMCE, Chile), respectively. We also thank Francisca Peña, Emilio Paredes, Iván Wainnright and Nelson Adrian in charge of Zoonoses and Vector Sections (Regional Ministerial Secretary, SEREMI) of Región de Arica y Parinacota, Región de Tarapacá, Región de Valparaíso and Región de O’Higgins, respectively. We specially thank Patricio Arroyo and Universidad Católica del Maule-Chile for allowing us to use their computers, and Ricardo Campos for providing location data.
This work was supported by research grants CONICYT-FONDECYT 1170367 and 1140521. Partial support was obtained from CONICYT-FONDECYT 1180940. The funding sources had no role in the design of the study, collection, analysis, interpretation of data or in writing the 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.
Ecoregions present in the study area, which includes Chile and Peru.
Complete dataset for Mepraia spinolai occurrences. Dataset S2. Complete dataset for Mepraia gajardoi occurrences.
List of articles used to obtain occurrence coordinates for Mepraia spinolai and Mepraia gajardoi.
Flow diagram of the methods followed in this study.
Correlation matrix for bioclimatic variables used to model Mepraia spinolai potential distribution. Table S3. Correlation matrix for bioclimatic variables used to model Meparaia gajardoi potential distribution.
About this article
Cite this article
Garrido, R., Bacigalupo, A., Peña-Gómez, F. et al. Potential impact of climate change on the geographical distribution of two wild vectors of Chagas disease in Chile: Mepraia spinolai and Mepraia gajardoi. Parasites Vectors 12, 478 (2019). https://doi.org/10.1186/s13071-019-3744-9
- Wild vectors
- Trypanosoma cruzi
- Climate change
- Species distribution models
- Climate suitability
- Future scenarios
- Kissing bugs