Skip to main content

Effect of larval density and substrate quality on the wing geometry of Stomoxys calcitrans L. (Diptera: Muscidae)



In insects, oviposition decisions may lead to egg deposition in substrates with different larval density and nutritional levels. Individuals developing in such substrates may present plasticity in their phenotype. Here, we investigated the effect of two factors related to oviposition decisions, namely larval density and substrate quality, on the wing size and wing shape of the stable fly, Stomoxys calcitrans L. (Diptera: Muscidae).


We reared S. calcitrans larvae at different densities (5, 15 and 25) and on different substrates (camel, cow, donkey and sheep dung). For each fly that emerged, we recorded body weight, and detached, slide-mounted and photographed the right wing. Next, we collected 15 landmarks on each photographed wing, and applied geometric morphometric analysis to assess variation in wing size and wing shape of S. calcitrans across the different larval densities and substrate types.


We observed that wing size and wing shape of S. calcitrans were affected by larval density and the nature of the developmental substrate. Flies reared in a group of 5 had larger wing centroid size, wing length, wing width, wing area and wing loading compared with those reared in a group of 25. Also, flies developed in donkey and sheep dung had larger wing centroid size, wing length, wing width, wing area and wing loading in comparison with those grown in camel and cow dung. Canonical variate analysis followed by discriminant analysis revealed significant wing shape variation in S. calcitrans across the different densities and substrates. Wing size had a significant but weak positive effect on wing shape.


This study demonstrates the high sensitivity of S. calcitrans wings to variation in larval density and developmental substrate, and that use of landmark-based geometric morphometric analysis could improve our understanding of how flies of veterinary importance respond to environmental variability.


In holometabolous insects, individual fitness mostly relies on oviposition decisions by gravid females. When, where and how mothers deposit their eggs can affect the performance and phenotype of their progeny [1]. Therefore, it is important for females to oviposit on a substrate that provides the best conditions for the next generation. However, a female may fail to make the seemingly optimal choice to oviposit on an appropriate substrate that could enhance offspring fitness. For instance, Heard [2] found that in the pitcher plant mosquito, Wyeomyia smithii Coquillett, although larval fitness is better in pitchers with fewer conspecific and more midge larvae, gravid females did not deposit more eggs in such pitchers. Instead, they laid more eggs in pitchers containing either midges or conspecific larvae. Wong et al. [3] found that larval survival and development of Aedes aegypti L. was poor in containers where gravid females laid more eggs. This imperfection in oviposition decisions generally leads to phenotypic variation [4, 5], in which individuals react to the inputs of their breeding substrate with a change in their form, state, movement, or rate of activity [6]. These inputs include environmental factors such as the dietary value of the substrate and the number of individuals sharing the same substrate (density) [7]. In evolutionary ecology, understanding how these factors influence organism phenotype is a fundamental concern because such flexibility can affect fitness, generate novelty, facilitate evolution, and structure ecological communities [8].

Insect wings are good indicators of population responses to changes that occur in their environment [9]. Lin et al. [10] demonstrated that variation in food nutrient content and density are key ecological factors related to the expression of condition-dependent, adaptive phenotypes such as wing polyphenisms. They found that in the brown planthopper, Nilaparvata lugens Stål, a serious rice pest, an increase of long-winged N. lugens in a population is related to higher glucose levels in host rice plants. Conversely, the appearance of the short-winged form of N. lugens is linked to a reduction in host glucose level. Dipteran species are particularly well suited for studying phenotypic changes induced by the environment because their wings are highly plastic and wing landmarks are homologous across various species [11]. Variation in food quality and population density are key factors associated with fly wing polyphenisms. In Drosophila buzzatii Patterson & Wheeler and Drosophila koepferae Soto et al. [12] detected significant differences in wing size and shape between flies that were reared on different cactus hosts. In Ae. aegypti, males and females have longer wings when developed in conditions of low larval density [13].

Changes in wing morphology are known to affect insect flight aerodynamics. Long and slender wings are optimal for long-duration flight, while short and broad wings are optimal for slow and agile flight [14]. Also, broad wing bases allow a wider range of speed and a narrow wing tip allows less costly, extensive flight [15]. As a consequence, wing morphology is closely related to several insect behavioural activities including food searching, location of breeding sites and sexual partners, and avoidance of natural enemies. In some mosquito species such as Anopheles gambiae Giles where wing size is positively correlated with body size, an increase in wing size augments the frequency of blood meals [16]. This then leads to an increase in the likelihood of pathogen dissemination [16]. In Aedes albopictus Skuse, there is a positive correlation between wing length, larval diet quality, and the number of eggs laid [17]. It has been demonstrated that males of the olive fruit fly, Bactrocera oleae Rossi, with large wings (characterised by a high vibration frequency) achieve higher mating success than males with smaller wings [18].

Fly wing vein networks are excellent models for statistical analysis of size and shape variation [19]. In recent years, landmark-based geometric morphometric analysis has been increasingly used to analyse insect wings to address intraspecific variation [20, 21], interspecific variation [22, 23], sexual dimorphism [24, 25], parasite detection [26, 27], laboratory strain separation [28] and phenotypic plasticity [12, 29, 30]. Geometric morphometric techniques are potent tools to assess the correlation between the size and shape of organisms and environmental variables. The approach uses coordinates of identified morphological “landmarks” to study the form of biological structures in two or three dimensions. It involves several statistical techniques that preserve shape information and detect even subtle morphological variations [31]. Moreover, geometric morphometric techniques are cheap, simple and fast [32]. Using geometric morphometric analysis, this study examined the changes that occur in the wings of the stable fly, Stomoxys calictrans L. (Diptera: Muscidae) reared on different substrates and over a range of larval densities.

Stomoxys calcitrans is a cosmopolitan haematophagous fly that mechanically transmits viruses (e.g. West Nile fever virus, Rift Valley fever virus), bacteria (e.g. Bacillus anthracis, Pasteurella multocida), protozoans (e.g. Trypanosoma evansi, Besnoitia besnoiti), and helminths (e.g. Habronema microstoma, Dirofilaria repens) to their hosts, which include cattle, camels, horses, dogs, and humans [33,34,35]. During outbreaks, S. calcitrans can reduce weight gain in cattle by up to 19%, and lead to a 40–60% reduction in milk yields [36, 37]. In the USA, Taylor et al. [38] estimated economic losses attributable to S. calcitrans infestation at around $2.211 billion per year. Gravid female S. calcitrans oviposit on vertebrate herbivore dung including that of camel, cow, donkey and sheep, with the latter two the most preferred [39]. It has already been demonstrated that the fitness of S. calcitrans immature stages (hatchability, developmental time, emergence time, larval and pupal weight) varies across these substrates due to differences in their physicochemical composition [39]. However, the way in which preferred and non-preferred substrates affect S. calcitrans wing size and wing shape remains unclear. Furthermore, not only the substrate nutrient quality but also larval density should be assessed. We hypothesised that larval density and vertebrate herbivore dung type on which S. calcitrans develop would affect wing size and shape.


Biological material

Stomoxys calcitrans flies were obtained from a single culture that had been established for approximately 8 months at the Duduville campus of the International Centre for Insect Physiology and Ecology (icipe) in Nairobi (1°13′12″S, 36°52′48″ E; c.1600 m above sea level). By sourcing experimental flies from a laboratory culture, we minimised potential variation between populations. Adults reared from rabbit faeces were kept in cages (75 × 60 × 45 cm) under conditions of 25 ± 5 °C and 65 ± 5% relative humidity with a photoperiod of 12L:12D. Flies were fed twice per day (8:00 and 16:00 h) on defibrinated bovine blood poured on moistened cotton.

Density experiment

Gravid female S. calcitrans from the established colony were allowed to oviposit on donkey dung placed in plastic containers (21.5 × 14.5 × 7.4 cm). Baleba et al. [39] demonstrated that this dung is best for S. calcitrans development. To assess the effect of density of S. calcitrans on wing size and shape, we reared S. calcitrans larvae at varying densities by gently transferring (using soft forceps) 5, 15 and 25 first-instar larvae to plastic cups (200 ml) filled with 25 g of donkey dung. We replicated this process several times to obtain 30 emerged females and 30 emerged males. After emergence, each individual was weighed, killed in 70% ethanol, and its right wing was gently removed from the thorax using a fine clamp. The removed wings were slide-mounted (dorsally placed between two microscope slides) to avoid deformation and to enhance accuracy during photography and landmark collections [40]. We photographed the wings at 16× magnification with a Leica DFC320 digital camera coupled to a Leica S6 microscope.

Substrate quality experiment

To test the effect of substrate on S. calcitrans wing size and shape, 10 first-instar larvae were transferred to and permitted to develop on 25 g of camel, cow, donkey or sheep dung. It has been shown that these substrates differ in their physico-chemical composition and affect development and adult body weight [39]. We replicated the process several times to obtain 30 emerged females and 30 emerged males. After emergence, each individual of S. calcitrans was weighed, killed in 70% ethanol, the right wing was removed and slide-mounted, and wings were photographed as described above.

Wing geometric morphometric analysis

To collect wing landmark coordinates, we opened the digital photographs in ImageJ software [41] and generated Cartesian coordinates for 15 wing landmarks (Fig. 1). To quantify measurement error relative to the landmark digitalisation, we collected landmarks for the wings of all 30 individuals reared from donkey dung three times. After executing the generalised Procrustes analysis to extract shape information from the data and eliminate differences in orientation, position and isometric size, we ran analysis of variance (ANOVA) and multivariate analysis of variance (MANOVA) tests to determine if wing size and the wing shape of S. calcitrans varied across the three landmark collections [42] (Table 1).

Fig. 1

Dorsal view of the right wing of S. calcitrans. Numbers indicate the location of 15 selected landmarks (described in Table 1)

Table 1 Description of the 15 anatomical landmarks used to characterise S. calcitrans wing geometry. Numbers relate to landmarks shown in Fig. 1

Four parameters were derived to describe wing size of S. calcitrans: (i) centroid size; (ii) wing length (distance between the 6th and 12th landmark); (iii) wing width (distance between the 3rd and 8th landmark) (Fig. 1); and (iv) wing area. The centroid size, also called the “configuration barycentre” is a global size (or multidimensional measurement) calculated as the square root of the sum of squared Euclidean distances between each landmark and the wing centroid. We computed centroid size using PAST software V.3.09 [43]. Also, based on the adult weight parameter, we calculated wing loading (in kg/m2) using the formula: wl = mass/wing area [44]. To assess the effect of larval density (5, 15 and 25) and substrate type (camel, cow, donkey and sheep dung) on the parameters described above, we ran analyses of variance (ANOVA) followed by post-hoc Student-Newman-Keuls (SNK) tests after checking the wing size parameters for normality using the Shapiro-Wilk test (P > 0.05). To identify correlations between centroid size, wing length, wing width, wing area, adult weight and wing loading, we performed separate principal components analysis (PCA) for larval density and substrate type. We used R version 3.5.1 software [45] to compute all statistical analyses.

To assess wing shape variation across the different densities and substrates, we imported the raw landmark Cartesian coordinates into MorphoJ software [46]. This software was first used to perform a generalised Procrustes analysis to extract shape information from the data and eliminate differences in orientation, position and isometric size. Afterwards, we ran separate multivariate analyses of variance (MANOVA) to compare wing shapes across the different larval densities (5, 15 and 25) and substrates (camel, cow, donkey and sheep dung). Using PAST software, we performed thin plate spline analysis to visualise wing shape deformations. We used canonical variate analysis combined with discriminant analysis to analyse the relative similarities and dissimilarities of the different wing groups. To determine the significance of pairwise differences in mean shapes, we performed permutation tests (10,000 rounds) with Mahalanobis distances and Procrustes distances. To assess the effect of wing size on wing shape (allometry), we fit a linear regression between the Procrustes coordinates and the centroid size, using a permutation test with 10,000 randomisations. For all of these analyses, we excluded the effect of sex because preliminary analyses indicated that wing shape of females and males did not differ.


Measurement error test

The landmarks measured repeatedly on the same individual wing (3 times) were not significantly different for both size (ANOVA, F(2,84) = 0.03; P = 0.97) and shape (MANOVA, F(48,2016) = 0.6; Pillai’s trace = 0.16; P = 0.99). Therefore, we ruled out error due to landmark digitalisation, and we considered that any differences found in the morphology of S. calcitrans wings resulted from the two factors manipulated in our study (larval density and substrate type).

Effect of larval density on the wing geometry of S. calcitrans

Wing size parameters

Larval density significantly affected wing centroid size (F(2,84) = 104.1, P < 0.0001), wing length (F(2,84) = 97.91, P < 0.0001), wing width (F(2,84) = 85.63, P < 0.0001), wing area (F(2,84) = 22.67, P < 0.0001), adult weight (F(2,84) = 4.51, P = 0.014) and wing loading (F(2,84) = 14.35, P < 0.0001) of S. calcitrans. We obtained the largest wing centroid size (Fig. 2a.i), wing length (Fig. 2a.ii), wing width (Fig. 2a.iii), wing area (Fig. 2a.iv), adult weight (Fig. 2a.v) and wing loading (Fig. in flies reared from a group of five larvae. The biplot from the principal component analysis separated flies reared in a group of five from those reared in a group of 25; flies reared in a group of 15 occupied an intermediate position (Fig. 2a). The two first dimensions accounted for 85.6% of the total wing size variation. Dimension 1 explained 60.6% of the total variation, with wing length as the major contributor. Dimension 2 accounted for 25% of the total variation, with wing loading as the major contributor. Except for an absence of correlation between adult weight and centroid size (r = 0.18, P = 0.094), and adult weight and wing area (r = 0.02, P = 0.83), all other parameters were significantly correlated either negatively or positively (Table 2). For instance, wing loading was positively correlated with adult weight (r = 0.50, P < 0.0001) and negatively correlated with wing width (r = − 0.43, P < 0.0001), wing length (r = − 0.46, P < 0.0001), centroid size (r = − 0.50, P < 0.0001), and wing area (r = − 0.60, P < 0.0001).

Fig. 2

Wing size and adult weight of S. calcitrans is significantly affected by larval density. a Boxplots depicting variation in wing centroid size (i), wing length (ii), wing width (iii), wing area (iv), adult weight (v), and wing loading (vi) across the different larval densities. Boxplot whiskers indicate ± 1.5 interquartile range limits. Boxplots with different letters show significant differences as grouped by ANOVA tests followed by SNK post-hoc tests (P < 0.05, n = 30). b Principal components biplot showing the similarities (or dissimilarities) existing among flies reared at different densities (5, 15 and 25)

Table 2 Correlation matrix between body weight and wing size parameters of S. calcitrans reared at different densities. Values above the diagonal represent P-values from Pearson correlation tests; values below the diagonal represent correlation coefficient. P-value in bold is not significant

Wing shape parameters

The wing shape of S. calcitrans significantly differed between the larval densities (MANOVA, F(52,1872) = 2.26; Pillai’s trace = 0.99; P = 0.0059). The thin plate spline (Fig. 3a) showed the variation in expansions and contractions in wing vein intersections (landmarks) of flies emerged from densities of 5, 15 and 25. For instance, in flies from densities of 5 and 15, the landmarks 4, 5 and 6 underwent expansion movement, while in the flies from a density of 25, the same landmarks contracted. Additionally, in flies from densities of 5 and 25, expansions in the landmarks 12, 13 and 14 were more pronounced (red coloured) compared to those of flies from a density of 15 (yellow coloured). The landmarks 7, 8, 9 and 10 expanded in wings emerged from a density of 15, while in fly wings from densities of 5 and 25, these landmarks contracted.

Fig. 3

Wing shape of S. calcitrans is significantly affected by larval density. a Thin plate spline deformation grids modelling the difference of S. calcitrans wing shape across the different larval densities. The number displayed on each grid represents the landmark positions. Yellow to orange-red colours indicate landmark expansions, light to dark-blue indicates landmark contraction. b Scatter plot showing the difference in wing shapes of S. calcitrans individuals reared at different larval densities along the first two canonical variate axes [CV1 (82.56%) and CV2 (17.44%)] with 90% confidence ellipses. c Discrimination histograms comparing the S. calcitrans wing shape between pairs of larval densities

Canonical variate analysis discriminated flies emerged from each density based on wing shape (Fig. 3b). The two first dimensions accounted for 100% of the total shape variation (CV1 = 82.56% and CV2 = 17.44%), and clustered wing shapes in three distinct groups based on the three larval densities. Pairwise comparisons using discriminant analysis with Mahalanobis distances revealed a highly significant difference in S. calcitrans wing shapes (Fig. 3c; Table 3; permutation test, 10000 replicates, P < 0.0001). When Procrustes distances were used, we found that wing shape of flies reared from densities of 15 and 25 were similar (P = 0.16). Regression of Procrustes coordinates on centroid size between densities was significant (permutation test with 10000 rounds, P = 0.008), with allometry explaining 2.97% of the total shape variation.

Table 3 Difference in the shape of right wings from S. calcitrans reared at a density 5, 15 and 25. P-values (above the diagonal); distances between populations (below the diagonal). P < 0.05 denotes a significant difference

Effect of developmental substrate on the wing geometry of S. calcitrans

Wing size parameters

Wing centroid size (Fig. 4a.i; F(3,122) = 39.24, P < 0.0001), wing length (Fig. 4a.ii; F(3,122) = 34.9, P < 0.0001), wing width (Fig. 4a.iii; F(3,122) = 46.98, P < 0.0001), wing area (Fig. 4a.iv; F(3,122) = 31.02, P<0.0001), adult weight (Fig. 4a.v; F(3,122) = 140.4, P < 0.0001) and the wing loading (Fig.; F(3,122) = 67.79, P < 0.0001) of S. calcitrans individuals reared from various animal dung differed significantly. All these parameters were highest in flies reared on donkey and sheep dung. The principal components analysis differentiated flies emerged from the different herbivore dung (Fig. 4b). The first dimension accounted for 55.5% of the variation in wing shape and were highly correlated with adult weight. The second dimension explained 19.4% of the total variation and was highly associated with wing length. Other than a lack of correlation between wing loading and wing length (r = 0.16, P = 0.08), all wing size parameters were positively correlated (Table 4).

Fig. 4

Wing size and adult weight of S. calcitrans significantly varies across different developmental substrates. a Boxplots depicting variation in wing centroid size (i), wing length (ii), wing width (iii), wing area (iv), adult weight (v), and wing loading (vi) across the different dung substrates. The limit of each boxplot whiskers represents the minimum and maximum of all the data. Boxplots with different letters depict significant differences as grouped by ANOVA tests followed by SNK post-hoc tests (P < 0.05, n = 30). b Principal components biplot showing the similarities (or dissimilarities) existing among flies reared from different dung types

Table 4 Correlation matrix between body weight and wing size parameters of S. calcitrans reared from different developmental substrates. Values above the diagonal represent P-values from Pearson correlation tests; values below the diagonal represent correlation coefficient. P-value in bold is not significant

Wing shape parameter

Developmental substrate significantly affected the wing shape of S. calcitrans (MANOVA, F(78,3172) = 4.07, Pillai’s trace = 1.17, P < 0.0001). This is clearly illustrated in the thin plate spline deformation grid (Fig. 5a). For instance, in flies emerged from camel dung (Fig. 5a.i), landmarks 12, 13, 14 and 15 underwent expansion movement, whereas in flies from cow (Fig. 5a.ii), donkey (Fig. 5a.iii) and sheep (Fig. 5a.iv) dung, the same landmarks contracted. Canonical variate analysis (Fig. 5b) and discriminant analysis (Fig. 5c) separated S. calcitrans wing shapes according to the dung in which flies developed. The first two dimensions of the canonical variate analysis explained 85.06% of the total S. calcitrans wing shape variation (Fig. 5b; CV1 = 66.06% and CV2 = 16%). All pairwise permutation tests performed with Mahalanobis distances revealed that the shape of S. calcitrans wings diverged significantly when reared from the different animal dung (Table 5; 10,000 rounds, P < 0.0001). With Procrustes distance estimators, we obtained a non-significant difference in wing shapes only in flies emerged from camel and cow dung (P = 0.2). In the allometry test, the centroid size had a significant effect on wing shape (10,000 rounds of permutation tests, P = 0.0015), with a variance prediction of 2.45%.

Fig. 5

Wing shape in S. calcitrans is significantly affected by developmental substrate. a Thin plate spline (TPS) deformation grids modelling the difference of S. calcitrans wing shape across the different larval densities. The number displayed on each grid represents the landmark positions. Yellow to orange-red colours indicate landmark expansions, light to dark-blue indicates landmark contraction. b Scatter plot showing the difference in wing shapes of S. calcitrans individuals reared from different developmental substrates along the first two canonical variate axes [CV1 (66.06%) and CV2 (19)] with 90% confidence ellipses. c Discrimination histograms comparing S. calcitrans wing shape between developmental substrate pairs

Table 5 Difference in the shape of right wings from S. calcitrans reared from various animal dung. P-values (above the diagonal); distances between populations (below the diagonal). P <  0.05 denotes significant difference


Our study showed that the size and shape of S. calcitrans wings exhibit a plastic response to larval density and the quality of the larval substrate. The study also demonstrated the power of the wing landmark-based geometric approach for studying phenotypic plasticity.

We showed that S. calcitrans wing size parameters (centroid size, length, width, area, and loading) are affected by larval density (5, 15 and 25) and substrate type (camel, cow, donkey and sheep dung). This indicates the effect of larval developmental conditions on adult wing size. Our study is consistent with previous work on Ae. aegypti that revealed the influence of larval density and substrate quality on wing size [47, 48]. The wing size variation obtained in our study may be due to variability in nutrients. Baleba et al. [39] previously determined that donkey and sheep dung, from which emerging flies had the largest wing size, had higher concentrations of specific micronutrients (nitrogen, phosphorous, potassium and zinc) in comparison with camel and cow dung. Furthermore, competition for limited nutrients by S. calcitrans larvae may also affect wing size. Dutra et al. [49] found that under high larval density, Wolbachia-uninfected Ae. aegypti presented reduced wing size (centroid size) and a lower body glucose concentration.

Wing sizes are closely related to the flight capacity in insects. Individuals with longer wings are better at flying compared to those with shorter wings [14]. Long wings favour wider variation in speed and long flight duration [50]. We obtained longer wings in S. calcitrans reared in a group of five or when reared from donkey and sheep dung. Long wings allow insects to fly at a great speed for a long period of time and cover a large area [51, 52]. For instance, released Ae. aegypti with larger wings are more successful in host-seeking and oviposition site location [53,54,55]. We predict that S. calcitrans developed under lower density conditions are more likely to have good flight capacity. Also, developmental substrates such as donkey and sheep dung may increase the efficiency of S. calcitrans flight. This result supports Baleba et al. [39], who found that dung types preferred by female S. calcitrans for oviposition were best for offspring growth and development. In this case, preferred substrates lead to potential adult fitness benefits associated with wing size. Low larval density leads to larger adults with long wings due to higher resource availability for growth and nutrient storage in the larval stage. Similarly, development of larvae in preferred dung types would lead to larger adults with larger wings. Such changes may influence dispersal, mating, and vector competency of S. calcitrans, rendering this fly more capable of spreading pathogens.

Our study also showed that larval density and substrate type affected S. calcitrans wing shape. The thin-plate spline analyses showed that most of the shape changes (landmark movements) occurred on the radial (landmarks 11, 12, 13, 14 and 15) and medial (2, 3, 7, 8, 9 and 10) portions of the S. calcitrans wing. Oguz et al. [56] observed the same variation in the radial portion of the wing of Phlebotomus tobbi Adler & Theodore. Pieterse et al. [59] found variation in wing landmarks located at the costal, sub-costal and radial veins of Bactrocera dorsalis (Hendel) and Ceratitis capitata (Wiedemann) reared from nectarine, plum, pear and apple. According to Wootton et al. [57] and Shimmi et al. [58], the radial and the medial portions of insect wings play a critical function in the aerodynamics of insect flight. Wootton [60] suggested such changes may influence the wing strength, beat pattern and ultimately the dispersal potential of a fly. Therefore, the wing shape deformation observed here may affect the flight performance of S. calcitrans, their ability to find a host for a blood meal, and consequently, vectorial capacity. Several studies with no emphasis on wing morphology have already demonstrated the indirect effect of larval density and food quality on vector competence. For instance, in Ae. albopictus, a greater dissemination rate of Sindbis virus by the adult is the consequence of high levels of competition experienced by the larvae [61]. In Anopheles stephensi Liston, larvae developed in a nutritious substrate are more likely to transmit the human malaria parasite, Plasmodium falciparum, than those developed in a substrate with a poor nutritional value [62]. The discriminant factors on which the differentiation between flies reared from different density or dung type was based were not free of some allometric effects. In other words, wing size contributed significantly to wing shape variation. However, in the case of both larval density and developmental substrate, less than 3% of variation in wing shape was attributed to size. Such low residual variation indicates that changes in the relative position of landmarks as wing size increases are minimal [63].


This study highlights the effect of larval density and developmental substrate of wing size and wing shape of S. calcitrans using the landmark-based geometric morphometric method. The method satisfactorily discriminated S. calcitrans emerged from different larval densities and substrates based on the size and the shape of their wings. While there was a significant effect of size variation on variation in shape, but this accounted for less than 3% of variation. Future studies on flight performance of S. calcitrans as well as their vectorial capacity in pathogen transmission when reared under different larval conditions are required. However, our results demonstrate a role for larval density and developmental substrate to influence wing size and to some extent wing shape, which might have a significant effect on flight and dispersal of adult S. calcitrans.



Analysis of variance


Multivariate analysis of variance






  1. 1.

    Mousseau T. The adaptive significance of maternal effects. Trends Ecol Evol. 1998;13:403–7.

    CAS  Article  Google Scholar 

  2. 2.

    Heard SB. Imperfect oviposition decisions by the pitcher plant mosquito (Wyeomyia smithii). Evol Ecol. 1994;8:493–502.

    Article  Google Scholar 

  3. 3.

    Wong J, Morrison AC, Stoddard ST, Astete H, Chu YY, Baseer I, et al. Linking oviposition site choice to offspring fitness in Aedes aegypti: consequences for targeted larval control of dengue vectors. PLoS Negl Trop Dis. 2012;6:e1632.

    Article  Google Scholar 

  4. 4.

    Jirakanjanakit N, Dujardin J-P. Discrimination of Aedes aegypti (Diptera: Culicidae) laboratory lines based on wing geometry. Southeast Asian J Trop Med Public Health. 2005;36:4.

    Google Scholar 

  5. 5.

    Wolf JB, Wade MJ. What are maternal effects (and what are they not)? Philos Trans R Soc Lond B Biol Sci. 2009;364:1107–15.

    Article  Google Scholar 

  6. 6.

    West-Eberhard MJ. Developmental plasticity and the origin of species differences. Proc Natl Acad Sci USA. 2005;102:6543–9.

    CAS  Article  Google Scholar 

  7. 7.

    Dogan M, Gunay F, Puggioli A, Balestrino F, Oncu C, Alten B, et al. Establishment of a satellite rearing facility to support the release of sterile Aedes albopictus males. I. Optimization of mass rearing parameters. Acta Trop. 2016;159:62–8.

    Article  Google Scholar 

  8. 8.

    Whitman DW, Agrawal AA. What is phenotypic plasticity and why is it important? In: Whitman DW, Ananthakrishnan TN, editors. Phenotypic plasticity of insects. Boca Raton: CRC Press; 2009. p. 1–63.

    Google Scholar 

  9. 9.

    Johansson F, SöDerquist M, Bokma F. Insect wing shape evolution: independent effects of migratory and mate guarding flight on dragonfly wings: dragonfly wing shape evolution. Biol J Linn Soc. 2009;97:362–72.

    Article  Google Scholar 

  10. 10.

    Lin X, Xu Y, Jiang J, Lavine M, Lavine LC. Host quality induces phenotypic plasticity in a wing polyphenic insect. Proc Natl Acad Sci USA. 2018;115:7563–8.

    CAS  Article  Google Scholar 

  11. 11.

    Fraimout A, Jacquemart P, Villarroel B, Aponte DJ, Decamps T, Herrel A, et al. Phenotypic plasticity of Drosophila suzukii wing to developmental temperature: implications for flight. J Exp Biol. 2018;221:jeb166868.

    Article  Google Scholar 

  12. 12.

    Soto IM, Carreira VP, Soto EM, Hasson E. Wing morphology and fluctuating asymmetry depend on the host plant in cactophilic Drosophila. J Evol Biol. 2008;21:598–609.

    CAS  Article  Google Scholar 

  13. 13.

    Jong Z-W, Kassim NFA, Naziri MA, Webb CE. The effect of inbreeding and larval feeding regime on immature development of Aedes albopictus. J Vector Ecol. 2017;42:105–12.

    Article  Google Scholar 

  14. 14.

    DeVries PJ, Penz CM, Hill RI. Vertical distribution, flight behaviour and evolution of wing morphology in Morpho butterflies: wing evolution in Morpho butterflies. J Anim Ecol. 2010;79:1077–85.

    CAS  Article  Google Scholar 

  15. 15.

    Betts CR, Wootton RJ. Wing shape and flight behaviour in butterflies (Lepidoptera: Papilionoidea and Hesperioidea): a preliminary analysis. J Evol Biol. 1988;138:271–88.

    Google Scholar 

  16. 16.

    Takken W, Klowden MJ, Chambers GM. Effect of body size on host seeking and blood meal utilization in Anopheles gambiae sensu stricto (Diptera: Culicidae): the disadvantage of being small. J Med Entomol. 1998;35:639–45.

    CAS  Article  Google Scholar 

  17. 17.

    Yoshioka M, Couret J, Kim F, McMillan J, Burkot TR, Dotson EM, et al. Diet and density dependent competition affect larval performance and oviposition site selection in the mosquito species Aedes albopictus (Diptera: Culicidae). Parasit Vectors. 2012;5:225.

    Article  Google Scholar 

  18. 18.

    Benelli G, Donati E, Romano D, Ragni G, Bonsignori G, Stefanini C, et al. Is bigger better? Male body size affects wing-borne courtship signals and mating success in the olive fruit fly, Bactrocera oleae (Diptera: Tephritidae): size and mating success in B. oleae. Insect Sci. 2016;23:869–80.

    Article  Google Scholar 

  19. 19.

    Trotta V, Duran Prieto J, Battaglia D, Fanti P. Plastic responses of some life history traits and cellular components of body size in Aphidius ervi as related to the age of its host Acyrthosiphon pisum: wing shape variation in A. ervi. Biol J Linn Soc. 2014;113:439–54.

    Article  Google Scholar 

  20. 20.

    Motoki MT, Suesdek L, Bergo ES, Sallum MAM. Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions. Infect Genet Evol. 2012;12:1246–52.

    Article  Google Scholar 

  21. 21.

    Gómez GF, Márquez EJ, Gutiérrez LA, Conn JE, Correa MM. Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation. Acta Trop. 2014;135:75–85.

    Article  Google Scholar 

  22. 22.

    Changbunjong T, Sumruayphol S, Weluwanarak T, Ruangsittichai J, Dujardin J-P. Landmark and outline-based geometric morphometrics analysis of three Stomoxys flies (Diptera: Muscidae). Folia Parasitol. 2016;63:037.

    Article  Google Scholar 

  23. 23.

    Wilke ABB, Christe R, Multini LC, Vidal PO, Wilk-da-Silva R, de Carvalho GC, et al. Morphometric wing characters as a tool for mosquito identification. PLoS ONE. 2016;11:0161643.

    Google Scholar 

  24. 24.

    Vidal PO, Suesdek L. Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti. Infect Genet Evol. 2012;12:591–6.

    Article  Google Scholar 

  25. 25.

    Virginio F, Oliveira Vidal P, Suesdek L. Wing sexual dimorphism of pathogen-vector culicids. Parasit Vectors. 2015;8:159.

    Article  Google Scholar 

  26. 26.

    Yeap H, Axford JK, Popovici J, Endersby NM, Iturbe-Ormaetxe I, Ritchie SA, et al. Assessing quality of life-shortening Wolbachia-infected Aedes aegypti mosquitoes in the field based on capture rates and morphometric assessments. Parasit Vectors. 2014;7:58.

    Article  Google Scholar 

  27. 27.

    Sendaydiego JP, Demayo CG. Describing variations in wing shapes of Anopheles flavirostris detected positive and negative of filaria using relative warp and Euclidean distance matrix analysis. Int J Mosq Res. 2015;2:09–13.

    Google Scholar 

  28. 28.

    Kitthawee S, Rungsri N. Differentiation in wing shape in the Bactrocera tau (Walker) complex on a single fruit species of Thailand. Scienceasia. 2011;37:308–13.

    Article  Google Scholar 

  29. 29.

    Sandoval Ramirez CM, Nieves Blanco EE, Gutiérrez Marin R, Jaimes Mendez DA, Rodríguez NO, Otálora-Luna F, et al. Morphometric analysis of the host effect on phenotypical variation of Belminus ferroae (Hemiptera: Triatominae). J Entomol. 2015;2015:1–12.

    Google Scholar 

  30. 30.

    Alves VM, Moura MO, de Carvalho CJB. Wing shape is influenced by environmental variability in Polietina orbitalis (Stein) (Diptera: Muscidae). Rev Bras Entomol. 2016;60:150–6.

    Article  Google Scholar 

  31. 31.

    Viscosi V, Cardini A. Leaf morphology, taxonomy and geometric morphometrics: a simplified protocol for beginners. PLoS ONE. 2011;6:e25630.

    CAS  Article  Google Scholar 

  32. 32.

    Lorenz C, Almeida F, Almeida-Lopes F, Louise C, Pereira SN, Petersen V, et al. Geometric morphometrics in mosquitoes: what has been measured? Infect Genet Evol. 2017;54:205–15.

    Article  Google Scholar 

  33. 33.

    Lehane MJ. The biology of blood-sucking in insects. New York: Cambridge University Press; 2005.

    Google Scholar 

  34. 34.

    Pitzer JB, Kaufman PE, Tenbroeck SH, Maruniak JE. Host blood meal identification by multiplex polymerase chain reaction for dispersal evidence of stable flies (Diptera: Muscidae) between livestock facilities. J Med Entomol. 2011;48:53–60.

    CAS  Article  Google Scholar 

  35. 35.

    Baldacchino F, Muenworn V, Desquesnes M, Desoli F, Charoenviriyaphap T, Duvallet G. Transmission of pathogens by Stomoxys flies (Diptera, Muscidae): a review. Parasite. 2013;20:26.

    Article  Google Scholar 

  36. 36.

    Carn VM. The role of dipterous insects in the mechanical transmission of animal viruses. BVA. 1996;152:377–93.

    CAS  Google Scholar 

  37. 37.

    Walker AR. Disease caused by arthropods. In: Sewell MMH, Brocklesby DW, editors. Handbook on animal diseases in the tropics. 4th ed. London: Bailliere Tindall; 1990. p. 1–385.

    Google Scholar 

  38. 38.

    Taylor DB, Moon RD, Mark DR. Economic impact of stable flies (Diptera: Muscidae) on dairy and beef cattle production. J Med Entomol. 2012;49:198–209.

    Article  Google Scholar 

  39. 39.

    Baleba SBS, Torto B, Masiga D, Weldon CW, Getahun MN. Egg-laying decisions based on olfactory cues enhance offspring fitness in Stomoxys calcitrans L (Diptera: Muscidae). Sci Rep. 2019;9:3850.

    Article  Google Scholar 

  40. 40.

    Suesdek L, Lorenz C. Evaluation of chemical preparation on insect wing shape for geometric morphometrics. Am J Trop Med Hyg. 2013;89:928–31.

    Article  Google Scholar 

  41. 41.

    Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.

    CAS  Article  Google Scholar 

  42. 42.

    Klingenberg CP, Barluenga M, Meyer A. Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry. Evolution. 2002;56:1909–20.

    Article  Google Scholar 

  43. 43.

    Hammer O, Harper DAT, Ryan PD. PAST: Paleontological Statistics Software: package for education and data analysis. Palaeontol Electron. 2001;4:4–9.

    Google Scholar 

  44. 44.

    Ribak G, Barkan S, Soroker V. The aerodynamics of flight in an insect flight-mill. PLoS One. 2017;12:e0186441.

    Article  Google Scholar 

  45. 45.

    R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.

  46. 46.

    Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics: computer program note. Mol Ecol Resour. 2011;11:353–7.

    Article  Google Scholar 

  47. 47.

    Jirakanjanakit N, Dujardin J-P. Discrimination of Aedes aegypti (Diptera: Culicidae) laboratory lines based on wing geometry. Southeast Asian J Trop Med Public Health. 2005;36:4.

    Google Scholar 

  48. 48.

    Stephens CR, Juliano SA. Wing shape as an indicator of larval rearing conditions for Aedes albopictus and Aedes aegypti (Diptera: Culicidae). J Med Entomol. 2012;49:927–38.

    CAS  Article  Google Scholar 

  49. 49.

    Yeap HL, Hoffmann AA, Ross PA, Endersby NM. Larval competition extends developmental time and decreases adult size of wMelPop Wolbachia-infected Aedes aegypti. Am J Trop Med Hyg. 2014;91:198–205.

    Article  Google Scholar 

  50. 50.

    Dutra HLC, Lopes da Silva V, da Rocha Fernandes M, Logullo C, Maciel-de-Freitas R, Moreira LA. The influence of larval competition on Brazilian Wolbachia-infected Aedes aegypti mosquitoes. Parasit Vectors. 2016;9:282.

  51. 51.

    Hoffmann AA, Ratna E, Sgrò CM, Barton M, Blacket M, Hallas R, et al. Antagonistic selection between adult thorax and wing size in field released Drosophila melanogaster independent of thermal conditions. J Evol Biol. 2007;20:2219–27.

    CAS  Article  Google Scholar 

  52. 52.

    Berwaerts K, Dyck HV, Aerts P. Does flight morphology relate to flight performance? An experimental test with the butterfly Pararge aegeria. Funct Ecol. 2002;16:484–91.

    Article  Google Scholar 

  53. 53.

    Davis AK, Holden MT. Measuring intraspecific variation in flight-related morphology of monarch butterflies (Danaus plexippus): which sex has the best flying gear? J Insects. 2015;2015:1–6.

    CAS  Article  Google Scholar 

  54. 54.

    Kay BH, Muir LE. Aedes aegypti survival and dispersal estimated by mark-release-recapture in northern Australia. Am J Trop Med Hyg. 1998;58:277–82.

    Article  Google Scholar 

  55. 55.

    Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, Clark GG, et al. Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg. 2005;72:209–20.

    Article  Google Scholar 

  56. 56.

    Oguz G, Kasap OE, Alten B. Wing morphology variations in a natural population of Phlebotomus tobbi Adler and Theodor 1930. J Vector Ecol. 2017;42:223–32.

    Article  Google Scholar 

  57. 57.

    Wootton RJ, Herbert RC, Young PG, Evans KE. Approaches to the structural modelling of insect wings. Philos Trans R Soc Lond Biol Sci B. 2003;358:1577–87.

    CAS  Article  Google Scholar 

  58. 58.

    Shimmi O, Matsuda S, Hatakeyama M. Insights into the molecular mechanisms underlying diversified wing venation among insects. Philos Trans R Soc Lond Biol Sci B. 2014;281:20140264.

    Google Scholar 

  59. 59.

    Pieterse W, Benítez HA, Addison P. The use of geometric morphometric analysis to illustrate the shape change induced by different fruit hosts on the wing shape of Bactrocera dorsalis and Ceratitis capitata (Diptera: Tephritidae). Zool Anz. 2017;269:110–6.

    Article  Google Scholar 

  60. 60.

    Wootton RJ. Support and deformability in insect wings. J Zool. 1981;193:447–68.

    Article  Google Scholar 

  61. 61.

    Alto BW, Lounibos LP, Higgs S, Juliano SA. Larval competition differentially affects arbovirus infection in Aedes mosquitoes. Ecology. 2005;86:3279–88.

    Article  Google Scholar 

  62. 62.

    Shapiro LLM, Murdock CC, Jacobs GR, Thomas RJ, Thomas MB. Larval food quantity affects the capacity of adult mosquitoes to transmit human malaria. Proc R Soc Biol Sci B. 2016;283:20160298.

    Article  Google Scholar 

  63. 63.

    Dujardin J-P. Morphometrics applied to medical entomology. Infect Genet Evol. 2008;8:875–90.

    Article  Google Scholar 

Download references


We thank Dr Joshua Njelembo Mbewe for the insightful advises and comments provided during the implementation of this work. We are also grateful to Pascal Atalor who generously assisted in the photography of wings.


We thank Deutscher Akademischer Austauschdienst (DAAD) for providing SBSB with a doctoral scholarship through the ARPPIS-DAAD scholarship programme. This work was supported by the IBCARP camel, grant no. DCI-FOOD/2014/ 346-739 - mainly by the European Union and Max Planck Institute for chemical ecology - icipe partner group to MNG. We also gratefully acknowledge the financial support for this research by the following organizations and agencies: Swedish International Development Cooperation Agency (Sida); UK Department for International Development (DFID); the Swiss Agency for Development and Cooperation (SDC); and the Kenyan Government. The views expressed herein do not necessarily reflect the official opinion of the donors.

Availability of data and materials

Data supporting the conclusions of this article are included within the article. All data collected during this study were archived in the Dryad data repository (

Authorsʼ contributions

SBSB: conceived the research idea, designed the experiment, collected and analysed the data, wrote the manuscript. DM: supervision and manuscript preparation. BT: supervision and manuscript preparation. CWW: supervision, contributed to experimental design, manuscript preparation and English editing. MNG: experimental design, supervision and manuscript preparation. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Author information



Corresponding author

Correspondence to Steve B. S. Baleba.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Baleba, S.B.S., Masiga, D., Torto, B. et al. Effect of larval density and substrate quality on the wing geometry of Stomoxys calcitrans L. (Diptera: Muscidae). Parasites Vectors 12, 222 (2019).

Download citation


  • Stomoxys calcitrans
  • Wing morphology
  • Larval density
  • Developmental substrate
  • Geometric morphometrics
  • Phenotypic plasticity