Does Triatoma brasiliensis occupy the same environmental niche space as Triatoma melanica?
- Rita de Cássia Moreira de Souza†1Email author,
- Gabriel H Campolina-Silva†1,
- Claudia Mendonça Bezerra2,
- Liléia Diotaiuti1 and
- David E Gorla3
© Souza et al. 2015
Received: 22 May 2015
Accepted: 30 June 2015
Published: 10 July 2015
Triatomines (Hemiptera, Reduviidae) are vectors of Trypanosoma cruzi, the causative agent of Chagas disease, one of the most important vector-borne diseases in Latin America. This study compares the environmental niche spaces of Triatoma brasiliensis and Triatoma melanica using ecological niche modelling and reports findings on DNA barcoding and wing geometric morphometrics as tools for the identification of these species.
We compared the geographic distribution of the species using generalized linear models fitted to elevation and current data on land surface temperature, vegetation cover and rainfall recorded by earth observation satellites for northeastern Brazil. Additionally, we evaluated nucleotide sequence data from the barcode region of the mitochondrial cytochrome c oxidase subunit I (CO1) and wing geometric morphometrics as taxonomic identification tools for T. brasiliensis and T. melanica.
The ecological niche models show that the environmental spaces currently occupied by T. brasiliensis and T. melanica are similar although not equivalent, and associated with the caatinga ecosystem. The CO1 sequence analyses based on pair wise genetic distance matrix calculated using Kimura 2-Parameter (K2P) evolutionary model, clearly separate the two species, supporting the barcoding gap. Wing size and shape analyses based on seven landmarks of 72 field specimens confirmed consistent differences between T. brasiliensis and T. melanica.
Our results suggest that the separation of the two species should be attributed to a factor that does not include the current environmental conditions. However, as the caatinga is a biome that has existed in the area for at least the last 18,000 years, past conditions might have had an influence in the speciation process. The DNA Barcoding approach may be extended to these members of the subfamily Triatominae.
KeywordsTriatominae Geometric Morphometrics DNA barcode Chagas disease Geographic distribution Ecological niche modelling
The epidemiological importance of the Triatominae, a subfamily of blood-sucking bugs, derives from its capacity to transmit Trypanosoma cruzi, the etiological agent of Chagas disease, one of the most important vector borne diseases affecting around 8 million people in Latin America . At present, 148 triatomine species are recognized, whose taxonomic knowledge reflects a balance between the use of traditional approaches and a wide variety of evolutionary interpretations [2–8].
Triatoma brasiliensis sensu lato Neiva, 1911 is one of the most important species for Chagas disease transmission and the main vector of T. cruzi in semi-arid areas of Northeastern Brazil, colonizing both natural and artificial ecotopes [9, 10]. This species shows remarkable chromatic variation across its geographic distribution, and indeed, it was concluded that T. brasiliensis is a complex of species.
Multidisciplinary studies have indicated that the Triatoma brasiliensis complex is a monophyletic group [11, 12], comprising four species, one of which includes two distinct subspecies. The relevant taxa in the group are T. brasiliensis brasiliensis Neiva, 1911; T. brasiliensis macromelasoma Galvão, 1956; Triatoma juazeirensis Costa & Felix, 2007; Triatoma melanica Neiva & Lent, 1941 and Triatoma sherlocki Papa et al. 2002. Each member of the group can be identified by external morphological characteristics and a taxonomic key was recently published .
Correct taxonomic identification of these insects is important because triatomine bugs behave differently in terms of their adaptation to artificial environments. Some species are primarily found in sylvatic ecotopes, where they live associated with small nest-building mammals or birds. Others invade and can even colonize the peridomestic environment and/or the human dwellings, exposing domestic mammals and humans to the infection of T. cruzi [13–17]. Knowledge about the environmental space occupied by the triatomine species is an important factor for the development and/or improvement of vector control and surveillance strategies .
Currently, complementary methods to the classical taxonomy have been used in the comparative study of Triatominae species, such as geometric morphometric and molecular techniques. Wing shape is increasingly utilized in species identification and characterization  and in distinguishing sibling species [20–22].
Wings are excellent structures for studying morphological variation because the intersections of the wing veins provide many well-defined landmarks suitable for morphometrics, and the metric properties of the wing provide precise quantitative information for the identification of species complexes [23, 24] and within-species variations . In triatomines, the study of morphological variation by geometric morphometrics has been helpful to solve taxonomic problems [20, 26, 27], to analyze the house reinfestation phenomenon [28, 29] and to distinguish cryptic species [23, 30, 31].
Variations that are not expressed phenotypically ascribes importance to molecular studies for the comprehension of systematic issues. In this sense, the DNA sequencing has been widely used for identification, determination of phylogenetic relationships and inferences. The DNA Barcoding method is a system of universal taxonomic identification, fast with free access, based on the variation of a part of the sequence of the mitochondrial gene cytochrome c oxidase subunit I (CO1) . This method remains in constant evaluation and, among other groups, the use of the CO1 sequences has shown promise for species-level identification in insects, including cryptic species of triatomines [33–36] and other insects [37–41].
After the confirmation of the species composition of the T. brasiliensis complex, we started exploration about the causes of the speciation process within the complex. Guided by the general question on speciation causes, this study aims at verifying the hypothesis that T. brasiliensis and T. melanica occupy different environmental niche spaces as a consequence of the speciation process within the T. brasiliensis complex. For the first time, we report findings on DNA barcoding as a tool for the identification of these species, and to complement, we also report the analysis of wing geometric morphometrics.
Triatoma brasiliensis and Triatoma melanica collection sites
Cachoeira do Júlio (CE)
João Costa (PI) 
Independencia (CE) 
Monteiro (PB) 
Mãe d’Água (PB) 
Lagoa Grande (PE) 
Serra Talhada (PE) 
São José (PB) 
Santa Cruz (PB) 
São Francisco (PB) 
Jaguaruana (CE) 
Morada Nova (CE) 
Caicó (RN) 
Santa Fé (CE)
Merejo do Angico (CE)
São Bento do Incra (CE)
Morada Nova do Pedro (CE)
São Cristóvão (CE)
Várzea do Touro (CE)
Jasmin do Aluísio (CE)
Açude Novo do Satiro (CE)
Benfica do Incra (CE)
Morada Nova do Tomaz (CE)
Cachoeira dos Pedrosas (CE)
Belo Horizonte do Alfredo (CE)
Monte Azul (MG)
Porteirinha (MG) 
Centroid size (CS), an isometric estimator derived from coordinate data, was computed as a measure of overall wing size. It is defined as the square root of the sum of the squared distances between the centre of the configuration of landmarks and each individual landmark . The centroid size values of different species and sexes were compared by multiple linear regression using the “lm function” of R. The presence of atypical points (outliers) was verified and a 95 % confidence interval for the estimated regression coefficients was constructed. In this analysis, males and females were examined separately to check for the presence of sexual dimorphism.
Generalized Procrustes Analysis (GPA) superimposition algorithm  was used to produce shape variables. Here, both non-uniform (‘Partial Warps’) and the uniform component were used as shape variables, which describe respectively the local and global variation as deviations from a consensus shape . In order to explore the wing shape variation among species, a multivariate discriminant analysis was performed on shape variables and the main results were shown by a factorial map of the first two discriminant factors. Mahalanobis distances were also computed from the same variables and used to illustrate morphological divergence among groups (species and sexes) by a dendrogram constructed according to the UPGMA algorithm. Statistical significance of these distances was assessed by non-parametric analyses based on permutation tests (1000 runs), and corrected by the Bonferroni method for P ≤ 0.05.
The residual relationship between shape and size variables was evaluated using a multivariate regression analysis of the first two shape discriminant factors against the centroid size values. This procedure allowed measurement of the contribution of size to shape variation (allometry), where statistical significance was assessed according to the non-parametric method based on permutation (1000 runs) .
We used the modules COO, TET, MOG, COV and PAD from the CLIC package (freely available at http://mome-clic.com/) to collect landmarks, to compute centroid size and shape variables, as well as to evaluate the residual relationship between shape and size variables, and to generate discriminant models and Mahalanobis distances. For comparing overall wing size among groups and to elaborate quantile-plots of CS values, the lm function of R was used. The UPGMA dendrogram based on Mahalanobis distances was constructed using the neighbor module from the PHYLIP package (by Joe Felsenstein, available at http://evolution.genetics.washington.edu/phylip.html).
Genomic DNA was extracted from individual specimen legs (stored at −70 °C) using the Wizard Genomic DNA Purification System (Promega, Madison, WI), re-suspended in 30 μl of DNA Rehydration Solution (10 mM Tris pH 8.0, 1 mM EDTA pH 8.0) and then stored at −80 °C until amplification by standard polymerase chain reaction (PCR). A fragment of 658 bp in the 5′ end of the CO1 gene was amplified using the primers LCO 1490 (5′-GGTCAACAAATCATAAAGATATTGG-3′) and HCO 2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA-3′) , following the thermal cycling conditions: an initial denaturation step at 94 °C for 3 min, 30 cycles at 94 °C for 1 min, 45 °C for 1 min and 72 °C for 1 min, with a final extension at 72 °C for 10 min. The 20 μl PCR reaction mixes included 10 μl of ultrapure water, 2 μl of 10x PCR buffer, 0,8 μl of each primer (10 μM), 2 μl of dNTP (2,5 mM), 2 μl of MgCl2 (50 mM), 0,4 μl of Taq polymerase (5 U/μl) and 2 μl of the DNA template (10 ng/μl). The PCR products were visualized on a 2 % agarose gel and selected for direct sequencing, where both DNA strands were sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Inc., Foster City, California, USA), following the manufacturer’s protocol on an ABI 3730XL Automated DNA Analyzer (Applied Biosystems) available at the Oswaldo Cruz Institute, FIOCRUZ (Rio de Janeiro, Brazil).
All sequences were aligned with the Clustal W and edited using the Geneious R6 software (Biomatters Ltd, Auckland, New Zealand®). The consensus sequences analysed in this study are deposited on GenBank under the accession numbers KJ580486 to KJ580495. Afterwards, we used MEGA v. 5.2  to carry the sequence analyses. Intraspecific and interspecific genetic distances (sequence divergence) were calculated using the Kimura 2-Parameter (K2P) nucleotide substitution model . In this approach, we also quantified the number of variables sites among sequences, as well as the nucleotide frequencies. A neighbor-joining (NJ) tree was constructed based on the K2P distances for representation of relationships among the analysed sequences, and the reliability of each node was estimated using the bootstrap method with 2,000 replicates. The tree was rooted using sequences of Triatoma sordida Stål, 1851 from Minas Gerais, Brazil (G. H. Campolina-Silva, unpublished data) as an outgroup.
All barcoded specimens have been deposited in the Coleção de Vetores da Doença de Chagas (COLVEC), Centro de Pesquisas René Rachou/FIOCRUZ.
Ecological niche modelling
Ecological niche modelling methods were used for the study of the geographic distribution of biological species in general and triatomines in particular [17, 58, 59], using environmental variables as proxies of species niche dimensions. Two main approaches were used as environmental data sources: the Worldclim database  and remotely sensed variables recorded by earth observation satellites.
Previously published studies on geographic distribution of the T. brasiliensis species complex [10, 12, 17] used occurrence coordinates of specimens mainly collected during the 1980s and 1990s, before the species identification within the T. brasiliensis complex. To avoid confusion of occurrence points of possibly misclassified specimens, the present study of the geographic distribution of T. brasiliensis and T. melanica was based on an estimation of the environmental space from occurrence points where each species was collected either by our study group or taken from recent bibliographic records (15 for T. melanica and 140 for T. brasiliensis). The number of occurrence points in this study is smaller than the set analysed in previous studies (for example Costa et al. 2002) , as those studies used several occurrence points that in our consideration were initially identified as T. brasiliensis, but afterwards were synonymized to other specific names within the group complex. We used a set of occurrence points that have no possibility of taxonomic misidentification as T. brasiliensis brasiliensis. We consider that this conservative selection of occurrence points should show the distribution area of the species without errors. No absence points were used, but rather a random selection of 1000 background points. No biotic or abiotic barriers are known or can be imagined for the dispersion of the studied species, so that the accessible area for them is the whole study area. Although flying and/or walking are the most frequent mechanisms for active (and short ranged) dispersion in triatomines, passive dispersion (especially through human migrations) is the main mechanism for long range dispersion, as a number of very early studies demonstrated and recently confirmed using genetic markers .
A great variety of methods exist to estimate species distribution models. The maximum entropy is one of the most frequently used, because of the simplicity of use of the Maxent software. However, based on the simplicity criteria , the possibility of testing explicit assumptions, and on the equivalence between the results obtained using Maxent and Poisson point process models , we used a generalized linear model with a binomial link through the glm function of the R package.
A major objective of studies examining niche overlap is to determine if two species occupy the same niche. Even if two or more niches are identical, there will be some differences in the data purely by chance. To rule out detecting two niches as different, that actually only differ due to sampling variation, niche comparisons between two species must be done statistically to determine whether the same probability distribution describes the niche of two species, or whether there is evidence of some difference . The estimated distribution models for the T. melanica and T. brasiliensis were compared using the Schoener index of niche overlap D, an index varying between 0–1 (no niche overlap and identical niche, respectively), and used to perform niche equivalency and similarity tests . We followed the framework proposed by Broennimann et al. (2012)  (based on Warren et al. 2008) , and specifically an ordination method (PCA-env) to estimate niche overlap, as the method disentangles the dependence of species occurrence from different climatic conditions (correcting for relative availability of environments) and from environmental data resolution. For the niche equivalency test all occurrences were pooled and randomly split into two datasets, maintaining the number of occurrences as in the original datasets, and the niche overlap statistic D was calculated. The process was repeated 100 times and a histogram of simulated values was constructed. If the observed value of D falls within the density of 95 % of the simulated values, the null hypothesis of niche equivalency cannot be rejected. For the niche similarity test we randomly shift the entire observed density of occurrences in one range and calculate the overlap of the simulated niche with the observed niche in the other range. The test of niche similarity was also based on 100 repetitions. If the observed overlap is greater than 95 % of the simulated values, the entity occupies environments in both of its ranges that are more similar to each other than expected by chance. All calculations were carried out using an R script produced by Broennimann et al. (2012) , and appropriately modified.
Environmental variables included in this study were day and night land surface temperature (LST) and vegetation index (NDVI) estimated from raster images produced by the MODIS sensor during the period 2009 – 2013 (downloaded as granules h12v09, h13v09, h13v10, h13v11, h14v09, h14v10, h14v11 from the http://e4ftl01.cr.usgs.gov/MOLT/ site); rainfall estimated by the Tropical Rainfall Measuring Mission (TRMM) during the period 2005 – 2013 (downloaded from the Goddard Earth Science site http://mirador.gsfc.nasa.gov/collections/TRMM_3B43__007.shtml) and elevation (produced by the Shuttle RT Mission and downloaded as a digital elevation model (DEM) from the CGIAR Consortium for Spatial Information at http://www.cgiar-csi.org/data/srtm-90m-digital-elevation-database-v4-1). Temporal resolution was 8, 16 and 30 days for LST, NDVI and rainfall, respectively. The spatial resolution at the data sources were 1 km for LST, 250 metres for NDVI, 0.5° (nominally 25 kms at the Equator) for rainfall and 90 metres for the DEM. All imagery was reprojected to a common geographic projection with WGS84 datum and resampled to 250 metres spatial resolution. As some points of triatomine collection were closer than 250 m apart, the number of occurrence points considered for T. melanica decreased to 8 and 132 for T. brasiliensis. A total of 200 images of LST (day and night) dates, 100 of NDVI dates and 108 monthly rainfall time series were processed calculating average and standard deviation of the time series for the period 2009–2013. The ecological niche model was thus based on 9 environmental predictors. A separate model for each species was calculated using the glm function of R version 3.0.2. Variables were selected to minimize the AIC (Akaike Information Criteria) using the stepAIC package. Colinearity among predictor variables was evaluated using the variance inflation factor (using vif function of the R car package) and a variable remained in the model if vif <10. Area under the receiver operation curve was calculated for model evaluation estimating partial AUC using the pROC package of R  that take into account the criticisms raised by Lobo et al. 2007 on the use of AUC . Cross-validation (through the cv.glm function) was used to measure the robustness of model estimation, using K = 100 for T. brasiliensis and K = 8 for T. melanica.
DNA Barcoding analysis
Summary of the variation along of the CO1 sequences analyzed
Pairwise intra and interspecific genetic distances, based on the Kimura 2-Parameter model
1. T.me (KJ580495)
2. T.me (KJ580494)
3. T.me (KJ580493)
4. T.me (KJ580492)
5. T.br (KJ580490)
6. T.br (KJ580486)
7. T.br (KJ580488)
8. T.br (KJ580491)
9. T.br (KJ580487)
10. T.br (KJ580489)
Generalised linear model fit
−1.96e + 02***
−6.527e + 02***
The model for T. melanica (n = 8) showed that elevation, annual averages of LSTd, NDVI and variability of accumulated rainfall were the significant variables (Table 4). The occurrence chance of this species is positively related with elevation and averages of LSTd and NDVI, and negatively related with the variability of accumulated rainfall. Model fit was high, with pAUC = 98 % (CI95 of 97.2 – 99.4) and low cross validation error = 0.0085.
Niche comparison using principal component analysis showed that the two first principal components explain 58 % of the variation included in the 9 environmental variables used for the study. Niche overlap was estimated as D = 0.291. This value is larger than expected by chance, thus leading to acceptance of niche similarity (P > 0.05), although leading to rejection of niche equivalency (P < 0.05).
Based on similarities in morphology, geographical distribution profile, epidemiological significance and phylogenetic relationships, a review  analyzing the systematics and evolution of the Triatominae subfamily, list the probable existence of 8 species complexes within the genus Triatoma Laporte, 1832. Among these, the T. brasiliensis complex has been continuously evaluated and is currently composed by T. melanica, T. juazeirensis, T. sherlocki and T. brasiliensis, the latter having two subspecies, T. brasiliensis brasiliensis and T. brasiliensis macromelasoma .
Among the five members of this complex, the subspecies T. b. brasiliensis has been considered as the main concern in terms of Chagas disease transmission, because this subspecies is the most widespread and shows the highest rate of domestic capture and of natural infection by T. cruzi [9, 10]. On the other hand, to date, T. melanica has been found almost exclusively in wild environments, and collected specimens have been shown to be restricted to municipalities in the northern region of Minas Gerais and the southern region of Bahia .
Analysis involving Triatominae mtDNA sequences are regularly used in order to elucidate phylogenetic relationships among closely related Triatominae species. In this respect, the most widely used markers include fragments of gene encoding proteins such as cytochrome b [11, 34, 36, 72, 73] and subunits of ribosomal RNA (12S and 16S) [33, 74–76]. The present study shows the effectiveness of DNA barcoding in identifying species of the T. brasiliensis complex and reinforces its contribution to classic taxonomy. The analyzed CO1 fragment was sufficient to discriminate T. brasiliensis and T. melanica, showing that these species are separated by a large genetic distance. The discrimination of these species can be well represented by the NJ tree with bootstrap values of 100 % for the two identified clades. This shows that our analysis of pair wise distance-based K2P within the T. brasiliensis complex supported a clear barcoding gap between the intra- and interspecific variation .
In a recent study, Justi et al. conclude that DNA barcoding is not applicable for the identification of Southern American Triatoma species, which may have diverged recently, because they observed at least one intraspecific genetic distance greater than interspecific distance in assessed species . The authors mention that to be considered appropriate to identify species, intraspecific distances must always be lower than interspecific ones. Unlike the findings of Justi et al., the results of this study show that the CO1 barcode region is a robust marker for differentiating this species complex, and it may constitute a valuable tool for both epidemiologic studies and for Chagas disease control, including subspecies differentiation, as observed by Vendrami et al. . Furthermore, this methodology can be especially useful for identifying immature stages, whose traditional taxonomic method is made impossible in smaller nymphs, especially considering the species complexes that exhibit wide morphological and chromatic variation.
The barcode sequences result found is congruent with the geometric morphometric analysis. Centroid size variation was as significant in the interspecific analysis as it was in the separation between males and females. Size variation related to gender, with females of both species larger than males, was expected with the presence of sexual dimorphism being a known characteristic in triatomines . The morphometric analyses of wing was able to discriminate the two members of the complex (DF1 = 94 %), allowing assignment without error of each single wing to its corresponding taxon, and producing a pattern of relationships in agreement within known genetic distances between populations, as observed by Costa and colleagues . As a result of there being no allometric trends for wing shape comparisons between T. brasiliensis and T. melanica, it is suggested that the change in this structure’s shape should mainly reflect the variation of genetic origin . Furthermore, the recorded sexual dimorphism suggests stability and adaptation of both species to the sylvatic ecotope, without evidence of exposure to environmental stress conditions, which may lead to reduction in size of females, as proposed by Dujardin et al. 1999 .
The ecological niche model using elevation and current data on land surface temperature, vegetation cover and rainfall showed that elevation, the average of LSTd, accumulated rainfall and variability of LSTd and LSTn are the environmental variables that best describe the geographic distribution of T. brasiliensis, whereas elevation, the average of LST, NDVI and variability of accumulated rainfall are the ones best describing the distribution of T. melanica. The result for T. melanica is probably less reliable, because of the small sample size of occurrence points, and restricted to a small area. Conversely, the prediction map for T. brasiliensis is very reliable, and indicates a potential distribution particularly associated with the caatinga biome. Of particular importance for the objective of this study is the niche comparison between the species. Results show that niches are similar but not equivalent, suggesting that T. brasiliensis and T. melanica share similar environmental constraints but that a different set of variables within the environmental niche space restrict their distribution. It is interesting that the model for T. brasiliensis predicts occurrence where T. melanica was collected, and the model of T. melanica predicts occurrence (although loosely) where T. brasiliensis was collected. Both species are associated with the caatinga region and have rupicolous habits. No difference between the species at the microscale is expected within the macro landscape of the caatinga, but this is an aspect that should be verified under field conditions. In evolutionary terms, the results of this study suggest that the separation in two different species should be attributed to a factor that does not include the current environmental conditions. However, as the caatinga is a biome that has existed in the area for at least the last 18,000 years [80, 81], past conditions might have had an influence in the speciation process.
The potential geographic distribution maps obtained differed with previously published studies [17, 82]. The differences are attributed to a) our use of a different set of occurrence points of specimens, based on criteria of minimum chance of species misidentification, a weakness that old records, previous to the species identification within the T. brasiliensis complex might possess, b) our use of current environmental conditions of temperature, vegetation and rainfall (measured by sensors onboard earth observation satellites, reflecting the land cover situation within the last 5 years of an area that shows an important landscape modification within the last 20 years), instead of the Bioclim database  that reflect past conditions (averaged between 1950–2000), through an interpolation process of ground meteorological stations, that in the area were extremely sparse during the considered period and might produce low reliability for the temperature and rainfall variables for the area, and c) our approach used for the whole study area as potentially accessible by the species, without imposing political boundary limits to map predictions.
Despite the successful vector control programme of Triatoma infestans in most parts of the southern cone countries of South America, native vectors persistently reinfest insecticide-treated households, and sylvatic triatomines maintain disease transmission without colonizing human dwellings. In these scenarios, fine-scale vector studies are essential to define epidemiological risk patterns and clarify the involvement of the little-known triatomine taxa in disease transmission. These eco-epidemiological investigations, as well as the planning and monitoring of control interventions, rely on accurate taxonomic judgments. The problems of cryptic speciation and phenotypic plasticity illustrate this need.
Although DNA barcoding is a straightforward approach, it is applicable for identifying species of the T. brasiliensis complex and contributes to classic taxonomy. The analyses of wing shape was able to differentiate the two members of the complex, and showed that the species are not under stressful conditions. As a result of there being no allometric trends in the wing shape comparison between T. brasiliensis and T. melanica, it is suggested that the change in this structure’s shape should mainly reflect the variation of genetic origin. In evolutionary terms, the result of the environmental niche space occupied by the studied species suggests that their separation as two different entities should be attributed to a factor that does not include the current environmental conditions. However, as the caatinga is a biome that has existed in the area for at least the last 18,000 years, past conditions might have had an influence in the speciation process.
This study makes up part of the results reached by the Brazilian Network for Molecular Identification of Biodiversity (Br-Bol), financed by the National Council for Scientific and Technological Development (CNPq) through public notice 50/2010 PPBio. We acknowledge the Ciências sem Fronteiras (CNPq) programme that supported the working periods of DEG at FIOCRUZ-MG. We would like to thank Dr. Karina Alessandra Morelli, who was responsible for the Molecular Identification Platform of Parasites and Vectors, for her help and clarification; Raíssa Nogueira Brito and Juliana Assis Geraldo for their support in the analyses; the PCDCh health workers for their help in capturing the triatomines.
- Rassi Jr A, Rassi A, Marin-Neto JA. Chagas disease. Lancet. 2010;375:1388–402.PubMedView ArticleGoogle Scholar
- Schofield CJ, Galvão C. Classification, evolution, and species groups within the Triatominae. Acta Trop. 2009;110:88–100.PubMedView ArticleGoogle Scholar
- Frias-Lasserre D. A new species and karyotype variation in the bordering distribution of Mepraia spinolai (Porter) and Mepraia gajardoi Frias et al. (Hemiptera: Reduviidae: Triatominae) in Chile and its parapatric model of speciation. Neotrop Entomol. 2010;39:572–83.PubMedView ArticleGoogle Scholar
- Rosa JA, Rocha CS, Gardim S, Pinto MC, Mendonça VJ, Ferreira Filho JCR, et al. Description of Rhodnius montenegrensis n. sp. (Hemiptera: Reduviidae: Triatominae) from the state of Rondônia, Brazil. Zootaxa. 2012;3478:62–76.Google Scholar
- Abad-Franch F, Pavan MG, Jaramillo ON, Palomeque FS, Dale C, Chaverra D, et al. Rhodnius barretti, a new species of Triatominae (Hemiptera: Reduviidae) from western Amazonia. Mem Inst Oswaldo Cruz. 2013;108:92–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Gonçalves TC, Teves-Neves SC, Santos-Mallet JR, Carbajal-de-la-Fuente AL, Lopes CM. Triatoma jatai sp. nov. in the state of Tocantins, Brazil (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2013;108:429–37.PubMed CentralPubMedView ArticleGoogle Scholar
- Jurberg J, Cunha V, Cailleaux S, Raigorodschi R, Lima MS, Rocha DS, et al. Triatoma pintodiasi sp. nov. do subcomplexo T. rubrovaria (Hemiptera: Reduviidae: Triatominae). Pan-Amaz Saude. 2013;4:43–56.View ArticleGoogle Scholar
- Poinar Jr GO. Panstrongylus hispaniolae sp. n. (Hemiptera: Reduviidae: Triatominae), a new fossil triatomine in Dominican amber, with evidence of gut flagellates. Palaeodiversity. 2013;6:1–8.Google Scholar
- Silveira AC, Vinhaes MC. Elimination of vector-borne transmission of Chagas disease. Mem Inst Oswaldo Cruz. 1999;94:405–11.PubMedView ArticleGoogle Scholar
- Costa J, Almeida CE, Dotson EM, Lins A, Vinhaes M, Silveira AC, et al. The epidemiologic importance of Triatoma brasiliensis as a Chagas disease vector in Brazil: a revision of domiciliary captures during 1993–1999. Mem Inst Oswaldo Cruz. 2003;98:443–9.PubMedView ArticleGoogle Scholar
- Monteiro FA, Donnelly MJ, Beard CB, Costa J. Nested clade and phylogeographic analyses of the Chagas disease vector Triatoma brasiliensis in Northeast Brazil. Mol Phylogenet Evol. 2004;32:46–56.PubMedView ArticleGoogle Scholar
- Costa J, Correia NC, Neiva VL, Goncalves TC, Felix M. Revalidation and redescription of Triatoma brasiliensis macromelasoma Galvão, 1956 and an identification key for the Triatoma brasiliensis complex (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2013;108:785–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Morel CM, Lazdins J. Chagas disease. Nat Rev Microbiol. 2003;1:14–5.PubMedView ArticleGoogle Scholar
- Schofield CJ, Jannin J, Salvatella R. The future of Chagas disease control. Trends Parasitol. 2006;22:583–8.PubMedView ArticleGoogle Scholar
- Lima MM, Sarquis O. Is Rhodnius nasutus (Hemiptera; Reduviidae) changing its habitat as a consequence of human activity? Parasitol Res. 2008;102:797–800.PubMedView ArticleGoogle Scholar
- Silveira AC, Dias JCP. O controle da transmissão vetorial. Rev Soc Bras Med Trop. 2011;44:52–63.PubMedView ArticleGoogle Scholar
- Gurgel-Gonçalves R, Galvão C, Costa J, Peterson AT. Geographic distribution of Chagas disease vectors in Brazil based on ecological niche modeling. J Trop Med. 2012:ID705326. doi: 10.1155/2012/705326. Epub 2012 Feb 27Google Scholar
- Abad-Franch F, Monteiro FA. Molecular research and the control of Chagas disease vectors. An Acad Bras Cienc. 2005;77:437–54.PubMedView ArticleGoogle Scholar
- Henry A, Thongsripong P, Fonseca-Gonzalez I, Jaramillo NO, Dujardin JP. Wing shape of dengue vectors from around the world. Inf Gen Evol. 2010;10:207–14.View ArticleGoogle Scholar
- Baylac M, Villemant C, Simbolotti G. Combining geometric morphometrics with pattern recognition for the investigation of species complexes. Biol J Linn Soc. 2003;80:89–98.View ArticleGoogle Scholar
- Matias A, de la Rive JX, Torrez M, Dujardin JP. Rhodnius robustus in Bolivia identified by its wings. Mem Inst Oswaldo Cruz. 2001;96:947–50.PubMedView ArticleGoogle Scholar
- Abad-Franch F, Monteiro FA, Jaramillo N, Gurgel-Gonçalves R, Dias FBS, Diotaiuti L. Ecology, evolution, and the long-term surveillance of vector-borne Chagas disease: A multi-scale appraisal of the tribe Rhodniini (Triatominae). Acta Trop. 2009;110:159–77.PubMedView ArticleGoogle Scholar
- Villegas J, Feliciangeli MD, Dujardin JP. Wing shape divergence between Rhodnius prolixus from Cojedes (Venezuela) and Rhodnius robustus from Merida (Venezuela). Infect Genet Evol. 2002;2:121–8.PubMedView ArticleGoogle Scholar
- Calle D, Quiñones ML, Erazo HF, Jaramillo N. Morphometric Discrimination of Females of Five Species of Anopheles of the Subgenus Nyssorhynchus from Southern and Northwest Colombia. Mem Inst Oswaldo Cruz. 2002;97:1191–5.View ArticleGoogle Scholar
- Dujardin JP, Le Pont F, Baylac M. Geographic versus interspecific differentiation of sandflies: a landmark data analysis. Bull Entomol Res. 2003;93:87–90.PubMedView ArticleGoogle Scholar
- Gumiel M, Catala S, Noireau F, De Arias AR, Garcia A, Dujardin JP. Wing geometry in Triatoma infestans (Klug) and T. melanosoma Martinez, Olmedo & Carcavallo (Hemiptera: Reduviidae). Syst Entomol. 2003;28:173–9.View ArticleGoogle Scholar
- Vargas E, Espitia C, Patino C, Pinto N, Aguilera G, Jaramillo C, et al. Genetic structure of Triatoma venosa (Hemiptera:Reduviidae): molecular and morphometric evidence. Mem Inst Oswaldo Cruz. 2006;101:39–45.PubMedView ArticleGoogle Scholar
- Dujardin JP, Beard CB, Ryckman R. The relevance of wing geometry in entomological surveillance of Triatominae, vectors of Chagas disease. Inf Gen Evol. 2007;7:161–7.View ArticleGoogle Scholar
- Feliciangeli MD, Sanchez-Martin M, Marrero R, Davies C, Dujardin JP. Morphometric evidence for a possible role of Rhodnius prolixus from palm trees in house re-infestation in the State of Barinas (Venezuela). Acta Trop. 2007;101:169–77.PubMedView ArticleGoogle Scholar
- Gurgel-Gonçalves R, Abad-Franch F, Ferreira JB, Santana DB, Cuba CA. Is Rhodnius prolixus (Triatominae) invading houses in central Brazil? Acta Trop. 2008;107:90–8.PubMedView ArticleGoogle Scholar
- Gurgel-Gonçalves R, Ferreira JB, Rosa AF, Bar ME, Galvão C. Geometric morphometrics and ecological niche modelling for delimitation of near-sibling triatomine species. Med Vet Entomol. 2011;25:84–93.PubMedView ArticleGoogle Scholar
- Hebert PD, Cywinska A, Ball SL, De Waard JR. Biological identifications through DNA barcodes. Proc Biol Sci. 2003;270:313–21.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia BA, Powell JR. Phylogeny of species of Triatoma (Hemiptera: Reduviidae) based on mitochondrial DNA sequences. J Med Entomol. 1998;35:232–8.PubMedView ArticleGoogle Scholar
- Pfeiler E, Bitler BG, Ramsey JM, Palacios-Cardiel C, Markow TA. Genetic variation, population structure, and phylogenetic relationships of Triatoma rubida and T. recurva (Hemiptera: Reduviidae: Triatominae) from the Sonoran Desert, insect vectors of the Chagas’ disease parasite Trypanosoma cruzi. Mol Phylogenet Evol. 2006;41:209–21.PubMedView ArticleGoogle Scholar
- Cortez MR, Emperaire L, Piccinali RV, Gurtler RE, Torrico F, Jansen AM, et al. Sylvatic Triatoma infestans (Reduviidae, Triatominae) in the Andean valleys of Bolivia. Acta Trop. 2007;102:47–54.PubMedView ArticleGoogle Scholar
- Gardim S, Almeida CE, Takiya DM, Oliveira J, Araújo RF, Cicarelli RMB, et al. Multiple mitochondrial genes of some sylvatic Brazilian Triatoma: Non-monophyly of the T. brasiliensis subcomplex and the need for a generic revision in the Triatomini. Inf Gen Evol. 2014;23:74–9.View ArticleGoogle Scholar
- Hajibabei M, Janzen DH, Burns JM, Hallwachs W, Hebert PD. DNA barcodes distinguish species of tropical Lepidoptera. Proc Natl Acad Sci U S A. 2006;103:968–71.View ArticleGoogle Scholar
- Cywinska A, Hunter FF, Hebert PD. Identifying Canadian mosquito species through DNA barcodes. Med Vet Entomol. 2006;20:413–24.PubMedView ArticleGoogle Scholar
- Park DS, Foottit R, Maw E, Hebert PD. Barcoding bugs: DNA-based identification of the true bugs (Insecta: Hemiptera: Heteroptera). PLoS One. 2006;6, e18749.View ArticleGoogle Scholar
- Scarpassa VM, Alencar RB. Molecular taxonomy of the two Leishmania vectors Lutzomyia umbratilis and Lutzomyia anduzei (Diptera: Psychodidae) from the Brazilian Amazon. Parasit Vectors. 2013;6:258.PubMed CentralPubMedView ArticleGoogle Scholar
- Rebijith KB, Asokan R, Kumar NK, Krishna V, Chaitanya BN.Google Scholar
- Costa J, Peterson AT, Dujardin JP. Morphological evidence suggests homoploid hybridization as a possible mode of speciation in the Triatominae (Hemiptera, Heteroptera, Reduviidae). Inf Gen Evol. 2009;9:263–70.View ArticleGoogle Scholar
- Batista VSP, Fernandes FA, Cordeiro-Estrela P, Sarquis O, Lima MM. Ecotope effect in Triatoma brasiliensis (Hemiptera: Reduviidae) suggests phenotypic plasticity rather than adaptation. Med Vet Entomol. 2013;27(3):247–54.PubMedView ArticleGoogle Scholar
- Sarquis O, Carvalho-Costa FA, Oliveira LS, Duarte R, D’Andrea PS, Oliveira TG, et al. Ecology of Triatoma brasiliensis in northeastern Brazil: seasonal distribution, feeding resources, and Trypanosoma cruzi infection in a sylvatic population. J Vec Ecol. 2010;35(2):385.View ArticleGoogle Scholar
- Velença-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.View ArticleGoogle Scholar
- Borges EC, Dujardin JP, Schofield CJ, Romanha AJ, Diotaiuti L. Dynamics between sylvatic, peridomestic and domestic populations of Triatoma brasiliensis (Hemiptera: Reduviidae) in Ceará State, Northeastern Brazil. Acta Tropica. 2005;93:119–26.PubMedView ArticleGoogle Scholar
- Sarquis O, Carvalho-Costa FA, Toma HK, Georg I, Burgoa MR, Lima MM. Eco-epidemiology of Chagas disease in northeastern Brazil: Triatoma brasiliensis, T. pseudomaculata and Rhodnius nasutus in the sylvatic, peridomestic and domestic environments. Parasitol Res. 2012;110:1481–5.PubMedView ArticleGoogle Scholar
- Câmara ACJ, Varela-Freire AA, Valadares HMS, Macedo AM, D’Ávila DA, et al. Genetic analyses of Trypanosoma cruzi isolates from naturally infected triatomines and humans in northeastern Brazil. Acta Tropica. 2010;115:205–11.PubMedView ArticleGoogle Scholar
- Vendrami DP, Ceretti-Junior W, Obara MT, Marrelli MT. Mitochondrial PCR-RFLP Assay to Distinguish Triatoma brasiliensis macromelasoma from Triatoma brasiliensis brasiliensis Subspecies (Hemiptera: Reduviidae). J Trop Med. 2013;ID 305198:7.Google Scholar
- Correia N, Almeida CE, Lima-Neiva V, Gumiel M, Dornak LL, Lima MM, et al. Cross-mating experiments detect reproductive compatibility between Triatoma sherlocki and other members of the Triatoma brasiliensis species complex. Acta Trop. 2013;128:162–7.PubMedView ArticleGoogle Scholar
- Bookstein FL. Morphometric Tools for Landmark Data: Geometry and Biology. Cambridge UK: Cambridge University Press; 1991.Google Scholar
- Rohlf FJ. Morphometric spaces, shape components and the effects of linear transformations. In Advances in Morphometrics. Plenum Publication. Edited by Marcus LF, Corti M, Loy A, Naylor GJP, Slice D. NATO ASI, Series A Life Sciences. N Y. 1996;117–129.Google Scholar
- Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Geometric Morphometrics for Biologists: A Primer. Edited by Elsevier. New York: Academic; 2004.Google Scholar
- Good P. Permutation Tests: A Practical Guide to Resampling Methods for Testing Hypotheses. Edited by Springer. Universidade de Michigan: Academic Press. Hardin JW: Springer series in Statistics, 2ª Ed; 2000.Google Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.PubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.PubMedView ArticleGoogle Scholar
- Gorla DE. Variables ambientales registradas por sensores remotos como indicadores de la distribucíon geográfica de Triatoma infestans (Heteroptera: Reduviidae). Ecol Austral. 2002;2:117–27.Google Scholar
- de La Fuente AL C, Porcasi X, Noireau F, Diotaiuti L, Gorla D, et al. The association between the geographic distribution of Triatoma pseudomaculata and Triatoma wygodzinskyi (Hemiptera: Reduviidae) with environmental variables recorded by remote sensors. Inf Gen Evol.Google Scholar
- 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.View ArticleGoogle Scholar
- Costa J, Peterson A, Beard C. Ecologic niche modeling and differentiation of populations of Triatoma brasiliensis Neiva, 1911, the most important Chagas’ disease vector in northeastern Brazil (Hemiptera, Reduviidae, Triatominae). Am J Trop Med Hyg. 2002;67:516–20.PubMedGoogle Scholar
- Piccinali RV, Canale DM, Sandoval AE, Cardinal MV, Jensen O, Kitron U, et al. Triatoma infestans bugs in Southern Patagonia, Argentina Emerging. Inf Diseases. 2010;16:887–9.View ArticleGoogle Scholar
- Merow C, Smith MJ, Edwards TC, Guisan A, McMahon SM, Normand S, et al. What do we gain from simplicity versus complexity in species distribution models? Ecography. 2014;37:1267–81.View ArticleGoogle Scholar
- Renner IW, Warton DI. Equivalence of MAXENT and Poisson point process models for species distribution modeling in ecology. Biometrics. 2013;69:274–81.PubMedView ArticleGoogle Scholar
- Geange SW, Pledger S, Burns KC, Shima JS. A unified analysis of niche overlap incorporating data of different types. Meth Ecol Evol. 2011;2:175–84.View ArticleGoogle Scholar
- Schoener TW. No synchronous spatial overlap of lizards in patchy habitats. Ecology. 1970;51:408–18.View ArticleGoogle Scholar
- Broennimann O, Fitzpatrick MC, Pearman PB, Petitpierre B, Pellissier L, Yoccoz NG, et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecol Biogeography. 2012;21:481–97.View ArticleGoogle Scholar
- Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution. 2008;62:2868–83.PubMedView ArticleGoogle Scholar
- Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.PubMed CentralPubMedView ArticleGoogle Scholar
- Lobo JM, Jiménez-Valverde A, Real R. AUC: a misleading measure of the performance of predictive distribution models. Glob Ecol Biogeo. 2008;17:145–51.View ArticleGoogle Scholar
- Hebert PD, Stoeckle MY, Zemlak TS, Francis CM. Identification of Birds through DNA Barcodes. PLoS Biol. 2004;2, e312.PubMed CentralPubMedView ArticleGoogle Scholar
- Monteiro FA, Perez R, Panzera F, Dujardin JP, Galvão C, Rocha DS, et al. Mitochondrial DNA Variation on Triatoma infestans populations and its implication on the specific status of T. melanosoma. Mem Inst Oswaldo Cruz. 1999;94:851–954.View ArticleGoogle Scholar
- Monteiro FA, Barret TV, Fitzpatrick S, Cordon-Rosales C, Feliciangeli D, Beard CB. Molecular phylogeography of the Amazonian Chagas disease vectors Rhodnius prolixus and R. robustus. Mol Ecol. 2003;12:997–1006.PubMedView ArticleGoogle Scholar
- Garcia BA, Moriyama EN, Powell JR. Mitochondrial DNA sequences of triatomines (Hemiptera: Reduviidae): phylogenetic relationships. J Med Entomol. 2011;38:675–83.View ArticleGoogle Scholar
- Hypša V, Tietz DF, Zrzavy J, Rego ROM, Galvão C, Jurberg J. Phylogeny and biogeography of Triatominae (Hemiptera: Reduviidae): molecular evidence of a New World origin of the Asiatic clade. Mol Phylogenet Evol. 2002;23:447–57.PubMedView ArticleGoogle Scholar
- Sainz AC, Mauro LV, Moriyama EN, Garcia BA. Phylogeny of triatomine vectors of Trypanosoma cruzi suggested by mitochondrial DNA sequences. Genetica. 2004;121:229–40.PubMedView ArticleGoogle Scholar
- Justi SA, Dale C, Galvão C. DNA barcoding does not separate South American Triatoma (Hemiptera: Reduviidae), Chagas Disease vectors. Parasit Vectors. 2014;7:519.PubMed CentralPubMedView ArticleGoogle Scholar
- Lent H, Wygodzinsky P. Revision of the Triatominae (Hemiptera, Reduviidae), and their significance as vectors of Chagas’ Disease. Bull Amer Mus Nat Hist. 1979;163:123–520.Google Scholar
- Dujardin JP, Steindel M, Chavez T, Machane M, Schofield CJ. Changes in the Sexual Dimorphism of Triatominae in the Transition from Natural to Artificial Habitats. Mem Inst Oswaldo Cruz. 1999;94(4):565–9.PubMedView ArticleGoogle Scholar
- Adams JM. Global land environments since the last interglacial. Oak Ridge National Laboratory, TN, USA. [http://www.esd.ornl.gov/ern/qen/nerc.html].
- Adams JM, Faure H: Review and Atlas of Palaeovegetation: Preliminary land ecosystem maps of the world of the since the Last Glacial Maximum. Oak Ridge National Laboratory, TN, USA. [http://www.esd.ornl.gov/ern/qen/adams1.html].
- 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. Parasites Vectors. 2014;7:238.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.