
==== Front
Mem Inst Oswaldo Cruz
Mem Inst Oswaldo Cruz
mioc
Memórias do Instituto Oswaldo Cruz
0074-0276
1678-8060
Instituto Oswaldo Cruz, Ministério da Saúde

10.1590/0074-02760230100
01126
Research Article
Climate change-induced degradation of expert range maps drawn for kissing bugs (Hemiptera: Reduviidae) and long-standing current and future sampling gaps across the Americas
https://orcid.org/0000-0002-3589-9699
Shirey Vaughn processed the occurrence data performed the analysis interpreted the results wrote the manuscript together 1 2 +
Rabinovich Jorge came up with the original idea for analysis provided the range map data interpreted the results wrote the manuscript together 3
1 University of Southern California, Department of Biological Sciences, Marine and Environmental Biology Section, Los Angeles, CA, United States
2 Georgetown University, Department of Biology, Washington, DC, United States
3 Universidad Nacional de La Plata, Centro de Estudios Parasitológicos y de Vectores, La Plata, Argentina
+ Corresponding author: vshirey@usc.edu
The authors declare no conflicts of interest for this research

23 9 2024
2024
119 e23010016 6 2023
23 7 2024
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License
BACKGROUND

Kissing bugs are the vectors of Trypanosoma cruzi, the etiological agent of Chagas disease (CD). Despite their epidemiological relevance, kissing bug species are under sampled in terms of their diversity and it is unclear what biases exist in available kissing bug data. Under climate change, range maps for kissing bugs may become less accurate as species shift their ranges to track climatic tolerance.

OBJECTIVES

Quantify inventory completeness in available kissing bug data. Assess how well range maps are at conveying information about current distributions and potential future distributions subject to shift under climate change. Intersect forecasted changes in kissing bug distributions with contemporary sampling gaps to identify regions for future sampling of the group. Identify whether a phylogenetic signal is present in expert range knowledge as more closely related species may be similarly well or lesser understood.

METHODS

We used species distribution models (SDM), specifically constructed from Bayesian additive regression trees, with Bioclim variables, to forecast kissing bug distributions into 2100 and intersect these with current sampling gaps to identify priority regions for sampling. Expert range maps were assessed by the agreement between the expert map and SDM generated occurrence probability. We used classical hypothesis testing methods as well as tests of phylogenetic signal to meet our objectives.

FINDINGS

Expert range maps vary in their quality of depicting current kissing bug distributions. Most expert range maps decline in their ability to convey information about kissing bug occurrence over time, especially in under sampled areas. We found limited evidence for a phylogenetic signal in expert range map performance.

MAIN CONCLUSIONS

Expert range maps are not a perfect account of species distributions and may degrade in their ability to accurately convey distribution knowledge under future climates. We identify regions where future sampling of kissing bugs will be crucial for completing biodiversity inventories.

Key words:

climate change
bias
biodiversity knowledge
kissing bugs
range maps
Triatominae
VS# 1937959 VS is supported by a United States National Science Foundation Graduate Research Fellowship (# 1937959), the Department of Biology at Georgetown University; and a David H Smith Postdoctoral Conservation Research Fellowship while conducting this work
==== Body
pmcKissing bugs (Hemiptera: Reduviidae: Triatominae) are a largely hematophagous group of insects that are all potential vectors for Trypanosoma cruzi (Trypanosomatidae: Trypanosoma), a parasite that causes Chagas disease (CD), a potentially life-threatening condition that can cause long-term heart disease and damage to other organs. The parasite is typically passed into the human bloodstream after the insect completes a blood meal, and defecates on the hosts’ skin, passing the T. cruzi parasite from the feces into the host’s circulatory system. 1 Importantly, with increasing habitat fragmentation and climate change, it is projected that several kissing bug species may become epidemiologically more relevant, both in the endemic range of many species and in new areas that may become suitable to kissing bugs in the future under changing environmental conditions. 2 , 3 , 4 , 5 , 6

Being important disease vectors, kissing bug species have been previously examined through the lens of species distribution modeling (SDM) at both local and broader geographic scales. 7 , 8 , 9 , 10 , 11 Kissing bugs are largely Neotropical in their distribution with some species extending their range into the southern Nearctic. Their presence seems to be largely driven by climatic tolerances 5 , 6 , 12 - 17 as well as by the presence of suitable hosts from which to take a blood meal. 18 , 19 , 20 , 21 , 22 For example, several previous works have demonstrated that before the advent of SDMs, the ecological factors considered most important for influencing the distribution of kissing bugs were temperature and humidity. 23 Since then, it was recognized that some species endure more or less variation while others are harmed by the same fluctuations which was corroborated by SDMs analyses. 6 , 17 , 24 , 25 Temperature in particular is important for thermal preference 26 , 27 , 28 as well as host finding, feeding, reproduction, and development. 29 , 30 , 31 , 32

The early detection of the dispersal and/or range shifts of triatomine species to new areas (especially those populated by humans) is critical for assessing potential future public health threats. 3 , 33 , 34 , 35 , 36 Opportunistically sampled data such as those from museum collections can provide a historical basis for triatomine species ranges; additionally, community science is playing an increasing role in triatomine species early detection, 37 , 38 , 39 as has been the case for invasive insect species elsewhere. 40 Thus, it is important to establish a baseline assessment of where the knowledge of triatomine occurrence is most complete, and to examine which regions are under-sampled in order to target those areas in future sampling.

Given the importance of kissing bugs to public health, we sought to assess the current state of our knowledge of kissing bug occurrence from publicly available datasets. These datasets include DataTri 41 a curated database of triatomine occurrence records, and records from the Global Biodiversity Information Facility (GBIF), 42 which include “Research Grade” records from community science platforms such as iNaturalist as well as Integrated Digitized Biocollections (iDigBio) and Symbiota Collections of Arthropods Network (SCAN). 43 Previous work on inventory completeness and the assessment of bias in insect data has demonstrated that regions experiencing drastic climate change are likely to be under sampled 44 as sampling efforts typically focus on regions of high human population density or at the interface of anthropogenic infrastructure (e.g., roads and recreational trails). 44 , 45 , 46 Recent research has shown that these biases not only occur in triatomine sampling schemes, 47 but also increase in magnitude over time. 46

We expect kissing bugs to exhibit similar trends, mainly because ― being vectors of a human disease ― most of the sampling efforts are made at or near domiciliary rural sites 48 and because their wild habitats are difficult to sample. 49 Here we focus on how expertly drawn range maps compared to maps produced by species distribution models. In particular, we highlight how range maps may degrade in their ability to convey accurate information about species occurrence, especially in light of forecasted climate change. Finally, we performed an analysis of several kissing bug species, their probability of occurrence over time and gaps in species inventories to identify priority regions for the future sampling of this group.

MATERIALS AND METHODS

Species range maps and occurrence data - Species range maps were obtained from the triatomine geographical Atlas by Carcavallo et al. 50 These maps cover the breadth of potential kissing bug distributions throughout the neotropics and adjacent regions including the southern United States. Those maps were drawn by hand by the senior author (Dr Rodolfo Carcavallo) and there was no explanation for the methodology or procedures used in their drawing. Those range maps were scanned and digitized, recording the areas of presence at a scale of 0.1 x 0.1 coordinates degrees; so, our working data were a set of latitude/longitude coordinates representing those range maps, both graphically and as data in a spreadsheet. The author of the Atlas was an excellent “expert” in triatomines, collected tirelessly in all of the Americas personally, and was an inexhaustible collector of bibliography (his atlas is based not only in his personal experience but also on hundreds of locations cited by many other researchers). There were 115 triatomine species included in this original atlas; however, due to taxonomic revisions since the creation of the Atlas, there are now 112 valid species in the atlas. 51 , 52 , 53

Confirmed occurrence records for triatomine species were gathered from DataTri, 41 the Global Biodiversity Information Facility (GBIF), 42 the Integrated Digitized Biocollections (iDigBio), and the Symbiota Collections of Arthropods Network (SCAN) (see Supplementary data for a list of collections accessed). These occurrence records were then filtered to include only species for which we had range map information and in which the point occurrence intersected with the range map for that species for our inventory analysis. We retained all records for our SDM and expert score assessments. Additionally, we only included records of kissing bugs observed/collected from 1910-2021, further partitioning these into two distinct datasets for our analysis. All subsequent described analyses were conducted in R v. 4.2.1. 54

Inventory completeness analysis - Expected richness was calculated for 100×100 km square cells by layering expert range maps and counting the number of overlapping ranges per square cell: observed richness was calculated by counting the number of unique species from the filtered occurrence data for each grid cell. Thus, if expert ranges for three kissing bug species overlapped in a given cell, the expected richness was three. The ratio was calculated using the following equation:

Inventory completeness=(Observed richness from ocurrences) (Expected richness from range maps)

A ratio of observed richness to expected richness was calculated from these two values and fell between zero (no records of kissing bugs for which we had range maps) and one (complete recording of expected kissing bug richness based on range maps). This approach for inventory completeness has been performed elsewhere. 44

Species distribution models - We ran species distribution models using Bayesian Additive Regression Trees (BARTs) via the package “embarcadero”. 55 Like other machine-learning methods, BARTs compute a binary representation (using a logit-link function where applicable) of habitat suitability for a species (see Carlson 55 for more detail on the algorithm). In addition, BARTs have the added benefit of reduced overfitting problems that can be common in decision tree approaches. 56 We assessed the performance of our SDMs by examining the average area under the receiver-operator curve (AUC) across five top models per species. AUC indicates the overall performance of the model with respect to true and false positive rates. 57 Values above 0.5 indicate the model is performing better than random chance with respect to predictive capacity. The top models for each iteration of the BART procedure were selected using the highest training true skill statistic (TSS). 58 The AUC and TSS scores for all discovered models can be found in the Supplementary data.

We required species to have at least 25 unique occurrence records (from the period prior to 1999) to be considered for distribution modeling. Unfortunately, this eliminated Triatoma brasiliensis from our pre-1999 analysis; however, given its epidemiological importance, we included it in the post-1999 analysis (where it meets our minimum requirement). Random background, or pseudoabsence, points were generated such that the number of pseudoabsence points was equal to the number of occurrence records used for modeling. Five identically specified models were run per species to obtain an average probability of occurrence. 59 We used the Bioclim dataset for environmental predictors in our model. 60 In the past, kissing bug distributions have been modeled using the full Bioclim dataset, with the variables BIO4 (Temperature Seasonality), BIO5 (Maximum Temperature of the Warmest Month), BIO6 (Minimum Temperature of the Coldest Month), BIO13 (Precipitation of the Wettest Month), BIO14 (Precipitation of the Driest Month), and BIO15 (Precipitation Seasonality) found to be highly predictive of kissing bug occurrence in the past. 6 We used correlation plots to identify similar variables for inclusion in our present study to reduce issues with multicollinearity. In the case of purely predictive approaches, issues such as multicollinearity are not as relevant; 61 , 62 however, we have provided correlation plots highlighting the correlation between BIOs for each species in our analysis in the Supplementary data. Studies that aim to infer the specific, causal influence of environmental factors on distributions should more explicitly consider multicollinearity among other signals of potential confounders. Finally, we partitioned our occurrence dataset into data that were collected before 1999 (to remain congruent with when the expert range maps were drawn); and data from across all periods (to examine how well models conditioned on data collected both pre- and post-1999 agreed with expert range maps). We used a 2.5-minute resolution of the Bioclim dataset across all analyses. Full specifications of the models as well as AUC diagnostics can be found in Supplementary data (Table I) (pre-1999 model diagnostics) and Supplementary data (Table II) (all occurrence model diagnostics). We also include all model files as an additional Supplementary data to this work. We used a minimum bounding box which fully encompassed both the occurrence records and expert range maps with a 500-kilometer buffer zone as the calibration area for the model. We used this minimum bounding box to constrain the analysis and reduce the potential that more global areas of calibration might have on predictions, 63 , 64 including the buffer zone to allow for projections to potentially suitable environmental spaces that kissing bugs could plausibly occur but have not been explicitly sampled.

Expert score analysis - Following the construction of our SDM, we used the package, “expertscore” to compute our metrics for expert range maps. 65 Expert scores close to one indicate a high congruence between expert range maps and species distribution models while values close to zero indicate that the expert map is no more predictive of species distribution than a null model map. Expert scores less than zero indicate that the null model map has greater predictive accuracy than the expert range map. In our case, the null model map was the minimum quadrilateral polygon that encompasses both the expert range map and the point occurrence data with a buffer of 500 km. Following the calculation of expert scores, we compared performance of the pre-1999 and all occurrence record analysis to each other and from the current time period into 2100 using a two-sided Wilcox test with a two-sided alternative hypothesis. We conducted this test because our expert score data did not meet assumptions of normality of residuals when conducting a t-test. These comparisons were made across datasets to assess the influence that additional occurrences may have had on the stability of expert range maps constructed from occurrences pre-1999.

Using our distribution models, we forecasted each of the species’ ranges into the periods 2041-2060, 2061-2080, and 2081-2100 under RCP 8.5 (GCM ACCESS-ESM1-5, 2.5-minute resolution) to assess the rate at which expertly drawn range maps become less/more reliable over time according to their expert scores. Model transfer to future climatic conditions was performed using the “predict” function in the R package “raster” 66 with our BART models that were trained using current climatic conditions and occurrences. This allowed us to project species responses to the current environment onto forecasted environmental conditions. We also intersected the average differences in SDM predicted occurrence probabilities with our inventory completeness analysis to assess how many areas forecasted to increase in overall mean kissing bug occurrence probability (across all species) overlapped with regions that currently have little to no publicly available kissing bug data (i.e., large sampling gaps). In doing this we generated a map which shows where kissing bug probability of occurrence is likely to increase intersected with where the current knowledge gaps (represented by inventory completeness) occur. We split these metrics into four equally sized partitions representing the shift in kissing bug occurrence probabilities (low to high) and the completeness of contemporary inventories (complete to incomplete). These categories were used to produce a bivariate map highlighting where future sampling might be prioritized.

Finally, we aimed to test if expert scores and degradation/improvement of these scores over time exhibited a phylogenetic signal across a tree of kissing bugs obtained from Ceccarelli et al. 67 Several of the species in our analysis were not mapped onto an existing phylogenetic tree of kissing bug species and so we pruned those tips from our tree and subsequent analyses. This resulted in 27 species for the pre-1999 scores and 28 species for the scores using all available occurrence data. We tested for phylogenetic signal among expert scores and the magnitude/direction of expert score shifts across our tree using Blomberg’s K, K-star, 68 and Pagel’s λ. 69 Blomberg’s K and K-star values closer to or greater than 1.0 are indicative variance in the trait value being distributed among clades while values less than 1.0 are indicative of the variance being distributed within clades (i.e., closely related species do not resemble one another). 68 Pagel’s λ values close to 1.0 indicate a high correlation between species traits equal to Brownian motion. 69 We used the package “phylosignal” to calculate these metrics. 70

RESULTS

Records for kissing bugs are scarce, especially in remote and sparsely populated regions, particularly in the Amazon River Basin and northern Mexico (Figs 1-2). With respect to inventory completeness, several regions of low completeness emerged including southern Brazil, Uruguay, Peru, Ecuador, Colombia, Venezuela, Guyana, Suriname, and French Guiana. Notably, areas in close proximity to major cities such as Santiago, Chile; Córdoba, Argentina; Mexico City, Mexico; and San Diego/San Antonio, United States demonstrates a high density of occurrence records. In contrast, however, many regions have no single digitized occurrence record for a kissing bug species (e.g., much of the Amazon River basin). Overall, our SDMs performed well in their ability to predict current kissing bug distributions according to AUC metrics ENT#091;Supplementary data (Tables I, II)ENT#093;.

Fig. 1: point occurrence data and years of observation (inset line plot) colored by each record’s ‘basisOfRecord’ field (of either museum specimens or human observations). A random sampling of 50% of the occurrences (n = 13,641) are illustrated on the map to avoid overplotting. Green = observation; Orange = collected specimens.

Fig. 2: kissing bug inventory completeness at the 100×100 km spatial resolution (darker colors indicate greater inventory completeness). The inset map indicates expected richness based on overlapping expert range maps (darker colors indicate greater expected kissing bug species richness).

Expertly drawn range maps (constructed pre-1999) varied in their level of agreement with species distribution models constructed from occurrence data in the same time period (Figs 3A, B, C, 4, Supplementary data - Table III). When we included additional occurrence records collected post-1999, the agreement between expertly drawn range maps and SDM improved slightly overall ENT#091;Fig. 4, Supplementary data (Table IV)ENT#093;. In total, we obtained enough information to model 32 (pre-1999) and 33 (all occurrence) species of kissing bugs. We used the full dataset of either 32 or 33 species for all subsequent analyses. Among all species modeled in both the pre-1999 and all occurrence record frameworks, Panstrongylus geniculatus expressed the least agreement (0.24 and 0.29) with its expertly drawn range map. P. megistus and T. guasayana expressed the best agreement in the pre-1999 analysis (0.91 and 0.89) and all occurrence analysis (0.95 and 0.92) respectively. The average contemporary expert map score was 0.68 ENT#091;standard deviation (SD) +/- 0.18ENT#093; for the pre-1999 models and 0.73 (SD +/- 0.17) for the models using all occurrence records (Fig. 4). The distribution of expert scores across species from both datasets were not significantly different according to the Wilcox test (W = -423, p-value = 0.171).

Fig. 3: agreement between expert range maps (drawn in 1999 and shown here as red polygons) and species distribution model (SDM) output produced from only occurrence records collected before 1999 (A ,B, C) and occurrence records collected in all time periods (D, E, F). Only three species are highlighted here. Species and corresponding scores, from left to right pairs of panels, are (A, D) Triatoma delpontei (0.53, 0.60); (B, E) Rhodnius pallescens (0.54, 0.58); and (C, F) Rhodnius prolixus (0.51, 0.52).

Fig. 4: degradation of expert ranges maps over time under climate change scenario RCP 8.5 into 2041-2060, 2061-2080, and 2081-2100 (Global Circulation Model: ACCESS-ESM1-5, 2.5-minute resolution) based on species distribution projections from (A) the only pre-1999 occurrence record model and (B) the full occurrence record model. Species are summarized into two groups, declining expert scores over time (red) and increasing export scores over time (blue). The shaded area represents on standard deviation of variation around the mean trend line.

Over time, most expertly drawn range maps degrade in their ability to convey accurate information about kissing bug species distributions under a “business-as-usual” carbon emissions scenario (RCP 8.5) (Fig. 4). On average, expert range maps scores from now through the 2100s shifted by -0.046 (SD +/- 0.11) and -0.054 (SD +/- 0.13) for the pre-1999 models and all occurrence models respectively (Fig. 4). There was no significant difference between the shift in expert score across the pre-1999 and all occurrence models (W = 553, p-value = 0.75). Expert range maps for T. lecticularia (-0.35) had the greatest decline in being able to convey information about kissing bug occurrence into the 2081-2100 time period from the pre-1999 model. When examining declines in the ability of expert range maps to convey accurate information about species using all available occurrence data, T. lecticularia (-0.36) and T. brasiliensis (-0.35) emerged as range maps with sharply declining accuracy into 2081-2100. Notably, several species range maps became more accurate over time including R. pallescens (+0.26) and T. rubrofasciata (+0.12) (pre-1999 models); and R. pallescens (+0.25) and E. cuspidatus (+0.16) (all occurrence models).

When examining the mean occurrence probability shift into 2081-2100 for all kissing bugs over our time period of inference, regions experiencing notable sample gaps today are also regions forecasted to increase in the occurrence probability of kissing bugs on average (Fig. 5). Notably, regions already experiencing sampling gaps with respect to inventory completeness such as the eastern Amazon River Basin, parts of Venezuela, central Mexico, the Argentina/Paraguay/Bolivia border, and the Mexico/United States border were identified as key regions where average kissing bug occurrence is forecasted to increase and notable contemporary sampling gaps occur.

Fig. 5: potential regions in which to better sample kissing bug occurrence across the Americas. Targeted regions of sampling (green axis) between the (A) pre-1999 occurrence and (B) all occurrence models mostly overlap. Please note that this figure does not include measures of disease risk, simply projected distributions, and existing sampling gaps. Grey-green scale indicates higher probability of future kissing bug occurrence in 2080-2100 (green indicates higher probability) while the purple-grey scale indicates inventory completeness (purple indicates higher inventory completeness).

Across the kissing bug tree, we found generally low support for expert scores/change in expert scores being related strongly to phylogeny. In our pre-1999 analysis, Blomberg’s K and K-star were estimated to be smaller than 1.0 (0.26 ENT#091;p = 0.50ENT#093; and 0.31 ENT#091;p = 0.49ENT#093; respectively). Pagel’s for this set of expert scores was estimated to be 0.31 (p = 0.15). The estimated shift in expert range map score also exhibited low phylogenetic signal for the pre-1999 analysis with Bloomberg’s K at 0.399 (p = 0.16), but K-star was estimated at 0.58 (p = 0.029) indicating stronger phylogenetic signal in the shift in expert score from now into 2100. Pagel’s λ for the shift in expert score metric in this analysis was estimated to be 6.2 x 10-5 (p = 1.0). In our analysis of expert score using all available occurrence information Blomberg’s K and K-star were estimated to be 0.31 (p = 0.33) and 0.33 (p = 0.49) for the contemporary score respectively. Pagel’s λ was estimated to be 0.35 (p = 0.03). With respect to the change in expert score from now into 2100, Blomberg’s K and K-star were estimated to be 0.36 (p = 0.27) and 0.54 (p = 0.09) respectively. Pagel’s λ was estimated to be 4.5 x 10-5 (p = 1.0). A full table including all scores can be found in Supplementary data (Table V). This pattern of detecting both low phylogenetic and slightly elevated phylogenetic signal across metrics may be, in part, due to the smaller size of the phylogeny we used in this analysis. 71

DISCUSSION

Despite their importance for public health, available information on kissing bug occurrence is shockingly sparse, mainly from regions that include much of the Amazon River Basin and parts of Mexico, Central and Northern South America (Figs 1-2). Even in regions with relatively high densities of occurrence records such as those in south-eastern Brazil and parts of Venezuela, inventory completeness for the group is low (Fig. 2), and this is even more surprising given that those areas have a high-density rural as well as metropolitan centers. In concordance with previous research, our expert range maps indicate high species richness for kissing bugs in Brazil 7 and across northern South America. 6 Specifically, Gurgel-Gonçalves et al. 7 in particular demonstrated high species richness along the Atlantic coast of Brazil which we also find from overlaying expert range maps (Fig. 2). We did not find high richness in the states of Ceará, Rio Grande do Norte, Paraíba, Pernambuco, and related regions, but this may be due to the range maps having been constructed from fewer data points pre-1999.

Our phylogenetic analysis of expert scores revealed low phylogenetic signal among contemporary scores (and, for that matter, scores across other projected time periods) ENT#091;Supplementary data (Table V)ENT#093;. We assumed that as all species belong taxonomically to one subfamily, the ability to accurately convey a range map may be similar for closely related species (e.g., closely related species may be similarly detectable or share a similar niche space). 72 , 73 , 74 Similarly, we expected that there may be a phylogenetic signal in the change in expert score over time given that closely related species may respond similarly to changing climatic conditions. This was not the case in our analysis as Blomberg’s K and Pagels λ estimates were quite low ENT#091;Supplementary data (Table V)ENT#093;. Testing expert scores and their forecasts shifts over time using a larger sample of kissing bug species and a tree with additional tip reflecting those species may help resolve this apparently conflicting evidence for phylogenetic signal in scoring.

Some of the same locations with notable sampling gaps (Fig. 2) correspond with where kissing bugs are expected to increase in their occurrence probability into 2081-2100 under a “business-as-usual” climate change scenario (Fig. 5). This may make these regions, especially the eastern Amazon River Basin, parts of Venezuela, central Mexico, the Argentina/Paraguay/Bolivia border, and the Mexico/United States border, important to sample in the future both to confirm the results of our species distribution forecasts but also to detect emerging and important vector species on the move. Baseline expert range map scores may also be impacted by misidentifications of historical material, as may be the case with the Rhodnius species group, where the expert range map indicates a wide distribution of R. proxilus (Fig. 3), but many specimens may have actually been R. neglectus based on a more recent morphometric analysis. 75

Several triatomine species have been recorded to have been expanding both in range and in habitat type: e.g., P. geniculatus, considered until recently a sylvatic triatomine 50 that fed almost exclusively on the nine-banded armadillo (Dasypus novemcinctus), is now being found in human habitations in Colombia, 11 , 76 Brazil, 77 , 78 and Venezuela. 79 However, to determine if the cause of such distribution/habitat changes can be attributed to climatic or non-climatic factors is extremely difficult, for usually several factors are concomitant in determining these geographic/habitat changes. Catalá 80 has shown that in T. infestans the most attractive habitats for dispersing bugs would be those at short distance, with high CO2 emission and strong IR radiation, indicative of host presence (within the domestic habitat goat corrals may be the most attractive habitat to disperse); additionally, dispersal would be favored in periods of low atmospheric water saturation when IR perception is highest. 80 Furthermore, Baines et al. 81 have suggested that phenotype-by-environment interactions strongly influence dispersal of populations. On the other hand, transient changes in dispersal are common in many species undergoing range expansion, and may have major population and biogeographic consequences. 82 In the case of triatomines, inter-specific competition seems to be one of those non-climatic factors; e.g., (a) T. rubrovaria in Brazil is found in the wild in rocky habitats, but between 1975 and 1997, a growing domiciliary and peridomiciliary invasion of T. rubrovaria has taken place since the control of T. infestans, 83 an effect determined by inter-specific competition, 84 (b) Abrahan et al. 85 found that, after insecticide spraying, the sylvatic triatomines T. guasayana, T. eratyrusiformis, T. garciabesi, and T. platensis, not targeted by insecticide spraying, were captured simultaneously within peri-domestic areas and showed higher house invasion pressure than T. infestans, and (c) in Venezuela, P. geniculatus has been recorded to be in a process of domiciliation as a result of the control of R. prolixus, suggesting a competition for resources; 79 maybe this range dispersal and habitat change of P. geniculatus explains why this species showed in our analysis the least agreement (0.23 and 0.28) with the expertly drawn range map. On the other hand, T. infestans in itself is an example of the complexity of factors affecting the geographic/habitat occupation of triatomines, mainly because it has expanded its range into rapidly developing cities of Latin America, 86 but also restricting its range due to the intensive use of insecticide campaigns to control this vector species. Reductions in the geographic area of these important vector species have been quantified in the literature after such vector control programs have been implemented. For example, Ribeiro Jr et al. 87 demonstrated a reduced area of occurrence for P. megistus and T. infestans but increases in T. sordida and T. pseudomaculata after control programs were implemented in Bahia, Brazil. 87 Further examples of the influence of control campaigns on kissing bug distributions are: (a) In 1997 Uruguay was the first country to receive the International Evaluation Commission certification for achieving the interruption of T. cruzi transmission, through the total elimination of T. infestans populations, 88 and a similar certification was given to Brazil in 1999; 89 these two cases of human interventions impinge on the disappearance of T. infestans from very large regions; (b) however, there is another human intervention that acts opposite to the point above: the passive dispersal of T. infestans carried by people (and their domestic animals); 90 further, Abrahan et al. 91 consider that passive dispersal is one of the most frequent ways of spreading for triatomines over large areas. Passive dispersal of Rhodnius spp. by birds has also been observed. 92 To make things more complex, Richer et al. 93 showed by means of the detection of restricted gene flow between close but distinct sylvatic sites, that wild T. infestans does not disperse by flying at high altitude (2,750 meters above sea level). Recently, genomic techniques have helped to better understand the dispersal and habitat adaptation in triatomines; e.g., Hernandez-Castro et al. 94 have shown that R. ecuadoriensis shares outlier loci consistent with local adaptation to the domestic setting, which mapped to genes involved with embryogenesis and saliva production; and in the case of T. infestans Panzera et al. 95 showed that the ribosomal patterns are associated with a particular geographic distribution, and that chromosomal markers allowed to detect the existence of a hybrid zone occupied by individuals derived from crosses between two chromosomal groups.

Importantly, we want to highlight the potential risks of inferring the direct transferability of our SDMs for kissing bugs through time. Non-stationarity is likely present in many biological populations. 96 In other words, the effect of a climatic factor such as temperature may not be consistent across a species’ range due to underlying biological variation. This problem may be exacerbated in species with fairly large distributions that likely experience and have locally adapted to a variety of climatic conditions across their range. 96 Additional uncertainty in our future projections may also come from the fact that we used a single RCP scenario, assuming a “business-as-usual” carbon emissions trajectory and one global circulation model. This scenario and its interaction with global circulation models may not entirely convey future conditions should human development or unforeseen climate effects take place, especially those farther out from the present. Future work in this sphere could examine the impact that these global models would have on inference regarding species ranges and expert range maps, especially for more fine-scale analyses than presented here. 97 Further, we did not include information in our distribution models about land-cover or human/livestock population densities, instead opting to focus on climate change, but these factors may also contribute strongly to the distribution of kissing bugs. 1 , 19 , 20 , 21 When considering these forecasted areas of sampling, it will be important to account for geographical spaces that kissing bugs could or could not occupy in relation to the environmental spaces predicted by our distribution models, for example in spaces that are hard to reach through dispersal. Finer scale study should explicitly incorporate this distinction. 98 Additionally, SDMs are likely to capture the realized rather than fundamental niche of a species 99 which may mean that our future projections are merely projections of the realized niche space which could change over time due to a myriad of factors. The temporal transfer of niche models to future conditions should be treated with caution as truly novel environments could indeed be suitable for species occurrence in addition to errors of environmental omission based on the extent of the calibration area. 100 , 101 Risks of over- and underpredicting are possible considering that projections to future scenarios may include falsely suitable regions and decisions made in modelling may lead to partial niche characterization. Finally, although closely related kissing bug species may have similar niches, they may not respond to changing climates in the same or similar way. Further study should aim to disentangle the roles of geography and evolution in driving responses to future climate change in the region. Thus, we encourage that our future projections are considered in the context of these omissions and modeling decisions.

Finally, we enthusiastically encourage the continued collection and monitoring of kissing bug distribution data which will undoubtable serve as a useful validation tool for this and future analyses. Early surveillance programs including those in Brazil 83 as well as the importance of community-science programs for the early detection of kissing bug occurrence are critical for educating and preventing potential outbreaks of CD. 102 , 103 Such programs will likely become more popular and better at identifying novel occurrences of kissing bugs, especially with the advent of artificial intelligence tools like computer vision, which are already being developed to catalogue and identify kissing bug species from cellphone and other photography. 103 , 104 , 105 , 106 , 107 , 108

Kissing bugs are an increasingly important group of vector insect species in the Americas. Here, we have demonstrated that our current understanding of their distributions via expert range maps is quite variable at the species level and that, on average, the ability of these expert range maps to convey accurate information about kissing bug distributions will likely decline under a “business-as-usual” carbon dioxide emissions scenario. Further, regions that are currently under-sampled for kissing bugs may also be regions that are increasing in climatic suitability for many species. With increasing human development and habitat fragmentation, interactions between kissing bugs and their hosts (both human and non-human) is quite complex: e.g., Ceballos et al. 109 showed that massive deforestation around villages or selective extraction of older trees in the Dry Chaco in Argentina, has led to reductions in opossum abundance jointly with increases in foxes and skunks, leading to a dramatic decrease of T. cruzi infection in wild reservoir hosts, but may be on the rise in urban and suburban habitats, despite the State and community triatomine control activities. As an important vector of CD, we should make a concerted effort to accurately and publicly document the occurrence of this important insect group well into the future.

ACKNOWLEDGEMENTS

To the countless people involved in curating, digitizing, and observing kissing bugs for their contributions to the data used in this work.

Financial support:VS is supported by a United States National Science Foundation Graduate Research Fellowship (# 1937959), the Department of Biology at Georgetown University; and a David H Smith Postdoctoral Conservation Research Fellowship while conducting this work.
==== Refs
REFERENCES

1 Soares AC Sant Anna MRV Gontijo NF Araújo RN Pessoa GCD Koerich LB Guarneri A Lorenzo M eds Triatominae - The biology of Chagas disease vectors Features of interaction between triatomines and vertebrates based on bug feeding parameters 2021 Switzerland Springer 239 264
2 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 Parasite Vectors 2014 7 238 238
3 Garza M Feria Arroyo TP CasillSupplas 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
4 Medone P Ceccarelli S Parham PE Figuera A Rabinovich J The impact of climate change on the geographic distribution of two vectors of Chagas disease implications for the force of infection Philos Trans R Soc Lond B Biol Sci 2015 370 1665 20130560 20130560 25688019
5 Garrido R Bacigalupo A Peña-Gómez F Bustamante RO Cattan PE Gorla DE Potential impact of climate change on the geographical distribution of two wild vectors of Chagas disease in Chile Mepraia spinolai and Mepraia gajardoi Parasite Vectors 2019 12 478 478
6 Eberhard FE Cunze S Kochmann J Klimpel S Modelling the climatic suitability of Chagas disease vectors on a global scale Elife 2020 9 e52072 32374263
7 Gurgel-Gonçalves R Galvao C Costa J Peterson AT Geographic distribution of Chagas disease vectors in Brazil based on ecological niche modeling J Trop Med 2012 2012 705326 705326 22523500
8 Parra-Henao G Suárez-Escudero LC González-Car S Potential distribution of Chagas disease vectors (Hemiptera, Reduviidae, Triatominae) in Colombia, based on ecological niche modeling J Trop Med 2016 2016 1439090 1439090 28115946
9 Moo-Llanes DA Montes de Oca-Aguilar AC.Rodríguez-Rojas JJ Pattern of climate connectivity and equivalent niche of Triatominae species of the Phyllosoma complex Med Vet Entomol 2020 34 4 440 451 32697402
10 Campos-Soto R Diaz-Campusano G Rives-Blanchard N Cianferoni F Torres-Perez F Biogeographic origin and phylogenetic relationships of Mepraia (Hemiptera, Reduviidae) on islands of northern Chile PLoS One 2020 15 6 e0234056 32525913
11 Zuluaga S Mejía P Velez-Mira A Quintero J Triana-Chavez O Cantillo-Barraza O Updated geographical distribution and natural infection of Panstrongylus geniculatus (Latreille, 1811) in Antioquia department, Colombia Parasite Epidemiol Control 2021 15 e00226 34646953
12 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 26321302
13 de Souza RCM Campolina-Silva GH Bezerra CM Diotaiuti L Gorla DE Does Triatoma brasiliensis occupy the same environmental niche space as Triatoma melanica Parasite Vectors 2015 8 1 1 14
14 de la Vega GJ.Medone P.Ceccarelli S.Rabinovich J.Schilman PE Geographical distribution, climatic variability and thermo-tolerance of Chagas disease vectors Ecography 2015 38 8 851 860
15 Gorla D Noireau F Geographic distribution of Triatominae vectors in America Elsevier 2017
16 Eduardo AA Santos LABO Reboucas MC Martinez PA Patterns of vector species richness and species composition as drivers of Chagas disease occurrence in Brazil Int J Environ Health Res 2018 28 6 590 598 30063379
17 Bender A Python A Lindsay SW Golding N Moyes CL Modelling geospatial distributions of the triatomine vectors of Trypanosoma cruzi in Latin America PLoS Negl Trop Dis 2019 14 8 0008411 0008411
18 Miles MA de Souza AA Povoa MM Chagas's disease in the Amazon basin III. Ecotopes of ten triatomine bug species (Hemiptera: Reduviidae) from the vicinity of Belém, Pará State, Brazi J Med Entomol 1981 18 4 266 278 6790704
19 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 2 303 307 11296839
20 Peterson AT Sánchez-Cordero V Beard CB Ramsey JM Ecologic niche modeling and potential reservoirs for Chagas disease, Mexico Emerg Infect Dis 2002 8 7 662 667 12095431
21 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 4 479 494 21739038
22 Testai R Ferreira de Siqueira M.Rocha DSB.Roque ALR.Jansen AM.Xavier SCDC Space-environment relationship in the identification of potential areas of expansion of Trypanosoma cruzi infection in Didelphis aurita in the Atlantic Rainforest PLoS One 2023 18 7 e0288595 37506103
23 De Lucena DT Ecología dos Triatomineos do Brasil Rev Bras Malariol Doenças Trop 1959 11 4 577 633 14467139
24 Diniz-Filho JAF Ceccarelli S Hasperué W Rabinovich J Geographical patterns of Triatominae (Heteroptera Reduviidae) richness and distribution in the Western Hemisphere Insect Conserv Divers 2013 6 6 704 714
25 De La Vega GJ.Schilman PE Ecological and physiological thermal niches to understand distribution of Chagas disease vectors in Latin America Med Vet Entomol 2018 32 1 1 13 28857300
26 Lazzari CR Temperature preference in Triatoma infestans (Hemiptera Reduviidae) Bull Entomol Res 1991 81 273 276
27 Schilman PE Lazzari CR Temperature preference in Rhodnius prolixus, effects and possible consequences Acta Trop 2004 90 115 122 14739030
28 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 17 9086704
29 Clark N The effect of temperature and humidity upon eggs of the bug, Rhodnius prolixus (Heteroptera, Reduviidae) J Anim Ecol 1953 1 82 87
30 Lazzari CR Nuñez JA The response to radiant heat and the estimation of the temperature of distant sources in Triatoma infestans J Insect Physiol 1989 35 525 529
31 Fresquet N Lazzari CR Daily variation of the response to heat in Rhodnius prolixus The roles of light and temperature as synchronisers J Insect Physiol 2014 70 36 40 25200474
32 Okasha AYK Effects of high temperature in Rhodnius prolixus (Stal) Nature 1964 204 1221 1222
33 Ceccarelli S Rabinovich JE Global climate change effects on Venezuela´s vulnerability to Chagas disease is linked to the geographic distribution of five Triatomine species J Med Entomol 2015 52 1333 1343 26336258
34 Badel-Mogollón J Rodríguez-Figueroa L Parra-Hena G Análisis espacio-temporal de las condiciones biofísicas y ecológicas de Triatoma dimidiata (Hemiptera: Reduviidae: Triatominae) en la región nororiental de los Andes de Colom Biomédica 2017
35 Ayala S Alvarado S Cáceres D Zulantay I Canals M Estimando el efecto del cambio climático sobre el riesgo de la enfermedad de Chagas en Chile por medio del número reproductiv Rev Med Chil 2019 147 6 683 692 31859820
36 Shi Y Wei Y Feng X Liu J Jiang Z Ou F Distribution, genetic characteristics and public health implications of Triatoma rubrofasciata, the vector of Chagas disease in Guangxi, China Parasit Vectors 2020 13 33 33 31959216
37 Cecere MC Rodríguez-Planes LI Vazquez-Prokopec G Kitron U Gurtler RE Community-based surveillance and control of Chagas disease vectors in remote rural areas of the Argentine Chaco a five-year follow-up Acta Trop 2019 191 108 115 30593817
38 Abrahan L Cavallo MJ Amelotti I Impact of involving the community in entomological surveillance of Triatoma infestans (Klug, 1834) (Hemiptera, Triatominae) vectorial control Parasit Vectors 2021 14 98 98 33546756
39 Pennington PM Rivera EP de Urioste-Stone SM Aguilar T Juárez JG Hendrichs J Pereira R Vreysen MJB Area-wide integrated pest management: development and field application A successful community-based pilot programme to control insect vectors of Chagas disease in rural Guatemala 2021 Boca Raton CRC Press 709 727
40 Larson ER Graham BM Achury R Coon JJ Daniels MK Gambrell DK From eDNA to citizen science emerging tools for the early detection of invasive species Front Ecol Environ 2020 18 4 194 202
41 Ceccarelli S Balsalobre A Vicente ME Curtis-Robles R Hamer SA Ayala Landa JM American triatomine species occurrences updates and novelties in the DataTri database GigaByte 2022 2022 gigabyte62 gigabyte62 10.46471/gigabyte.62 36824504
42 GBIF GBIF occurrence download 2022 10.15468/dl.3fjwjq
43 Heinrich PL Gilbert E Cobb NS Franz N Symbiota collections of arthropods network (SCAN): a data portal built to visualize, manipulate, and export species occurrences 2015 https://openknowledge.nau.edu/id/eprint/2258/
44 Shirey V Belitz MW Barve V Guralnick R A complete inventory of North American butterfly occurrence data narrowing data gaps, but increasing bias Ecography 2021 44 4 537 547
45 Girardello M Chapman A Dennis R Kaila L Borges PA Santangeli A Gaps in butterfly inventory data a global analysis Biol Conserv 2019 236 289 295
46 Bowler DE Callaghan CT Bhandari N Henle K Barth MB Koppitz C Temporal trends in the spatial bias of species occurrence records Ecography 2022 8 e06219
47 Valenca-Barbosa C Lima MM Sarquis O Bezerra CM Abad-Franch F Modeling disease vector occurrence when detection is imperfect II drivers of site-occupancy by synanthropic Triatoma brasiliensis in the Brazilian northeast PLoS Negl Trop Dis 2014 8 5 e2861 24811125
48 Ribeiro-Jr G Abad-Franch F de Sousa OMF dos Santos CGS Fonseca EOL dos Santos RF TriatoScore an entomological-risk score for Chagas disease vector control-surveillance Parasit Vectors 2021 14 492 492 34563255
49 Abad-Franch F Gurgel-Gonçalves R The ecology and natural history of wild Triatominae in the Americas Springer 2021
50 Carcavallo RU Curto de Casas SI Sherlock IA Galíndez Girón I Jurberg J Galvão C Geographical distribution and alti-latitudinal dispersion of Triatominae. In: Carcavallo RU, Galíndez Girón I, Jurberg J, Lent H, editors. Atlas of Chagas' disease vectors in the Americas Fiocruz 1999
51 Páez-Rondón O Otálora-Luna F Aldana E Revalidation of synonymy between Nesotriatoma flavida and N bruneri (Hemiptera, Reduviidae, Triatominae) J Arthropod Borne Dis 2017 11 4 446 452 29367921
52 Monteiro FA Weirauch C Felix M Lazoski C Abad-Franch F Evolution, systematics, and biogeography of the Triatominae, vectors of Chagas disease Adv Parasitol 2018 99 265 344 29530308
53 Oliveira Correia JPS Gil-Santana HR Dale C Galvão C Triatoma guazu is a junior synonym of Triatoma williami Insects 2022 13 591 591 35886767
54 R Core Team R: a language and environment for statistical computing 2022 Vienna R Foundation for Statistical Computing
55 Carlson CJ Embarcadero species distribution modeling with Bayesian additive regression trees in R Methods Ecol Evol 2020 11 7 850 858
56 Chipman HA George EI McCulloch RE BART Bayesian additive regression trees Ann Appl Stat 2010 4 1 266 298
57 Fielding AH Bell JF A review of methods for the assessment of prediction errors in conservation presence/absence models Environ Conserv 1997 24 1 38 49
58 Allouche O Tsoar A Kadmon R Assessing the accuracy of species distribution models prevalence, kappa and the true skill statistic (TSS) J Appl Ecol 2006 43 6 1223 1232
59 Barbet-Massin M Jiguet F Albert CH Thuiller W Selecting pseudo-absences for species distribution models How, where and how many? Methods Ecol Evol 2012 3 327 338
60 Fick SE Hijmans RJ WorldClim2 new 1km spatial resolution climate surfaces for global land areas Int J Climatol 2017 37 4302 4315
61 Wheelwright S Makridakis S Hyndman RJ Forecasting: methods and applications 1998 3 New York John Wiley & Sons
62 Shmueli G To explain or predict Stat Sci 2010 25 3 289 310
63 VanDerWal J Shoo LP Graham C Williams SE Selecting pseudo-absence data for presence-only distribution modeling how far should you stray from what you know? Ecol Model 2009 220 589 594
64 Barve N Barve V Jiménez-Valverde A Lira-Noriega A Maher SP Peterson AT The crucial role of accessible area in ecological niche modeling and species distribution Ecol Model 2011 222 1810 1819
65 Mainali KP Hefley TJ Ries L Fagan WF Matching expert range maps with species distribution model predictions Conserv Biol 2020 34 1292 1304 32115748
66 Hijmans RJ raster: geographic data analysis and modeling. R package version 2023 https://CRAN.R-project.org/package=raster
67 Ceccarelli S Justi SA Rabinovich JE Diniz JAF Filho Villalobos F Phylogenetic structure of geographical co-occurrence among New World Triatominae species, vectors of Chagas disease J Biogeogr 2020 47 6 1218 1231
68 Blomberg SP Garland T Jr Ives AR Testing for phylogenetic signal in comparative data behavioral traits are more labile Evolution 2003 57 4 717 745 12778543
69 Pagel M Inferring the historical patterns of biological evolution Nature 1999 401 6756 877 884 10553904
70 Keck F Rimet F Bouchez A Franc A phylosignal: an R package to measure.test.and explore the phylogenetic signal Ecol Evol 2016 6 9 2774 2780
71 Münkemüller T Lavergne S Bzeznik B Dray S Jombart T Schiffers K How to measure and test phylogenetic signal Methods Ecol Evol 2012 3 4 743 756
72 Weins J Graham CH Niche conservatism integrating evolution, ecology, and conservation biology Annu Rev Ecol Evol Syst 2005 36 1 519 539
73 Weins J Ackerly DD Allen AP Buckley LB Niche conservatism as an emerging principle in ecology and conservation biology Ecol Lett 2010 13 10 1310 1324 20649638
74 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 10 e3266 25356550
75 Gurgel-Gonçalves R Abad-Franch F Ferreira JBC Santana DB Cuba CAC Is Rhodinius proxilus (Triatominae) invading houses in central Brazil Acta Trop 2008 107 2 90 98 18550022
76 Reyes M Esteban L Torres FA Flórez M Agudelo JC Angulo VM Intrusión de Pastrongylus geniculatus y Rhodnius pallescens a viviendas y áreas sociales en un barrio de Bucaramanga, Santander, Colombia Biomédica 2011
77 Da Costa Valente V Potential for domestication of Panstrongylus geniculatus (Latreille, 1811) (Liemiptera, Reduviidae, Triatominae) in the municipality of Muaná, Marajó Island, State of Pará Mem Inst Oswaldo Cruz 1999
78 Ceretti-Junior W Vendrami DP de MO Matos-Junior Rimoldi-Ribeiro A Alvarez JV Marques S Occurrences of triatomines (Hemiptera: Reduviidae) and first reports of Panstrongylus geniculatus in urban environments in the city of São Paulo, Brazil Rev Inst Med Trop 2018
79 Feliciangeli MD Carrasco H Patterson JS Suarez B Martínez C Medina M Mixed domestic infestation by Rhodnius prolixus Stäl, 1859 and Panstrongylus geniculatus Latreille, 1811, vector incrimination, and seroprevalence for Trypanosoma cruzi among inhabitants in El Guamito, Lara State, Venezuela Am J Trop Med Hyg 2004 71 4 501 505 15516649
80 Catalá SS The infra-red (IR) landscape of Triatoma infestans An hypothesis about the role of IR radiation as a cue for Triatominae dispersal Infect Genet Evol 2011 11 1891 1898 21856443
81 Baines CB Ferzoco IMC McCauley SJ Phenotype-by-environment interactions influence dispersal J Anim Ecol 2020 88 8 1263 1274
82 Simmons AD Thomas CD Changes in dispersal during species' range expansions Am Nat 2004 164 3 378 395 15478092
83 Almeida CE Vinhaes MC Silveira AC de Almeide JR Silveira AC Costa J Monitoring the domiciliary and peridomiciliary invasion process of Triatoma rubrovaria in the State of Rio Grande do Sul, Brazil Mem Inst Oswaldo Cruz 2000 95 6 761 768 11080758
84 Berenger JM Pages F Les Triatominae une domestication qui se généralise (Triatominae: Growing Trend to Domesticity) Med Trop (Mars) 2007 67 217 222 17784669
85 Abrahan L Gorla D Catalá S Active dispersal of Triatoma infestans and other triatomines in the Argentinean arid Chaco before and after vector control interventions J Vector Ecol 2016 41 1 90 96 27232129
86 Khatchikian CE Foley EA Barbu CM Hwang J Ancca-Juárez J Borrini-Mayori K Population structure of the Chagas disease vector Triatoma infestans in an urban environment PLoS Negl Trop Dis 2015 9 2 e0003425 25646757
87 Ribeiro G Jr Araújo RF Carvalho CMM Cunha GM Lanza FC Miranda DLP Triatomine fauna in the state of Bahia, Brazil what changed after 40 years of the vector-control program? Rev Soc Bras Med Trop 2022 55 1 9
88 Salvatella R Rosa R La interrupción en Uruguay de la transmisión vectorial de Trypanosoma cruzi, agente de la enfermedad de Chagas, por control de Triatoma infestans Rev Patol Trop 2000 29 2 213 231
89 INCOSUR Evaluaciones nacionales en 2000 y 2001: Paraguay, Brasil, Argentina, Uruguay y Bolivia INCOSUR Chagas 2001
90 Cortez MR Monteiro FA Noireau F New insights on the spread of Triatoma infestans from Bolivia-implications for Chagas disease emergence in the Southern Cone Infect Genet Evol 2010 10 350 353 20060504
91 Abrahan LB Gorla DE Catalá SS Dispersal of Triatoma infestans and other Triatominae species in the arid Chaco of Argentina - Flying, walking or passive carriage The importance of walking females Mem Inst Oswaldo Cruz 2011 106 2 232 239 21537686
92 Gurgel-Gonçalves R Cuba CAC Predicting the potential geographical distribution of Rhodnius neglectus (Hemiptera, Reduviidae) based on ecological niche modeling J Med Entomol 2009 46 4 952 960 19645302
93 Richer W Kengne P Rojas Cortez MR Perrineau MM Cohuet A Fontenille D Active dispersal by wild Triatoma infestans in the Bolivian Andes Trop Med Int Health 2007 12 6 759 764 17550473
94 Hernandez-Castro LE Villacis AG Jacobs A Cheaib B Day CC Ocaña-Mayorga S Population genomics and geographic dispersal in Chagas disease vectors Landscape drivers and evidence of possible adaptation to the domestic setting PLoS Genet 2022 18 2 e1010019 35120121
95 Panzera F Ferreiro MJ Pita S Calleros L Pérez R Basmadjián Y Evolutionary and dispersal history of Triatoma infestans, main vector of Chagas disease, by chromosomal markers Infect Genet Evol 2014 27 105 113 25017654
96 Rousseau JS Betts MG Factors influencing transferability in species distribution models Ecography 2022 2022 7 e06060
97 Weins J Stralberg D Jongsomjit D Howell CA Snyder MA Niches, models, and climate change assessing the assumptions and uncertainties Proc Natl Acad Sci USA 2009 106 2 19729 19736 19822750
98 Brown JL Carnaval AC A tale of two niches methods, concepts, and evolution Front Biogeogr 2019 11 4 e44158
99 Veloz DS Williams JW Blois JS He F Otto-Bliesner B Liu Z No-analog climates and shifting realized niches during the late quaternary implications for 21st-century predictions by species distribution models Glob Chang Biol 2011 18 5 1698 1713
100 Owens HL Campbell LP Dornak L Saupe EE Barve N Soberónet J al Constraints on interpretation of ecological niche models by limited environmental range on calibration areas Ecol Model 2013 263 10 18
101 Feng X Park SD Liang Y Pandey R Pape SM Collinearity in ecological niche modeling confusions and challenges Ecol Evol 2019 9 18 10365 10376 31624555
102 Curtis-Robles R Wozniak EJ Auckland LD Hamer GL Hamer SA Combining public health education and disease ecology research using citizen science to assess Chagas disease entomological risk in Texas PLoS Negl Trop Dis 2015 9 12 e0004235 26658425
103 Delgado-Noguera LA Hernández-Pereira CE Ramírez JD Hernández C Velasquez-Ortíz N Clavijo J Tele-entomology and tele-parasitology a citizen science-based approach for surveillance and control of Chagas disease in Venezuela Parasite Epidemiol Control 2022 19 e00273 36118050
104 Khalighifar A Komp E Ramsey JM Gurgel-Gonçalves R Peterson AT Deep learning algorithms improve automated identification of Chagas disease vectors J Med Entomol 2019 56 5 1401 1410
105 Abdelghani BA Banitaan S Maleki M Mazen A Kissing bugs identification using convolutional neural network IEEE Access 2021 9 140539 140548
106 Cochero J Pattori L Balsalobre A Ceccarelli S Marti G A convolutional neural network to recognize Chagas disease vectors using mobile phone images Ecol Inform 2022 68 101587 101587
107 de Miranda VL de Souza EP Bambil D Khalighifar A Peterson AT de Oliveira Nascimento FA Cellphone picture-based, genus-level automated identification of Chagas disease vectors effects of picture orientation on the performance of five machine-learning algorithms Ecol Inform 2024 79 102430 102430
108 Gurgel-Gonçalves R de Miranda VL Khalighifar A Peterson AT Shooting in the dark automatic identification of disease vectors without taxonomic expert supervision Ecol Inform 2023 75 102029 102029
109 Ceballos LA Cardinal MV Vazquez-Prokopec GM Lauricella MA Orozco MM Cortinas R Long-term reduction of Trypanosoma cruzi infection in sylvatic mammals following deforestation and sustained surveillance in northwestern Argentina Acta Trop 2006 98 3 286 296 16839513
