Research article | Open | Published:
Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup
BMC Evolutionary Biologyvolume 9, Article number: 110 (2009)
Sexual dimorphism of body size has been the subject of numerous studies, but few have examined sexual shape dimorphism (SShD) and its evolution. Allometry, the shape change associated with size variation, has been suggested to be a main component of SShD. Yet little is known about the relative importance of the allometric and non-allometric components for the evolution of SShD.
We investigated sexual dimorphism in wing shape in the nine species of the Drosophila melanogaster subgroup. We used geometric morphometrics to characterise wing shape and found significant SShD in all nine species. The amount of shape difference and the diversity of the shape changes evolved across the group. However, mapping the divergence of SShD onto the phylogeny of the Drosophila melanogaster subgroup indicated that there is little phylogenetic signal. Finally, allometry accounted for a substantial part of SShD, but did not explain the bulk of evolutionary divergence in SShD because allometry itself was found to be evolutionarily plastic.
SShD in the Drosophila wing can evolve rapidly and therefore shows only weak phylogenetic structure. The variable contribution of allometric and non-allometric components to the evolutionary divergence of SShD and the evolutionary plasticity of allometry suggest that SShD and allometry are influenced by a complex interaction of processes.
Sexual dimorphism is one of the most striking and widespread sources of phenotypic variation in animals and plants and has therefore attracted considerable interest in evolutionary biology. The evolution of sex dimorphism has been extensively studied, but most studies have concerned dimorphism of size [1–4]. In contrast, sexual shape dimorphism (SShD) has been much less investigated.
Of those studies that considered SShD, most have discussed it as a diagnostic-trait for diverse purposes, such as sex identification or the analysis of ontogeny [5–14]. Only relatively few investigations have specifically considered the evolution of SShD, covering a wide range of study systems including the skull in primates [15–18], body proportions of lizards [19, 20], newts  or flies , the head shape of Chironomus larvae  and Lycium flowers . With few exceptions, these studies used only small numbers of species and no explicit phylogenetic framework for analyzing evolutionary change.
A factor that many studies have identified as playing a particularly important role for SShD is allometry, the relationship between size and shape [25–29]. Particularly in primatology, SShD has often been explained by allometric variation of shape and ontogenetic scaling [6, 17, 27], which raises the possibility that SShD may evolve primarily as a by-product of the evolution of sexual size dimorphism. To address this question, it is useful to separate an allometric component of SShD, for which size dimorphism accounts, from a non-allometric component of residual SShD from other sources. For instance, several studies in primates found that allometric scaling tends to account for much of SShD, but that non-allometric shape variation also contributes a significant part of SShD [16–18]. Similarly, allometric scaling accounts for most of SShD in some species of lizards , but other species show sexual dimorphism in body proportions without sexual size dimorphism . Altogether, these studies suggest that the role of allometry in SShD is more complex than previously assumed and therefore requires systematic investigation.
Here we address this issue by investigating the evolution of SShD and the role of allometry of the wing in the nine species of the Drosophila melanogaster subgroup. Drosophila is a long-standing model organism in evolutionary biology, and the melanogaster subgroup has been particularly well studied [31, 32]. The Drosophila wing is a model trait in evolutionary studies, which has a remarkable potential to evolve under selection [33, 34]. Because the intersections between wing veins provide landmarks that can be located precisely, shape variation is easy to quantify with the landmark-based methods of geometric morphometrics  and the Drosophila wing has therefore been widely used for morphometric studies [36–40]. SShD in the wing of Drosophila melanogaster has been found to be quite constant over latitudinal clines of wing shape on different continents, suggesting that SShD may be evolutionarily constrained . Allometry appears to be one of the evolutionary processes possibly involved in this SShD [36, 41], although this question has not been addressed directly. This study uses geometric morphometrics to examine SShD of the wing shape in all nine species of the melanogaster subgroup, to test for the presence of a phylogenetic signal in SShD and to study the role of allometry in its evolution.
Variation of wing size and shape
Size and shape of the wings were characterized by a set of 15 landmarks (Figure 1) and analyzed with the methods of geometric morphometrics. To quantify wing size, we used centroid size, which is a measure of the spread of landmarks around their centre of gravity . An analysis of variance indicates that both species and sex have significant effects on centroid size (Table 1). Moreover, the significant species × sex interaction indicates a divergence of sexual size dimorphism among species.
Wing shape was extracted from the landmark coordinate data with a generalized Procrustes fit . A multivariate analysis of variance (MANOVA) was used to test the effects of species, sex, size and their interactions on wing shape (Table 2). Species of the melanogaster subgroup differ in their wing shapes (significant main effect of species). The MANOVA results also indicate that there is sexual shape dimorphism; however, the marginal significance of the main effect of sex and the highly significant species × sex interaction suggest that sex dimorphism has diverged considerably among species. Allometry contributed to differences in wing shape (significant main effect of centroid size), but there appear to be differences in the allometric patterns among species and sexes (significant interactions of centroid size × sex, centroid size × species and centroid size × species × sex).
Variation in sexual shape dimorphism
To quantify the amount of SShD, we computed the Procrustes tangent distance  between the mean shapes of the males and females of each species. On the whole, these shape changes were fairly subtle. The magnitudes of total SShD differed among species, with D. orena displaying the least shape dimorphism and D. mauritiana being the most dimorphic species in the subgroup (Figure 2; numbers next to the diagrams for total SShD).
To visualize the shape changes associated with SShD, we graphed the shape changes from male to female average shapes for all nine species (Figure 2, left). Here we also describe the shape features associated with SShD as changes from the male to female shape. A common feature of SShD in most species is an overall narrowing of the distal part of the wing blade from male to female, particularly in the region between the L3 and L4 veins (landmarks 13 and 14; weaker in D. erecta) and to some extent also the distal ends of the L2 and L5 veins (landmarks 12 and 15). This general narrowing was accompanied by various species-specific shape changes from male to female: a posterior shift of landmark 5 (D. orena, D. santomea, and D. sechellia) and a distal shift of the posterior crossvein (D. mauritiana, D. teissieri and D. yakuba), of the anterior crossvein (D. erecta) or of both crossveins (D. melanogaster and D. simulans).
The SShD of different species were distinct from each other (Table 3): pairwise distances between the SShD vectors of different species ranged from 0.0078 (D. teissieri versus D. yakuba) to 0.0140 (D. mauritiana versus D. melanogaster). The differences between SShD of different species were therefore of a similar magnitude as the SShD in each species (Figure 2), indicating a considerable diversification. Bootstrap tests  against the null hypothesis of equal SShD vectors were all nominally significant (all P < 0.05; Table 3) and the false discovery rate  is controlled at the 5% level (i.e. it is expected that fewer than 5% of the significant tests are false positives).
The variation of SShD in the melanogaster subgroup was displayed graphically by a principal component analysis (Figure 3). None of the principal components (PCs) accounted for a very large proportion of the variation (25.56% for PC1), indicating that the variation of SShD diverged in many directions of the shape space. The confidence ellipses for the SShD of most species were fairly large (Figure 3), which suggests a considerable degree of uncertainty about the SShD for those species. The overlap of the confidence ellipses in the PC plots, despite the significant separation of all species from each other, points to the high-dimensional nature of the variation.
To reconstruct the evolutionary history of SShD, we mapped the SShD onto a phylogeny of the melanogaster subgroup (Figure 4) assembled from published sources [44–46]. The values of SShD corresponding to the internal nodes of the phylogeny were estimated by squared-change parsimony [47, 48] and the tree was graphed in the principal component plots [24, 49–51] so that reconstructed evolutionary trajectories of SShD can be visualized (Figure 3).
In the resulting graph (Figure 3), many terminal branches of the tree appeared to be fairly long, whereas the internal branches tended to be relatively short. This reflected the fact that some closely related species showed clearly different SShD (e.g. the pairs of sister species D. simulans and D. mauritiana, D. santomea and D. yakuba). Moreover, there were many sharp changes of direction in the reconstructed trajectories of evolutionary change and branches appeared to cross in all three PC plots in Figure 3, indicating the presence of homoplasy. Yet, the divergence of SShD also seemed to contain some phylogenetic structure, as the shape dimorphisms of the species belonging to the same complexes were often found in the same areas of the PC plots.
To test for the presence of a phylogenetic signal in SShD, we used a permutation procedure that is a multivariate generalization of the method previously used for scalar traits such as body size . This test simulated the null hypothesis of the complete absence of a phylogenetic signal by randomly permuting SShD vectors among the terminal nodes of the phylogenetic tree and evaluated the resulting tree length for each permutation. The tree length obtained by mapping of the original data onto the phylogenetic tree by weighted squared-change parsimony was 0.000392 (in units of squared Procrustes distance). The permutation procedure resulted in an equal or shorter tree length in nearly half the permutations (proportion 0.48, which corresponds to the P-value of the test) and therefore provides no evidence for a phylogenetic signal in SShD.
Allometry in SShD
To assess the role of allometry for SShD, we decomposed the total SShD for each species into allometric and non-allometric components (Figure 5). We used multivariate regression of shape on centroid size  as the general approach to analyze allometry. To obtain a common estimate for the allometry in both sexes, we ran a pooled within-sex regression in each species. The allometric component of SShD is the part predicted by the sexual size dimorphism (and is thus in the same direction in the shape tangent space as the regression vector), whereas the non-allometric component is the difference between the total SShD and the allometric component (Figure 5).
Because this approach assumes that the allometric regression vectors in males and females are parallel, we conducted a bootstrap test  of the null hypothesis that the regression vectors for males and females are the same. There was a significant and substantial difference between regression vectors for males and females in D. melanogaster (sum of squared differences 0.023; P = 0.022) and smaller differences in D. santomea (0.0045; P = 0.0094), D. orena (0.00096; P = 0.037), D. teissieri (0.0050; P = 0.023). The differences in the other species were small (0.0011 – 0060) and not statistically significant (P-values 0.12 – 0.55). Therefore, the separation of allometric and non-allometric components of SShD for D. melanogaster is likely to be unreliable, and should be interpreted with some caution in the remaining species as well.
The magnitudes of the allometric and non-allometric components of SShD tended not to add up to the magnitude of the total SShD (Figure 2), as would be expected if all are in the same direction of shape tangent space. Instead, the magnitudes of the allometric and non-allometric components were variable and often of a similar magnitude as the total SShD. In some species, the magnitudes of both components exceeded that of the total SShD (D. melanogaster, D. orena, D. teissieri), but these were among the species where the allometric regressions differed significantly between the sexes, so that this result may be related to poor estimates of a common allometric regression. In all species, however, the magnitudes of the total SShD and its components suggested that the direction of total SShD and the allometric component (and therefore also the non-allometric component) differ considerably.
This expectation was confirmed by inspection of the shape changes associated with SShD and its components (Figure 2). The allometric component of SShD was associated with a wide range of shape changes, which differed considerably from those of total SShD in most species. Allometric SShD tended to consist of a slight expansion of the proximal part of the wing and various changes in the distal part of the wing, which often included a distal shift of landmark 12. In contrast, non-allometric SShD did not include a distal shift of landmark 12 or only a slight one. Non-allometric SShD also showed variety of other shape changes, but there were no other shape features that were shared by most species and there was no apparent relation to the relatedness of species.
This study has yielded three main results concerning the evolution of sex dimorphism in Drosophila wing shape. First, we have shown that SShD has diverged among the nine species of the Drosophila melanogaster subgroup. Second, the evolution of SShD is little constrained by the phylogenetic structure of the subgroup. Third, allometry is an important component of sexual shape dimorphism in the wing but is not its main evolutionary driver. Here we present the implications of our results for the evolution of the Drosophila wing and sexual shape dimorphism.
Our study found substantial evolution of SShD in the Drosophila melanogaster subgroup. The MANOVA of shape revealed a clearly significant sex × species interaction (Table 2), which indicates that the effect of sex depends on the species. Similarly, most of the pairwise comparisons of SShD showed highly significant differences between species (Table 3). With a range of 0.0078 to 0.014 (in units of Procrustes distance, Table 3), the magnitudes of pairwise differences between the SShD of different species are therefore of a similar to the magnitudes of SShD, which ranged from 0.0077 to 0.017 (same units, Figure 2). This result is in marked contrast with the finding of Gilchrist et al. , who found SShD to be remarkably conserved across latitudinal clines of wing shape that had evolved separately on three continents. Those authors attributed the relative evolutionary conservatism of SShD to constraints by high genetic correlations between male and female wing traits , which would be expected to impede the evolution of differences between the sexes. As Gilchrist et al.  pointed out, such constraints by genetic correlations can be broken over longer-term evolution. Our results suggest that this may have happened in the evolution of the melanogaster subgroup. At the evolutionary scale of the entire clade, rather than a single species, SShD appears to be a trait that readily evolves.
The shape changes involved in SShD consist of the same features that are widely observed in variation within or between populations of Drosophila [36, 38, 40], such as proximal or distal shifts of the crossveins along the longitudinal veins or variation in the positions at which the L2 and L5 veins meet the wing margin. For SShD in different Drosophila species, shape features like these seem to be combined in various ways. The PCA of SShD variation among species suggested that variation is distributed over many dimensions rather than being concentrated primarily in one or a few dimensions (the PC1 accounted for 25.6% of the variation and the subsequent eigenvalues decreased gradually). Moreover, the scatter of SShD in shape space is irregular with no dominant patterns or strongly preferred directions (Figure 3). This suggests that evolutionary divergence may flexibly combine a set of wing shape features to make up the SShD for each species.
This impression of flexibility is further reinforced by the observation that there was little phylogenetic structure in the divergence of SShD. There appeared to be no association between relatedness of species and the difference of their SShD (Table 3). The mapping of SShD onto the phylogeny by squared-change parsimony indicated that changes tended to be small for internal branches and large for terminal branches (Figure 3), suggesting that much of the evolutionary divergence of SShD is separately derived for each species. This overall impression was confirmed by the permutation test that yielded a clearly non-significant result, consistent with the null hypothesis of the lack of any phylogenetic structure in SShD.
Allometry has been generally recognized to be a key factor for SShD [6, 13, 16–18, 27] and previous studies have demonstrated allometry in Drosophila wings [36, 41, 55]. Our study confirms that allometry is an important component in the evolution of SShD in the melanogaster subgroup; yet, the results also suggest that allometry is not a stringent constraint on the evolution of SShD. In all nine species, both the allometric and non-allometric components of SShD are fairly large (at least half the magnitude of the total SShD, Figure 2). On the one hand, this indicates that allometry can account for a considerable part of SShD, but on the other hand, it also shows that the directions of the total SShD and the allometric regression vector differed considerably from each other and that a substantial part of SShD therefore cannot be explained by allometry. Clearly, factors other than allometry are contributing substantially to SShD in the melanogaster subgroup. This is in agreement with studies in other organisms [13, 16–18] that also have found non-allometric components of SShD.
Not only is there a substantial contribution by non-allometric SShD, which suggests that allometry is not a severe constraint on the evolution of SShD, but the allometries themselves seem to be variable. The shape features associated with allometry differ from species to species (Figure 2) and the bootstrap test indicated that for some species even the allometries of the two sexes are not the same. This suggests that allometry itself can evolve, which may further erode its role as an evolutionary constraint.
Overall, SShD appears to be a trait that can readily respond to selection or may evolve by drift. The absence of a detectable phylogenetic signal may be the consequence of such selection that has overridden phylogenetic structure in SShD or, alternatively, it may be the result of evolution by random drift (e.g., in association with speciation and founder events). It is plausible that both sexual selection and natural selection affect the evolution of SShD. For instance, the convergent evolution of an unusual sexual dimorphism of body shape in two distantly related species of flies appears to be due to sexual selection from male-male interactions , whereas dimorphism of body proportions in lizards has been related to ecological factors [19, 20]. In Drosophila, the two most likely selective pressures influencing the evolution of SShD in the wing are flight and the production of courtship song. Because females fly to find oviposition sites whereas males search mating opportunities, flight requirements and optimal wing shapes may differ between sexes and result in different regimes of natural selection. Such a difference has been reported for a midge, where different flight behaviour in males and females is associated with sexual shape dimorphism in the wings . Sexual selection may affect SShD in Drosophila wings because males use the wings for producing courtship songs with species-specific characteristics . Experimental alteration of sexual selection can produce changes in courtship song , which may be related to changes in wing shape. Because the species of the melanogaster subgroup live in a range of different environments [59, 60], flight behaviour and mating systems differ between species and are likely to yield divergent regimes of natural and sexual selection that may affect SShD in the wings. Because little is known about natural or sexual selection on wing shape, any hypotheses about their effects on SShD must remain very tenuous until more information is available.
This study characterised SShD in all nine species of the Drosophila melanogaster subgroup. We showed a significant divergence of SShD in this clade. There was no detectable phylogenetic signal in the differences, and SShD has both allometric and non-allometric components. Overall, there was little evidence for evolutionary constraints and our results suggested that SShD can readily evolve in response to selection. Particularly because little is known about the natural and sexual selection on wing shape in Drosophila, the mechanisms involved in the evolution of SShD are far from being understood.
Samples and measurements
We used flies from long-term laboratory populations of each of the species known to belong to the Drosophila melanogaster subgroup: D. erecta (49 females; 51 males), D. mauritiana (47 females; 43 males), D. melanogaster (49 females; 53 males), D. orena (49 females; 50 males), D. santomea (50 females; 50 males), D. sechellia (50 females; 34 males), D. simulans (47 females; 37 males), D. teissieri (52 females; 53 males), D. yakuba (52 females; 52 males). The right wing of each fly was mounted on a microscope slide in a 70% lactic acid/30% ethanol solution and photographed in standard conditions with camera mounted on a Leica microscope.
Wing shape was assessed with 15 landmarks located at vein intersections with other veins and margins (Figure 1). The landmarks were digitised from the images with ImageJ  and a plug-in specifically written for this purpose.
Morphometric and statistical analyses
To quantify wing size, centroid size  was computed from the raw coordinates of the landmarks. A generalised Procrustes superimposition  was performed on the landmark coordinates to extract the information on wing shape (Procrustes coordinates). All morphometric and statistical analyses were performed using SAS 9.01 software , including routines for the Procrustes fit as well as the permutation and bootstrap procedures.
Because sex dimorphism is computed from the mean sizes or shapes of the sexes within each species (unlike, e.g., studies of fluctuating asymmetry ) and because our sample sizes are relatively large, measurement error is not a serious concern. Moreover, the methods like those used to acquire coordinate data produce measurement errors that are much smaller than the variation among individuals , so that the effect of measurement error on the estimates of sex dimorphism or their statistical variability is negligible.
The effects of species and sex on centroid size were tested with an ANOVA (type III sums of squares). The effects of species, sex and centroid size on shape were tested with a MANOVA (type III sums of squares and cross-products).
SShD for each species was computed as the mean vector for female shape minus the mean vector for male shape. Because male and female samples are independent within each species, covariance matrices for SShD were computed as the sum of the covariance matrices for the mean shapes of the males and females.
To compare SShD between pairs of species, we used a bootstrap test  against the null hypothesis that the species share the same SShD. To simulate this null hypothesis, the data for each pair of species were first modified so that they had the same SShD. Bootstrap samples were then repeatedly drawn from the modified data, SShD was computed in both species, and the magnitude of the difference in SShD was then compared to the value obtained from the original data. For each pairwise combination of species, 10,000 rounds of bootstrap resampling were used.
The principal component analysis of SShD was based on the covariance matrix of the SShD vectors for the nine species, and therefore maintains the metric of the Procrustes tangent distances . The 95% confidence ellipses  for PC scores were computed from the covariance matrices for SShD of the respective species.
Phylogenetic signal in SShD
The topology of the phylogenetic tree of the Drosophila melanogaster subgroup was assembled as a consensus from the highly supported phylogenetic information based on the entire genomes of five species of the melanogaster subgroup  together with single-sequence phylogenies of the entire subgroup [44, 45] and information from several other studies focussing on specific species complexes [65–69]. Points of disagreement concerned the rooting of the tree, for which the whole-genome analyses  provide decisive support, and the sequence of splits between D. simulans, D. sechellia and D. mauritiana, where several sources of evidence favour an earlier divergence of D. sechellia [69–71]. Branch lengths for the resulting consensus tree (Figure 4) were compiled from published estimates of divergence times [66, 72–75].
To reconstruct the SShD corresponding to the internal nodes of the phylogeny, we used squared-change parsimony, weighted by branch lengths [47, 48]. To visualize evolutionary trajectories, we computed the PC scores corresponding to the reconstructed SShD values for the internal nodes of the phylogeny and included them in the PC plots of Figure 3 [49–51].
To test for a phylogenetic signal in SShD against the null hypothesis of the absence of any phylogenetic signal, we used a permutation procedure that is a multivariate extension of a method developed by Laurin . The test uses a permutation approach  to simulate the null hypothesis that there is no phylogenetic signal. The nine observed SShD vectors were reassigned to the terminal nodes of the phylogenetic tree in all 362,880 possible permutations of the original order. For each permutation, the new values were mapped onto the phylogeny by weighted squared-change parsimony and the tree length (the sum of the squared Procrustes distances between ancestral and descendant SShD for all branches of the tree) was computed. If the phylogenetic signal in the data is strong, it is expected that random permutation of SShD will usually result in a greater tree length than that of the original data. Therefore, the test computes the P-value as the proportion of permutations that result in a tree length that is equal to or less than the one observed for the original data (more detail in Klingenberg and Gidaszewski, unpublished).
Allometric and non-allometric components of SShD
Total SShD was decomposed into allometric and non-allometric components of variation in each species (Figure 5). We characterized allometry with a multivariate regression of shape (Procrustes coordinates) on centroid size . To obtain a single allometric regression for each species, we used pooled within-sex regression of shape on centroid size.
To test the assumption that the allometric regression is the same in both sexes of each species, we conducted bootstrap tests  with the null hypothesis of a common allometric regression between sexes. To simulate this null hypothesis, separate regressions of shape on centroid size were conducted in each sex and the residuals from these regressions were computed. A pooled within-sex regression of shape on centroid size was then used to obtain a common regression vector. The predicted values from this common regression were added to the residuals of the two separate within-sex regressions to generate a modified dataset for which the null hypothesis of a common allometry in both sexes holds true, but which otherwise corresponds to the original samples. Bootstrap samples for males and females were randomly drawn, with replacement, from the modified dataset and the multivariate regressions of shape on centroid size were computed. To quantify the difference between the regression vectors, the sum of squared differences between corresponding elements was computed, and this value was compared to the difference obtained from the allometric regressions with the original data. The bootstrap resampling procedure was repeated 10,000 times for each species and the achieved significance level was computed as the proportion of the rounds of bootstrap resampling for which this difference was equal to or greater than the difference obtained for the original data.
The allometric component of SShD was computed as the vector of coefficients of allometric regression times the difference of centroid size between males and females (Figure 5). We computed the non-allometric component of SShD as the total SShD minus the allometric SShD. The magnitudes of total SShD and the allometric and non-allometric components are the lengths of the corresponding vectors, expressed in units of Procrustes distance.
Lande R: Sexual dimorphism, sexual selection and adaptation in polygenic characters. Evolution. 1980, 34 (2): 292-305. 10.2307/2407393.
Lande R, Arnold SJ: Evolution of mating preference and sexual dimorphism. Journal of Theoretical Biology. 1985, 117: 651-664. 10.1016/S0022-5193(85)80245-9.
Hedrick AV, Temeles EJ: The evolution of sexual dimorphism in animals – Hypotheses and tests. Trends Ecol Evol. 1989, 4 (5): 136-138. 10.1016/0169-5347(89)90212-7.
Fairbairn DJ, Blanckenhorn W, Szekely T: Sex, size and gender roles: Evolutionary studies of sexual size dimorphism. 2007, Oxford: Oxford University Press
Ploch I: Taxonomic interpretation and sexual dimorphism in the early cretaceous (Valanginian) ammonite Valanginites nucleus (Roemer, 1841). Acta Geologica Polonica. 2003, 53 (3): 201-208.
O'Higgins P, Johnson DR, Moore WJ, Flinn RM: The variability of sexual dimorphism in the hominoid skull. Experientia. 1990, 46: 670-672. 10.1007/BF01939929.
Taylor AB: Size and shape dimorphism in great ape mandibles and implications for fossil species recognition. American Journal of Physical Anthropology. 2006, 129: 82-98. 10.1002/ajpa.20266.
Hingst-Zaher E, Marcus LF, Cerqueira R: Application of geometric morphometrics to the study of postnatal size and shape changes in the skull of Calomys expulsus. Hystrix. 2000, 11 (1): 99-113.
Valenzuela N, Adams DC, Bowden RM, Gauger AC: Geometric morphometrics sex estimation for hatchling turtles: a powerful alternative for detecting subtle sexual shape dimorphism. Copeia. 2004, 4: 735-742. 10.1643/CH-03-248R1.
Oettlé AC, Pretorius E, Steyn M: Geometric morphometric analysis of mandibular ramus flexure. American Journal of Physical Anthropology. 2005, 128: 623-629. 10.1002/ajpa.20207.
Pretorius E: Using geometric morphometrics to investigate wing dimorphism in males and females of Hymenoptera – a case study based on the genus Tachysphex Kohl (Hymenoptera : Sphecidae : Larrinae). Australian Journal of Entomology. 2005, 44: 113-121. 10.1111/j.1440-6055.2005.00464.x.
Pretorius E, Steyn M, Scholtz Y: Investigation into the usability of geometric morphometric analysis in assessment of sexual dimorphism. American Journal of Physical Anthropology. 2006, 129: 64-70. 10.1002/ajpa.20251.
Bruner E, Costantini D, Fanfani A, Giacomo DO: Morphological variation and sexual dimorphism of the cephalic scales in Lacerta bilineata. Acta Zoologica (Stockholm). 2005, 86: 245-254. 10.1111/j.1463-6395.2005.00206.x.
Wood CG, Lynch JM: Sexual dimorphism in the craniofacial skeleton of modern humans. Advances in Morphometrics. Edited by: Marcus LF, Corti M, Loy A, Naylor GJP, Slice DE. 1996, New York: Plenum Press, 407-414.
O'Higgins P, Dryden IL: Sexual dimorphism in hominoids: further studies of craniofacial shape differences in Pan, Gorilla and Pongo. Journal of Human Evolution. 1993, 24: 183-205. 10.1006/jhev.1993.1014.
O'Higgins P, Collard M: Sexual dimorphism and facial growth in papionin monkeys. Journal of Zoology, London. 2002, 257: 255-272. 10.1017/S0952836902000857.
Leigh SR: Cranial ontogeny of Papio baboons (Papio hamadryas). American Journal of Physical Anthropology. 2006, 130: 71-84. 10.1002/ajpa.20319.
Schaefer K, Mitteroecker P, Gunz P, Bernhard M, Bookstein FL: Craniofacial sexual dimorphism patterns and allometry among extant hominids. Annals of Anatomy. 2004, 471-478. 10.1016/S0940-9602(04)80086-4.
Butler MA, Losos JB: Multivariate sexual dimorphism, sexual selection, and adaptation in Greater Antillean Anolis lizards. Ecological Monographs. 2002, 72 (4): 541-559.
Butler MA, Sawyer SA, Losos JB: Sexual dimorphism and adaptive radiation in Anolis lizards. Nature (Lond). 2007, 447: 202-205. 10.1038/nature05774.
Malmgren JC, Thollesson M: Sexual size and shape dimorphism in two species of newts, Triturus cristatus and T. vulgaris (Caudata: Salamandridae). Journal of Zoology, London. 1999, 249: 127-136. 10.1111/j.1469-7998.1999.tb00750.x.
Bonduriansky R: Convergent evolution of sexual shape dimorphism in Diptera. Journal of Morphology. 2006, 267: 602-611. 10.1002/jmor.10426.
Atchley WR: Components of sexual dimorphism in Chironomus larvae (Diptera: Chironomidae). American Naturalist. 1971, 105 (945): 455-466. 10.1086/282737.
Miller JS, Venable DL: Floral morphometrics and the evolution of sexual dimorphism in Lycium (Solanacae). Evolution. 2003, 57 (1): 74-86.
Mosimann JE: Size allometry: size and shape variables with characterizations of the lognormal and generalized gamma distributions. J Am Stat Assoc. 1970, 65: 930-945. 10.2307/2284599.
Gould SJ: Allometry and size in ontogeny and phylogeny. Biological Reviews. 1966, 41: 587-640. 10.1111/j.1469-185X.1966.tb01624.x.
Gould SJ: Allometry in Primates, with emphasis on scaling and the evolution of the brain. Approaches to Primate Paleobiology. Edited by: Szalay. 1975, Basel: Karger, 5: 244-292.
Klingenberg CP: Multivariate allometry. Advances in Morphometrics. Edited by: Marcus LF, Corti M, Loy A, Naylor GJP, Slice DE. 1996, New York: Plenum Press
Klingenberg CP: Heterochrony and allometry: the analysis of evolutionary change in ontogeny. Biological Reviews of the Cambridge Philosophical Society. 1998, 73 (1): 79-123. 10.1017/S000632319800512X.
Schwarzkopf L: Sexual dimorphism in body shape without sexual dimorphism in body size in water skinks (Eulamprus quoyii). Herpetologica. 2005, 61 (2): 116-123. 10.1655/04-66.
Ashburner M, Golic KG, Hawley RS: Drosophila : a laboratory handbook. 2005, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press
Powell JR: Progress and prospects in evolutionary biology: the Drosophila model. 1997, New-York, Oxford: Oxford University Press
Hansen TF, Houle D: Evolvability, stabilizing selection, and the problem of stasis. Phenotypic integration – studying the ecology and evolution of complex phenotypes. Edited by: Pigliucci M, Preston K. 2004, New York: Oxford University Press
Weber KE: Increased selection response in larger populations. I. Selection for wing-tip height in Drosophila melanogaster at three populations sizes. Genetics. 1990, 125: 579-584.
Dryden IL, Mardia KV: Statistical shape analysis. 1998, Chichester: Wiley
Gilchrist AS, Azevedo RBR, Partridge L, O'Higgins P: Adaptation and constraint in the evolution of Drosophila melanogaster. Evolution and Development. 2000, 2 (2): 114-124. 10.1046/j.1525-142x.2000.00041.x.
Palsson A, Gibson G: Association between nucleotide variation in Egfr and wing shape in Drosophila melanogaster. Genetics. 2004, 167: 1187-1198. 10.1534/genetics.103.021766.
Klingenberg CP, Zaklan SD: Morphological integration between developmental compartments in the Drosophila wing. Evolution. 2000, 54: 1273-1285.
Breuker CJ, Patterson JS, Klingenberg CP: A single basis for developmental buffering of Drosophila wing shape. PLoS ONE. 2006, 1: e7-10.1371/journal.pone.0000007.
Dworkin I, Gibson G: Epidermal growth factor receptor and transforming growth factor-β signaling contributes to variation for wing shape in Drosophila melanogaster. Genetics. 2006, 173: 1417-1431. 10.1534/genetics.105.053868.
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 (12): 2773-2784.
Efron B, Tibshirani RJ: An introduction to the bootstrap. 1993, New York: Chapman & Hall
Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
Cariou M-L, Silvain J-F, Daubin V, Da Lage J-L, Lachaise D: Divergence between Drosophila santomea and allopatric or sympatric populations of D. yakuba using paralogous amylase genes and migration scenarios along the Cameroon volcanic line. Molecular Ecology. 2001, 10: 649-660. 10.1046/j.1365-294x.2001.01225.x.
Lachaise D, Harry M, Solignac M, Lemeunier F, Bénassi V, Cariou M-L: Evolutionary novelties in islands: Drosophila santomea, a new melanogaster sister species from Sao Tome. Proceedings of the Royal Society of London B, Biological Sciences. 2000, 267: 1487-1495. 10.1098/rspb.2000.1169.
Drosophila 12 Genomes Consortium: Evolution of genes and genomes on the Drosophila phylogeny. Nature (Lond). 2007, 450: 203-218. 10.1038/nature06341.
Rohlf JF: Comparative methods for the analysis of continuous variables: geometric interpretations. Evolution. 2001, 55 (11): 2143-2160.
Maddison WP: Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree. Systematic Biology. 1991, 40: 304-314.
Klingenberg CP, Ekau W: A combined morphometric and phylogenetic analysis of an ecomorphological trend: pelagization in Antarctic fishes (Perciformes: Nototheniidae). Biol J Linn Soc. 1996, 59 (2): 143-177. 10.1111/j.1095-8312.1996.tb01459.x.
Rohlf FJ: Geometric morphometrics and phylogeny. Morphology, shape, and phylogeny. Edited by: MacLeod N, Forey PL. 2002, London: Taylor & Francis, 175-193.
Nicola PA, Monteiro LR, Pessôa LM, Von Zuben FJ, Rohlf FJ, Dos Reis SF: Congruence of hierachical, localized variation in cranial shape and molecular phylogenetic structure in spiny rats, genus Trinomys (Rodentia: Echimyidae). Biol J Linn Soc. 2003, 80: 385-396. 10.1046/j.1095-8312.2003.00245.x.
Laurin M: The evolution of body size, Cope's rule and the origin of amniotes. Systematic Biology. 2004, 53: 594-622. 10.1080/10635150490445706.
Monteiro LR: Multivariate regression models and geometric morphometrics: the search for causal factors in the analysis of shape. Syst Biol. 1999, 48: 192-199. 10.1080/106351599260526.
Cowley DE, Atchley WR, Rutledge JJ: Quantitative genetics of Drosophila melanogaster. I. Sexual dimorphism in genetic parameters for wing traits. Genetics. 1986, 114: 549-566.
Baylac M, Penin X: Wing static allometry in Drosophila simulans males (diptera, drosophilidae) and its relationships with developmental compartments. Acta Zoologica Academia Scientiarum Hungaricae. 1998, 44 (1–2): 97-112.
McLachlan AJ: Sexual dimorphism in Midges: strategies for flight in the rain-pool dweller Chironomus imicola (Diptera: Chironomidae). Journal of Animal Ecology. 1986, 55 (1): 261-267. 10.2307/4706.
Cowling DE, Burnet B: Courtship songs and genetic control of their acoustic characteristics in sibling species of the Drosophila melanogaster subgroup. Animal Behaviour. 1981, 29: 924-935. 10.1016/S0003-3472(81)80030-9.
Snook RR, Robertson A, Crudgington HS, Ritchie MG: Experimental manipulation of sexual selection and the evolution of coutrship song in Drosophila pseudoobscura. Behav Genet. 2005, 35: 245-255. 10.1007/s10519-005-3217-0.
Lachaise D, Cariou M-L, David JR, Lemeunier F, Tsacas L, Ashburner M: Historical biogeography of the Drosophila melanogaster species subgroup. Evolutionary Biology. 1988, 22: 159-225.
Lachaise D, Capy P, Joly D, Lemeunier F, David JR: Nine relatives from one African ancestor: Population biology and evolution of the Drosophila melanogaster subgroup species. The evolution of population biology. Edited by: Singh RS, Uyenoyama MK. 2003, Cambridge: Cambridge University Press
Rasband WS: ImageJ. 1997, Bethesda, Maryland, USA: U. S. National Institutes of Health
SAS Institute I: SAS System 9.1. Cary, NC, USA
Klingenberg CP, McIntyre GS: Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with Procrustes methods. Evolution. 1998, 52: 1363-1375. 10.2307/2411306.
Sokal RR, Rohlf JF: Biometry: the principles and practice of statistics in biological research. 2000, New-York: W. H. Freeman and company
Parsch J: Selective constraints on intron evolution in Drosophila. Genetics. 2003, 165: 1843-1851.
Kliman RM, Andolfatto P, Coyne JA, Depaulis F, Kreitman M, J BA, McCarter J, Wakeley J, Hey J: The population genetics of the origin and divergence of the Drosophila simulans complex species. Genetics. 2000, 156: 1913-1931.
Ting C-T, Tsaur S-C, Wu C-I: The phylogeny of closely related species as revealed by the genealogy of a speciation gene. Proceedings of the National Academy of Sciences. 2000, 97: 5313-5316. 10.1073/pnas.090541597.
Harr B, Weiss S, David JR, Brem G, Schlötterer C: A microsatellite-based multilocus phylogeny of the Drosophila melanogaster species complex. Current Biology. 1998, 8: 1183-1186. 10.1016/S0960-9822(07)00490-3.
Harr B, Schlötterer C: Patterns of microsatellite variability in the Drosophila melanogaster complex. Genetica (Dordr). 2004, 120: 71-77. 10.1023/B:GENE.0000017631.00820.49.
Harr B, Weiss S, David JR, Brem G, Schlötterer C: A microsatellite-based multilocus phylogeny of the Drosophila melanogaster species complex. Curr Biol. 1998, 8: 1183-1186. 10.1016/S0960-9822(07)00490-3.
Kliman RM, Andolfatto P, Coyne JA, Depaulis F, Kreitman M, Berry AJ, McCarter J, Wakeley J, Hey J: The population genetics of the origin and divergence of the Drosophila simulans complex species. Genetics. 2000, 156: 1913-1931.
Llopart A, Elwyn S, Lachaise D, Coyne JA: Genetics of a difference of pigmentation between Drosophila yakuba and Drosophila santomea. Evolution. 2002, 56: 2262-2277.
Russo CAM, Takezaki N, Nei M: Molecular phylogeny and divergence times of drosophilid species. Mol Biol Evol. 1995, 12: 391-404.
Tamura K, Subramanian S, Kumar S: Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocks. Mol Biol Evol. 2004, 21: 36-44. 10.1093/molbev/msg236.
Li Y-J, Satta Y, Takahata N: Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method. Genes and Genetic Systems. 1999, 74: 117-127. 10.1266/ggs.74.117.
Manly BFJ: Randomization, bootstrap and Monte Carlo methods in biology. 2007, Boca Raton, FL: Chapman & Hall/CRC
We thank Jean-François Ferveur, Daniel Lachaise and Laurence Monti for providing the fly samples. We are grateful to three anonymous reviewers who provided comments that helped us to improve the manuscript. NAG was supported by a NERC studentship and a bursary from the Faculty of Life Sciences at The University of Manchester.
The project was initiated by CPK and MB. MB arranged for samples of several species. NAG prepared the wings, collected the data and performed the analyses. CPK supervised the project and programmed ImageJ and SAS routines used in the analyses. NAG and CPK wrote the paper.