Characterization of melanic and non-melanic forms in domestic and peridomestic populations of Triatoma infestans (Hemiptera: Reduviidae)
Parasites & Vectors volume 13, Article number: 47 (2020)
Melanic (dark) morphs have been barely reported in peridomestic and sylvatic conditions for Triatoma infestans, the most important vector of Chagas disease in the Southern Cone of South America. Adults with dark and small yellow markings on the connexivum were collected after manual searches conducted by technical personnel in 62 domiciliary units in Cruz del Eje, Córdoba Province, Argentina. The last community-wide insecticide spraying campaign before the study had been conducted three years earlier. We investigated if there was a measurable color morph variation (melanic and non-melanic) in wings and connexivum; we determined infestation, distribution of melanic and non-melanic forms, and correspondence of colorimetric variation with variations in morphology (wing size and shape and body length), development (wing fluctuating asymmetry), physiology (nutritional status) or behaviour (flight initiation).
Forty-nine females, 54 males and 217 nymphs were collected in 24 domiciliary units. House infestation and colonization were 53% and 47%, respectively. Most of the T. infestans individuals (83.2%) were collected in chicken coops; intradomicile infestation was recorded in only one case. The chromatic cluster analysis showed two well-defined groups: melanic and non-melanic. The melanic group included 17 (35%) females and 25 (46%) males. Peridomestic infestation was lower for melanic than for non-melanic adults. Melanic morphs were collected in houses from several localities. Sexual dimorphisms were confirmed by morphometric measurements. Body length was large in melanic adults (P < 0.01 only for males). Differences between groups were significant for wing size and shape, but not for weight or weight/body length ratio. Melanic females and males showed significantly higher fluctuating asymmetry (FA) indices than their non-melanic counterparts.
This is the second report of melanic forms of T. infestans in domestic and peridomestic habitats in the Dry Chaco region of Argentina. Although non-melanic adults exhibited a higher infestation rate, melanic adults were widespread in the area and were collected in the infested domicile and in most types of peridomestic annexes. Differences in morphometric variables between groups might be due to different ecological adaptations. The higher FA levels observed in melanic individuals suggest a higher developmental instability and a selective advantage of non-melanic individuals in domestic and peridomestic habitats.
Natural populations often exhibit great phenotypic variation. Differences in visible traits among organisms, such as body color, are the most conspicuous evidence of morphological variability in nature. One of the simplest and most common examples of such conspicuous variation is melanism, the occurrence of variants that are mostly or completely dark in pigmentation. This type of variation may involve discrete melanic and non-melanic phenotypes or continuously varying pigmentation . Melanism is very common in insects, with melanic or dark morphs exhibiting an unusually high concentration of cuticular melanin . Melanin pigments and their precursors are important structural and protective components of the cuticle . Diverse selective pressures, including visual and non-visual cues, might be associated with melanism, (reviewed in ). Theory posits that individuals with dark color patterns would be at an advantage under low ambient temperature and high UV radiation conditions and should be less tolerant to desiccation [4,5,6,7]. Melanism also has correlated and/or pleiotropic effects, including morphological, developmental, physiological, behavioural and/or neurological traits (reviewed in ).
Triatoma infestans (Hemiptera: Reduviidae: Triatominae) is the most important vector of Trypanosoma cruzi, the etiological agent of Chagas disease, in the Southern Cone of South America . Insecticide spraying reduced the geographical range and abundance of T. infestans but did not interrupt vector transmission of T. cruzi in the Gran Chaco , a 1.3 million-km2 ecoregion extending across Paraguay, Bolivia and Argentina, including the northwestern extreme of Córdoba Province, Argentina . The northwest of Córdoba Province has been historically an endemic Chagas disease region. The latest epidemiological reports from this disease in Argentina indicate Córdoba as one of the provinces with the highest risk of vector and congenital transmission of T. cruzi .
Although color pattern quantification has been performed in other insect species (e.g. the local mimicry polymorphism of Heliconius butterflies, ), in Triatominae species, color variation has not yet been quantified. Melanic (dark) morphs of T. infestans were reported for the first time in sylvatic foci of the Bolivian Chaco in an isolated dry forest very far from human settlements [13, 14]. Morphologically, these specimens were very similar to domestic T. infestans collected from the Bolivian Chaco, except for their overall darker coloration with small yellow markings on the connexivum . Melanic forms were then reported near chicken coops in the humid forest of northeastern Argentina (an area not included in the Gran Chaco) , in fallen trees among shrubs close to an indigenous community in the Paraguayan Chaco  and in hollow tree trunks harboring parrot nests in the Argentine Dry Chaco . Recently, domestic colonies of the dark morph were reported in the Chaco Province from Argentina . Several studies using different methods confirmed that individuals of the dark morph were chromatic variants of T. infestans differing in color, head and wing morphometry, antennal sensilla patterns, chromosome C-banding, nuclear rDNA sequences, genome size, and mitochondrial cytochrome B and CO1 genes [13, 17, 19,20,21,22,23,24].
Several T. infestans adults collected in domestic and peridomestic habitats in Cruz del Eje, Córdoba Province, Argentina, in 2012 displayed an evident visual variation in color, with some individuals being darker than others and having small yellow markings on the connexivum. To our knowledge, there are only one record of melanic morphs in domestic and peridomestic populations of T. infestans from the Gran Chaco of Argentina . One of the most important current challenges of interrupting Chagas disease transmission is to know the origin of reinfestant insects after the application of vector control actions. Recent studies have shown that T. infestans populations that remain after spraying campaigns are highly connected to sylvatic individuals that could be involved in restoring the reinfestation process . These data show the importance of sylvatic T. infestans populations in recolonization of treated areas. Considering that melanism could be the result of an adaptation to certain environmental conditions, melanic individuals might be more adapted than non-melanic individuals to live in human-made structures, such as peridomiciles, which are less isolated from the surrounding environment and more exposed to harsh weather.
The main goal of this study was to analyze this particular chromatic polymorphism of T. infestans, its relationship with infestation parameters and its correlation with other variables. Particularly, we aimed to (i) investigate if there is measurable variation in color morphs (melanic and non-melanic) in adults collected in domestic and peridomestic habitats; (ii) determine infestation, frequency and distribution of melanic and non-melanic forms; and (iii) determine if colorimetric variation corresponds to other variations in morphology (wing geometric morphometry, total body length and total body area), development (wing developmental instability), physiology (nutritional status) and behaviour (flight initiation). Wing shape analyses performed to assess population structuring in T. infestans showed patterns compatible with independent geographical origins [25, 26] or the presence of different groups or variants [21, 27]. A different nutritional status between melanic and non-melanic forms suggests dissimilar feeding capacity and ability [28, 29]. Variations in the total weight/total body length (W/L) ratio were related to different flight initiation probabilities [30,31,32,33]. Differences in wing developmental instability may indicate differences in ecological adaptation to the environmental conditions of the habitat where the insects develop [34, 35]. We expected to find differences in morphometric traits between melanic and non-melanic groups as proxies of differences in their ecological adaptations.
All the insects included in this study were collected in a rural area (approximately 60 × 50 km) of Cruz del Eje department, Córdoba Province, northwestern Argentina, in December 2012 (Fig. 1). This area is part of the Argentine Dry Chaco, characterized by a subtropical dry climate with a summer season from October to March. The last community-wide insecticide spraying campaign before this study had been carried out by vector control personnel three years earlier. The present cross-sectional study was conducted in domiciliary units, which consisted of a domicile with human resting places and the associated peridomicile, usually chicken coops, goat or pig corrals and, less frequently, rabbit hutches or storerooms. Personnel of the National Vector Control Programme inspected 62 houses for the presence of T. infestans. After inspection, technicians sprayed all the infested houses with cypermethrin.
Triatomines were searched for inside domiciles and peridomestic structures using the hour-person technique . For this purpose, one person searched for the presence of triatomines at each site within each house during 20 min. All collected insects (103 adults and 217 nymphs) were taken to the laboratory under refrigerated conditions and species, sex and stage were identified . All T. infestans adults were evaluated for the presence of T. cruzi in feces using a conventional optical microscope; none of them were positive for the presence of the parasite.
Weight (W) of adults was recorded on an electronic scale (precision ± 0.001 g; Mettler, Denver, USA) the day after collection. Dorsal view digital photographs of all adults were taken under a white background with a reference scale under the same illumination conditions and camera position. After insect dissection, right and left wings of each individual were mounted on graph paper. Digital photographs of wings were taken with a Moticam 2 (CMOS, Richmond, Canada) camera connected to a stereomicroscope (Stemi 2000-C; Zeiss, Oberkochen, Germany) under 6 × magnification.
Infestation by T. infestans
A house was recorded as colonized when at least one T. infestans nymph was found in the evaluated site (intradomicile or peridomicile structure). We considered infestation or colonization in the intradomicile (IDI, IDC) or in the peridomicile (PDI, PDC) when T. infestans adults and/or nymphs were recorded inside the domicile or in the peridomicile structures, respectively. House infestation was calculated as the percentage of houses infested with adults only, with adults and nymphs, or with nymphs only of the total evaluated houses. House colonization was calculated as the percentage of colonized houses of the total of evaluated houses. Infestation was also estimated for melanic and non-melanic forms.
Color was quantified from image data of wings and connexivum (Fig. 2a). Consistent comparisons of color from images requires homologous alignment of the portion of tissue to be analyzed and color-based segmentation of the images . We performed homologous alignment of a portion of right wing and connexivum by extracting 8.5 × 20 mm rectangles from the wing and 1.5 × 18 mm rectangles from the right side of the connexivum (Fig. 2b). We chose the wing to quantify color because we observed the presence of a yellow spot on the upper third of the wings only in the clearer individuals. The colorimetric analysis was performed using the software Image Color Summarizer 0.76 (http://mkweb.bcgsc.ca/color-summarizer/). We obtained the average value of each red green blue (RGB) component, as three separate variables, for the portions of wing and connexivum of each individual (six variables in total). A K-means clustering analysis was used to define the number of colorimetric groups present in the sampled T. infestans populations . Euclidean distance measure was selected for the cluster analysis. The K-means clustering algorithm was run a priori considering up to 10 clusters. The optimal number of clusters for the data set was determined using the elbow and average silhouette methods . The elbow method looks at the percentage of explained variance, the sum of squared errors, as a function of the number of clusters. The chosen number of clusters is the one that does not improve with the addition of a new cluster. The average silhouette method considers that the optimal number of clusters is the one that maximizes the mean silhouette value, a measure of how similar an object is to its own cluster compared to other clusters .
Body data collection
Images of a dorsal view of each individual were processed using the UTHSCSA ImageTool for Windows ver. 3.00. Total body length (L) was measured from the base of the head to the base of the body. Body area (mm2) was calculated excluding the head, pronotum and scutellum. The relationship between W and L was also estimated as the W/L ratio.
Wing data collection
We used a landmark-based geometric morphometry approach to collect 10 type-I landmarks positioned at wing vein intersection, as described elsewhere . Landmarks were collected using TPSdig 2.31 . Comparison of wing size and shape between melanic and non-melanic groups for each sex involved the right wing.
To assess fluctuating asymmetry (FA), left wings were first reflected to their mirror images. Separate landmark configurations were digitized twice in both wings of each individual to estimate measurement error (ME) . Asymmetry can be partitioned in directional and non-directional asymmetry. Directional asymmetry (DA) occurs whenever there is, on average, a greater development of a character on one side of the plane of symmetry than on the other. DA is presumably unrelated to developmental stability . Non-directional asymmetry measures FA.
K-means clustering analysis was done using the software InfoStat . The elbow and silhouette methods used to determine the optimal number of clusters were performed with factoextra R package (https://cran.r-project.org/web/packages/factoextra/index.htmll). After assigning individuals to melanic and non-melanic groups via K-means clustering, comparisons between females and males within each group were significant in most of the performed analyses; therefore, males and females were treated separately in the comparisons between groups.
After checking that assumptions for parametric analyses were met, W, L, W/L ratio, and body area were compared between melanic and non-melanic groups for each sex using a t-test. Wing size was measured using the centroid size (CS) variable. This measure is a single variable of size that integrates different axes of growth, and is measured as the square root of the sum of the squared distances between the centre of the configuration of landmarks and each individual landmark . Wing size (CS) was compared between groups for each sex through a one-way analysis of variance (ANOVA). These analyses were done with the software InfoStat .
Wing shape was compared using a Procrustes approach, by computing differences of landmark coordinates after a full Procrustes fit superimposition of both sexes and groups. Comparisons were performed using discriminant function analysis (DFA), first between sexes and then between melanic and non-melanic groups for each sex. Mahalanobis distances between pairs of species were calculated and their significance was evaluated using a non-parametric test based on permutations (1000 runs). The percentage of phenotypic similarity between pairs of species was calculated using the cross-check test of discriminant analysis. The relationship between CS and shape within each group (allometry) was estimated using a multivariate regression between the Procrustes coordinates (dependent variables) and the CS (independent variable). Morphometric analysis was done using MorphoJ 1.05f .
To assess FA, configurations for melanic and non-melanic females and males were superimposed using the least-squares Procrustes method . A two-way mixed ANOVA with side and individual as fixed and random factors, respectively was performed for each group within each sex to assess the occurrence of DA and FA . A two-way mixed ANOVA for wing size and a Procrustes ANOVA for wing shape were performed in each group to estimate the occurrence of FA in size and shape, respectively . FA index for size and shape was estimated as the mean square (MS) of the interaction between side and individual of the Procrustes ANOVA. These indices were corrected for measurement error (MS/ME × 103).
To analyze allometry a multivariate regression was computed between the Procrustes coordinates (dependent variables) and CS (independent variable) for melanic and non-melanic individuals within each sex. MorphoJ 1.05f  was used for FA analysis.
Colorimetric analysis of melanic and non-melanic forms
An example of the colorimetric analysis procedure is presented in Fig. 2c. The six variables derived from the RGB threshold analysis were used in the K-means clustering analysis. The number of optimal clusters, evaluated with the elbow and silhouette method, was two and three groups, respectively (Fig. 3a, Additional file 1: Figure S1a). A PCA showing the distribution of the individuals in two clusters suggested relatively well-defined melanic and non-melanic groups, with the non-melanic group being more variable than the melanic one (Fig. 3b). The first PCA axis (explaining 49.54% of the total variance) was mainly associated with wing RGB channel colors and the second axis (explaining 45.77% of the total variance), was mainly associated with connexivum RGB colors (Table 1). When three clusters were graphed in a biplot for the first two PCA axes, the melanic group remained the same, whereas the non-melanic group was split into two groups (Fig. 3b, Additional file 1: Figure S1b). For this reason, we chose the criterion of two clusters (melanic and non-melanic groups) for the subsequent analyses.
House infestation and distribution, and frequency of melanic and non-melanic forms
All the studied ecotopes harbored relatively small populations, with low abundance of adults (between 1 and 18) and nymphs (between 1 and 36) of T. infestans (Table 2). House colonization rate was 46.77% (29/62), and house infestation rate was 53% (33/62). Only T. infestans adults occurred in 5 of the 33 infested houses. An established colony was collected in the intradomicile in a single domiciliary unit (IDI = 4.17%, IDC = 4.17%) (Table 2). For peridomestic habitats, PDI was 37.10 % and PDC 43.54 %. One house was infested in more than one peridomestic structure (Table 2). Most (83.2%) T. infestans individuals (161 nymphs and 81 adults) were collected in chicken coops, whereas the remaining ones (17.8%, 37 nymphs and 16 adults) were collected in goat corrals, pig corrals and storerooms. A rabbit hutch harbored only one male of T. infestans (Table 2).
Melanic individuals were present in 13 (54%) of the 24 infested houses. Of the 103 collected adults, 17 (35%) females and 25 (46%) males were assigned to the melanic morph. Houses with presence of melanic individuals were distributed in all sampled localities, except in the only sampled house from Puesto Torrado (Table 2). Melanic and non-melanic adults were collected together in the same habitat in 7 of the 13 houses with melanic individuals (Table 2). In four of the five houses that exhibited only melanic individuals, nymphs were also collected, indicating ongoing colonized sites (Table 2). In 12 of the 24 houses (50%) only non-melanic individuals were collected. Of these 12 houses, nymphs were not collected in only three (25%) (Table 2). In the only house that exhibited domiciliary infestation, all collected adults were melanic (Table 2). Adult peridomestic infestation (APDI) was lower for melanic than for non-melanic adults (APDImelanic = 19.35, APDInon-melanic = 32.26). Considering only infested houses, 45.83% exhibited melanic and 83.33%, non-melanic adults.
Linear morphometric analysis and relationship with nutritional status
Sexual dimorphism was confirmed for L, W and W/L (one-way ANOVA between sexes, P < 0.001 in all cases). The one-way ANOVAs for L and total body area did not show significant differences between melanic and non-melanic females (t(45) = 1.38, P = 0.1736; t(45) = 1.18, P = 0.1619 for body length and body area, respectively), although melanic females tended to have greater body length and area (Fig. 4a, b). For males, L and body area showed significant differences between groups (t(48) = 3.27, P = 0.0020; t(48) = 2.72, P = 0.0091, respectively). Melanic males had greater body area than non-melanic ones (Fig 4a, b). For the entire adult sample, W/L ratio varied between 7.10–17.41 mg/mm. For females, W/L varied between 9.09–14.74 mg/mm for melanic individuals and between 7.12–17.41 mg/mm for non-melanic individuals, whereas for males, W/L varied between 7.11–13.33 for melanic individuals and between 7.10–12.35 mg/mm for non-melanic individuals. W and W/L did not show significant differences between colorimetric groups, either for females or for males (P > 0.05 in all cases). For the houses where both colorimetric groups coexisted and each group was represented by more than one adult (i.e. houses ID No. 19, 805, 809 and 812), comparisons between W/L were made. No significant differences were found in any of the four houses (P > 0.05 in all cases).
Geometric morphometric analysis
Comparisons of wing size (CS) between groups showed significant differences both for females and males (Females: t(46) = 2.01, P = 0.050; males: t(50) = 2.97, P = 0.005), with melanic females and males exhibiting larger wings than their non-melanic counterparts (Fig. 4c).
The first axis of the DFA for wing shape variation between melanic and non-melanic females explained the total variation (100%). Mahalanobis distances between groups were highly significant (Mahalanobis distance: 1.953, P = 0.0008). Non-melanic females showed the lowest misclassification error (90.91% and 84.21% of correct classification from cross-validation for non-melanic and melanic females, respectively). For males, the first axis of the DFA also explained the total variation. Mahalanobis distances were significant between groups (Mahalanobis distance: 1.496, P = 0.0054). Non-melanic males showed the lowest misclassification error (80.77% and 71.43% of correct classification from cross-validation for non-melanic and melanic males, respectively).
Allometric analysis performed through multivariate regressions between the Procrustes coordinates and the CS for melanic and non-melanic individuals within each sex were non-significant (P > 0.05 in all cases).
Fluctuating asymmetry analysis
Significant evidence of FA in wing size was revealed by the two-way mixed ANOVAs for melanic and non-melanic females and males (Table 3). Melanic females and males exhibited significantly higher wing size FA indices than their non-melanic counterparts (Fig. 5a). Procrustes ANOVA for wing shape FA also showed significant FA occurrence for both melanic and non-melanic females and males (Table 3). As for wing size, melanic females and males showed higher FA than non-melanic individuals (Fig. 5b).
Our results showed the existence of a measurable colorimetric variation among the sampled adults of T. infestans collected mainly in peridomestic habitats. To our knowledge, this is the first time that color variation is quantified accurately in a species of the Triatominae. The technique allowed us to group the individuals into two color morphs (i.e. melanic and non-melanic groups) and emerges as a promising tool to study this kind of variation in triatomines, as previously shown in other insect species . Fine-scale quantification of phenotypic color variation has potential applications in the study of populations, sibling species, hybrids and closely related taxa [1, 3, 4].
This is the second report of melanic forms of T. infestans in domestic and peridomestic habitats in the Dry Chaco region of Argentina. Our results showed that melanic adults did not have an aggregated distribution, but were scattered throughout the study area; moreover, in 7 of the 13 infested houses where they were detected, they were collected along with non-melanic individuals. Taking into account that the last community-wide insecticide spraying campaign conducted before the study was three years earlier, house infestation and colonization (53% and 47%, respectively) were relatively high, indicating a good reestablishment of T. infestans colonies. This inference is based on the comparison with other endemic areas of the Argentine Chaco. For example, house infestation was 37% in a rural area located in the southern region of Los Llanos (La Rioja Province, Argentina) after decades of non-systematic vector surveillance and control . In the locality of Pampa del Indio (Chaco Province, Argentina) house infestation was 31.9%, a much lower value than that expected from lack of recent insecticide spraying . In an area in the Monte Desert ecoregion of mid-western Argentina, where the last insecticide spraying campaign conducted before the study was between 2 and 10 years earlier, depending on house accessibility, house infestation was 21.2% .
Heritability (h2) of melanization has been estimated in a few insect species (h2 = 0.73 in Philaenus spumarius , 0.77 in Drosophila immigrans , 0.17–0.37 in Spodoptera littoralis [53, 54] and 0.61 in Gryllus firmus ), indicating that this trait can respond rapidly to selection. In many species, melanism is genetically determined by only one or a few loci (e.g. Lepidoptera [56, 57], Diptera , Coleoptera [59, 60], Orthoptera [61, 62]). Moreover, it has been observed that the pattern of inheritance of melanism follows Mendel’s simple laws of segregation and an independent assortment of alleles in some species [63, 64], like the triatomine species Rhodnius nasutus . In the latter case, the melanic phenotype also seems to be encoded by recessive alleles. The fact that in our study area melanic T. infestans adults were widely spread, with relatively high representation in the area at the house infestation level (45.83%) compared with non-melanic adults (83.33%), suggests that this trait may have been present in these populations for a sufficient number of generations to be well established and represented in this area. Chromatic variants have been reported in other Triatominae species. In Triatoma rubrovaria, four discrete chromatic morphotypes were identified in the same geographical area in the State of Rio Grande do Sul, Brazil . In this case, chromatic variation was represented by patterns of discrete chromatic variants of the pronotum and the frequency of each morphotype was variable according to different ecotopes or municipalities in the same area. Chromatic phenotypes were studied in Triatoma brasiliensis macromelasoma in various locations between the distributions of T. b. brasiliensis and T. juazeirensis. Thirteen different color patterns were characterized, with nine of them being intermediate phenotypes . For Rhodnius nasutus, the occasional appearance of adults with atypical dark coloring was observed in the F3 generation of a colony from Sobral, State of Ceará, Brazil. In this case, dark and typical morphs were discrete color variants easily identifiable in both nymphs and adults .
Environmental factors are considered determinant of variation within structures or traits. Insects generally grow to smaller sizes at high temperatures  and on low quality diets . However, how these environmental signals alter the developmental programme to produce smaller individuals is still unknown. Our results showed that melanic females and males exhibited larger body size, body area and wings than the non-melanic ones (although this tendency was not significant for body size or body area in females). In the Triatominae, evidence of size reduction was reported for wings of domestic specimens of T. infestans collected from two sites in Cochabamba, Bolivia (Laguna Angostura and Jamach’uma) compared with sylvatic dark morphs collected 1 km away from these domestic foci . A reduction in body size was also observed in Panstrongylus rufotuberculatus (Triatominae) in the transition from sylvatic to domestic habitats  and from natural to artificial (laboratory) habitats for Panstrongylus geniculatus, Triatoma brasiliensis and Triatoma flavida [72, 73]. Considering that in our study area there is no geographical or even physical separation (in seven houses melanic and non-melanic individuals were collected in the same sites) between collected melanic and non-melanic individuals, differences in wing size may be attributed to differences in ecological adaptation of both colorimetric groups to the same habitats. It has been postulated that melanism is an important trait related to thermoregulation, with darker individuals probably heating up faster and reaching higher body temperatures than light-colored ones due to the higher absorption of solar radiation . Thermal capacity is strongly related to body size; larger bodied individuals would exhibit lower cooling rates (higher heat conservation) but also lower heating rates (lower capacity to gain heat) than smaller bodied ones.
Shape variation in T. infestans is considered frequently driven by genetic variations . Wing shape differences between groups showed that these populations are structured, suggesting different origins for both morphs. The fact that the area of study had been treated with insecticide three years before this study suggests that T. infestans from both colorimetric groups may have become from different sources (non-treated nearby areas and/or residual or untreated foci [76,77,78,79]). A strong spatial structuring at a microgeographical scale using either phenotypic or genotypic markers after community-wide or selective residual spraying with insecticides has been reported for T. infestans populations in other areas of the Dry Chaco [80,81,82]. Considering that starved individuals would disperse to locate eventual blood sources, high population structuring might be expected to be associated with high nutritional status . Nutritional status is, among others, a factor that modulates flight dispersal in the Triatominae . For T. infestans, many authors use the W/L ratio as an indicator of flight initiation [30,31,32,33]. Our results showed that this ratio did not differ between groups. For T. infestans, adults with W/L < 8 mg/mm have high probability of flight initiation [30, 31]. The W/L ratio for the adults analyzed in this study varied between 7.10–17.41 mg/mm. Eight (8.7%) of the adults exhibited a W/L ratio < 8 mg/mm, of which 6 (75%) were non-melanic and 7 (87.5%) were males. Despite these observations, within each house, results showed no significant differences in W/L between melanic and non-melanic groups, confirming that flight initiation is not associated with a particular colorimetric group.
Our results showed that melanic females and males exhibited significantly higher FA indices than non-melanic ones (24.6% and 69.9% higher for wing size and shape for females, respectively, and 72.7% and 34.8% higher for wing size and shape for males, respectively). This apparent difference in the magnitude of FA between colorimetric groups was also observed in T. infestans collected in different ecotopes, seasons and years, and before and after a community-wide pyrethroid spraying campaign in rural areas of Argentina [84, 85]. In these studies, different developmental instability suggests congruence of FA with habitat/host combinations and/or timing of flight dispersal for the former example  and differential survival of adults with more symmetric wings for the latter study . In the context of the present study, differences in wing developmental instability suggest a different influence of environmental conditions during development between morphs. It is well known that there is a wide range of phenotypic effects correlated with melanism, including developmental time, behavioural traits and different tolerance to environmental characteristics, such as temperature, UV radiation and relative humidity (see Table 2c in ). For example, in Manduca sexta, melanic forms are less tolerant to desiccation than typical non-melanic individuals . Environmental characteristics of the different domestic and peridomestic habitats where individuals develop might have a differential effect on the level of instability of melanic and non-melanic individuals. Another possibility that may explain the differences found in FA indices could be associated with the general hypothesis that melanization is costly and will result in trade-offs [1, 54]. This hypothesis is based on two premises: first, tyrosine, the main compound of the pro-phenol oxidase system, can only be obtained from ingested food; secondly, melanin is a nitrogen-rich compound that may require substantial nitrogen or protein investment for its synthesis . Several traits described in the literature as associated with melanism are likely to have negative consequences on fitness . In this sense, melanic individuals should be less adapted to domestic and peridomestic habitats where non-melanic ones seem to be better adapted, considering the recorded levels of FA. While the physical attributes of melanism may have obvious selective advantages (e.g. thermal regulation), the effect of energetic costs on fitness-related traits such as fecundity depends on how such costs interact with these traits. Moreover, dispersive capacity of melanic individuals may be mechanically compromised due to their less symmetrical wings, considering flight performance and wing function .
Our results showed that although chromatic cluster analysis clearly defined a melanic and a non-melanic group, within the non-melanic group, there are also two well-differentiated groups by the PC2 (cluster 2 and 3, Additional file 1: Figure S1). These two groups showed significant differences in total body length and total body area for females and significant differences in wing shape for males. Fluctuating asymmetry analysis showed differences between groups (4.6% and 2.4% for wing size and shape for females, respectively, and 3.2% and 5.8% for wing size and shape for males, respectively). These results indicate that differences between groups were not always consistent between sexes and do not allow us to draw a clear conclusion about whether the color variants within the non-melanic groups can be attributed to different origins of reinfestation or are only a common colorimetric variants exhibited by a population.
To the best of our knowledge, this study is a first approach to the melanic and non-melanic characterization in T. infestans morphs associated with domestic and peridomestic structures. Understanding other pleiotropic effects associated with melanism, identification of responsible genes of melanism and selective responses of melanic alleles in nature, among other factors, would help elucidate the evolutionary history, the epidemiological importance and the possible origins of melanic forms in domestic and peridomestic T. infestans populations.
Two colorimetric groups of T. infestans were recorded in a rural area from Cruz del Eje, Córdoba Province, Argentina, with a high infestation rate (53%). Melanic adults were widespread throughout the study area, although less represented than non-melanic ones. Linear and geometric morphometric analyses showed a greater body size for melanic females and males than for non-melanic ones, suggesting possible differences in ecological adaptation. Wing shape analysis revealed differences between groups, which may be related to their different origins. Melanic and non-melanic individuals exhibited good nutritional status and low probability of flight initiation (W/L > 8 mg/mm). Wing size and wing shape FA results are in agreement with lower adaptation of melanic individuals to domestic and peridomestic environments or a cost associated with melanism. This study provides evidences of a well-established T. infestans chromatic polymorphism in the area and possible pleiotropic effects of melanism associated with environmental conditions.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional file. Raw data are available from the corresponding author upon reasonable request.
adult peridomestic infestation
infestation in the intradomicile
colonization in the intradomicile
total body length
infestation in the peridomicile
colonization in the intradomicile
red, green, blue
total weight/total body length ratio
True JR. Insect melanism: the molecules matter. Trends Ecol Evol. 2003;18:640–7.
Majerus MEN. Melanism: evolution in action. 1st ed. New York: Oxford University Press; 1998.
Dubovskiy IM, Whitten MMA, Kryukov VY, Yaroslavtseva ON, Grizanova EV, Greig C, et al. More than a colour change: insect melanism, disease resistance and fecundity. Proc R Soc B. 2013;280:20130584.
Watt WB. Adaptive significance of pigment polymorphism in Colias butterflies. I: variation of melanin in relation to thermoregulation. Evolution. 1968;22:437–58.
Safranek L, Riddiford LM. The biology of the black larval mutant of the tobacco hornworm. Manduca sexta. J Insect Physiol. 1975;21:1931–8.
Huey R, Kingsolver J. Evolution of thermal sensitivity of ectotherm performance. Trends Ecol Evol. 1998;4:131–5.
de Jong PW, Gussekloo SWS, Brakefield PM. Differences in thermal balance, body temperature and activity between non-melanic and melanic two-spot ladybird beetles (Adalia bipunctata) under controlled conditions. J Exp Biol. 1996;199:2655–66.
Schofield CJ, Jannin J, Salvatella R. The future of Chagas disease control. Trends Parasitol. 2006;22:583–8.
WHO. Chagas disease (American trypanosomiasis). Geneva: World Health Organization; 2017. http://www.who.int/mediacentre/factsheets/fs340/en/. Accessed 26 Jul 2019.
Bucher EH, Huszar PC. Sustainable management of the Gran Chaco of South America: ecological promise and economic constraints. J Environ Manage. 1999;57:94–108.
Spillmann C, Burrone S, Coto H. Análisis de la situación epidemiológica de Chagas en Argentina: avances en el control, 2012. Rev Argent Salud Pública. 2013;4:40–4.
Le Poul Y, Whibley A, Chouteau M, Prunier F, Llaurens V, Joron M. Evolution of dominance mechanisms at a butterfly mimicry supergene. Nat Commun. 2014;5:5644.
Noireau F, Flores R, Gutierrez T, Dujardin JP. Detection of sylvatic dark morphs of Triatoma infestans in the Bolivian Chaco. Mem Inst Oswaldo Cruz. 1997;92:583–4.
Noireau F, Flores R, Gutierrez T, Abad-Franch F, Flores E, Vargas F. Natural ecotopes of Triatoma infestans dark morph and other wild triatomines in the Bolivian Chaco. Trans R Soc Trop Med Hyg. 2000;94:23–7.
Martínez A, Olmedo RA, Carcavallo RU. Una nueva subespecie argentina de Triatoma infestans. Chagas. 1987;4:479–80.
Yeo M, Acosta N, Llewellyn M, Sánchez H, Adamson S, Miles GA, et al. Origins of Chagas disease: Didelphis species are natural hosts of Trypanosoma cruzi I and armadillos hosts of Trypanosoma cruzi II, including hybrids. Int J Parasitol. 2005;35:225–33.
Ceballos LA, Piccinali RV, Berkunsky I, Kitron U, Gürtler RE. First finding of melanic sylvatic Triatoma infestans (Hemiptera: Reduviidae) in the Argentine Chaco. J Med Entomol. 2009;46:1195–202.
Lobbia PA, Alvarez R, Picollo MI, Mougabure-Cueto G. First record of domestic colonies of the dark chromatic variant of Triatoma infestans (Hemiptera: Reduviidae). Rev Soc Entomol Arg. 2019;78:33–7.
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.
Noireau F, Bastrenta B, Catalá S, Dujardin JP, Panzera F, Torres M, et al. Sylvatic population of Triatoma infestans from the Bolivian Chaco: from field collection to characterization. Mem Inst Oswaldo Cruz. 2000;95:119–22.
Gumiel M, Catalá S, Noireau F, Arias AR, García A, Dujardin JP. Wing geometry in Triatoma infestans (Klug) and T. melanosoma Martinez, Olmedo and Carcavallo (Hemiptera: Reduviidae). Syst Entomol. 2003;28:173–80.
Bargues MD, Klisiowicz DR, Panzera F, Noireau F, Marcilla A, Pérez R, et al. Origin and phylogeography of the Chagas disease main vector Triatoma infestans based on nuclear rDNA sequences and genome size. Infect Genet Evol. 2006;6:46–62.
Piccinali RV, Marcet PL, Ceballos LA, Kitron UD, Gurtler RE, Dotson M. Genetic variability, phylogenetic relationships and gene flow in Triatoma infestans dark morphs from the Argentinean Chaco. Infect Genet Evol. 2011;11:895–903.
Ceballos LA, Piccinali RV, Marcet PL, Vazquez-Prokopec GM, Cardinal MV, Schachter-Broide J, et al. Hidden sylvatic foci of the main vector of Chagas disease Triatoma infestans: threats to the vector elimination campaign? PLoS Negl Trop Dis. 2011;5:e1365.
De Rojas Arias A, de la Carbajal Fuente AL, Gómez A, Cecere MC, Rolón M, Vega Gómez MC, et al. Morphometric wings similarity among sylvatic and domestic populations of Triatoma infestans (Hemiptera: Reduviidae) from the Gran Chaco region of Paraguay. Am J Trop Med Hyg. 2017;97:481–8.
Gaspe MS, Schachter-Brodie J, Gurevitz JM, Kitron U, Gürtler RE, Dujardin JP. Microgeographic spatial structuring of Triatoma infestans (Hemiptera: Reduviidae) populations using wing geometric morphometry in the Argentine Chaco. J Med Entomol. 2012;49:504–14.
Schachter-Broide J, Dujardin J-P, Kitron U, Gürtler RE. Spatial structuring of Triatoma infestans (Hemiptera, Reduviidae) populations from northwestern Argentina using wing geometric morphometry. J Med Entomol. 2004;41:643–9.
Carey JR. Insect biodemography. Annu Rev Entomol. 2001;46:79–110.
Lehane MJ. The biology of blood-sucking in insects. 2nd ed. New York: Cambridge University Press; 2005.
Lehane MJ, Schofield CJ. Flight initiation in Triatoma infestans (Klug) (Hemiptera: Reduviidae). Bull Entomol Res. 1982;72:497–510.
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.
Ceballos LA, Vazquez-Prokopec GM, Cecere MC, Marcet PL, Gurtler RE. Feeding rates, nutritional status and flight dispersal potential of peridomestic populations of Triatoma infestans in rural northwestern Argentina. Acta Trop. 2005;95:149–59.
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.
Bijlsma R, Loeschcke V. Environmental stress, adaptation and evolution: an overview. J Evol Biol. 1997;18:744–9.
Møller AP, Swaddle JP. Asymmetry, developmental stability and evolution. Oxford: Oxford University Press; 1997.
Schofield CJ. A comparison of sampling techniques for domestic populations of Triatominae. Trans R Soc Trop Med Hyg. 1978;72:449–55.
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.
Van Belleghem SV, Papa R, Ortiz-Zuazaga H, Hendrickx F, Jiggins CD, Owen McMillan W. patternize: an R package for quantifying colour pattern variation. Methods Ecol Evol. 2018;9:390–8.
Hartigan JA, Wong MA. Algorithm AS 136: A K-means clustering algorithm. J Royal Stat Soc Series C. 1979;28:100–8.
Kaufman L, Rousseeuw PJ. Finding groups in data. 1st ed. New York: Wiley; 2005.
Rohlf FJ. TPSDig2 2.31. New York: Stony Brook University; 2017.
Palmer AR. Fluctuating asymmetry analyses: a primer. In: Markow TA, editor. Developmental instability: its origins and evolutionary implications. Dordrecht: Kluwer Academic Publishers; 1994. p. 335–64.
Di Rienzo JA, Casanoves F, Balzarini MG, González L, Tablada M, Robledo CW. InfoStat versión 2018. Grupo InfoStat. Córdoba: FCA, Universidad Nacional de Córdoba; 2018.
Bookstein FL. Morphometric tools for landmark data. 1st ed. Cambridge: Cambridge University Press; 1991.
Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011;11:353–7.
Klingenberg CP, McIntyre GS. Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with Procrustes methods. Evolution. 1998;52:1363–75.
Palmer AR, Strobeck C. Fluctuating asymmetry: measurement, analysis, pattern. Annu Rev Ecol Syst. 1986;17:391–421.
Gorla D, Porcasi X, Hrellac H, Catalá S. Spatial stratification of house infestation by Triatoma infestans in La Rioja, Argentina. Trans R Soc Trop Med Hyg. 2009;80:405–9.
Gaspe MS, Provecho YM, Cardinal MV, Fernández MP, Gürtler RE. Ecological and sociodemographic determinants of house infestation by Triatoma infestans in indigenous communities of the Argentine Chaco. PLoS Negl Trop Dis. 2015;9:e000361.
de la Carbajal Fuente AL, Provecho YM, Fernández MP, Cardinal MV, Lencina P, Spillman C, et al. The eco-epidemiology of Triatoma infestans in the temperate monte desert ecoregion of mid-western Argentina. Mem Inst Oswaldo Cruz. 2017;112:698–708.
Yurtsever S. The high heritability estimates of ventral darkness variation in Philaenus spumarius (Homoptera: Cercopidae). Israel J Zool. 2000;46:249–59.
Singh S, Ramniwas S, Parkash R. Fitness consequences of body melanization in Drosophila immigrans from montane habitats. Entomol Res. 2009;39:182–91.
Cotter SC, Myatt JP, Benskin CMH, Wilson K. Selection for cuticular melanism reveals immune function and life-history trade-offs in Spodoptera littoralis. J Evol Biol. 2008;21:1744–54.
Lee KP, Wilson K. Melanism in a larval Lepidoptera: repeatability and heritability of a dynamic trait. Ecol Entomol. 2006;31:196–205.
Roff DA, Fairbairn DJ. The costs of being dark: the genetic basis of melanism and its association with fitness-related traits in the sand cricket. J Evol Biol. 2013;26:1406–16.
Grant BS. Allelic melanism in American and British peppered moths. J Hered. 2004;95:97–102.
Suzuki Y, Nijhout H. Genetic basis of adaptive evolution of a polyphenism by genetic accommodation. J Evol Biol. 2007;21:56–66.
Wittkopp PJ, Beldade P. Development and evolution of insect pigmentation: genetic mechanisms and the potential consequences of pleiotropy. Semin Cell Dev Biol. 2009;20:65–71.
de Jong PW, Brakefield PM, Geerinck BP. The effect of female mating history on sperm precedence in the two-spot adybird, Adalia bipunctata (Coleoptera, Coccinellidae). Behav Ecol. 1998;9:559–65.
Michie LJ, Mallard F, Majerus MEN, Jiggins FM. Melanic through nature or nurture: genetic polymorphism and phenotypic plasticity in Harmonia axyridis. J Evol Biol. 2010;23:1699–707.
Tanaka S. Hormonal deficiency causing albinism in Locusta migratoria. Zool Sci. 1993;10:467–71.
Maeno K, Tanaka S. Genetic and hormonal control of melanization in reddish-brown and albino mutants in the desert locust Schistocerca gregaria. Physiol Entomol. 2010;35:2–8.
Roff DA, Fairbairn DJ. The evolution of alternative morphologies: fitness and wing morphology in male sand crickets. Evolution. 1993;47:1572–84.
McRobie H, Thomas A, Kelly J. The genetic basis of melanism in the gray squirrel (Sciurus carolinensis). J Hered. 2009;100:709–14.
Braga Stehling Dias F, 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. 1859;2014(47):637–41.
Almeida CE, Pacheco RS, Noireau F, Costa J. Triatoma rubrovaria (Blanchard, 1843) (Hemiptera: Reduviidae) I: isoenzymatic and chromatic patterns of five populations from the State of Rio Grande do Sul. Brazil. Mem Inst Oswaldo Cruz. 2002;97:829–34.
Costa J, Bargues MD, Lima Neiva V, Lawrence GG, Gumiel M, Oliveira G, et al. Phenotypic variability confirmed by nuclear ribosomal DNA suggests a possible natural hybrid zone of Triatoma brasiliensis species complex. Infect Genet Evol. 2016;37:77–87.
Atkinson D. Temperature and organism size—a biological law for ectotherms. Adv Ecol Res. 1994;25:1–58.
Chapman RF. The insects: structure and function. 4th ed. Cambridge: Cambridge University Press; 2012.
Dujardin JP, Bermudez H, Casini C, Schofield CJ, Tibayrenc M. Metric differences between sylvatic and Triatoma infestans (Heteroptera: Reduviidae) in Bolivia. J Med Entomol. 1997;34:544–51.
Dujardin JP, Forgues G, Torrez M, Martinez E, Córdoba C, Gianella A. Morphometrics of domestic Panstrongylus rufotuberculatus in Bolivia. Ann Trop Med Parasitol. 1998;92:219–28.
Jaramillo-O N, Castillo D, Wolff M. Geometric morphometric differences between Panstrongylus geniculatus from field and laboratory. Mem Inst Oswaldo Cruz. 2002;97:667–73.
Borges EC, Dujardin JP, Schofield CJ, Romanha AJ, Diotaiuti L. Dynamics between sylvatic, peridomestic and domestic populations of Triatoma brasiliensis (Hemiptera: Reduviidae) in Ceará State, northeastern Brazil. Acta Trop. 2005;93:119–26.
Gates DM. Biophysical ecology. 1st ed. New York: Springer; 1980.
Dujardin JP, Costa J, Bustamante D, Jaramillo N, Catalá S. Deciphering morphology in Triatominae: the evolutionary signals. Acta Trop. 2009;110:101–11.
Germano MD, Santo Orihuela P, Roca Acevedo G, Toloza AC, Vassena C, Picollo MI, et al. Scientific evidence of three different insecticide-resistant profiles in Triatoma infestans (Hemiptera: Reduviidae) populations from Argentina and Bolivia. J Med Entomol. 2012;49:1355–60.
Gurevitz JM, Gaspe MS, Enriquez GF, Provecho YM, Kitron U, Gürtler RE. Intensified surveillance and insecticide-based control of the Chagas disease vector Triatoma infestans in the Argentinean Chaco. PLoS Negl Trop Dis. 2013;7:e2158.
Picollo MI, Vassena C, Orihuela PS, Barrios S, Zaidemberg M, Zerba E. High resistance to pyrethroid insecticides associated with ineffective field treatments in Triatoma infestans (Hemiptera: Reduviidae) from northern Argentina. J Med Entomol. 2005;42:637–42.
de Rojas Arias A, Abad-Franch F, Acosta N, López E, González N, Zerba E, et al. Post-control surveillance of Triatoma infestans and Triatoma sordida with chemically-baited sticky traps. PLoS Negl Trop Dis. 2012;6:e1822.
Gaspe MS, Gurevitz JM, Gürtler RE, Dujardin J-P. Origins of house reinfestation with Triatoma infestans after insecticide spraying in the Argentine Chaco using wing geometric morphometry. Infect Genet Evol. 2013;17:93–100.
de Pérez Rosas AR, Segura EL, Fichera L, García BA. Macrogeographic and microgeographic genetic structure of the Chagasʼ disease vector Triatoma infestans (Hemiptera: Reduviidae) from Catamarca, Argentina. Genetica. 2008;133:247–60.
Marcet PL, Mora MS, Cutrera AP, Jones L, Gürtler RE, Kitron U, et al. Genetic structure of Triatoma infestans populations in rural villages of Santiago del Estero, northern Argentina. Infect Genet Evol. 2008;8:835–46.
Hernández ML, Abrahan LB, Dujardin JP, Gorla DE, Catalá SS. 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 Zoonot Dis. 2011;11:503–13.
Nattero J, Cecere MC, Gürtler RE. Temporal variations of fluctuating asymmetry in wing size and shape of Triatoma infestans populations from northwest Argentina. Infect Genet Evol. 2015;56:133–42.
Nattero J, Piccinali RV, Gaspe MS, Gürtler RE. Fluctuating asymmetry and exposure to pyrethroid insecticides in Triatoma infestans populations in northeastern Argentina. Infect Genet Evol. 2019;74:103925.
Gonzalez-Santoyo I, Córdoba-Aguilar A. Phenoloxidase: a key component of the insect immune system. Entomol Exp Appl. 2012;142:1–16.
McLachlan AJ. Fluctuating asymmetry in flies, what does it mean? Symmetry. 2010;2:1099–107.
We thank the technicians of the National Programme of Chagas for technical and logistical support in the field. We thank Ana Graciela López, Carola Soria and Florencia Carnicero for their contributions in laboratory. We thank the two anonymous reviewers whose suggestions helped improve and clarify the paper. JN, ALCF and RVP are members of the CONICET Researcher’s Career. MC profits from a doctoral scholarship from CONICET.
This study was supported by SECYT (National University of Córdoba, Grant No. 05/I554), CONICET (PIP No 11220110101102) and the Ministry of Science and Technology of Córdoba.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
K-means clustering analysis. a Plot of the curve of the average silhouette for 10 K-values (silhouette method). The optimal number of clusters suggested by this method (three) is indicated with a dashed line. b Biplot derived from the first two axes of a principal component analysis showing the spatial distribution of the three clusters derived from a K-means clustering analysis (K = 3).
About this article
Cite this article
Nattero, J., Carbajal de la Fuente, A.L., Piccinali, R.V. et al. Characterization of melanic and non-melanic forms in domestic and peridomestic populations of Triatoma infestans (Hemiptera: Reduviidae). Parasites Vectors 13, 47 (2020). https://doi.org/10.1186/s13071-020-3912-y