Chagas vectors Panstrongylus chinai (Del Ponte, 1929) and Panstrongylus howardi (Neiva, 1911): chromatic forms or true species?
Parasites & Vectors volume 13, Article number: 226 (2020)
Chagas disease is a parasitic infection transmitted by “kissing bugs” (Hemiptera: Reduviidae: Triatominae) that has a huge economic impact in Latin American countries. The vector species with the upmost epidemiological importance in Ecuador are Rhodnius ecuadoriensis (Lent & Leon, 1958) and Triatoma dimidiata (Latreille, 1811). However, other species such as Panstrongylus howardi (Neiva, 1911) and Panstrongylus chinai (Del Ponte, 1929) act as secondary vectors due to their growing adaptation to domestic structures and their ability to transmit the parasite to humans. The latter two taxa are distributed in two different regions, they are allopatric and differ mainly by their general color. Their relative morphological similarity led some authors to suspect that P. chinai is a melanic form of P. howardi.
The present study explored this question using different approaches: antennal phenotype; geometric morphometrics of heads, wings and eggs; cytogenetics; molecular genetics; experimental crosses; and ecological niche modeling.
The antennal morphology, geometric morphometrics of head and wing shape and cytogenetic analysis were unable to show distinct differences between the two taxa. However, geometric morphometrics of the eggs, molecular genetics, ecological niche modeling and experimental crosses including chromosomal analyses of the F1 hybrids, in addition to their coloration and current distribution support the hypothesis that P. chinai and P. howardi are separate species.
Based on the evidence provided here, P. howardi and P. chinai should not be synonymized. They represent two valid, closely related species.
Chagas disease constitutes a serious problem in Latin America. Today, six to seven million people worldwide, mostly in Latin America, are infected with Trypanosoma cruzi, the causative agent of this disease . The population most commonly at risk is typically of low socioeconomic status . The main way of transmission is through the contact with infected feces of insects of the subfamily Triatominae, so vector control is the most useful method to prevent Chagas disease in Latin America .
Currently, 16 triatomine species have been reported in Ecuador, with Rhodnius ecuadoriensis and Triatoma dimidiata being the main Chagas disease vectors [4, 5]. Six species of the genus Panstrongylus have also been reported in Ecuador: Panstrongylus chinai; P. howardi; P. geniculatus (Latreille,1811); P. lignarius/P. herreri (Walker, 1873); and P. rufotuberculatus (Champion, 1899) [6, 7]. The epidemiological importance of several of these species, considered to act as secondary vectors, has been dramatically increased by their invasion into human structures . Some species of Panstrongylus have a high capacity as Chagas disease vectors because of their longevity, rapid response to the presence of a host, the large volume of blood ingested, and frequent defecation during the feeding process .
Panstrongylus chinai is distributed in Loja and El Oro provinces, in the southern Andean region of Ecuador. This species has also been reported in the north of Peru, on the border with Loja, in different wild habitats and colonizing human dwellings . In Ecuador, the species is found in domestic and peridomestic environments such as chicken coops, pigeon and rat nests. Inside houses, it is found mainly under beds or in holes in the walls. In Loja (Ecuador) and some areas of Peru, P. chinai has been categorized as a potential vector of Chagas disease requiring continuous entomological surveillance .
Panstrongylus howardi is an endemic species for Ecuador, restricted to the Manabí Province (central coast region) [6, 7]. This species has been found only once in wild conditions  but is frequently collected in peridomestic habitats associated with rodent nests located between brick piles [6, 7, 11]. Abundant colonies can be also found in wood piles, as well as in the “piñuelas” plant Aechmea magdalenae, again in association with nesting places of rodents or marsupials . Panstrongylus howardi shows high levels of infection with Trypanosoma parasites (T. cruzi and T. rangeli), which suggests the existence of active transmission between this vector species and its associated vertebrate hosts .
Panstrongylus chinai and P. howardi have been regarded as belonging to the same evolutionary clade [3, 13]. Their morphological similarity has led some authors to suspect P. chinai to be a melanic form of P. howardi [13, 14]. The present study aims to explore this question through an integrative approach including quantitative morphology, chromosome analyses, molecular genetics, mating behavior and ecological arguments.
Integrative taxonomy has received considerable attention in the last years as a multidisciplinary approach that aims at delimiting the units of life’s diversity from multiple and complementary perspectives [15, 16]. It is useful to evaluate the species limits and is of great help in cases of recent speciation processes, cryptic diversity, or simply intraspecific variation (phenotypic plasticity) [17,18,19]. Historically, taxonomists have been concerned with classifying organisms into groups based on shared traits, and then further classifying those groups into the categories of the taxonomic hierarchy, from kingdom to species. In contrast, modern systematic biologists, despite the fact that they still use data taking the same basic form of similarities and differences among organisms, are increasingly devoting their efforts to testing hypotheses about lineage boundaries and phylogenetic relationships . A review of methods used in taxonomic studies related to the Triatominae concluded that an integrative approach was necessary, putting forward the concept of “evolutionary lineage” . The evolutionary concept allows the discussion about infraspecific populations and decisions about the utility of giving them subspecific ranks . It was applied successfully to various questions related to speciation in the subfamily Triatominae [23,24,25,26]. Current taxonomy of the Triatominae has been examined in the light of the main species concepts (morphological, biological, Hennigian, phylogenetic and evolutionary concepts) by Bargues et al. [19, 21].
Morphological variation in the Triatominae is clearly modulated by ecological factors. Features like color patterns or pilosity can change within a single species when populations are under disruptive selection or collected in a wide geographical range, a common observation which has been attributed to apparently weak canalization mechanisms . Developmental canalization mechanisms maintain the stability of a given phenotype in the face of environmental disturbance. The existence of distinct forms, morphs or ecotypes in many species is indicative of a weak developmental canalization in the Triatominae . This has been quantified between populations or among different species by studies on the size, shape (head or wing), or on the number of sensorial receptors (sensilla) covering the antennae [27,28,29,30,31,32]. In the Triatominae, the morphological variation of eggs among species and populations has also been used in qualitative [33,34,35,36,37] and quantitative analyses [38,39,40,41,42]. In this study, we examined the geometric morphometrical variation of both adults and eggs of these two species.
Chromosomal studies have been applied extensively in this subfamily to differentiate species or determine common characteristics in evolutionarily close species. Chromosomal variations in the amount of constitutive heterochromatin and in the position of the ribosomal clusters are the main intraspecific and interspecific features for karyotypic differentiation in the Triatominae [43, 44].
Molecular tools applied to Ecuadorian triatomines have been poorly explored . The mitochondrial cytochrome b (cytb) gene has been found to be a very useful molecular marker to differentiate cryptic species or closely related species in different triatomine genera [46,47,48].
Experimental crosses (EC) had an important role in clarifying controversial taxonomic issues within the Triatominae [49,50,51,52,53,54,55]. Hybridization studies may be very helpful to formulate hypotheses concerning the origin and divergence of species . In case of allopatry, it is not possible to ascertain that the hybridization (or lack of hybridization) observed in the laboratory would be indeed the natural behavior of these taxa if they could meet in natural conditions. Reproductive isolation, even when observed experimentally, is one of the best criteria to evaluate the taxonomic status of morphologically or genetically close populations .
A complementary approach to the systematics of the Triatominae has been the recent application of ecological niche modeling (ENM) . The ENM is a method that predicts the potential geographical distribution of a species under certain environmental conditions, even in territories where the species may have not reached yet [57,58,59,60,61]. The known distribution of P. howardi and P. chinai is restricted to separate territories in Ecuador, so that it is relevant to explore their potential geographical distribution. Panstrongylus chinai occurs in an arid zone (Loja), located in a wide range of altitudes (175 to 2003 meters above sea level, masl) , while P. howardi is restricted to the more humid and much lower (up to 400 masl) territory of the Manabí Province [10, 63].
The present study tries to discern if P. howardi and P. chinai are valid species or are two populations of the same species with morphological variations, mainly color and size differences, as proposed by Abad-Franch & Aguilar  and Patterson et al. . In the Triatominae, it has been shown that morphological divergence may occur rapidly, before the establishment of reproductive barriers . If this is the case, P. chinai and P. howardi should not, or barely, depart from what is expected between populations of the same species when compared through diverse biological analyses. For this reason, to demonstrate if P. chinai and P. howardi are chromatic forms or true species, the following approaches were used: comparative assessment of the antennal phenotype (APH); geometric morphometrics (GM); cytogenetics (CYT); molecular genetics (MG); experimental crosses (EC) including chromosomal analyses of the hybrids; and ecological niche modeling (ENM).
Study area and triatomine collection
In Ecuador, specimens under study are exclusively found in the Central Coastal region (P. howardi, Manabí Province) and in the southern Andean Region (P. chinai, Loja Province) (Fig. 1). The triatomines were collected in 11 communities of Manabí and 17 localities of Loja during the years 2004–2011 (Additional file 1: Table S1). Both species were collected in domestic and peridomestic habitats in Loja and Manabí provinces as described by Grijalva et al.  and Villacís et al. , respectively. The domestic searches included all the rooms inside each house. The peridomestic habitat was defined as all surrounding structures and potential microhabitats within 20 m of the house (including chicken coops and pigeon nests; accumulations of stones, tiles, wood, bricks, piles of agricultural products, storage buildings; and other structures near the houses).
Manabí Province is located along the Central Coast of Ecuador at an altitude ranging from 0 to 400 masl, and receives an average annual rainfall of 563 mm . The walls and floors of houses in this region are principally constructed with bamboo cane “caña guadúa” (Guadua angustifolia) or wood, the roofs are commonly made of the cade palm frowns, and in less proportion made of zinc .
The Province of Loja, in the southern Andean region of Ecuador, is characterized by a mix of hilly and mountainous topography, with an altitude ranging from 120 to 3800 masl. The region includes inter-Andean temperate valleys and has an average rainfall of 400 mm . Houses are typically made of adobe walls, the floors with dirt, and the roofs are made with ceramic tile [65, 67, 68].
Analysis of antennal phenotype
The antennae of 44 adult P. chinai (15 females and 29 males) and 23 P. howardi (8 females and 15 males) were analyzed following the methodology of Abrahan et al. . One antenna was cut from each individual at the level of the scape (segment 1), and slide-mounted in glycerine. By optical microscopy at 400× (Olympus BX41) with a lucid camera (Olympus U-DA 9E12246; Olympus corporation, Tokyo, Japan), the number of sensilla of different types was counted over the whole ventral surface of the three distal segments of the antenna, i.e. the pedicel (Pedi), and the two flagellar segments: flagellum 1 (Flag1) and flagellum 2 (Flag2). For this study, we considered the mechanoreceptors (bristles, BR) and three types of chemoreceptors: thin-walled trichoidea (TH); thick-walled trichoidea (TK); and basiconics (BA) (nomenclature following Catalá & Schofield  (Fig. 2). For each segment, the number of each kind of sensilla per ventral face (antennal phenotype) was used for statistical comparisons between sexes and species.
Morphometric analyses of adults and eggs followed two main steps: (i) extraction of size and shape variables; and (ii) discrimination between species or groups based on shape variables.
Wings and heads: a total of 79 right hemelytra and heads of 18 females and 25 males of P. chinai, and of 15 females and 21 males of P. howardi, were prepared for geometric morphometric analyses. The wings and heads were deposited on microscope slides. The wings were cautiously placed on a coverslip, and the heads were positioned parallel to the focal plane, using an entomological pin. A digital image of each organ was obtained using a MiScope®-MIP (http://www.zarbeco.com). Nine landmarks were identified for each wing (Fig. 3a) and 16 landmarks were identified for each head (Fig. 3b). For the head, the coordinates of corresponding landmarks on the left and right sides were averaged, thus reducing the number of landmarks to 8, prior to statistical analyses. Collection of landmarks (on heads and wings) and subsequent statistical analysis were performed using the various modules of the CLIC package (version 99), written by one of us (JPD) and freely available at http://xyom-clic.eu.
Eggs: the 78 P. chinai eggs came from 48 females (6 localities) and the 75 P. howardi eggs came from 8 females (one locality) (Additional file 1: Table S1). All of them were viable eggs, obtained from virgin, laboratory females coming from field-collected specimens. To ensure a reproducible protocol of image capture (MiScope-MIP), the 153 viable eggs were photographed one by one at the same developmental time (25 days of development) and exactly the same position (ventral) on a small platform (Fig. 4), as described in Santillán-Guayasamín et al. . For morphometric comparisons, we only considered the complete contour of the egg, thus including the operculum (Fig. 3c), using the outline-based approach [71, 72]. As for the landmark-based analyses of head and wings, pseudolandmarks were collected and analyzed using the CLIC package.
Testes were removed from freshly killed adults, fixed in an ethanol-acetic acid mixture (3:1) and stored at − 20 °C. Chromosome studies were made on 2 male specimens of P. chinai from Loja province and 5 specimens of P. howardi from Manabí Province (Additional file 1: Table S1). Meiotic processes of two male F1 hybrids resulting from the cross between a male of P. chinai (Loja) with a female of P. howardi (Manabí) and the reverse cross were also analyzed. Chromosomal techniques were applied using squashed gonad preparations made on slides in a drop of 50% acetic acid, freezing them in liquid nitrogen and removing the coverslip with a razor blade. The C-banding technique was performed as reported by Panzera et al.  in order to observe the distribution and behavior of C-heterochromatin during mitosis and meiosis. Fluorescence in situ hybridization (FISH) assay was carried out according to Panzera et al.  in order to determine the chromosomal location of the 45S ribosomal DNA clusters. Slides were analyzed using a Nikon Eclipse 80i epifluorescence microscope (Nikon, Kanagawa, Japan). Images were obtained with a Nikon DS-5Mc-U2 digital, cooled camera (Nikon). For analyses of cytogenetic images, we used the software IPP plus, Nikon Nis Elements 3.1 Advanced Research, and Adobe Photoshop.
For cytochrome b (cytb) marker, newly generated sequences for 26 P. chinai and 33 P. howardi deposited in the GenBank database were employed for MG analyses (Table 1, Additional file 2: Table S2). Additionally, sequences for three Panstrongylus species, i.e. P. megistus (GenBank: KC249230, KC249231) P. tupynambai (GenBank: KC249233, KC249234) and P. lutzi (GenBank: KC249262) and three Triatoma species, i.e. T. tibiamaculata (GenBank: KC249296, KC249297), T. brasiliensis (GenBank: AY336524) and T. infestans (GenBank: AY062165) were used for comparative purposes. The two species studied here represent 20 haplotypes: 9 of P. chinai and 11 of P. howardi.
Collected insects from six communities in Loja and one community in Manabí (Additional file 1: Table S1) were maintained under suitable conditions at the CISeAL (Quito, Ecuador). This facility is equipped with a “dual chamber incubator” where the original microhabitat temperature and humidity conditions were the same as those described in Villacís et al.  and Santillán-Guayasamín et al. . Fifth-instar nymphs (NV) of P. chinai and P. howardi were separated because virgin females were needed for crossing. Once the NV emerged to adults, both females and males were placed in plastic vials (9 × 5 × 4 cm) fitted with fan-folded Whatman filter paper to facilitate movement of specimens and absorb excess humidity. Blood meals were offered fortnightly using immobilized pigeons for 30 min. This experiment was conducted using protocol 15-H-034 approved by American Association for Laboratory Animal Science - IACUC.
Nine interspecific crosses were conducted and inspected daily looking for copulation posture and eggs. Five experimental crosses were performed for one female (♀) of P. chinai with two males (♂) of P. howardi, labeled “♀c × ♂h”; and four crosses were performed for one female of P. howardi with two males of P. chinai, labeled “♀h × ♂c”. Once the adults F1 emerged, we constituted nine couples in an attempt to produce F2 offspring. We also conducted control crosses for each species, 24 crosses for P. chinai and 4 for P. howardi. These parental crosses were composed of two females and three males. The eggs of these crosses were used in the morphometric analyses.
Ecological niche modeling (ENM)
We used data collected in Loja and Manabí provinces (Additional file 1: Table S1). Panstrongylus chinai specimens came from 41 communities in Loja Province, and P. howardi specimens from 27 communities in Manabí Province. Each locality was georeferenced in decimal coordinates (decimal coordinates have a reasonable error range ± 8–9 m). The ENM was designed to estimate the potential distribution area of each species [75, 76]. Given the suspicion that P. chinai and P. howardi could be local populations of the same species , we used the ENM method to determine if their potential distribution area might show some overlap. Due to the lack of relevant data from Peru, where P. chinai occupies a small territory at the border with Ecuador (Loja), we decided to define the entire coastal area up to 2000 masl of Ecuador, excluding Peru, as our calibrated area. We eliminated 4 of the 19 bioclimatic variables to avoid “distortions” to the distribution model . To determine the potential distribution areas for P. howardi and P. chinai, we used an environmental dataset from the WorldClim data archive , selecting a set of 10 “bioclimatic” data layers of 30 arc-seconds (~ 1 km). These data layers were selected after determining that they do not have a climatic correlation with each other. Each data layer’s contribution was high in a jacknife analysis of the model.
Bioclim variables included total monthly precipitation, maximum and minimum annual temperatures and other variables derived from monthly rainfall and temperature. These variables represent annual trends, seasonality and extreme or limiting environmental factors .
Phenotype: sensilla distribution on antennae
The number and type of antennal sensilla [bristles (BR) thin-walled trichoidea (TH) thick-walled trichoidea (TK) basiconic (BA)] in each segment (Pedi, Flag1, Flag2) were recorded and descriptive statistics were calculated. Homogeneity of variance was checked by the Levene test and normality was checked using the Kolmogorov-Smirnov test. Homoscedastic variables were analyzed using analysis of variance (one-way ANOVA), whereas heteroscedastic variables were analyzed using the Kruskal–Wallis non-parametric test. Univariate analyses of receptor densities were performed using the SPSS package (Statistical Package for Social Sciences) for Windows, version 19.0 (SPSS Inc., Chicago, Illinois). Based on the significant variables, discriminant analyses (one for each sex) between taxa and corresponding cross-validated reclassifications were performed with the CLIC package (version 99; https://xyom-clic.eu).
Size of wings and heads. Mean and variance of size were estimated in each sex and compared. Also, using a non-parametric ANOVA based on random sampling with replacement (bootstrap), the amount of sexual size dimorphism (SSD) of P. chinai was compared with that of P. howardi. In each sex of each species, the linear correlation coefficient (r) was computed between head and wing size dimensions of each individual.
Shape of wings and heads. Shape variables were the partial warps computed after the Procrustes superimposition of each configuration on the consensus configuration . Shape divergence between taxa was estimated in each sex by the Mahalanobis distance computed from these shape variables. In the same way as for the size variation, the observed Mahalanobis distances were subjected to non-parametric tests.
In addition to the comparison of mean shapes, we compared the variances of shape or metric disparity (DM), as described by Zelditch et al. . We performed shape-based validated reclassifications  between species in each sex, using either head or wing shape. As for size, we estimated the shape difference between sexes (“SShD”, for sexual shape dimorphism) and we verified that its importance was similar in both taxa. Using the same non-parametric test as for shape variation between species (see above), we tested in each species the significance of the Mahalanobis distance between sexes. To compare SShD between taxa, a two-factor MANCOVA was used to test shape variation against two main effects (species and sex) and against their interaction. Since shape difference is not only represented by the magnitude of shape change but also by its direction in the morphospace, we adopted the general approach for the statistical comparison of multivariate vectors of phenotypic change as described by Collyer & Adams , and applied to the Triatominae by Caro-Riaño et al. .
To explore the co-variation of head and wing shape, we computed the Escoufier’s “Rv” coefficient . It is a measure of correlation between two sets of variables, working directly on variance and covariance rather than on correlation matrices . The statistical significance of the Rv coefficient was assessed by randomizing 1000 times the covariance structure between blocks without affecting variance-covariance structure within the blocks, computing each time a pseudo-coefficient and counting the number of times where randomly generated Rv could be equal to or larger than the observed one.
GM, adults: size and shape. A multivariate regression was computed with size as independent variable and shape descriptors (partial warps) as dependent variables, and a non-parametric test of statistical significance was applied as described by Good . A multivariate test of significance for “different slopes” versus “common slope” model was computed according to the Wilk’s test, for heads and for wings, in both sexes. It is based on the comparison of two matrices derived from residuals of the regression of shape variables against size. The contribution of size to shape divergence between either sexes or species was computed as the coefficient of determination when regressing the (unique) discriminant factor on size.
GM, eggs: size and shape variation. Egg size was estimated as the perimeter of the egg contour and compared between groups using non-parametric tests. For shape variable definition, we exclusively used the elliptic Fourier analysis (EFA) . To accurately represent a closed curve, many harmonics may be needed, each with four coefficients, so that the number of variables could be too numerous relative to the number of specimens. The normalized coefficients (NEF) were thus submitted to a principal components analysis (PCA), and the principal components (PCs) were the final shape variables. The Mahalanobis distances between them were used to test for validated reclassification. The contribution of size variation to shape-based discrimination was estimated through the determination coefficient between the first (and unique) discriminant factor and the estimator of size. As for heads and wings, we compared the variance of size and of shape between taxa, we tested the hypothesis of a common allometric model, as well as the contribution of size to shape-based discrimination. All morphometric analyses (size and shape of adults and eggs) were performed using the CLIC package (version 99).
Molecular genetics: alignment, model selection and phylogenetic analyses
Alignments were performed using MAFFT . The GTR+I nucleotide substitution model was determined as the best-fitted model by model generator , under the Bayesian information criterion (BIC). The maximum likelihood (ML) tree was obtained by PhyML , with a nodal support of 1000 bootstrap pseudo-replicates. The phylogenetic trees were visualized and edited using FigTree , and ggtree . Intra- and interspecific sequence divergence for the cytb gene was assessed with Kimura 2-parameter (K2P) distance model in MEGA v6 .
Mean and standard deviation of pre-oviposition time were computed for each taxon and for hybrids. We also estimated the percentage of successful crosses, average number of eggs obtained by female, developing time and percentage of viable eggs. The latter was calculated as the ratio of total hatched eggs over total laid eggs. These analyses were performed using the SPSS package, version 19.0.
Ecological niche modeling
We used MaxEnt 3.4.1  to determine the potential distribution of each species; the software is useful for estimating a distribution across geographical space [95, 96], based on the presence or absence of some species. It used bioclimatic layers belonging to the Bioclim world climate database. The results of ecological niche modeling (ENM) analyses were exported to the ArcGIS (Geographic Information System) Software by Esri (Environmental Systems Research Institute) for Windows, version 10.3. (Release 10. Redlands, California). Because the number of points from presence was low, we utilized 10,000 points as background. To measure the robustness of model estimation, we used a cross-validation procedure, and a random seed with 10 replicates. To evaluate the validity of the model we used the “area under curve” (AUC estimation. Finally, we also estimated the index of niche overlap “D” using the software ENM tools for Windows, version 1.4.4 .
Table 2 shows the means and standard deviations for the number of antennal sensilla in both sexes and in each species without discriminating by habitat. Significant differences were detected between taxa for TH from all examined segments of the antenna (Pedi, Flag1 and Flag2). For males only, the BR of the pedicel also showed significant differences. For the two species, no sexual dimorphism in the antennal phenotype could be observed. Multivariate analyses (one for each sex) using as input the significant variables (4 for males, 3 for females) did not reveal significant Mahalanobis distances (P > 0.05) between the two species. Cross-check classification resulted in relatively low scores (from 56% to 68%).
Geometric morphometrics: adults
For head and wing shape (Fig. 5a, b), a significant difference was generally found between taxa for males and females (Table 3). The reclassification produced better results when based on wing shape relative to head shape, with scores generally better for males (85%, 92%, 76%, 80%) than for females (93%, 66%, 73%, 77%). The variance of shape (metric disparity, MD) did not show statistically significant differences between taxa, except for the variance of head shape, larger in male P. howardi (P = 0.012). Within each taxon, the variance of head shape (MD range of 10.1–16.8) was approximately two times the variance of wing shape (MD range of 6.8–8.6).
Size variation was illustrated for wings, heads and eggs by quantile boxes (Fig. 6). In both taxa, sexual dimorphism of size was not detected for heads, and was nearly significant for wings (Table 4). For the head, shape variation between sexes was found to be significant in both species. For the wing, shape variation was significant between sexes in P. chinai only. However, there was no significant difference in the amount of sexual dimorphism between taxa for head or wings (Table 4). Actually, the Procrustes distances between sexes were very similar between taxa for heads (0.036 for P. chinai and 0.031 for P. howard) and for wings (0.017 and 0.011, respectively).
Within each species and in each sex, for heads and for wings, the non-parametric multivariate regression of shape against size variation was not significant. As expected, for heads or for wings, both taxa shared a common allometric axis, suggesting similar patterns of growth (details not shown). A much stronger contribution of size to head shape divergence (52–58%) was found, relative to the contribution of size to wing shape difference (2–18%). Head shape divergence was dominated by size variation, and the reclassification of specimens based on head shape variation was less satisfactory than that based on wings (Table 3).
The covariation of head and wing shape was very similar between taxa in both sexes (0.19 for female P. chinai and 0.22 for female P. howardi, 0.29 and 0.35 for corresponding males). Contrasting with this similarity, head size was differently correlated with the wing size according to the species (Table 5). The correlation of size variation between head and wing was high (0.81 in males and 0.84 in females) and significant (P < 0.001) in P. chinai but much lower in P. howardi (0.25 in males and 0.52 in females), where it was only significant for females.
Geometric morphometrics: eggs
While P. howardi showed a larger size than P. chinai for both heads and wings, eggs in this species were significantly smaller. Egg shape variation between taxa was very significant (Table 6), with outline divergence visible without image amplification (Fig. 5c). The variance of size and shape showed significantly higher values for P. chinai, and, contrary to the adult morphometric variation, these values were different between taxa. Also contrary to the adults, the common allometric model was rejected (P = 0.040, details not shown). Finally, the reclassification scores based on the contour shape among taxa were excellent: 97% (76/78) for P. chinai and 94% (71/75) for P. howardi (Table 6). The Mahalanobis distance (6.27) was highly significant (P < 0.0001), and the influence of size on this shape-based distinction was 19%.
The comparative analyses of P. chinai and P. howardi from Ecuador studied here with previously published data for P. chinai from Peru , showed that the two species have the same chromosomal characteristics. C-banding revealed they have the same chromosome diploid number constituted by 23 chromosomes (2n = 20 autosomes + X1X2Y) in males and similar amount of heterochromatic regions. The ten autosomal bivalents exhibited a C-positive heterochromatic block in both chromosomal ends (Fig. 7a, b). The Y sex chromosome was totally heterochromatic while the X1 and X2 sex chromosomes were the smallest of the complement and presented an intermediate staining (Fig. 7b, c). According to the FISH technique, the 45S ribosomal DNA cluster was located in the largest autosomal pair in both taxa (Fig. 7d).
The meiotic analysis of two male adults (F1) obtained from the two interspecific crosses showed important anomalies in the chromosomal behavior. In pachytene, association of different bivalents forming chains, univalents (consequence of the lack of pairing between homologous chromosomes) and chromosome fragments are observed, all of them representing products of alterations in chromosomal pairing (Fig. 7e). During late diplotene, associated bivalents and univalents were observed (Fig. 7f). The proportion of bivalents and univalents varied not only between the two analyzed individuals, but also between cells of the same specimen. These abnormalities produced metaphases I (Fig. 7g) or metaphases II (Fig. 7h) with variable numbers of autosomes and sex chromosomes. All metaphases I (Fig. 7g) and metaphases II (Fig. 7h) were aneuploid, with deficiency or excess of chromosomes. Normal first and second metaphases were not observed.
Intra-species genetic K2P distances for P. chinai (26 individuals) and P. howardi (33 individuals) varied between 0–0.94% (average 0.32%) and 0–1.9% (average 0.28%), respectively (Additional file 2: Table S2). Genetic distance between P. chinai and P. howardi was 7.82% on average (range 6.25–8.73%). Lastly, P. chinai and P. howardi distance with the other three Panstrongylus species ranged between 17.93–24.11%. In the ML tree (Fig. 8), all nodes showed high support (> 70) and strongly supported the monophyly of P. chinai and P. howardi (bootstrap value 100).
Copulation and viable eggs were observed in seven out of the nine interspecific couples (F1 generation). Most of the couples were successful, 4 out of 5 of the ♀c × ♂h crosses and 3 out of 4 of the ♀h × ♂c crosses. All adult F1 hybrids of each sex showed a general color typical of the P. howardi phenotype, while they showed a general size and shape closer to the P. chinai phenotype (Fig. 9). In spite of the daily examination of the crosses between F1 hybrids, no copulation was observed. Eggs appeared in 6 out of the 9 couples, all of them empty, with no visible embryo; these eggs promptly shrank and never hatched.
The pre-oviposition period was slightly lower when the females were crossed with males of the same species (parental crosses) (Table 7). The average number of eggs was the greatest for P. howardi, but their viability was two times lower than for other crosses (Table 7). No significant differences were observed in the time of development between the parental crosses with the interspecific crosses.
Ecological niche modeling
The ENM model for P. chinai adequately predicted the current distribution of this species but also predicted its occurrence in the north of the province of Loja, between Loja, Azuay and El Oro provinces (Fig. 10a). Figure 10b shows the distribution of P. howardi in Manabí Province, with a geographical expansion of this species to Esmeraldas, Guayas, El Oro and some parts of Loja. Figure 10c provides the predictive models of P. chinai and P. howardi together, showing the overlapping potential areas of P. howardi with P. chinai in Loja Province. Thus, Loja and Manabí apparently share environmental characteristics suitable for the colonization by both species. For P. chinai, the model fit had an AUC of 98%, while for P. howardi the AUC was 96%. The index of the niche overlap was D = 0.140 with P > 0.05, suggesting similarity rather than equivalence.
The distribution model for P. chinai showed that the mean diurnal range (maximum temperature − minimum temperature, BIO2; Table 8) significantly contributed (77.7%) to the presence of this species, together with the annual mean temperature (BIO1, 7.6%) and the temperature seasonality (BIO4, 5.8%; Table 8). The model fit had an AUC of 98%. The distribution model for P. howardi showed that BIO12 (annual precipitation) (Table 8) had the highest contribution (61.3%), followed by the temperature seasonality (BIO4, 17.3%) and the annual mean temperature (BIO1, 14.2%; Table 8).
Species of the Triatominae display a high degree of morphological plasticity . Due to this, entomologists face situations where morphological differences exist between groups (mainly size or/and color variation) without clear genetic support, and others where genetic variation seems important without morphological discriminating traits, or even without established reproductive isolation . Thus, morphological or genetic variation should not lead per se to the erection of a new species. In the spirit of integrative taxonomy, we followed the recommendations of Bargues et al. [19, 21]: “a multidisciplinary approach would converge to the hypothesis of an evolutionary lineage, or, conversely, to the idea of an intraspecific variation”. In the speciation process, the different biological properties examined here are not expected to evolve at the same time and in a regular order [20, 99], so that opposite trends could be expected. We observed indeed incongruous results preventing us to modify the current taxonomic status of P. chinai and P. howardi.
Most of the specimens (94%) were collected in domestic and peridomestic environments; however, most P. chinai (95%) were found inside houses while most P. howardi (90%) were collected in peridomestic habitats. In addition to a behavioral difference between the two species, this observation probably reflects the difference in habitats between Loja and Manabí. Houses in Loja are made of adobe, where the bugs can hide , while houses in Manabí have cane walls not suitable for bug infestation [11, 63]. Different rates of trypanosomatid infection were reported for these two taxa, with T. cruzi prevalence of 57.5% for P. howardi  and 12.8% for P. chinai .
According to a cladistic analysis based on morphological characters, Lent & Wygodzinsky  suggested that P. chinai and P. howardi constitute a separate clade of two closely related species in the genus Panstrongylus. What justified their distinct species status was mainly the clear-cut difference in color. Many species in the Triatominae may be recognized thanks to their color patterns. In the case of P. chinai and P. howardi, the color difference could suggest the existence of dark or melanic forms as already observed within some other species, such as T. infestans , R. stali  and R. nasutus . Such melanization is not necessarily suggestive of distinct species [29, 102].
The antennal phenotype has proven to be a useful tool for the taxonomic comparison among species of the Triatominae . Between geographically separate populations, significant differences may also arise, as it was the case for wild and domestic populations of T. infestans , R. ecuadoriensis  and six other Triatoma species . The differences disclosed here between the antennae of P. howardi and P. chinai were lower than those found between populations of R. ecuadoriensis .
We detected a lower morphometric variation of the head and wings than expected for distinct species. Altogether, the amount of metric divergence could be in conformity with local populations of the same species. Size difference is a common finding between local populations. Despite living in a warmer climate, the Manabí specimens tended to be larger, which aligns with Bergmann’s rule. This rule also applies to intraspecific populations as a species criterion by Marcondes et al. . However, this argument may not be as relevant amongst populations located close to the equatorial line . For instance, the Bergmann’s rule was contradicted for populations of R. ecuadoriensis collected in the same provinces, Loja and Manabí . The clear-cut difference in the correlations of head and wing size is not an argument suggesting different species; it has been observed between populations of the same species from different habitats . Shape differences are expected between species, but they are also common between populations, either due to isolation or adaptation to different environments, or both. Shape changes per se cannot distinguish the variation due to simple genetic drift between isolated subpopulations from that due to evolutionary divergence between species. However, the relatively low reclassification scores, the lack of difference in the magnitude of sexual dimorphism, the apparent lack of difference in the covariation of head-wing shape, and a common allometric axis do not advocate for very different organisms.
Metric differences of the eggs between taxa, especially shape differences, were clearly more pronounced than those of head or wing. By comparing laboratory lines reared under the same conditions, fed on the same blood source at the same frequency, using eggs at similar developmental stages and positions, we could reduce the possible environmental influence on shape variation. However, the metric differences between eggs of both taxa could be attributed in part to the relatively restricted spatial and temporal sampling of one of them. The eggs of P. howardi were collected from a single community, on the same day. As a consequence, the P. howardi sample could provide a truncated representation of the species variability, as suggested by its lower variance of size and shape (Table 6). However, another reason for the clear-cut difference found between P. howardi and P. chinai could be a true evolutionary difference.
C-heterochromatin distribution and the chromosome position of ribosomal clusters are the principal cytogenetic traits to aboard the analysis of karyotype differentiation in the Triatominae [42, 43]. In several triatomine groups, such as infestans and sordida subcomplexes, closely related species exhibit striking differentiation for these cytogenetic markers. However, in other subcomplexes such as phyllosoma and brasiliensis, valid species did not present cytogenetic differentiation . Therefore, the chromosomal identity found here between P. chinai and P. howardi does not mean that they should be the same species (Fig. 7a–d). In fact, the chromosomal analysis of the hybrids showed significant alterations in the chromosome pairing of the bivalents which reveals deep chromosomal differences in the parental species not detected with the techniques applied here (Fig. 7e–h). In this way, the meiotic process of the hybrids does not end correctly so that viable gametes are not formed and, as a consequence, these hybrid individuals are infertile. Thus, the meiotic analysis of F1 hybrids strongly suggests that P. chinai and P. howardi are distinct species.
The cytb phylogenetic analyses strongly suggest that P. chinai and P. howardi are included in the same clade, representing sister species (Fig. 8). Intraspecific average distances were 0.32% and 0.28% for P. chinai and P. howardi, respectively. These distances are within the expected range for differentiation within a species, which indicates that the species identification for each individual was correct with this marker. The average distance value K2P (7.8%) between P. chinai and P. howardi is slightly higher than that reported to separate other closely related triatomine species (e.g. 7.5% for the cytb marker) . Thus, our molecular data support considering P. chinai and P. howardi as two separate sister species, in agreement with the suggestion by Barnabe et al. .
The hypothesis of a single species was not supported by our experimental crosses and the chromosomal analyses of the F1 hybrids. The main results of our crossing experiments were: (i) the ease to obtain a viable F1 offspring between the two species studied; and (ii) the failure to obtain an F2 offspring due to the lack of copulation between the F1 hybrids. Inter-specific hybrid generation (adult F1), both naturally and experimentally, is a common phenomenon in the Triatominae, and occurs between both closely and evolutionarily distant species of the genera Meccus, Rhodnius and Triatoma . The ease of obtaining inter-specific F1 adults reveals that several pre-zygotic isolation mechanisms can be easily avoided by experimental hybridization. For example, the morphology of the genital structures does not represent an important barrier to mating that prevents interspecific crosses in some Triatominae. However, the great majority of these hybrids are infertile, unable to produce viable offspring . In the Triatominae, interspecific hybrids can be totally fertile or have different degrees of infertility . In the species studied here, the lack of copulation among the F1 hybrids and their inability to produce viable gametes (chromosome studies in males) reveal that different pre-zygotic isolation mechanisms are acting which lead to a total infertility of the F1 hybrids. In conclusion, our experimental crosses and the chromosome studies of the F1 hybrids strongly suggest that both Panstrongylus taxa are separate species.
The ENM correctly predicted the current distribution of P. chinai and P. howardi, but also predicted other locations with a partial overlap. This overlap, restricted to a few sites, suggests that both species share similar not equivalent niches. The ENM analyses suggest that both taxa could be different species.
Overall, antennal morphology, geometric morphometrics of the head and wing shape, and cytogenetic analysis did not indicate distinct differences between P. howardi and P. chinai. Conversely, geometric morphometrics of the eggs, ENM analyses, molecular phylogeny, and crossing experiments (including chromosome analyses of the hybrids), in addition to the color pattern and current distribution, support the hypothesis that these bugs are separate species. We conclude that P. howardi and P. chinai should not be synonymized, they represent two valid closely related species.
Availability of data and materials
All data generated or analyzed during this study are included within the article and its additional files. The newly generated sequences were submitted to the GenBank database under the accession numbers JX400933-JX401001.
area under curve
index of niche overlap
elliptic Fourier analysis
ecological niche modelling
fluorescence in situ hybridization
meters above sea level
normalized elliptic Fourier coefficients
partial least square
sexual size dimorphism
sexual shape dimorphism
WHO. Chagas disease (American trypanosomiasis). Geneva: World Health Organization; 2019. http://www.who.int/mediacentre/factsheets/fs340/en/index.html. 2019. Accessed 20 Jan 2019.
WHO. Control of Chagas disease. Second report of the WHO expert committee. WHO technical report series. Geneva: World Health Organization; 2002.
Lent H, Wygodzinsky P. Revision of the Triatominae (Hemiptera: Reduviidae), and their significance as vectors of Chagas disease. Bull Am Mus Nat Hist. 1979;163:123–520.
Grijalva MJ, Villacís AG. Presence of Rhodnius ecuadoriensis in sylvatic habitats in the southern highlands (Loja Province) of Ecuador. J Med Entomol. 2009;46:708–11.
Villacís AG, Grijalva MJ, Catalá SS. Phenotypic variability of Rhodnius ecuadoriensis populations at the Ecuadorian central and southern Andean region. J Med Entomol. 2010;47:1034–43.
Abad-Franch F, Paucar A, Carpio C, Cuba CA, Aguilar HM, Miles MA. Biogeography of Triatominae (Hemiptera: Reduviidae) in Ecuador: implications for the design of control strategies. Mem Inst Oswaldo Cruz. 2001;96:611–20.
Villacís AG, Ocaña-Mayorga S, Lascano MS, Yumiseva CA, Baus EG, Grijalva MJ. Abundance, natural infection with trypanosomes, and food source of an endemic species of Triatomine, Panstrongylus howardi (Neiva, 1911), on the Ecuadorian Central Coast. Am J Trop Med Hyg. 2015;92:187–92.
dos Santos CM, Jurberg J, Galvão C, da Silva D, Rodríguez J. Estudo morfométrico do gênero Panstrongylus Berg, 1879 (Hemiptera, Reduviidae, Triatominae). Mem Inst Oswaldo Cruz. 2003;98:939–44.
Mosquera KD, Villacís AG, Grijalva MJ. Life cycle, feeding, and defecation patterns of Panstrongylus chinai (Hemiptera: Reduviidae: Triatominae) under laboratory conditions. J Med Entomol. 2016;53:776–81.
Suarez-Davalos V, Dangles O, Villacís AG, Grijalva MJ. Microdistribution of sylvatic triatomine populations in Central-Coastal Ecuador. J Med Entomol. 2010;47:80–8.
Grijalva MJ, Villacís AG, Ocaña-Mayorga S, Yumiseva CA, Baus EG. Limitations of selective deltamethrin application for triatomine control in central coastal Ecuador. Parasit Vectors. 2011;4:20.
Pinto CM, Ocaña-Mayorga S, Lascano MS, Grijalva MJ. Infection by trypanosomes in marsupials and rodents associated with human dwellings in Ecuador. J Parasitol. 2006;92:1251–5.
Patterson JS, Barbosa SE, Feliciangeli MD. On the genus Panstrongylus Berg, 1879: evolution, ecology and epidemiological significance. Acta Trop. 2009;110:187–99.
Abad-Franch F, Aguilar VHM. Control de la Enfermedad de Chagas en el Ecuador OPS/OMS (Publicación auspiciada por Ayuda Popular Noruega, Catholic Relieve Services, COOPI, Médicos Sin Fronteras y Oxfam) Quito, Ecuador. 2003. http://www.opsecu.org/publicaciones/OPS.doc. Accessed 20 Nov 2004.
Wake MH. What is “Integrative Biology ”? Integr Comp Biol. 2003;43:239–41.
Dayrat B. Towards integrative taxonomy. Biol J Linn Soc Lond. 2005;85:407–15.
Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, Winker K, Ingram KK, Das I. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2006;22:148–55.
Abad-Franch F, Pavan MG, Jaramillo-O N, Palomeque FS, Dale C, Chaverra D, Monteiro FA. Rhodnius barretti, a new species of Triatominae (Hemiptera: Reduviidae) from western Amazonia. Mem Inst Oswaldo Cruz. 2013;108:92–9.
Bargues MD, Schofield CJ, Dujardin JP. Classification and phylogeny of the Triatominae. In: Telleria J, Tibayrenc M, editors. American trypanosomiasis Chagas disease: one hundred years of research. Amsterdam: Elsevier; 2010. p. 117–48.
de Queiroz KA. Unified concept of species and its consequences for the future of taxonomy. Proc Calif Acad Sci. 2005;56:196–215.
Bargues MD, Schofield C, Dujardin JP. Classification and systematics of the Triatominae. In: Telleria J, Tibayrenc M, editors. American trypanosomiasis (Chagas disease): one hundred years of research. 2nd ed. Amsterdam: Elsevier; 2017. p. 113–43.
Calleros L, Panzera F, Bargues MD, Monteiro FA, Klisiowicz DR, Zuriaga MA, Mas-Coma S, Pérez R. Systematics of Mepraia (Hemiptera-Reduviidae): cytogenetic and molecular variation. Infect Genet Evol. 2010;10:221–8.
da Rosa JA, Rocha CS, Gardim S, Pinto MC, Mendonça VJ, Filho CJRF, et al. Description of Rhodnius montenegrensis n. sp. (Hemiptera: Reduviidae: Triatominae) from the state of Rondônia. Brazil. Zootaxa. 2012;3478:62–6.
Campos-Soto R, Panzera F, Pita S, Lages C, Solari A, Botto-Mahan C. Experimental crosses between Mepraia gajardoi and M. spinolai and hybrid chromosome analyses reveal the occurrence of several isolation mechanisms. Infect Genet Evol. 2016;45:205–12.
Mendonça VJ, Alevi KCC, Pinotti H, Gurgel-Gonçalves R, Pita S, Guerra AL, Panzera F, de Araújo RF, Azeredo-Oliveira MTV, da Rosa JA. Revalidation of Triatoma bahiensis Sherlock & Serafim, 1967 (Hemiptera: Reduviidae) and phylogeny of the T. brasiliensis species complex. Zootaxa. 2016;4107:239–54.
Dujardin JP, Costa J, Bustamante D, Jaramillo N, Catalá SS. Deciphering morphology in Triatominae: the evolutionary signals. Acta Trop. 2009;110:101–11.
Matias A, De la Riva JX, Torrez M, Dujardin JP. Rhodnius robustus in Bolivia identified by its wings. Mem Inst Oswaldo Cruz. 2001;96:947–50.
Villegas J, Feliciangeli M, Dujardin J. Wing shape divergence between Rhodnius prolixus from Cojedes (Venezuela) and Rhodnius robustus from Mérida (Venezuela). Infect Genet Evol. 2002;2:121–8.
Gumiel M, Catalá S, Noireau F, Rojas de Arias A, García A, Dujardin JP. Wing geometry in Triatoma infestans (Klug) and T melanosoma Martinez, Olmedo & Carcavallo (Hemiptera: Reduviidae). Syst Entomol. 2003;28:173–80.
Hernández ML, Abrahan L, Moreno M, Gorla D, Catalá S. Phenotypic variability associated to genomic changes in the main vector of Chagas disease in South America. Acta Trop. 2008;106:60–7.
Carbajal de la Fuente AL, Jaramillo N, Barata JMS, Noireau F, Diotaiuti L. Misidentification of two Brazilian triatomes, Triatoma arthurneivai and Triatoma wygodzinskyi, revealed by geometric morphometrics. Med Vet Entomol. 2011;25:178–83.
Barata JMS. Aspectos morfológicos de ovos de Triatominae. Rev Saude Publica. 1981;15:490–542.
Barata JMS. Macroscopic and exochorial structures of Triatominae eggs (Hemiptera, Reduviidae). In: Carcavallo RU, Girón IG, Jurberg J, Lent H, editors. Atlas of Chagas’ disease vectors in the Americas. Rio de Janeiro: Editora Fiocruz; 1998. p. 409–29.
Obara MT, Barata JMS, da Silva NN, Ceretti Júnior W, Urbinatti PR, da Rosa JA, Jurberg J, Galvão C. Estudo de ovos de quatro espécies do gênero Meccus (Hemiptera, Reduviidae, Triatominae), vetores da doença de Chagas. Mem Inst Oswaldo Cruz. 2007;102:13–9.
dos Santos CM, Jurberg J, Galvão C, da Rosa JA, Júnior WC, Barata JMS, Obara MT. Comparative descriptions of eggs from three species of Rhodnius (Hemiptera: Reduviidae: Triatominae). Mem Inst Oswaldo Cruz. 2009;104:1012–8.
Rivas N, Sánchez ME, Martínez-Ibarra A, Camacho AD, Tovar-Soto A, Alejandre-Aguilar R. Morphological study of eggs from five Mexican species and two morphotypes in the genus Triatoma (Laporte, 1832). J Vector Ecol. 2013;38:90–6.
Monte Gonçalves TC, Jurberg J, Costa J, de Souza W. Estudio morfológico comparativo de ovos e ninfas de Triatoma maculata (Erichson, 1848) e Triatoma pseudomaculata Correa & Espinola, 1964 (Hemiptera, Reduviidae, Triatominae). Mem Inst Oswaldo Cruz. 1985;80:263–76.
Costa J, Barth OM, Marchon-Silva V, de Almeida CE, Freitas-Sibajev M, Goreti R, Panzera F. Morphological studies on the Triatoma brasiliensis Neiva, 1911 (Hemiptera, Reduviidae, Triatominae). Genital structures and eggs of different chromatic forms. Mem Inst Oswaldo Cruz. 1997;92:493–8.
Barbosa SE, Dujardin JP, Soares RPP, Pires HHR, Margonari C, Romanha AJ, et al. Interpopulation variability among Panstrongylus megistus (Hemiptera: Reduviidae) from Brazil. J Med Entomol. 2003;40:411–20.
Obara MT, da Rosa JA, da Silva NN, Ceretti Júnior W, Urbinatti PR, Barata JMS, et al. Estudo morfológico e histológico dos ovos de seis espécies do gênero Triatoma (Hemiptera: Reduviidae). Neotrop Entomol. 2007;36:798–806.
Santillán-Guayasamín S, Villacís AG, Grijalva MJ, Dujardin JP. The modern morphometric approach to identify eggs of Triatominae. Parasit Vectors. 2017;10:55.
Panzera F, Pérez R, Panzera Y, Ferrandis I, Ferreiro MJ, Calleros L. Cytogenetics and genome evolution in the subfamily Triatominae (Hemiptera, Reduviidae). Cytogenet Genome Res. 2010;128:77–87.
Panzera Y, Pita S, Ferreiro MJ, Ferrandis I, Lages C, Pérez R, Silva AE, Guerra M, Panzera F. High dynamics of rDNA cluster location in kissing bug holocentric chromosomes (Triatominae, Heteroptera). Cytogenet Genome Res. 2012;138:56–7.
Monteiro FA, Wesson DM, Dotson EM, Schofield CJ, Beard CB. Phylogeny and molecular taxonomy of the Rhodniini derived from mitochondrial and nuclear DNA sequences. Am J Trop Med Hyg. 2000;62:460–5.
Villacís AG, Marcet PL, Yumiseva CA, Dotson EM, Tibayrenc M, Brenière SF, Grijalva MJ. Pioneer study of population genetics of Rhodnius ecuadoriensis (Hemiptera: Reduviidae) from the central coast and southern Andean regions of Ecuador. Infect Genet Evol. 2017;53:116–27.
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.
Abad-Franch F, Monteiro FA. Molecular research and the control of Chagas disease vectors. An Acad Bras Cienc. 2005;77:437–50.
Mas-Coma S, Bargues MD. Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA. Acta Trop. 2006;110:112–36.
Galíndez-Girón I, Barazarte R, Márquez J, Oviedo M, Márquez Y, Morón L, Carcavallo RU. Relaciones reproductivas entre Rhodnius prolixus Stål y Rhodnius robustus Larrousse (Hemiptera, Reduviidae, Triatominae) bajo condiciones de laboratorio. Entomol Vectores. 1994;1:3–13.
Pérez R, Hernández M, Quintero O, Scvortzoff E, Canale D, Méndez L, Cohanoff C, Martino M, Panzera F. Cytogenetic analysis of experimental hybrids in species of Triatominae (Hemiptera: Reduviidae). Genetica. 2005;125:261–70.
Martínez-Ibarra JA, Ventura-Rodríguez LV, Meillón-Isáis K, Barajas-Martínez HM, Alejandre-Aguilar P, Lupercio-Coronel R, Rocha-Chávez G, Nogueda-Torres B. Biological and genetic aspects of crosses between species of the Phyllosoma complex (Hemiptera: Reduviidae Triatominae). Mem Inst Oswaldo Cruz. 2008;103:236–43.
Costa J, Peterson AT, Dujardin JP. Morphological evidence suggests homoploid hybridization as a possible mode of speciation in the Triatominae (Hemiptera, Heteroptera, Reduviidae). Infect Genet Evol. 2009;9:263–70.
Martínez-Hernández F, Martínez-Ibarra JA, Catalá S, Villalobos G, De la Torre P, Laclette JP, Alejandra-Aguilar R, Espinoza B. Natural crossbreeding between sympatric species of the Phyllosoma complex (Insecta: Hemiptera: Reduviidae) indicates the existence of only one species with morphologic and genetic variations. Am J Trop Med Hyg. 2010;82:74–82.
Correia N, Almeida CE, Lima-Neiva V, Gumiel M, Lima MM, Medeiros LMO, et al. Crossing experiments confirm Triatoma sherlocki as a member of the Triatoma brasiliensis species complex. Acta Trop. 2013;128:162–7.
Mendonça VJ, Alevi KC, Medeiros LM, Nascimento JD, de Azeredo-Oliveira MT, da Rosa JA, et al. Cytogenetic and morphologic approaches of hybrids from experimental crosses between Triatoma lenti Sherlock & Serafim, 1967 and T. sherlocki Papa et al., 2002 (Hemiptera: Reduviidae). Infect Genet Evol. 2014;26:123–31.
Souza NA, Andrade-Coelho CA, Vigoder FM, Ward RD, Peixoto AA. Reproductive isolation between sympatric and allopatric Brazilian populations of Lutzomyia longipalpis s.l. (Diptera: Psychodidae). Mem Inst Oswaldo Cruz. 2008;103:216–9.
Gorla DE. Variables ambientales registradas por sensores remotos como indicadores de la distribución geográfica de Triatoma infestans. Austral Ecol. 2002;12:117–27.
Peterson AT, Sánchez-Cordero V, Beard B, Ramsey JM. Ecological niche modeling and potential reservoirs for Chagas disease, México. Emerg Infect Dis. 2002;8:662–7.
Sandoval-Ruiz CA, Zumaquero-Rios JL, Rojas-Soto OR. Predicting geographic and ecological distributions of triatomine species in the southern Mexican state of Puebla using ecological niche modeling. J Med Entomol. 2008;45:540–6.
Peterson AT. Ecologic niche modeling and spatial patterns of disease transmission. Emerg Infect Dis. 2006;12:1822–6.
Warren DL, Seifert SN. Ecological niche modeling in Maxent: the importance of model complexity and performance of model selection criteria. Ecol Appl. 2011;21:335–42.
Grijalva MJ, Villacís AG, Ocaña-Mayorga S, Yumiseva CA, Moncayo AL, Baus EG. Comprehensive survey of domiciliary triatomine species capable of transmitting Chagas disease in southern Ecuador. PLoS Negl Trop Dis. 2015;9:e0004142.
Grijalva MJ, Villacís AG, Moncayo AL, Ocaña-Mayorga S, Yumiseva CA, Baus EG. Distribution of triatomine species in domestic and peridomestic environments in central coastal Ecuador. PLoS Negl Trop Dis. 2017;11:e0005970.
Dujardin JP, Panzera P, Schofield CJ. Triatominae as a model of morphological plasticity under ecological pressure. Mem Inst Oswaldo Cruz. 1999;94:223–8.
Grijalva MJ, Palomeque-Rodríguez F, Costales JA, Davila S, Arcos-Teran L. High household infestation rates by synanthropic vectors of Chagas disease in southern Ecuador. J Med Entomol. 2005;42:68–74.
INAMHI: Servicio Metereológico: características generales del clima en la provincia de Manabí y Loja-Ecuador. 2020. http://www.serviciometeorologico.gob.ec. Accessed 18 Mar 2020.
Black CL, Ocaña-Mayorga S, Riner D, Costales JA, Lascano MS, Davila S, et al. Household risk factors for Trypanosoma cruzi seropositivity in two geographic regions of Ecuador. J Parasitol. 2007;93:12–6.
Nieto-Sánchez C, Baus EG, Guerrero D, Grijalva MJ. Positive deviance study to inform a Chagas disease control program in southern Ecuador. Mem Inst Oswaldo Cruz. 2015;110:299–309.
Abrahan L, Hernández L, Gorla D, Catalá S. Phenotypic diversity of Triatoma infestans at the microgeographic level on the Gran Chaco of Argentina and the Andean valleys of Bolivia. J Med Entomol. 2008;45:660–6.
Catalá S, Schofield C. Antennal sensilla of Rhodnius. Int J Morphol. 1994;219:193–204.
Lu KH. Harmonic analysis of the human face. Biometrics. 1965;21:491–505.
Dujardin JP, Kaba D, Solano P, Dupraz M, McCoy KD, Jaramillo N. Outline-based morphometrics, an overlooked method in arthropod studies? Infect Genet Evol. 2014;28:704–14.
Villacís AG, Arcos-Teran L, Grijalva MJ. Life cycle, feeding and defecation patterns of Rhodnius ecuadoriensis (Lent & León 1958) (Hemiptera: Reduviidae: Triatominae) under laboratory conditions. Mem Inst Oswaldo Cruz. 2008;103:690–5.
Santillán-Guayasamín S, Villacís AG, Grijalva MJ, Dujardin JP. Triatominae: does the shape change of non-viable eggs compromise species recognition? Parasit Vectors. 2018;11:543.
Elith J, Leathwick J. Species distribution models: ecological explanation and prediction across space and time. Annu Rev Ecol Evol Syst. 2009;40:677–97.
Franklin J. Mapping species distributions. Spatial inference and prediction. Cambridge: Cambridge University Press; 2009.
Escobar LE, Lira-Noriega A, Medina-Vogel G, Peterson AT. Potential for spread of the white-nose fungus (Pseudogymnoascus destructans) in the Americas: use of Maxent and NicheA to assure strict model transference. Geospat Health. 2014;9:221–9.
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.
Nix HA. A biogeographic analysis of Australian elapid snakes. In: Longmore R, editor. Atlas of elapid snakes of Australia. Australian Flora and Fauna Series No. 7. Canberra: Australian Government Publishing Service; 1986. p. 4–15.
Rohlf FJ. Rotational fit (Procrustes) methods. In: Rohlf F, Bookstein F, editors. Proceedings of the Michigan Morphometrics Workshop. Special publication number 2. Ann Arbor: The University of Michigan Museum of Zoology; 1990. p. 380.
Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Geometric morphometrics for biologists: a primer. New York: Elsevier Academic Press; 2004.
Manly BFJ. Multivariate statistical methods: a primer. 3rd ed. Boca Raton: Chapman & Hall; 2004.
Collyer M, Adams DC. Analysis of two-state multivariate phenotypic change in ecological studies. Ecology. 2007;88:683–92.
Caro-Riaño H, Jaramillo N, Dujardin JP. Growth changes in Rhodnius pallescens under simulated domestic and sylvatic conditions. Infect Genet Evol. 2009;9:162–8.
Escoufier Y. Le traitement des variables vectorielles. Biometrics. 1973;29:751–60.
Claude J. Morphometrics with R. New York: Springer; 2008.
Good P. Permutation tests: a practical guide to resampling methods for testing hypotheses. New York: Springer; 2000.
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.
Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McInerney JO. Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006;6:29.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Rambaut, A. FigTree Version 1.4.2. 2014. http://tree.bio.ed.ac.uk/software/figtree/. Accessed 15 Jan 2019.
Yu G, Smith DK, Zhu H, Guan Y, Lam TTY. Ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution. 2008;62:2868–83.
Phillips SJ, Anderson RP, Schapire RE. Maximum, entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.
Phillips SJ, Dudik M. Modeling of species distributions with Maxent: new extensions and comprehensive evaluation. Ecography. 2008;31:161–75.
Pérez R, Hernández M, Caraccio M, Rose V, Valente SA, Valente VC, Moreno J, Angulo V, Ramírez CM, Roldán J, Vargas F, Wolff M, Panzera F. Chromosomal evolution trends of the genus Panstrongylus (Hemiptera, Reduviidae), vectors of Chagas disease. Infect Genet Evol. 2002;2:47–56.
Dujardin JP, Schofield CJ, Panzera F. Los vectores de la enfermedad de Chagas. Investigaciones taxonómicas, biológicas y genéticas. Tome 25, fasc. 3. Brussels: Académie Royale des Sciences d’Outre-Mer, Classe des Sciences naturelles et médicales. 2002.
de Queiroz K. The general lineage concept of species, species criteria, and the process of speciation. In: Howard DJ, Berlocher SH, editors. Endless forms: species and speciation. Oxford: Oxford University Press; 1998. p. 57–75.
Noireau F, Flores R, Gutiérrez T, Dujardin JP. Detection of silvatic dark morphs of Triatoma infestans in the Bolivian Chaco. Mem Inst Oswaldo Cruz. 1997;92:583–4.
Dias FBS, Jaramillo-O N, Diotaiuti L. Description and characterization of the melanic morphotype of Rhodnius nasutus Stål, 1859 (Hemiptera: Reduviidae: Triatominae). Rev Soc Bras Med Trop. 2014;47:637–41.
Monteiro FA, Pérez R, Panzera F, Dujardin JP, Galvão C, Rocha D, et al. Mitochondrial DNA variation of Triatoma infestans populations and its implication on the specific status of T. melanosoma. Mem Inst Oswaldo Cruz. 1999;94:229–38.
Rodríguez RJ, Catalá S, González OF, Nodarse AJF. Antennal phenotype of two Cuban species of flavida: Triatoma flavida and Triatoma bruneri (Reduviidae: Triatominae). Rev Cubana Med Trop. 2009;61:226–31.
Carbajal de la Fuente AL, Catalá S. Relationship between antennal sensilla pattern and habitat in six species of Triatominae. Mem Inst Oswaldo Cruz. 2002;97:1121–5.
Marcondes B, Leuch-Lozovei C, Galati A, Eunice AB, Taniguchi H. The usefulness of Bergmannʼs rule for the distinction of members of Lutzomyia intermedia species complex (Diptera, Psychodidae, Phlebotominae). Mem Inst Oswaldo Cruz. 1998;93:363–4.
Brehm G, Fiedler K. Bergmann’s rule does not apply to geometrid moths along an elevational gradient in an Andean montane rain forest. Glob Ecol Biogeogr. 2004;13:7–14.
Hernández ML, Abrahan L, Dujardin JP, Gorla D, Catalá S. Phenotypic variability and population structure of peridomestic Triatoma infestans in rural areas of the arid Chaco (western Argentina): spatial influence of macro and microhabitats. Vector Borne Zoonotic Dis. 2011;11:503–13.
Barnabé C, Grijalva MJ, Santillán-Guayasamín S, Yumiseva CA, Waleckx E, Brenière SF, Villacís AG. Genetic data support speciation between Panstrongylus howardi and Panstrongylus chinai, vectors of Chagas disease in Ecuador. Infect Genet Evol. 2020;78:104103.
Special thanks to the inhabitants of the visited communities and the personnel of the National Chagas Control Program of the Ministry of Health of Ecuador who participated in the collection of the triatomines. Technical advice was provided by Miryam Rivera, María José Navarrete and Omar Torres-Carvajal. Technical assistance was provided by Ignacio Pacheco. Chromosomal analyses were supported by the project No. 160 “Comisión Sectorial de Investigación Científica” (CSIC, Udelar, Uruguay). English editing was provided by Lori Lammert.
Financial support was received from the Pontifical Catholic University of Ecuador (K13063 and K13023); the National Institutes of Health, Fogarty International Center, Global Infectious Disease Training Grant (1D43TW008261-01A1) and the National Institute of Allergy and Infectious Diseases, Division of Microbiology and Infectious Diseases, Academic Research Enhancement Award (1R15AI077896-01); Children’s Heartlink; the European Commission Framework Programme 7 Project “Comparative epidemiology of genetic lineages of Trypanosoma cruzi” ChagasEpiNet, Contract No. 223034, and Institut de Recherche pour le Développement (IRD).
Ethics approval and consent to participate
This study was performed within the guidelines established by the American Association for Laboratory Animal Science (IACUC), using protocol 15-H-034-IACUC, title “Life-cycle, feeding and defecation patterns of kissing bugs (Hemiptera: Reduviidae: Triatominae) under laboratory conditions”, approved by Ohio University Institutional Animal Care and Use Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Geographical location and altitude of studied communities in Loja and Manabí provinces, Ecuador.
Genetic 2k-p distances among Panstrongylus chinai (26 individuals), P. howardi (33 individuals), P. megistus (2 individuals), P. tupynambai (2 individuals), P. lutzi (1 individual), Triatoma tibiamaculata (2 individuals), T. infestans (1 individual) and T. brasiliensis (1 individual).
About this article
Cite this article
Villacís, A.G., Dujardin, JP., Panzera, F. et al. Chagas vectors Panstrongylus chinai (Del Ponte, 1929) and Panstrongylus howardi (Neiva, 1911): chromatic forms or true species?. Parasites Vectors 13, 226 (2020). https://doi.org/10.1186/s13071-020-04097-z