Skip to main content

Phenotypic variability in traits related to flight dispersal in the wing dimorphic species Triatoma guasayana



Triatoma guasayana is considered an emerging vector of Chagas disease in the Southern Cone of South America. The presence of a triatomine population with brachypterous individuals, in which both wings are reduced, has recently been reported for this species. The aim of the present study was to determine if flight-related traits varied across populations, if these traits could explain differences in flight capacity across populations and if flight-related traits are associated with geographic and/or climatic variation.


The study involved 66 male T. guasayana specimens from 10 triatomine populations. Digital images of wing, head and pronotum were used to estimate linear and geometric morphometric variables. Variations in size and shape were analysed using one-way analysis of variance and canonical variate analysis (CVA), respectively. Mantel tests were applied to analyse the relationship between morphometric and geographic distances, and the association between size measurements was analysed using Pearson’s correlation. We explored covariation between size and shape variables using partial least square analyses (PLS). The association of geographic and climatic variables with size measurements was tested using linear regression analyses. We performed PLS analyses for shape measurements.


Wing size differed significantly across triatomine populations. The CVA showed that wing shape of the brachypterous population is well discriminated from that of the other populations. The Mantel test showed a positive and significant association between wing shape and geographic distances. The heads of the brachypterous population were significantly larger than those of the other populations. Similar to wing shape, the head shape of the brachypterous population was well discriminated from those of the other populations. Pronotum width did not show significant differences across populations. Geographic and climatic factors were associated with size and shape of both the wing and head, but not with pronotum width.


Most of the traits related to flight dispersal varied across populations. Wing shape and head shape were found to be better markers for differentiated morphological variation across populations. Head measurements also varied in accordance with this condition. Geographic and climatic variables were associated with most of the flight-related traits.

Graphical Abstract


In many organisms, environmental heterogeneity and ecological gradients can generate phenotypic variation [1]. Phenotypic patterns of morphological variations within insect species might reflect different dispersal patterns and habitat availability. Intraspecific morphometric analyses should focus on samples from a broad geographical range to understand variations in morphometry of individuals that occur in diverse environments. Additionally, to understand how environmental heterogeneity contributes to changes in phenotypic variation, knowledge on whether some morpho-functional mechanisms generate this pattern of variation is desirable [2].

The subfamily Triatominae (Hemiptera: Reduviidae) includes hematophagous insects that act as vectors of Chagas disease, one of the most important parasitic diseases in Latin America [3]. Currently, there are 157 living species [4] of triatomines, which differ in their epidemiological importance and in aspects of their biology [5]. In recent years, several studies have reported the occurrence of sylvatic triatomines dispersing actively to domestic environments in Argentina, Brazil and Paraguay [6,7,8]. Although sylvatic species are considered to be of secondary epidemiological importance, they are responsible for maintaining the sylvatic transmission cycle of Trypanosoma cruzi [9].

Triatoma guasayana is a sylvatic species that invades peridomestic and domestic habitats. It is distributed in the dry western Chaco region of Argentina, Bolivia and Paraguay. Some laboratory and field studies have shown that this triatomine species has vectorial capacity [8, 10,11,12,13]. In sylvatic ecotopes, T. guasayana has been found in dry cacti, bromeliads and fallen dry logs, and in peridomestic habitats, it is associated with chicken coops, goat or sheep corrals, piled materials and orchard fences [14,15,16]. Triatoma guasayana is traditionally considered to be part of the “sordida subcomplex” [17]. However, phylogenetic analysis indicated that T. guasayana is more closely related to the “rubrovaria subcomplex” than to members of the “sordida subcomplex” [18, 19]. According to Gorla and Noireau [20], T. guasayana, among other species, may be considered an emerging vector in the Southern Cone of South America.

It has been postulated that flight dispersal is the main mechanism of invasion, infestation and/or recolonization of houses by adult triatomines [21, 22]. Nutritional and reproductive status as well as population density and environmental conditions are known factors that modulate flight dispersal in member of the Triatominae [23,24,25,26]. There are interspecific differences in terms of dispersal capacity by flight; for example, the number of flying individuals was found to be higher in T. guasayana than in T. infestans and other secondary triatomines [12]. Flight dispersal may vary markedly even within the same species. Some species present polymorphisms in wing length and/or in the development of the flight muscles that affect the active dispersal capacity [27, 28]. Wing polymorphism implies variation in the size of the wings or in their presence/absence [27]. The diversity of wing morphotypes is usually associated with other morphological attributes, such as head shape, thorax size and olfactory capacity. In Mepraia spinolai, a species with wing dimorphism, macropterous individuals (with normal wing size) exhibited a greater development of the thorax than micropterous specimens [29]. For T. guasayana, the presence of a population with brachypterous individuals, in which both wings are reduced, has recently been reported [30]. In this study, changes in the shape of the head according to wing characteristics were also determined, with macropterous individuals found to have a shorter head and a greater distance between eyes [30]. These observations in M. spinolai and T. guasayana suggest that the varieties of wing morphotypes imply a diversity of body structures. Polymorphism in flight muscles has been reported and includes developmental changes that affect dispersal capacity [31]. In some triatomine species, a relationship between wing size and thorax development has been demonstrated. Hybrids between Triatoma sherlocki and Triatoma juazeirensis have an increased dispersal capacity that coincides with a wider thorax, despite having wings of intermediate size. A wider thorax may imply greater development of the wing muscles. Mepraia spinolai showed an increase in thorax size as a function of wing development [29].

Taking into account the existence of wing dimorphism in T. guasayana and that this condition determines flight functionality, the aim of this study was to understand: (i) if traits related to dispersal characteristics varied across populations; (ii) if size and/or shape variation of flight-related traits could explain differences in flight capacity across populations; and (iii) if flight-related traits are associated with geographic and/or climatic variation.



The study involved 66 males from 10 populations of T. guasayana (Table 1; Fig. 1). Only males were studied because the sample included more males than females, and because other studies suggest sexual dimorphism ([30] and J Nattero, unpublished data). For the Cochabamba population, from the 16 collected males, wings of five individuals could be included since the remaining individuals had not been properly preserved. All populations included in this study were collected from peridomiciliary structures. All adults from the Cochabamba population are brachypterous, a condition first reported by Hernández et al. [30] (Fig. 2a). For this population, the dispersal capacity in the field was evaluated with light traps, active searching and live-bait mouse traps (also known as Noireau traps). Although individuals were captured by active searching and with Norieau traps, no individuals were captured with light traps (J Espinoza, personal communication). A number of assays were carried out in the laboratory to verify the mode of dispersal; only dispersal by walking was observed (J Espinoza, personal communication). Atrophy of the alar muscles in the thoracic box was also observed in a few inspected individuals (ML Hernández, personal communication). We took digital images of the dorsal view of head, pronotum and right wing using digital cameras (Lumix DMC-ZS7, Panasonic Corporation, Kadoma, Osaka, Japan; and Nikon S9900, Nikon Corporation, Tokyo, Japan) attached to a stereomicroscope (Stemi SV-11; Carl Zeiss, Wetzlar, Germany) at 6× magnification. All images included a reference scale.

Table 1 Geographical location, coordinates and individual number of collected Triatoma guasayana individuals at each collection site
Fig. 1
figure 1

a Distribution area of Triatoma guasayana based on [63]. b Location of the study populations

Fig. 2
figure 2

a Specimens of macropterous (left) and brachypterous (right) Triatoma guasayana. be Linear measurements of the head (b) and pronotum (c) and landmark position for the head (d) and wing (e) in the T. guasayana populations included in this study. Head: AD, Anteocular distance; EW, maximum width across the eyes; MTW, maximum thorax width; PW, preocular width; WAT, width between antenniferous tubercle (WAT)

Geographic and climatic characterization

To obtain the geographic location for each population, we recorded latitude, longitude and altitude. The geographical coordinates of the 10 populations were imported into the Google Earth Engine platform to extract the climatic variables for each site. We used ERA5 (Copernicus Climate Data Store;!/home) to obtain the climate data with a spatial resolution of 11 × 11 km. Of the 50 climate variables available, in this work we only included data from annual rainfall and annual mean temperature. We also considered the average monthly dew point temperature and the average monthly temperature and estimated the relative humidity (RH), including these data as a third variable to characterize climate variation. RH was calculated as [32]:

$$RH =100 -5 x (T -Td)$$

where RH indicates relative humidity, T represents the average monthly temperature and Td is the average monthly dew point temperature. We only included these three environmental variables because they are determinants of flight dispersal in insects [33, 34].

Linear morphometry

Linear measurements of the head and pronotum were made using the free software package tpsDig2 version 2.31 ( The linear measurements of the heads included anteocular distance (AD), width between antenniferous tubercles (WAT), preocular width (PW) and maximum width across the eyes (EW) (Fig. 2b). For the pronotum, we measured the maximum thorax width (MTW) (Fig. 2c). All linear measurements showed a normal distribution.

Geometric morphometry

A landmark-based approach was applied to study geometric morphometrics of the heads and right forewings. We defined and collected five coplanar type II landmarks of the ventral view of the head, and eight type I landmarks of the wings using the software package tpsdig2 (Fig. 2d, e). We extracted data on head and wing shape with a generalized full Procrustes fit and a projection to shape tangent space [35]. We used Procrustes coordinates as shape variables. Centroid size (CS; i.e. the square root of the sum of squared distances from each landmark to the centroid of the configuration) was computed as a measure of wing size and head. All of these steps of the morphometric analysis were performed using the free software MorphoJ version 1.07a [36].

Statistical analysis

Linear and geometric variables of the head, pronotum and wings were analysed separately. Variations across populations for the single linear measurements of the wing, head and pronotum were explored via one-way analysis of variance (ANOVA) and Tukey’s post-hoc tests for each variable, using the software InfoStat version 2016 [37].

For shape measurements of the head and wings, we calculated Procrustes distances between pairs of populations and evaluated their significance via a non-parametric test based on permutations (1000 runs), using MorphoJ [36]. We then performed canonical variate analyses (CVA) and calculated the percentage of phenotypic similarity between pairs of populations using the cross-check test of discriminant analysis [38] in InfoStat software [37]. The Procrustes distances were represented in unrooted neighbour-joining (NJ) trees using the free software MEGA X 10.2.6 [39]. We calculated bootstrap values for 1000 replications following Ascarranuz et al. [40]. The relationship of geographic distances with linear and geometric morphometric measurements of head, pronotum and wing was analysed using Mantel tests with the free software PASSAGE 2 [41]. To construct the morphological distance matrices between populations, we used Procrustes distances for shape measurements and Euclidean distances for size measurements.

The association between wing, head and pronotum size was explored via a Pearson’s correlation, using InfoStat software [37]. We analysed the covariation between wing, head and pronotum shape and between size and shape using partial least square analyses (PLS) in MorphoJ [36].

To summarize the variation and describe the environmental characterization of the six geographic and climatic factors (i.e. latitude, longitude, altitude, annual rainfall, annual mean temperature and RH), we used a principal component analysis (PCA) with R package FactorMineR 2.4 [43]. The first two principal component (PC) scores described 95.04% of the total variance (59.72–35.69% for PC1 and PC2, respectively). In PC1, the factors with relevant coefficients (> 0.7) were latitude, altitude, temperature, precipitation and RH (Additional file 1: Table S1). In PC2, only longitude was an important factor (Additional file 1: Table S1). We explored the association between PC1 and PC2 and measurements of the head, wing and pronotum size across the 10 populations studied using a multiple linear regression analysis [43]. The covariation between geographic and climate variables and wing and head shape variation was explored using PLS.


Changes in wing size across populations

Significant differences in wing size were observed among populations (F(9, 42) = 2.99, P = 0.008) (Fig. 3). Tukey’s post-hoc tests showed significant differences between the Loreto population and two other populations: Cochabamba, the brachypterous population, and Independencia (P < 0.05) (Fig. 3). The Loreto population presented a smaller wing size than the Cochabamba and Independencia populations, whereas the latter two had the largest wings of the 10 populations. The Mantel test showed a negative and significant association between wing size and geographic distances (Z = 73,661.06, r = − 0.05, P = 0.04), indicating that the larger the differences in wing size, the smaller the geographic distance.

Fig. 3
figure 3

Box plot (median and standard deviation) of wing centroid size across populations of T. guasayana. Solid horizontal line represents the median value; dashed horizontal line, the median value; the whiskers, the standard deviation. Different uppercase letters above boxes indicate statistical difference at P < 0.05 (Tukey’s post-hoc tests)

Wing shape of the brachypterous population differs from that of the other populations

Wing shape differed among populations. The first two discriminant factors of CVA (CV1, CV2) explained 68.13% of the total variation (CV1: 42.60, CV2: 25.53%). The brachypterous population, Cochabamba, was found to be well discriminated from the other populations in the space of the first discriminant factor (CV1) (Fig. 4a). All individuals from this population, as well as those of the Campo Alegre, Cruz del Eje, Ischilín and Valle Fértil populations, were correctly assigned (0% of misclassified individuals) (Additional file 1: Table S2). In contrast, the San Martín population exhibited 60% of misclassified individuals (Additional file 1: Table S2). Procrustes distances between populations were significantly different (P < 0.05) for all of the comparisons that included the Cruz del Eje and Cochabamba populations (Additional file 1: Table S3). The NJ tree based on Procrustes distances is presented in Fig. 5a, which shows that the only bootstrap-supported clade was the one that included the Independencia and Cruz del Eje populations and the Ischilín and Campo Alegre populations (Fig. 5a). The Mantel test showed a positive and significant association between wing shape and geographic distances (Z = 2483.882, r = 0.44, P = 0.001), indicating that the smaller the differences in wing shape, the smaller the geographic distances.

Fig. 4
figure 4

Factorial maps in the plane of the two first axes of the canonical variate analysis for shape measurements of wings and heads for populations of T. guasayana. a Wings, b heads

Fig. 5
figure 5

Neighbour-joining trees based on Procrustes distances of shape measurements for wings (a) and heads (b) for populations of T. guasayana. Numbers close to the nodes are 1000 replicates of bootstrap values

Head size of the brachypterous population is larger than that of the other populations

The ANOVAs performed for the four head linear measurements showed significant differences in all cases (Table 2). Tukey’s post-hoc tests showed that the brachypterous population had the highest values for all four measurements, although these were significantly different from those of the other nine populations only for the anteocular distance (P < 0.05) (Table 2). In the Mantel test, head size showed a positive and significant correlation with geographic distances (Z = 83,069.24, r = 0.324, P = 0.001,) indicating that the smaller the differences in head size, the smaller the geographic distance.

Table 2 Mean values and analysis of variance results between populations for linear head measurements of T. guasayana

Head shape of the brachypterous population is different from that of the other populations

The first two discriminant factors accumulated 96.72% of the total variation (90.18 and 6.62% for the first and second axes, respectively). The space delimited by the first two discriminant factors showed three well-delimited groups, of which one included only the Cochabamba population, the brachypterous population, one included the Independencia, San Martín and Valle Fértil populations and the third included the remaining populations (Fig. 4b). All individuals of the Cochabamba and Ischilín populations were correctly assigned (0% of misclassified individuals) (Additional file 1: Table S2), while The population from 12 de Junio exhibited the highest number of misclassified individuals (85.71%) (Additional file 1: Table S2). Procrustes’ distances between populations were significantly different (P < 0.05) for all pairs of populations that included Cochabamba (Additional file 1: Table S3). The NJ tree based on Procrustes distances is presented in Fig. 5b. Bootstrap values supported different groups from different geographic areas; in particular, the brachypterous population was in a group that was separated from the rest (Fig. 5b). The Mantel test showed a positive and significant association between head shape and geographic distances (Z = 4499.254, r = 0.19, P = 0.005), indicating that the smaller the differences in head shape, the smaller the geographic distance.

Pronotum width is similar across populations

The ANOVA for the pronotum width did not show significant differences across populations (F(9, 51) = 1.78, P = 0.094). The Mantel test showed a negative and significant correlation with geographic distances (Z = 71,936.392, r = − 0.147, P = 0.024), indicating that the larger the differences in pronotum width, the smaller the geographic distance.

Wing, head and pronotum measurements are associated

Pearson’s correlations showed positive and significant associations between wing, head and pronotum size measurements (Table 3). The PLS analyses showed no covariation between wing and head shape, wing shape and head size, or between wing shape and pronotum width (Table 3). Head shape showed significant covariation with pronotum width and wing size (Table 3).

Table 3 Results of Pearson’s correlation between linear measurements of wing, head and pronotum, and partial least square analyses between shape measurements of wing and head and between size and shape measurements of wings, head and pronotum in T. guasayana populations

Linear and shape measurements are associated with geographic and climatic variables

Multiple linear regression analysis showed that head linear measurements were significantly associated with geographic and environmental factors with both PC1 and PC2 coefficient (with the exception of width between antenniferous tubercles with PC2) (Table 4). In contrast, pronotum width and wing size were not associated with geographic or climatic factors (Table 4). PLS analysis showed significant covariation of head and wing shape with all geographic and climatic variables (Table 5).

Table 4 Results of multiple linear regression analysis of morphometric traits and geographic and climatic variables in T. guasayana populations
Table 5 Results of the partial least square analyses of wing shape and head shape with geographic and climatic variables in T. guasayana populations


Triatoma guasayana exhibits a good dispersal capacity. In a study comparing this species with other triatominae species (T. infestans, T. platensis and T. eratyrusiformis) from the arid Chaco of Argentina, T. guasayana was found to have the highest average number of flying individuals [12]. Adult T. guasayana bugs may invade a wide variety of ecotopes and have a high potential for night flight dispersal during the hot dry months (August-November) [21, 44]. Since dispersal is linked to flight, the existence of wing dimorphism at the population level provides a suitable model to study not only evolutionary patterns, such as dispersal processes and population structure [45, 46], but also a range of ecological factors, including habitat structure and climate [47, 48].

Our results showed differences in wing size across T. guasayana populations. Although all triatomine males collected from Cochabamba showed a brachypterous condition, individuals from this population as well as those of the Independencia population exhibited a bigger wing size than the macropterous populations. Individuals from the brachypterous population had significantly greater body size than individuals from the two macropterous populations of Ischilín and Cruz del Eje (F(2,11) = 6.76, P = 0.012). This result suggests that the individuals from the Cochabamba population are larger and that their body is comparatively larger than their wings, showing the presence of brachyptery in this population. The same tendency has been observed in other insects, such as in carabids, with brachypterous individuals exhibiting a larger body size than their macropterous counterparts. This variation is attributed to differences in dispersal strategies, taking into account that brachypterous individuals disperse by running [49].

Regarding wing shape, our results showed differentiation between the brachypterous population and the other nine populations included in this study. This result suggests that differences in wing shape are associated with the brachyptery condition, and that this pattern could be influenced by the differences in flight capacity across populations. Differences in the shape but not in the size of a structure are considered to be determined by genetic processes [50]. Our results also showed that there is no covariation between wing shape and size measurements of the wing, head and pronotum. These results suggest a lack of functional integration (i.e. coordination among parts) of wing shape and size measurements across populations, and may be related to wing dimorphism. Studies combining biomechanics of flight and morphometric approaches could contribute to a better understanding of the functional implications of morphological changes associated with wing dimorphism [61].

Insects use a combination of sensory inputs, including cephalic sensilla, ocelli, antennal sensilla and eyes, to control flight activity [52,53,54]. Our results show that, in general terms, the brachypterous population presents a larger head size than the remaining populations. For T. guasayana, anterocular distance and distance between eyes have been suggested as the differences in the head that could be associated with flight capacity; it has also been proposed that shorter heads are related to changes in the position of the ocelli [30]. We found significantly longer head measurements for the brachypterous population, which agrees with the association suggested by [30]. Regarding head shape, our results showed that the head shape of the brachypterous population differed from the that of the other nine populations. The NJ tree for Procrustes distances showed similarities between populations that are geographically close.

Flight muscles are located in the pronotum. Our results showed that pronotum width was the only measured trait that did not differ between populations. There are alternative, non-mutually exclusive explanations that might account for this result. First, among insects, there is muscle polymorphism associated with the degree of muscle development. A greater pronotum width does not necessarily indicate a greater development of these muscles, as has been reported for other insect species [55]. Alary polymorphism is not necessarily accompanied with flight muscle development [30]. Second, there are functional or developmental constraints that might limit the potential for evolutionary change of pronotum width and reduce the potential of the population to respond to selection or to undergo non-adaptive evolution by drift [56].

The association of morphological distance with geographic distances between populations showed contrasting results. While wing shape and head size and shape exhibited positive associations, wing and pronotum size showed negative associations. Thus, these results might imply geographic separation as a key factor in morphological differentiation of genetic exchange between populations. A positive association would imply that geographically closer populations are more similar to each other. It would also indicate morphological spatial structuring for wings and heads between the populations of this species. This positive spatial association would adjust to a model of isolation by distance, in which differentiation increases with increasing distance. By contrast, the sizes of the wing and pronotum showed negative associations, implying that the greater the distance between populations, the smaller the difference in wing and pronotum size, respectively. These results should be interpreted with caution. In some cases, the correlation values were low, such as those for wing size and shape and for pronotum width.

The association of geographic and climatic variables with size and shape measurements showed that the former variables better explain head size, head shape and wing shape than the other measurements. For wing size and pronotum width, only longitude showed association. This means that, considering geographic variables, all morphological measures followed an east–west variation, and that head size, head shape and wing shape additionally followed a north–south variation and an altitudinal variation. Considering the broad-scale geographic variation of the sampled populations, variation in geographic variables also determines variation in climatic conditions. Wing and head shape exhibit significant phenotypic plasticity in response to developmental temperature, and this plastic response (hot-to-cold) mirrors a pattern of geographic differentiation, as might be expected under a synergistic scenario of adaptive plasticity and adaptive genetic divergence [57, 58]. Similar results were found in the yellow dung fly Scathophaga stercoraria (Diptera: Scathophagidae) from the New World, with clines of variation associated with wing shape but not with wing size [59]. There are also examples in Drosophila flies of no association between wing size and temperature variation (13 °C and 22 °C) [60].

In M. spinolai (sub Triatoma spinolai [61]) collected from various sites in northern Chile, adults from coastal populations were always wingless, while inland populations showed both wing and wingless adults. Winged males would be able to disperse, although this event does not guarantee that fliers will land in favourable sites for feeding and/or mating [61]. Brachyptery in T. guasayana is a heritable trait for at least three laboratory generations [35]. In other Hemiptera, the genetic base of brachytherapy has been found to be polygenic [62]. The morphological characteristics of the heads and wings observed in the Cochabamba population suggest a lack of dispersal by flight; therefore, the population would be genetically isolated. The differentiation observed in all of the characters would support this assumption, especially for those characters related to shape variation, which suggest genetic differentiation and, as mentioned above, this genetic differentiation between populations could be indicating different modes of dispersion. On the other hand, our results show that both head shape and wing shape correlate with geographic and climatic variables. Other insect species, such as Drosophila subobscura, or the fly Scathophaga stercoraria, also show a geographic pattern of morphological differentiation, including significant clinal variation in wing shape [52, 60].


This study indicates that wing shape and head shape are the better markers for differentiated morphological variation across populations. Head measurements also varied across populations. Pronotum size is not necessarily a good indicator of muscle development and was found not to be directly related to body size in this species. Future studies should include genetic variation across populations and explore whether the brachypterous population indeed has limited genetic exchange.

Availability of data and materials

The datasets supporting the conclusions of this article are included in the article and its additional files. Raw data are available from the corresponding author on reasonable request.


  1. West-Eberhard MJ. Developmental plasticity and evolution. Oxford: Oxford University Press; 2003.

    Book  Google Scholar 

  2. Thompson DB. Different spatial scales of natural selection and gene flow: the evolution of behavioral geographic variation and phenotypic plasticity. In: Foster S, Endler J, editors. Geographic diversification of behavior: an evolutionary perspective. Oxford: Oxford University Press; 1999. p. 33–51.

    Google Scholar 

  3. WHO. Chagas disease (also known as American trypanosomiasis). 2021. Accessed 27 May 2022.

  4. de Paiva VF, Belintani T, de Oliveira J, Galvão C, da Rosa JA. A review of the taxonomy and biology of Triatominae subspecies (Hemiptera: Reduviidae). Parasitol Res. 2022;121:499–512.

    Article  Google Scholar 

  5. Schofield CJ, Galvão C. Classification, evolution and species groups within the Triatominae. Acta Trop. 2009;110:88–100.

    Article  CAS  Google Scholar 

  6. Cavallo MJ, Amelotti I, Gorla DE. Invasion of rural houses by wild Triatominae in the arid Chaco. J Vector Ecol. 2016;41:97–102.

    Article  Google Scholar 

  7. Brito RN, Gorla DE, Diotaiuti L, Gomes ACF, Souza RCM, Abad-Franch F. Drivers of house invasion by sylvatic Chagas disease vectors in the Amazon-Cerrado transition: a multi-year, state-wide assessment of municipality-aggregated surveillance data. PLoS Negl Trop Dis. 2017;11:e0006035.

    Article  Google Scholar 

  8. Cardozo M, Fiad FG, Crocco LB, Gorla DE. Triatominae of the semi-arid Chaco in Central Argentina. Acta Trop. 2021;224:106158.

    Article  Google Scholar 

  9. Wisnivesky-Colli C, Schweigmann NJ, Alberti A, Pietrokovsky SM, Conti O, Montoya S, et al. Sylvatic American trypanosomiasis in Argentina. Trypanosoma cruzi infection in mammals from the Chaco forest in Santiago del Estero. Trans R Soc Trop Med Hyg. 1992;86:38–41.

    Article  CAS  Google Scholar 

  10. Rodríguez CS, Crocco LB, Nattero J. Competencia vectorial de Triatoma guasayana (Hemiptera: Reduviidae): patrón de alimentación y excreción. Rev Soc Entomol Arg. 2004;63:11–6.

    Google Scholar 

  11. Loza-Murguía M, Noireau F. Vectorial capacity of Triatoma guasayana (Wygodzinsky & Abalos) (Hemiptera: Reduviidae) compared with two other species of epidemic importance. Neotrop Entomol. 2010;39:799–809.

    Article  Google Scholar 

  12. Abrahan LB, Gorla DE, Catalá SS. Dispersal of Triatoma infestans and other triatominae species in the arid Chaco of Argentina—Flying, walking or passive carriage? The importance of walking females. Mem Inst Oswaldo Cruz. 2011;106:232–9.

    Article  Google Scholar 

  13. Lucero DE, Ribera W, Pizarro JC, Plaza C, Gordon LW, Peña R, et al. Sources of blood meals of sylvatic Triatoma guasayana near Zurima, Bolivia, assayed with qPCR and 12S cloning. PLoS Negl Trop Dis. 2014;8:e3365.

    Article  Google Scholar 

  14. Carcavallo RU, Canale DM, Martínez A. Habitats de triatominos argentinos y zonas ecológicas donde prevalecen. Chagas. 1988;5:8–17.

    Google Scholar 

  15. Canale DM, Cecere MC, Chuit R, Gurtler RE. Peridomestic distribution of Triatoma garciabesi and Triatoma guasayana in northwest Argentina. Med Vet Entomol. 2000;14:383–90.

    Article  CAS  Google Scholar 

  16. Rodríguez-Planes LI, Vazquez-Prokopec GM, Cecere MC, Canale DM, Gürtler RE. Selective insecticide applications directed against Triatoma infestans (Hemiptera: Reduviidae) affected a nontarget secondary vector of chagas disease Triatoma garciabesi. J Med Entomol. 2016;53:144–51.

    Article  Google Scholar 

  17. Dujardin JP, Schofield CJ, Panzera F. Les vecteurs de la maladie de Chagas: recherches taxonomiques, biologiques et génétiques. 2000. Brussels: Académie royale des sciences d'outre-mer.

  18. Almeida CE, Marcet PL, Gumiel M, Takiya DM, Cardozo-de-Almeida M, Pacheco RS, et al. Phylogenetic and phenotypic relationships among Triatoma carcavalloi (Hemiptera: Reduviidae: Triatominae) and related species collected in domiciles in Rio Grande do Sul State Brazil. J Vector Ecol. 2009;34:164–73.

    Article  Google Scholar 

  19. Justi SA, Russo CA, Mallet JRDS, Obara MT, Galvão C. Molecular phylogeny of Triatomini (Hemiptera: Reduviidae: Triatominae). Parasit Vectors. 2014;7:1–12.

    Article  Google Scholar 

  20. Gorla D, Noireau F. Geographic distribution of Triatominae vectors in America. In: Telleria J, Tibayrenc M, editors. American trypanosomiasis Chagas disease. Amsterdam: Elsevier; 2017. p. 197–221.

    Chapter  Google Scholar 

  21. Wisnivesky-Colli C, Gürtler RE, Solarz ND, Schweigmann NJ, Pietrokovsky SM, Alberti A, et al. Dispersive flight and house invasion by Triatoma guasayana and Triatoma sordida in Argentina. Mem Inst Oswaldo Cruz. 1993;88:27–32.

    Article  CAS  Google Scholar 

  22. Noireau F, Dujardin JP. Flight and nutritional status of sylvatic Triatoma sordida and Triatoma guasayana. Mem Inst Oswaldo Cruz. 2001;96:385–9.

    Article  CAS  Google Scholar 

  23. Lehane MJ, McEwen PK, Whitaker CJ, Schofield CJ. The role of temperature and nutritional status in flight initiation by Triatoma infestans. Acta Trop. 1992;52:27–38.

    Article  CAS  Google Scholar 

  24. Schofield CJ, Lehane M, McEwen P, Catalá SS, Gorla DE. Dispersive flight by Triatoma infestans under natural climatic conditions in Argentina. Med Vet Entomol. 1992;6:51–6.

    Article  CAS  Google Scholar 

  25. McEwen P, Lehane M. Factors influencing flight initiation in the triatomine bug Triatoma sordida (Hemiptera: Reduviidae). Insect Sci Applic. 1993;14:461–4.

    Article  Google Scholar 

  26. McEwen P, Lehane M, Whitaker C. The effect of adult population density on flight initiation in Triatoma infestans (Klug) (Hemiptera: Reduviidae). J Appl Entomol. 1993;116:321–5.

    Article  Google Scholar 

  27. Harrison RG. Dispersal polymorphisms in insects. Annu Rev Ecol Syst. 1980;11:95–118.

    Article  Google Scholar 

  28. Gurevitz JM, Kitron U, Gürtler RE. Flight muscle dimorphism and heterogeneity in flight initiation of field-collected Triatoma infestans (Hemiptera: Reduviidae). J Med Entomol. 2007;44:186–91.

    Article  Google Scholar 

  29. Hernández ML, Dujardin JP, Gorla DE, Catalá SS. Can body traits, other than wings, reflect the flight ability of Triatominae bugs? Rev Soc Bra Med Trop. 2015;48:682–91.

    Article  Google Scholar 

  30. Hernández ML, Espinoza J, Gomez M, Gorla D. Morphological changes associated with brachypterous Triatoma guasayana (Hemiptera, Reduviidae) and their relationship with flight. Int J Trop Insect Sci. 2020;40:413–21.

    Article  Google Scholar 

  31. Almeida CE, Oliveira HL, Correia N, Dornak LL, Gumiel M, Neiva VL, et al. Dispersion capacity of Triatoma sherlocki, Triatoma juazeirensis and laboratory-bred hybrids. Acta Trop. 2012;122:71–9.

    Article  Google Scholar 

  32. Lawrence MG. The relationship between relative humidity and the dewpoint temperature in moist air: a simple conversion and applications. Bull Am Meteorol Soc. 2005;86:225–34.

    Article  Google Scholar 

  33. Johnson CG. Migration and dispersal of insects by flight. London: Methuen; 1969.

    Google Scholar 

  34. Ekkens DB. Nocturnal flights of Triatoma (Hemiptera: Reduviidae) in Sabino Canyon, Arizona: I light collections. J Med Entomol. 1981;18:211–27.

    Article  Google Scholar 

  35. Dryden IL, Mardia KV. Statistical shape analysis. Chichester: Wiley; 1998.

    Google Scholar 

  36. Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011;11:353–7.

    Article  Google Scholar 

  37. Di Rienzo JA, Casanoves F, Balzarini MG, Gonzalez L, Tablada M, Robledo CW. InfoStat versión 2016 Universidad Nacional de Córdoba. Córdoba: FCA; 2016.

    Google Scholar 

  38. Lachenbruch PA. An almost unbiased method of obtaining confidence intervals for the probability of misclassification in discriminant analysis. Biometrics. 1967;23:639–45.

    Article  CAS  Google Scholar 

  39. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.

    Article  CAS  Google Scholar 

  40. Ascarrunz E, Claude J, Joyce WG. Estimating the phylogeny of geoemydid turtles (Cryptodira) from landmark data: an assessment of different methods. PeerJ. 2019;7:e7476.

    Article  Google Scholar 

  41. Rosenberg MS, Anderson CD. PASSaGE: pattern analysis, spatial statistics, and geographic exegesis. version 2. Methods Ecol Evol. 2011;2:229–32.

    Article  Google Scholar 

  42. Le S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.

    Article  Google Scholar 

  43. Bai Y, Dong JJ, Guan DL, Xie JY, Xu SQ. Geographic variation in wing size and shape of the grasshopper Trilophidia annulata (Orthoptera: Oedipodidae): morphological trait variations follow an ecogeographical rule. Sci Rep. 2016;6:1–15.

    Article  CAS  Google Scholar 

  44. Noireau F, Flores R, Vargas F. Trapping sylvatic Triatominae (Reduviidae) in hollow trees. Trans R Soc Trop Med Hyg. 1999;93:13–4.

    Article  CAS  Google Scholar 

  45. Hill JK, Thomas CD, Blakeley DS. Evolution of flight morphology in a butterfly that has recently expanded its geographic range. Oecologia. 1999;121:165–70.

    Article  CAS  Google Scholar 

  46. Hughes CL, Dytham C, Hill JK. Modelling and analysing evolution of dispersal in populations at expanding range boundaries. Ecol Entomol. 2007;32:437–45.

    Article  Google Scholar 

  47. Norberg U, Leimar O. Spatial and temporal variation in flight morphology in the butterfly Melitaea cinxia (Lepidoptera: Nymphalidae). Biol J Lin Soc. 2002;77:445–53.

    Article  Google Scholar 

  48. Sekar S. A meta-analysis of the traits affecting dispersal ability in butterflies: can wingspan be used as a proxy? J Anim Ecol. 2012;81:174–84.

    Article  Google Scholar 

  49. Gutierrez D, Rosa M. Patterns in the distribution, abundance and body size of carabid beetles (Coleoptera: Carabidae) in relation to dispersal ability. J Biogeogr. 1997;24:903–14.

    Article  Google Scholar 

  50. Dujardin JP, Costa J, Bustamante D, Jaramillo N, Catalá SS. Deciphering morphology in Triatominae: the evolutionary signals. Acta Trop. 2009;110:101–11.

    Article  CAS  Google Scholar 

  51. O’Higgins P, Cobb SN, Fitton LC, Gröning F, Phillips R, Liu J, et al. Combining geometric morphometrics and functional simulation: an emerging toolkit for virtual functional analyses. J Anat. 2011;218:3–15.

    Article  Google Scholar 

  52. Taylor CP. Contribution of compound eyes and ocelli to steering of locusts in flight: I. behavioural analysis. J Exp Biol. 1981;93:1–18.

  53. Moser JC, Reeve JD, Bento JMS, Della Lucia TM, Cameron RS, Heck NM. Eye size and behaviour of day-and night-flying leaf cutting ant alates. J Zool. 2004;264:69–75.

    Article  Google Scholar 

  54. Taylor GK, Krapp HG. Sensory systems and flight stability: what do insects measure and why? Insect Physiol. 2007;34:231–316.

  55. Suarez-Tovar CM, Sarmiento C. Beyond the wing planform: morphological differentiation between migratory and nonmigratory dragonfly species. J Evol Biol. 2016;29:690–703.

    Article  CAS  Google Scholar 

  56. Klingenberg CP. Developmental constraints, modules and evolvability. In: Hallgrímsson B, Hall BK, editors. Variation: a central concept in biology. Amsterdam: Academic Press; 2005. p. 219–47.

  57. Ghalambor CK, McKay JK, Carroll SP, Reznick DN. Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007;21:394–407.

    Article  Google Scholar 

  58. Schmid M, Guillaume F. The role of phenotypic plasticity on population differentiation. Heredity. 2017;119:214–25.

    Article  CAS  Google Scholar 

  59. Schäfer MA, Berger D, Rohner PT, Kjaersgaard A, Bauerfeind SS, Guillaume F, et al. Geographic clines in wing morphology relate to colonization history in new world but not old world populations of yellow dung flies. Evolution. 2018;72:1629–44.

    Article  Google Scholar 

  60. Santos M, Brites D, Laayouni H. Thermal evolution of pre-adult life history traits, geometric size and shape, and developmental stability in Drosophila subobscura. J Evol Biol. 2006;19:2006–21.

    Article  CAS  Google Scholar 

  61. Schofield CJ, Apt W, Sagua H, Panzera F, Dujardin JP. Alary polymorphism in Triatoma spinolai and its possible relationship with demographic strategy. Med Vet Entomol. 1998;12:30–8.

    Article  CAS  Google Scholar 

  62. Roff DA. The evolution of wing dimorphism in insects. Evolution. 1986;40:1009–20.

    Article  Google Scholar 

  63. Carcavallo RU, Curto de Casas SI, Sherlock IA, Galíndez Girón I, Jurberg J, Galvão C, et al. Geographical distribution and alti-latitudinal dispesion. In: Lent H, Carcavallo RU, Galíndez Girón I, Jurberg J, editors. Atlas of Chagas disease vectors in the Americas. Rio de Janeiro: Fiocruz Editorial; 1999. p. 747–92.

Download references


We thank Miriam Cardozo for help in the collection of biological material used in this study, Pablo Milla Carmona for sharing a script for bootstrap estimation and Jorgelina Brasca for language editing. GVG, FGF, PSC are doctoral fellows of CONICET. MLH, ALCF, RVP, PL, LA and JN are members of the CONICET Researcher’s Career.


This study was supported by grants from Agencia Nacional de Promoción Científica y Tecnológica (PICT 2016-2527 and PICT 2019-1884), CONICET (PIP N°1220200102527CO) and SECYT (Resol. 411/2018).

Author information

Authors and Affiliations



Conceived and designed the experiments: JN and CSR. Collected the insects: FGF, MLH, PSC, ARA PL, LA, MBG, JE, FC. Took size and shape measurements: GVG, CSR, FGF, MLH, JN. Analysed the data: GVG, CSR, FGF, JN. Wrote the manuscript: GVG, CSR, FGF, ALCF, RVP, JN. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Julieta Nattero.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Eigenvalues and percentage of total variance explained by a principal component analysis of geographic and environmental factors from the T. guasayana populations studied. Table S2. Reclassification of T. guasayana individuals with the number and percentage of assigned individuals derived from discriminant function analyses performed for wings and heads. Table S3. Procrustes distances between pairs of T. guasayana populations performed for shape variables of wing and heads.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gigena, G.V., Rodríguez, C.S., Fiad, F.G. et al. Phenotypic variability in traits related to flight dispersal in the wing dimorphic species Triatoma guasayana. Parasites Vectors 16, 8 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: