- Research article
- Open Access
Wing shape allometry and aerodynamics in calopterygid damselflies: a comparative approach
BMC Evolutionary Biology volume 13, Article number: 118 (2013)
Wing size and shape have important aerodynamic implications on flight performance. We explored how wing size was related to wing shape in territorial males of 37 taxa of the damselfly family Calopterygidae. Wing coloration was also included in the analyses because it is sexually and naturally selected and has been shown to be related to wing shape. We studied wing shape using both the non-dimensional radius of the second moment of wing area (RSM) and geometric morphometrics. Lower values of the RSM result in less energetically demanding flight and wider ranges of flight speed. We also re-analyzed previously published data on other damselflies and dragonflies.
The RSM showed a hump-shaped relationship with wing size. However, after correcting for phylogeny using independent contrast, this pattern changed to a negative linear relationship. The basal genus of the study family, Hetaerina, was mainly driving that change. The obtained patterns were specific for the study family and differed from other damselflies and dragonflies. The relationship between the RSM and wing shape measured by geometric morphometrics was linear, but relatively small changes along the RSM axis can result in large changes in wing shape. Our results also showed that wing coloration may have some effect on RSM.
We found that RSM showed a complex relationship with size in calopterygid damselflies, probably as a result of other selection pressures besides wing size per se. Wing coloration and specific behavior (e.g. courtship) are potential candidates for explaining the complexity. Univariate measures of wing shape such as RSM are more intuitive but lack the high resolution of other multivariate techniques such as geometric morphometrics. We suggest that the relationship between wing shape and size are taxa-specific and differ among closely-related insect groups.
Flight performance is a result of complex interactions between body morphology (i.e., wings, tails), behavior and the biological and physical environment [1–3]. In addition, body size has an important aerodynamic effect on flight performance [4, 5]. One of the main reasons is that Reynolds numbers change many orders of magnitude from small to large flying animals . At small body sizes, viscous effects due to small perturbations in the air are dissipated more rapidly, while at larger body sizes those small perturbations can result in stronger unsteadiness of the flow fields around wings . Another aspect of morphology that has a large impact on flight performance is wing shape [6–8]. For instance, in insects, long and slender wings are optimal at long duration flight, while short and broad wings are optimal at slow and agile flight [6, 7]. Moreover, a broad wing base allows a wider speed range  and a narrow wing tip allows less costly extensive flight . Other body traits such as the centers of body and wing mass are also very important in predicting flight performance [10, 11]. It is important to note that wing shape has an allometric component  and also a non-allometric component, since wings are complex structures and a simple proportional change of its shape with body size is unlikely [13–15]. For example, similar size clines may be achieved by different changes in wing proportions . Therefore, wing shape and wing size may function as independent components of wing morphology ( and references therein), although wing shape would be expected to improve flight performance for a given individual size . Hence, comparative studies on how wing shape changes with body size should be important for a better understanding of flight performance.
Apart from body size, other variables might also affect the optimal wing shape, for example flight behavior [2, 7, 18–21]. We should expect a similar wing shape within certain flight behaviors, and at the same time variation of this shape related to body size. In a comparative analysis between size and wing shape it is therefore important to take behavior into account or using species with the same behavior. Another variable that potentially could affect wing shape is the presence of ornaments on the wings that are used for sexual displays. These ornaments are typically positively sexually selected, either by male-male interactions and/or by female choice, and are condition dependent [22–25]. Wing shape may thus also play a role in ornament sexual signaling, since it might maximize both flight performance (e.g., flight style) and ornament display (e.g., size and shape of the ornament) [26–28]. For example, a large color patch at the tip of the wing should select for a broad wing tip shape  while a high agility (e.g., in territorial species) should select for a broad and short overall wing . Hence, the optimal shape might be somewhere in between these two. Therefore, the observed wing shape would differ from that predicted by aerodynamic theory per se.
In this study, we explore how wing shape is related to wing size (used as a proxy of body size), and discuss its implications on flight performance. For that purpose we compare males of damselfly taxa within the family Calopterygidae. These species display territorial behavior . Territorial flight is a common flight behavior in insects such as butterflies and damselflies [30, 31] and is associated with broad and short wings, which are expected to improve flight agility [6, 10, 32]. Many calopterygid males have wing coloration that is used for signaling in sexual selection processes [23, 33–39] and species discrimination [40–42], and it is also selected by bird predators [43, 44]. Although calopterygid males all have more or less short and broad wings, males still show a wide variation in wing morphology, which is related to wing size and coloration [28, 45, 46]. The present study differs from our previous works on wing shape in Calopterygidae [28, 45, 46] in that here we focus explicitly on the allometric effects on wing shape and the expected consequences on flight performance.
We studied wing shape using uni- and multivariate methods. First, we quantified wing shape using the non-dimensional radius of the second moment of wing area [10, 32, 47] (hereafter RSM). This parameter can be defined as wing area distribution along the wing axis and is proportional to the mean lift force and hence important for energetics of flight [10, 47]. Lower values involve a more basal distribution of the wing area and less energetically demanding flight with a wider range of available flight speeds [9, 10]. It has also been used as an estimate of flight agility , although this relationship might not be as straightforward. Second, we estimated wing shape using geometric morphometrics which allowed the analysis and visualization of shape in a very precise way . Since these two estimates of wing shape differ in how they capture shape and since wing shape affects flight performance, we also explored the relationship between them and with size.
In the present work, we specifically studied 37 taxa of male calopterygid damselflies that differ in wing shape, wing size and wing coloration patterns. We first explored the allometry of wing shape (using the RSM and geometric morphometrics) and the effect of wing coloration. We then compared wing shape captured by the RSM and by geometric morphometrics. We also investigated how the RSM and body length were related in a larger sample of both dragonflies and damselflies, by re-analyzing previously published data. Finally, we also performed independent contrasts for taking into account the evolutionary relationships of our study taxa.
RSM and size
Before correcting for phylogenetic effects, the RSM was related in a quadratic way to wing centroid size in both fore- and hindwings (Figure 1). Thus, the smallest and largest taxa showed the lowest values of the RSM. The quadratic term was significant both for fore- (linear regression: R2 = 0.013, F1,35 = 0.457, P = 0.503; quadratic regression: R2 = 0.424, F2,34 = 12.520, P < 0.001) and hindwings (linear regression: R2 = 0.151, F1,35 = 6.246, P = 0.017; quadratic regression: R2 = 0.435, F2,34 = 13.094, P < 0.001). Our study taxa were sorted in six different coloration groups according to the position, color and extension of wing coloration . The quadratic term was mainly driven by the coloration group 1, i.e. with a wing spot located at the wing base: the genera Hetaerina and Archineura (Figure 1). Using geometric morphometrics, we visualized how the RSM was related to wing shape and wing centroid size in our study taxa. For forewings, shape changed from a long slender wing with a broad base at small size to a wing with a broad tip and narrower base at intermediate sizes, then changing towards a longer wing with a broad base and a more slender tip at larger sizes (Figure 1). Hindwings showed a similar qualitative pattern. Thus, overall wing shape did not change in a merely proportional manner with wing centroid size (Figure 1).
The RSM showed a linear relationship with wing shape estimated from geometric morphometric methods as shape scores, both for fore- and hindwings (Figure 2). This relationship was supported by a MANCOVA on the shape components (obtained using geometric morphometrics) (fixed factor, fore- and hindwings: Wilk’s lambda = 0.157, F20,52 = 13.984, P < 0.001; covariate, RSM: Wilk’s lambda = 0.035, F20,52 = 72.692, P < 0.001). Nevertheless, for similar values of the RSM very different shapes can be obtained, partially as a result of differences in wing size (Figures 1 and 2).
We also tested the relationship between the RSM and body length in a larger size range of dragonflies and damselflies, by re-analyzing previously published data . The Odonata groups (calopterygid damselflies, non-calopteryid damselflies and dragonflies) differed significantly in the RSM, but body length did not show a significant effect in the model (Odonata group: F2,104 = 812.344, P < 0.001; Body length: F1,104 = 0.552, P = 0.459). Hence, we found no linear or hump-shaped relationship between the RSM and size when using a larger size range. Calopterygid damselflies showed intermediate values of the RSM with respect to non-calopterygid damselflies and dragonflies (Figure 3).
After accounting for phylogenetic non-independence among our study taxa, a visual inspection of the relationship between the RSM vs. wing centroid size revealed that a quadratic term was no longer present (Figure 4). Instead, a negative linear pattern was observed (forewings: R2 = 0.194, F1,35 = 8.450, P = 0.006; hindwings: R2 = 0.270, F1,35 = 12.954, P < 0.001). Therefore, the RSM tended to decrease with wing centroid size. We also note that at smaller sizes, the values of the RSM were more constrained, while at larger sizes there was more variation (Figure 4), i.e., at large size different wing shapes were observed.
RSM and wing coloration
In a previous study using 36 taxa and geometric morphometrics, we showed that the coloration groups differed in the non-allometric component of fore- and hindwing shape before correcting for phylogeny . Moreover, the non-allometric component of hindwing shape remained significant after correcting for the phylogenetic relationships . For comparative purposes, we ran a similar approach for the RSM. Before correcting for phylogeny, we did not find significant differences in the RSM among the coloration groups for forewings and marginally non-significant differences for hindwings (Table 1). A phylogenetic ANOVA on the non-allometric component of the RSM for hindwings revealed non-significant differences among the coloration groups (phylogenetic-P = 0.199).
Our study showed a quadratic relationship between the RSM and wing centroid size, implying that small and large species have lower values of this parameter compared to intermediate-sized species. This result was somewhat surprising since a linear relationship would be expected from aerodynamic theory , and in a previous study we found a linear relationship between wing shape (studied with geometric morphometrics) and wing centroid size . However, in the previous study  one of the largest species (genus Archineura) was not included. When we analyzed a larger size range of Odonata, we did not find the quadratic term, nor a linear relationship between the RSM and body length. These results suggest that our quadratic relationship is taxa-specific.
When phylogeny was taken into account, the RSM tended to decrease with wing centroid size. The change in this relationship after correcting for phylogenetic effects suggests a phylogenetic signal, i.e., that more closely related species tend to resemble each other more than a randomly chosen species sampled from the phylogeny. In fact, a significant phylogenetic signal of wing shape is present in this family of damselflies . We also found that at small size, the values of the RSM were more constrained than at large size. All together, these results suggest that shape changes along the size axis can be only partially related to the aerodynamic demands of flight at different sizes and more comparative studies in other groups are needed to explore if this pattern is general.
The linear negative relationship between the RSM and size was a result from small species showing a higher distribution of the wing area towards the wing apex, whereas large species showed a concentration of the wing area towards the wing base. A tendency to a higher concentration of wing area towards the wing base in the large species not only promotes less-energetically demanding flight [10, 47], it also results in a higher range of available flight speeds [49, 50], which would be beneficial for avoiding predators and/or increasing prey captures [51–53]. The wing shape shown at smaller sizes is expected to be more energetically costly in terms of flight energetics. A possible explanation to this costly design might be the need to display during territorial behaviors. For instance, butterfly Batesian mimics were shown to pay important aerodynamic costs for flying slower and with more agility, as distasteful mimics do, improving color signaling to predators .
The genus Hetaerina had the smallest wing centroid size and was a main driver in the quadratic relationship between the RSM and wing centroid size. Hetaerina is the basal genus in the whole Calopterygidae family (around 150 My)  (see also Methods) and this is likely to be the reason why after taking into account the phylogenetic relationships the quadratic term was no longer present. Hence, Hetaerina is probably on a different evolutionary pathway with regard to aerodynamic design. This is supported by the fact that this genus has not evolved a courtship flight, as in more recent genera such as Calopteryx. Hetaerina showed long, slender wings with a broad base and a slender tip. This shape was shared with the genus Archineura, which showed the largest wing centroid size and is a more recently diverged genus (less than 45 My) . Interestingly, both genera share the same wing color pattern. In fact, we only found marginally non-significant differences for the RSM of hindwings (controlling for size effects) which was driven by coloration group 1 (genera Hetaerina and Archineura). However, after correcting for phylogeny, that marginal difference was no longer present. Both genera might show some level of convergence in wing shape driven by a more effective display of the wing patch at the base of the wings. In this situation, a broader wing should provide a stronger signal of the color patch [27, 28]. Therefore, it is possible that the differences in the RSM among taxa are a consequence of a combination of wing size and wing coloration. However, since our results are correlative, it is still uncertain whether differences in wing shape are due to the presence of certain types of wing coloration or vice versa .
The lack of significant differences in the RSM of the hindwings among coloration groups after correcting for phylogeny contrasts with our previous study using geometric morphometrics, where we found significant differences . Thus, our results emphasize the need for using multivariate analyses of wing shape. In fact, although we showed that the RSM was linearly related to the shape scores, we obtained quite similar values of the RSM for very different wing shapes. For example, similar values of the RSM can be gained in a short and broad wing and a long and slender wing, given that the distribution of wing area is proportionally similar between both. Therefore, although the RSM has been broadly used before and its intuitive interpretation is very useful, the use of multivariate methods such as geometric morphometrics is much more precise and should, at least, be combined with other traditional measures of wing shape.
In this study we have shown that the relationships between wing shape and size are complex and taxon-specific, even within a group of species with similar flight behavior: territoriality. Moreover, the presence of wing coloration, a sexually and naturally selected trait, might influence the optimal wing shape within a trade-off between flight energetics and coloration display. The net selection pressures acting on flight performance, wing size, wing shape and body size probably differ among species, resulting in relationships that differ from the aerodynamic predictions alone. Experimental work in a comparative framework is needed to disentangle the role of wing morphology and size on flight performance.
Wing pictures of 338 males from 37 taxa of calopterygid damselflies (5–10 specimens per taxa) were collected from museum specimens (NCB Naturalis of Leiden and The Swedish Museum of Natural History, Stockholm), from colleagues or from our own samples (Table 2). Wings were either scanned in a flatbed scanner or photographed, together with a scale as a size reference. For the species of the genus Mnais Sélys, 1873, we used only territorial, colored morphs , so that only territorial taxa are included in our study.
Because wing shape is at least partially associated to wing coloration , taxa were grouped into different categories depending on the wing coloration position, color and extension (Table 2): 1) coloration only in the wing base; 2) yellow wing coloration covering most of the wing; 3) extensive dark wing coloration covering more than 85% of the wing area; 4) less extensive dark wing coloration, covering 30-85% and located either at the central part of the wing or at the wing apex; 5) very reduced wing coloration restricted to the wing tip and covering less than 20% of the wing area; and 6) no conspicuous coloration. Notice that in some species there are differences in coloration between fore- and hindwings (Table 2).
Wing shape analysis
We estimated shape using two different methods. First, we calculated the RSM. Second, we used geometric morphometrics techniques to study graphical changes in wing shape. The RSM has been used in the study of quasi-steady aerodynamics of hovering flight [10, 57]. This parameter can be defined as a quantitative measure of wing area distribution along the wing axis  and for a pair of wings is calculated as follows:
where ĉ is the normalized wing chord and is the non-dimensional radius equal to r/R, where r is the distance to the wing base on the chord c and R is wing length (Figure 5). This parameter is proportional to the mean lift force and thus lower values involve less energetically demanding flight [10, 47]. We note that the first three non-dimensional radii of moments of wing area are highly correlated [6, 10]. The RSM has been used before in dragonflies and damselflies [32, 50, 57] and thus we decided to use this specific radius for comparative purposes. We also note that in Odonata, fore- and hindwings do not form a common functional surface as in butterflies [49, 57]. Therefore, we calculated the RSM separately for fore- and for hindwings by following the protocol of Serrano-Meneses et al. , using ImageJ  and Microsoft Excel version 14.0 (Microsoft corp. 2010). We are aware that both pairs of wings may aerodynamically interact with each other, however, the aim of this study is rather descriptive and comparative.
We also used a geometric morphometrics approach to analyze wing shape variation. Although much more refined that the RSM, it does not provide an intuitive, quantitative estimate of wing shape as the RSM does, but it does provide an excellent visualization. In geometric morphometrics, shape is quantified from landmark coordinates after holding mathematically constant the effects of non-shape variation (position, orientation and scale) [48, 59, 60]. On the wing images, we digitized 12 landmarks and semi-landmarks (Figure 6) using tpsDig2 . Ten biologically homologous landmarks were located at the wing base or along the wing margin where some major veins terminate. To incorporate some aspects of the wing curvature, we also used two semi-landmarks (Figure 6). Only non-damaged wings were used. The landmarks and semi-landmarks were subjected to a Generalized Procrustes Analysis (GPA), where all specimens are translated to the origin, scaled to unit centroid size, and rotated to minimize the total sums-of-squares deviations of the landmark coordinates from all specimens to the average configuration . To minimize Procrustes distance between specimens during GPA, semi-landmarks were permitted to slide along their tangent directions [59, 63]. A single reference shape configuration (i.e., consensus wing) was obtained. The consensus wing was used for aligning all individual shape configurations and for computing the shape components (i.e., partial warps and the uniform component) in tpsRelw . Centroid size was computed for each wing and was logarithmic transformed in all analyses. Since centroid size had been previously shown to be highly correlated with body size in Calopterygidae [45, 65] we used it as a proxy for body size.
We used independent contrasts to correct for the phylogenetic non-independence of our study taxa. To obtain a phylogenetic tree that included all our study taxa, we re-analyzed previously published nucleotide sequences of the following genes: 18S, 5.8S, partial 28S rDNA, and of the spacers ITS1 and ITS2 [55, 66–69]. We used 70 taxa of the family Calopterygidae and three more taxa as outgroups (two non-calopterygids damselflies and one dragonfly) (see Supplementary Table 1 in ). We aligned the sequences using the ClustalW algorithm  in MEGA version 5 . Aligned sequences were analyzed by a Bayesian phylogenetic approach in the package BEAST version 1.7.1 , using a SRD06 model as the nucleotide substitution model, a relaxed molecular clock (uncorrelated lognormal) and a birth-death process as a tree prior. The Markov chain Monte Carlo sampling was run for 107 generations and logged every 1,000 generations. The consensus tree was pruned and was then used for the independent contrasts (Figure 7). The tree obtained was similar to previously published trees .
For all statistical tests, species were used as our level of replicate. Hence, we calculated for each species a mean value for each variable. We performed analyses both before and after accounting for the phylogenetic relationships among taxa. The statistical analyses were carried out in SPSS (IBM Corp.), unless otherwise stated.
We first performed a regression of the values of the RSM against wing centroid size. After visual inspection of the plot, we ran both a linear and a quadratic regression of the values of the RSM on wing centroid size. To visualize how wing shape was related to the variation in the values of the RSM and to wing centroid size, we computed thin-plate spline deformation grids in tpsSplin  using the shape components (partial warps and the uniform component). We did this by comparing the consensus wing shapes of certain taxa to the corresponding consensus wing shape, for both fore- and hindwings separately.
We also investigated how the RSM and the wing shape obtained by geometric morphometrics were related in our dataset, for fore- and hindwings separately. To do this, we first visualized this relationship by applying the approach developed by Drake and Klingenberg  which computes shape scores from multiple regressions of the Procrustes coordinates on a continuous variable (in our case, the RSM). The shape scores are the predicted shape variables in the regression including the residual variation in the direction of the shape space. Thus, the shape scores summarize all the landmark information into one single value. We ran this approach in R version 2.15 . To statistically support the obtained relationship, we also ran a MANCOVA using the shape components from geometric morphometrics (partial warps and the uniform component) as the dependent variables. The RSM was included as a covariate and wing (fore- and hindwing) was included as a fixed factor. Since the interaction term between the RSM and wing was non-significant, it was not included in the model.
To explore how the RSM was related to size over a larger size range in Odonata, we re-analyzed data published in Serrano-Meneses et al. , which contains 27 taxa of non-calopterygid damselflies (smaller than our study taxa), 54 taxa of dragonflies (covering a larger size range that our study taxa) and 27 taxa of calopterygid damselflies. These data were based on the mean of the values of the RSM between fore- and hindwings and size was measured as body length. We ran an ANCOVA to compare the values of the RSM among the three groups of Odonata, with body length as the covariate.
We also studied the relationship between the RSM and wing centroid size in our dataset using phylogenetic independent contrasts applying the method described by Felsenstein  for examining the correlated evolution of continuous traits. We did this for fore- and hindwings separately. We used the package PDAP:PDTREE version 1.16  for the Mesquite version 2.75 . After inspection of the plots of the absolute values of the standardized phylogenetically independent contrasts versus their standard deviations, we transformed the branch lengths of the tree to satisfy the requirements of the independent contrasts . All the obtained results were equivalent for both wings, so we only show the results for the exponentially transformed branch lengths. The plots of contrasts were standardized for both axes and also positivized for the X axis . After visual inspection of the relationships, least squared regressions were performed on the contrasts.
Finally, we inspected whether there were differences among the coloration groups in the values of the RSM. We have previously shown similar results for a dataset of 36 taxa of the same family, but using geometric morphometrics approaches . In the present study, we asked a similar question for comparative purposes between the two methods of measuring wing shape. We first explored the effect of wing coloration on wing shape without taking into account the phylogenetic effects. We ran a general linear model including the RSM as the response variable, wing centroid size and squared wing centroid size as covariates, and coloration group as a fixed factor, separately for fore- and hindwings. Non-significant interactions were removed one by one from the model. Second, we explored the effect of wing coloration on the RSM controlling for the phylogenetic non-independence. We removed the size effects on the RSM by performing a regression for each wing (fore- and hindwings). To correct for the quadratic relationship between the RSM and wing centroid size, we included in the model wing centroid size and squared wing centroid size as covariates. We saved the residuals from the regression, i.e., the non-allometric component of the RSM. We used this size-corrected RSM for examining the differences among the coloration groups from a phylogenetic perspective by running a phylogenetic ANOVA , for both fore- and hindwings separately. In this case, we used the R-package geiger , which returns a P-value based on Brownian motion simulations.
Norberg UM: How a long tail and changes in mass and wing shape affect the cost for flight in animals. Funct Ecol. 1995, 9: 48-54. 10.2307/2390089.
Berwaerts K, Van Dyck H, Aerts P: Does flight morphology relate to flight performance? An experimental test with the butterfly Pararge aegeria. Funct Ecol. 2002, 16: 484-491. 10.1046/j.1365-2435.2002.00650.x.
Clark CJ, Dudley R: Flight costs of long, sexually selected tails in hummingbirds. Proc R Soc B. 2009, 276: 2109-2115. 10.1098/rspb.2009.0090.
Sane SP: The aerodynamics of insect flight. J Exp Biol. 2003, 206: 4191-4208. 10.1242/jeb.00663.
Birch JM, Dickson WB, Dickinson MH: Force production and flow structure of the leading edge vortex on flapping wings at high and low Reynolds numbers. J Exp Biol. 2004, 207: 1063-1072. 10.1242/jeb.00848.
Betts CR, Wootton RJ: Wing shape and flight behaviour in butterflies (Lepidoptera: Papilionoidea and Hesperioidea): a preliminary analysis. J Exp Biol. 1988, 138: 271-288.
DeVries PJ, Penz CM, Hill RI: Vertical distribution, flight behaviour and evolution of wing morphology in Morpho butterflies. J Anim Ecol. 2010, 79: 1077-1085. 10.1111/j.1365-2656.2010.01710.x.
Förschler MI, Bairlein F: Morphological shifts of the external flight apparatus across the range of a passerine (Northern wheatear) with diverging migratory behaviour. PLoS One. 2011, 6: e18732-10.1371/journal.pone.0018732.
Wootton RJ, Newman DJS: Evolution, diversification, and mechanics of dragonfly wings. Dragonflies & damselflies. Model organisms for ecological and evolutionary research. Edited by: Córdoba-Aguilar A. 2008, New York: Oxford University Press, 261-274.
Ellington CP: The aerodynamics of hovering insect flight. II. Morphological parameters. Phil Trans R Soc Lond B. 1984, 305: 17-40. 10.1098/rstb.1984.0050.
Srygley RB: Locomotor mimicry in butterflies? The associations of positions of centres of mass among groups of mimetic, unprofitable prey. Phil Trans R Soc Lond B. 1994, 343: 145-155. 10.1098/rstb.1994.0017.
Stern DL, Emlen DJ: The developmental basis for allometry in insects. Development. 1999, 126: 1091-1101.
Debat V, Bégin M, Legout H, David JR: Allometric and nonallometric components of Drosophila wing shape respond differently to developmental temperature. Evolution. 2003, 57: 2773-2784.
Perrard A, Villemant C, Carpenter JM, Baylac M: Differences in caste dimorphism among three hornet species (Hymenoptera: Vespidae): forewing size, shape and allometry. J Evol Biol. 2012, 25: 1389-1398. 10.1111/j.1420-9101.2012.02527.x.
García Z, Sarmiento CE: Relationship between body size and flying-related structures in Neotropical social wasps (Polistinae, Vespidae, Hymenoptera). Zoomorphology. 2012, 131: 25-35. 10.1007/s00435-011-0142-z.
Huey RB, Gilchrist GW, Carlson ML, Berrigan D, Serra L: Rapid evolution of a geographic clines in size in an introduced fly. Science. 2000, 287: 308-309. 10.1126/science.287.5451.308.
Danforth BN: The evolution of hymenopteran wings - the importance of size. J Zool. 1989, 218: 247-276. 10.1111/j.1469-7998.1989.tb02536.x.
Wickman P-O: Sexual selection and butterfly design-a comparative study. Evolution. 1992, 46: 1525-1536. 10.2307/2409955.
Berwaerts K, Aerts P, Van Dyck H: On the sex-specific mechanisms of butterfly flight: flight performance relative to flight morphology, wing kinematics, and sex in Pararge aegeria. Biol J Linn Soc. 2006, 89: 675-687. 10.1111/j.1095-8312.2006.00699.x.
Fernández G, Lank DB: Variation in the wing morphology of western sandpipers (Calidris mauri) in relation to sex, age class, and annual cycle. Auk. 2007, 124: 1037-1046. 10.1642/0004-8038(2007)124[1037:VITWMO]2.0.CO;2.
Johansson F, Söderquist M, Bokma F: Insect wing shape evolution: independent effects of migratory and mate guarding flight on dragonfly wings. Biol J Linn Soc. 2009, 97: 362-372. 10.1111/j.1095-8312.2009.01211.x.
Andersson M: Sexual selection. 1994, Princeton: Princeton University Press
Contreras-Garduño J, Buzatto BA, Serrano-Meneses MA, Nájera-Cordero K, Córdoba-Aguilar A: The size of the red wing spot of the American rubyspot as a heightened condition-dependent ornament. Behav Ecol. 2008, 19: 724-732. 10.1093/beheco/arn026.
Legagneux P, Théry M, Guillemain M, Gomez D, Bretagnolle V: Condition dependence of iridescent wing flash-marks in two species of dabbling ducks. Behav Process. 2010, 83: 324-330. 10.1016/j.beproc.2010.01.017.
Rutowski RL, Nahm AC, Macedonia JM: Iridescent hindwing patches in the Pipevine Swallowtail: differences in dorsal and ventral surfaces relate to signal function and context. Funct Ecol. 2010, 24: 767-775. 10.1111/j.1365-2435.2010.01693.x.
Monteiro A, Brakefield PM, French V: The relationship between eyespot shape and wing shape in the butterfly Bicyclus anynana: A genetic and morphometrical approach. J Evol Biol. 1997, 10: 787-802. 10.1007/s000360050054.
Srygley RB: Locomotor mimicry in Heliconius butterflies: contrast analyses of flight morphology and kinematics. Phil Trans R Soc Lond B. 1999, 354: 203-214. 10.1098/rstb.1999.0372.
Outomuro D, Adams DC, Johansson F: The evolution of wing shape in ornamented-winged damselflies (Calopterygidae, Odonata). Evol Biol. 2013, 40: 300-309. 10.1007/s11692-012-9214-3.
Córdoba-Aguilar A, Cordero-Rivera A: Evolution and ecology of Calopterygidae (Zygoptera: Odonata): status of knowledge and research perspectives. Neotrop Entomol. 2005, 34: 861-879. 10.1590/S1519-566X2005000600001.
Kemp DJ, Wiklund C: Fighting without weaponry: a review of male-male contest competition in butterflies. Behav Ecol Sociobiol. 2001, 49: 429-442. 10.1007/s002650100318.
Marden JH, Waage JK: Escalated damselfly territorial contests are energetic wars of attrition. Anim Behav. 1990, 39: 954-959. 10.1016/S0003-3472(05)80960-1.
Serrano-Meneses MA, Córdoba-Aguilar A, Azpilicueta-Amorín M, González-Soriano E, Székely T: Sexual selection, sexual size dimorphism and Rensch’s rule in Odonata. J Evol Biol. 2008, 21: 1259-1273. 10.1111/j.1420-9101.2008.01567.x.
Grether GF: Intrasexual competition alone favors a sexually dimorphic ornament in the rubyspot damselfly Hetaerina americana. Evolution. 1996, 50: 1949-1957. 10.2307/2410753.
Hooper RE, Tsubaki Y, Siva-Jothy MT: Expression of a costly, plastic secondary sexual trait is correlated with age and condition in a damselfly with two male morphs. Physiol Entomol. 1999, 24: 364-369. 10.1046/j.1365-3032.1999.00152.x.
Siva-Jothy MT: Male wing pigmentation may affect reproductive success via female choice in a calopterygid damselfly (Zygoptera). Behaviour. 1999, 136: 1365-1377. 10.1163/156853999500776.
Rantala MJ, Koskimäki J, Taskinen J, Tynkkynen K, Suhonen J: Immunocompetence, developmental stability and wingspot size in the damselfly Calopteryx splendens L. Proc R Soc Lond B. 2000, 267: 2453-2457. 10.1098/rspb.2000.1305.
Siva-Jothy MT: A mechanistic link between parasite resistance and expression of a sexually selected trait in a damselfly. Proc R Soc Lond B. 2000, 267: 2523-2527. 10.1098/rspb.2000.1315.
Córdoba-Aguilar A: Wing pigmentation in territorial male damselflies, Calopteryx haemorrhoidalis: a possible relation to sexual selection. Anim Behav. 2002, 63: 759-766. 10.1006/anbe.2001.1974.
Svensson EI, Karlsson K, Friberg M, Eroukhmanoff F: Gender differences in species recognition and the evolution of asymmetric sexual isolation. Curr Biol. 2007, 17: 1943-1947. 10.1016/j.cub.2007.09.038.
Anderson CN, Grether GF: Interspecific aggression and character displacement of competitor recognition in Hetaerina damselflies. Proc R Soc Lond B. 2010, 277: 549-555. 10.1098/rspb.2009.1371.
Waage JK: Reproductive character displacement in Calopteryx (Odonata: Calopterygidae). Evolution. 1979, 33: 104-116. 10.2307/2407369.
Tynkkynen K, Rantala MJ, Suhonen J: Interspecific aggression and character displacement in the damselfly Calopteryx splendens. J Evol Biol. 2004, 17: 759-767. 10.1111/j.1420-9101.2004.00733.x.
Rantala MJ, Honkavaara J, Dunn DW, Suhonen J: Predation selects for increased immune function in male damselflies, Calopteryx splendens. Proc R Soc Lond B. 2011, 278: 1231-1238. 10.1098/rspb.2010.1680.
Svensson EI, Friberg M: Selective predation on wing morphology in sympatric damselflies. Am Nat. 2007, 170: 101-112. 10.1086/518181.
Outomuro D, Johansson F: The effects of latitude, body size, and sexual selection on wing shape in a damselfly. Biol J Linn Soc. 2011, 102: 263-274. 10.1111/j.1095-8312.2010.01591.x.
Outomuro D, Bokma F, Johansson F: Hind wing shape evolves faster than front wing shape in Calopteryx damselflies. Evol Biol. 2012, 39: 116-125. 10.1007/s11692-011-9145-4.
Weis-Fogh T: Quick estimates of flight fitness in hovering animals, including novel mechanisms for lift production. J Exp Biol. 1973, 59: 169-230.
Adams DC, Rohlf FJ, Slice DE: A field comes of age: Geometric morphometrics in the 21st century. Hystrix. in press
Wootton RJ: Functional morphology of insect wings. Annu Rev Entomol. 1992, 37: 113-140. 10.1146/annurev.en.37.010192.000553.
Wakeling JM: Odonatan wing and body morphologies. Odonatologica. 1997, 26: 35-52.
Hedenström A, Ellington CP, Wolf TJ: Wing wear, aerodynamics and flight energetics in bumblebees (Bombus terrestris): an experimental study. Funct Ecol. 2001, 15: 417-422. 10.1046/j.0269-8463.2001.00531.x.
Jantzen B, Eisner T: Hindwings are unnecessary for flight but essential for execution of normal evasive flight in Lepidoptera. PNAS. 2008, 105: 16636-16640. 10.1073/pnas.0807223105.
Combes SA, Crall JD, Mukherjee S: Dynamics of animal movement in an ecological context: dragonfly wing damage reduces flight performance and predation success. Biol Lett. 2010, 6: 426-429. 10.1098/rsbl.2009.0915.
Srygley RB: The aerodynamic costs of warning signals in palatable mimetic butterflies and their distasteful models. Proc R Soc Lond B. 2004, 271: 589-594. 10.1098/rspb.2003.2627.
Dumont HJ, Vanfleteren JR, De Jonckheere JF, Weekers PHH: Phylogenetic relationships, divergence time estimation, and global biogeographic patterns of Calopterygoid damselflies (Odonata, Zygoptera) inferred from ribosomal DNA sequences. Syst Biol. 2005, 54: 347-362. 10.1080/10635150590949869.
Tsubaki Y: The genetic polymorphism linked to mate-securing strategies in the male damselfly Mnais costalis Selys (Odonata: Calopterygidae). Popul Ecol. 2003, 45: 263-266. 10.1007/s10144-003-0162-8.
Wakeling JM, Ellington CP: Dragonfly flight. III. Lift and power requirements. J Exp Biol. 1997, 200: 583-600.
Rasband WS: ImageJ. 1997–2012, Bethesda: U. S. National Institutes of Health, http://imagej.nih.gov/ij/,
Bookstein FL: Morphometric tools for landmark data geometry and biology. 1991, Cambridge: Cambridge University Press
Rohlf FJ, Marcus LF: A revolution in morphometrics. Trends Ecol Evol. 1993, 8: 129-132. 10.1016/0169-5347(93)90024-J.
Rohlf FJ: tpsDig version 2.16. 2010, http://life.bio.sunysb.edu/morph/,
Rohlf FJ, Slice D: Extension of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990, 39: 40-59. 10.2307/2992207.
Gunz P, Mitteroecker P, Bookstein FL: Semilandmarks in three dimensions. Modern morphometrics in physical anthropology. Edited by: Slice DE. 2005, New York: Springer, 73-98.
Rohlf FJ: tpsRelw. Relative warps version 1.49. 2010, http://life.bio.sunysb.edu/morph/,
Sacchi R, Hardersen S: Wing length allometry in Odonata: differences between families in relation to migratory behaviour. Zoomorphology. 2013, 132: 23-32. 10.1007/s00435-012-0172-1.
Weekers PHH, De Jonckheere JF, Dumont HJ: Phylogenetic relationships inferred from ribosomal ITS sequences and biogeographic patterns in representatives of the genus Calopteryx (Insecta: Odonata) of the West Mediterranean and adjacent West European zone. Mol Phylogenet Evol. 2001, 20: 89-99. 10.1006/mpev.2001.0947.
Hayashi F, Dobata S, Futahashi R: Macro- and microscale distribution patterns of two closely related Japanese Mnais species inferred from nuclear ribosomal DNA, ITS sequences and morphology (Zygoptera: Calopterygidae). Odonatologica. 2004, 33: 399-412.
Dumont HJ, Vierstraete A, Vanfleteren JR: A molecular phylogeny of the Odonata (Insecta). Syst Entomol. 2010, 35: 6-18. 10.1111/j.1365-3113.2009.00489.x.
Guan Z, Han BP, Vierstraete A, Dumont HJ: Additions and refinements to the molecular phylogeny of the Calopteryginae Sl (Zygoptera: Calopterygidae). Odonatologica. 2012, 41: 17-24.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W - Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Drummond AJ, Rambaut A, Suchard M: BEAST version 1.7.1. Bayesian Evolutionary Analysis Sampling Trees. 2012, http://beast.bio.ed.ac.uk/Main_Page,
Rohlf FJ: tpsSplin. Thin-plate spline version 1.20. 2004, http://life.bio.sunysb.edu/morph/,
Drake AG, Klingenberg CP: The pace of morphological change: historical transformation of skull shape in St Bernard dogs. Proc R Soc Lond B. 2008, 275: 71-76. 10.1098/rspb.2007.1169.
R Development Core Team: R: A language and environment for statistical computing. Version 2.15.0. 2011, Vienna: R Foundation for Statistical Computing, http://www.r-project.org,
Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15. 10.1086/284325.
Midford PE, Garland T, Maddison WP: PDAP Package of Mesquite. Version 1.16. 2011
Maddison WP, Maddison DR: Mesquite: A modular system for evolutionary analysis. Version 2.75. 2011, http://mesquiteproject.org,
Garland T, Harvey PH, Ives AR: Procedures for the analysis of comparative data using phylogenetically independent contrasts. Syst Biol. 1992, 41: 18-32.
Garland T, Dickerman AW, Janis CM, Jones JA: Phylogenetic analysis of covariance by computer simulation. Syst Biol. 1993, 42: 265-292.
Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W: GEIGER: investigating evolutionary radiations. Bioinformatics. 2008, 24: 129-131. 10.1093/bioinformatics/btm538.
We are very grateful to M. Serrano-Meneses who kindly provided us with the protocol for calculating the RSM. We thank L. Faulks and W. Kunce for useful comments on the manuscript. We also thank K. D. B. Dijkstra for his support at the NCB Naturalis of Leiden and Gunvi Lindberg for her help at The Swedish Museum of Natural History in Stockholm. We also want to thank P. Brunelle, A. Córdoba-Aguilar, R. Futahashi, D. Halstead, M. Hämäläinen, I. Santoyo, G. Sims, Y. Tsubaki, H. Ubukata and X. Yu for their help in providing us with some of the taxa. This study has been supported by a postdoc position to D. Outomuro from the Spanish Ministry of Education. D. C. Adams was supported in part by NSF grant DEB-1118884 and F. Johansson was supported by The Swedish Research Council.
The authors declare that they have no competing interests.
DO coordinated sample collection, did data acquisition, analyzed data, interpreted the results and wrote the manuscript. DCA analyzed data, interpreted the results and revised the manuscript. FJ coordinated sample collection, interpreted the results and co-wrote the manuscript. All authors read and approved this paper.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Outomuro, D., Adams, D.C. & Johansson, F. Wing shape allometry and aerodynamics in calopterygid damselflies: a comparative approach. BMC Evol Biol 13, 118 (2013) doi:10.1186/1471-2148-13-118
- Geometric Morphometrics
- Non-dimensional Radius of the Second Moment of Wing Area
- Territorial Behavior
- Wing Coloration