Variation of wing shape in the Drosophila virilis species group (Diptera: Drosophilidae)

Wing shape variation was analysed with geometric morphometric methods in 17 laboratory strains, representing 11 closely related species (including two subspecies) of the Drosophila virilis group: D. virilis, D. lummei, D. novamexicana, D. americana americana, D. americana texana, D. montana, D. lacicola, D. flavomontana, D. borealis, D. littoralis, D. ezoana and D. kanekoi. Overall shape estimated using Procrustes coordinates of 14 landmarks was highly variable among strains and very similar in females and males. The landmarks in the distal part of the wing showed higher variation across strains than those in the proximal part. Procrustes distances between species were not consistent with phylogenetic distances previously suggested for the virilis group. Moreover, Procrustes distances between strains within species and within two major phylads (virilis and montana) were comparable with those between species and between phylads, respectively. The most different from other members of the group was the endemic D. kanekoi species, currently viewed as separate subphylad within the montana phylad. Allometric effects were found to be partly responsible for shape differences between the strains. Three most significant shape transformations were considered using the relative warp analysis and the strains were ordinated in accordance with transformation values. The pattern of relative warp scores could be easily interpreted only for the third warp explaining about 13% of shape variation. It separated the largest species, D. montana, D. ezoana and D. kanekoi, from other ones and was mainly associated with shape changes in the proximal region of the wing. The results of the present work suggest that wing shape in the virilis species group is not related to the speciation process. The observed proximal‐distal contrasts and allometric effects are in agreement with data of other studies, in which wing shape variation was analysed within Drosophila species.     

Principles and methods of geometric morphometrics

The basic concepts, notions and methods of geometric morphometrics (GM) are considered. This approach implies multivariate analysis of landmark coordinates located following certain rules on the surface of a morphological object. The aim of GM is to reveal differences between morphological objects by their shapes as such, the "size factor" being excluded. The GM is based on the concept of Kendall's space (KS) defined as a hypersphere with points distributed on its surface. These points are the shapes defined as aligned landmark configurations. KS is a non-Euclidian space, its metrics called Procrustes is defined by landmark configuration of a reference shape relative to which other shapes are aligned and compared. The differences among shapes are measured as Procrustes distances between respective points. For the linear methods of multivariate statistics to be applied to comparison of shapes, the respective points are projected onto the tangent plane (tangent space), the tangent point being defined by the reference. There are two principal methods of shape comparisons in GM: the Procrustes superimposition (a version of the least squares analysis) and thin-plate spline analysis. In the first case, Procrustes residuals are the outcome shape variables which remain after isometric alignment of the shapes being compared. Their summation over all landmarks yields Procrustes distances among these shapes. The Procrustes distances can be used in multivariate analyses just as the Euclidian distances. In the second case, the shapes are fitted to the references by stretching/compressing and shearing until complete identity of their landmark configurations. Eigenvectors of resulting bending energy matrix are defined as new shape variables, principal warps which yield another shape space with the origin defined by the reference. Projections of the shapes being compared onto principal warps yield partial warps, and their covariance matrix decomposition into eigenvectors yields relative warps which are similar to principal components (in particular, they are mutually orthogonal). Both partial and relative warps can be used in many multivariate statistic analyses as quantitative shape variables. Results of thin-plate spline analysis can be represented graphically by transformation grid which displays type, amount and localization of the shape differences. Basis rules of sample composition and landmark positioning to be used in GM are considered. At present, rigid (with minimal degrees of freedom) 2D morphological objects are most suitable for GM applications. It is important to recognize three type of real landmarks, and additionally semi-landmarks and "virtual" landmarks. Some procedures of thin-plate spline analysis are considered exemplified by some study cases, as well as applications of some standard multivariate methods to GM results. They make it possible to evaluate correlation between different shapes, as well as between a shape and some non-shape variables (linear measurements etc); to evaluate the differences among organisms by shape of a morphological structure; to identify landmarks which most accounted for both correlation and differences between the shapes. An annotated list of most popular softwares for GM is provided.     

Mandibular shape variation in mainland and insular hylobatids

Although hylobatids are the most speciose of the living apes, their morphological interspecies and intraspecies variation remains poorly understood. Here, we assess mandibular shape variation in two species of Hylobates, white-handed (Hylobates lar) and black-handed (Hylobates agilis) gibbons. Using 71 three-dimensional landmarks to quantify mandibular shape, interspecies and intraspecies variation and geographic patterns of mandibular shape are examined in a mixed sex sample of adult H. lar and H. agilis through generalized Procrustes analysis, Procrustes analysis of variance, and principal components analysis. We find that relative to H. agilis, H. lar exhibits a higher amount of variation in mandibular shape. Both species demonstrate similar allometric patterns in mandibular shape. We also highlight a geographic pattern in mandibular shape variation. Compared to mainland hylobatids, insular hylobatids have relatively lower, more posteriorly oriented, and anteroposteriorly wider mandibular condyles, with an increased distance between the condyles and the coronoid processes. This geographic pattern could reflect differences in functional demands on the mandible during mastication and/or could be driven by factors often associated with evolutionary pressures of island populations relative to mainland populations. The findings of this study highlight how little is known about Hylobates morphological variation and how important this is for using Hylobates to help interpret the primate fossil record. Understanding interspecific and intraspecific variation in extant primates is vital to interpreting variation in the primate fossil record.     

Are wing size,wing shape and asymmetry related to field fitness of Trichogramma egg parasitoids?

The effects of wing shape, wing size, and fluctuating asymmetry in these measures on the field fitness of T. nr. brassicae and T. pretiosum were investigated. Trichogramma wasps mass-reared on eggs of the factitious host Sitotroga cerealella were released in tomato paddocks and those females ovipositing on Helicoverpa spp. eggs were recaptured. Comparisons of the recaptured group with a sample from the release population were used to assess fitness. Wing data were obtained by positioning landmarks on mounted forewings. Size was then measured as the centroid size computed from landmark distances, while Procrustes analysis followed by principal component analysis was used to assess wing shape. Similar findings were obtained for both Trichogramma species: fitness of wasps was strongly related to wing size and some shape dimensions, but not to the asymmetries of these measures. Wasps which performed well in the field had larger wings and a different wing shape compared to wasps from the mass reared population. Both size and the shape dimensions were linearly associated with fitness although there was also some evidence for non-linear selection on shape. The results suggest that wing shape and wing size are reliable predictors of field fitness for these Trichogramma wasps.     

Evolutionary directional asymmetry and shape variation in Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae): an example using hind wings

The western corn rootworm Diabrotica virgifera virgifera LeConte is a pest of maize in the USA and Europe and especially a problem in particular regions of Croatia. In the present study, patterns of variation in hind wing shape were examined. The first objective was to examine the influence of soil type on 10 populations of D. v. virgifera sampled from three regions in Croatia that differed according to edaphic factors and climate. The second objective was to investigate the potential evolutionary presence of directional asymmetry on hind wings. Geometric morphometrics was used to examine these objectives by quantifying the morphological variation within and among individuals and populations. Overall, D. v. virgifera hind wing shape changed according to major soil type classifications in Croatia. The three hind wing morphotypes found varied because of basal radial vein differences, related to landmarks 1, 3, 7, and 14. The findings of the present study show that hind wing shape in D. v. virgifera can be used to differentiate populations based on edaphic factors and may have application as a monitoring tool in the integrated management of D. v. virgifera. In an evolutionary context, the presence of directional asymmetry in the hind wings of D. v. virgifera adds to the ever growing data on the evolution of insect wings. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 110–118.     

Fluctuating asymmetry in the honey bee,Apis mellifera: effects of ploidy and hybridization

We present a new measure of morphological asymmetry that avoids most of the statistical problems inherent in character-by-character analysis of size or shape. The method is an application of Procrustes analysis, which computes best-fitting super-positions of configurations of landmarks to the left and right sides of a single specimen. The Procrustes method combines subtle deviations in all aspects of the landmark configuration into one net asymmetry score. Directional asymmetry is separated from fluctuating asymmetry in a simple partition of a net sum-of-squares, and geometrical details of either component can be inspected by traditional methods of multivariate statistical analysis of landmarks. We demonstrate this method in a comparison of wing venation asymmetry in male (haploid) and female (diploid) honey bees (Apis mellifera). In addition we investigate the effects of ploidy and inter-subspecies hybridization on asymmetry and wing venation abnormalities, using the subspecies A. m. mellifera, A. m. carnica, and the hybrid strain “Nigra”. Results suggest that while the haploid males showed a higher frequency of wing venation abnormalities and greater total asymmetry than the diploid females, most of the asymmetry difference between males and females was in the form of directional, not fluctuating, asymmetry. Hybrid females had a higher frequency of wing venation abnormalities than females of either subspecies, but there were no significant differences in the mean level of asymmetry among females of A. m. mellifera, A. m. carnica and hybrid Nigra. Hybrid males had higher absolute frequency of wing venation abnormalities and asymmetry than males of either subspecies. However the mean frequency of venation abnormalities did not differ significantly between Nigra and A. m. carnica males, and mean asymmetries were not significantly different between Nigra and A. m. mellifera males. We discuss the relationship which is assumed to exist between developmental stability and fluctuating asymmetry in light of our result.     

Primate cranial diversity

Many studies in primate and human evolution focus on aspects of cranial morphology to address issues of systematics, phylogeny, and functional anatomy. However, broad analyses of cranial diversity within Primates as an Order are notably absent. In this study, we present a 3D geometric morphometric analysis of primate cranial morphology, providing a multivariate comparison of the major patterns of cranial shape change during primate evolution and quantitative assessments of cranial diversity among different clades. We digitized a set of 18 landmarks designed to capture overall cranial shape on male and female crania representing 66 genera of living primates. The landmark data were aligned using a Generalized Procrustes Analysis and then subjected to a principal components analysis to identify the major axes of cranial variation. Cranial diversity among clades was compared using multivariate measurements of variance. The first principal component axis reflects differences in cranial flexion, orbit size and orientation, and relative neurocranial volume. In general, it separates strepsirrhines from anthropoids. The second axis reflects differences in relative cranial height and snout length and primarily describes differences among anthropoids. Eulemur, Mandrillus, Pongo, and Homo are among the extremes in cranial shape. Anthropoids, catarrhines, and haplorhines show a higher variance than prosimians or strepsirrhines. Hominoids show the highest variance in cranial shape among extant primate clades, and much of this diversity is driven by the unique cranium of Homo sapiens. Am J Phys Anthropol 142:565–578, 2010. © 2010 Wiley‐Liss, Inc.     

A geometric morphometric approach to the analysis of skull shape in Triassic dicynodonts (Therapsida,Anomodontia) from South America

Dicynodont therapsids were a major component of the Permo-Triassic terrestrial ecosystems across Pangea and have been regarded as specialized herbivores. In South America, the group was represented by several taxa of the clade Kannemeyeriiformes spanning from the Middle to the Late Triassic. In order to evaluate if cranial differences among taxa are potentially related to differences in feeding function, we performed a geometric morphometric analysis on 28 South American dicynodont crania. We digitized 19 cranial landmarks and conducted generalized Procrustes analysis, principal component analysis (PCA), principal component analysis between groups (bg-PCA), and a branch weighted squared-change parsimony approach. Phylogenetic inertia was not a significant driver of cranial shape evolution in the group, whereas PCA and bg-PCA support that major morphological shape differences are concentrated in the preorbital region (relative length of the snout and width of the caniniform process), in the position of quadrate condyle in relation to the caniniform process, and in the increase in the intertemporal surface area. In this context, tusked Dinodontosaurus, “Kannemeyeria,” and Vinceria have relatively smaller adductor attachment areas and input moment arm than younger taxa lacking tusks, such as Ischigualastia, Stahleckeria, and Jachaleria. Differences in cranial morphology in later dicynodonts reflect modifications in feeding mechanics, probably due to changes in food resources (vegetation) in their habitats toward the end of the Triassic.     

Variation in elasmoid fish scale patterns is informative with regard to taxon and swimming mode

Variations in scales from nine regions on the flank of teleost fish were examined from the point of view of functional adaptation and with regard to which scales best differentiate species. Three teleost species were selected; two are from the genus Mugil, M. cephalus and M. curema, which are phylogenetically distant from the third, Dicentrarchus labrax. Scale form was described using seven landmarks, the coordinates of which were subjected to generalized Procrustes analysis followed by principal components analysis. Principal component scores were submitted to cross‐validated discriminant analysis to assess the utility of each scale in identifying species. The best discrimination (98%) was obtained with the scale from the central‐dorsal area. Scales from the anterior and central zones are relatively wide dorsoventrally and narrow anteroposteriorly. This appears to be related to the profile of the lateral body wall and with subcarangiform swimming. Scales from the posterior region are anteroposteriorly long and dorsoventrally narrow, this shape possibly being related to thrust. Despite the wide phylogenetic separation between mullets and D. labrax, the pattern of scale variation is similar. This may imply strong functional convergence, although studies of sister taxa with different swimming modes are required to confirm this. © 2009 The Linnean Society of London, Zoological Journal of the Linnean Society, 2009, 155 , 834–844.     

Evolution of shape in oligocene and mioceneNotocarinovalva (Ostracoda, Crustacea): A multivariate statistical study

Largely due to L. F. Bookstein, geometric morphometrics has been developed as an amalgamation of techniques drawn from mathematical statistics, non-Euclidean geometry and computer graphics, applied to labelled points (landmarks) and the biological images upon which they are registered. A tool of fundamental importance is the method of interpolation known as the thin-plate spline. Bookstein's sample-oriented procedure of relative warps, applied to data composed of coordinate pairs observed on eight landmarks on three samples of two species of fossil marine ostracods (bivalved microcrustacean),Notocarinovalva airella andN. yulecartensis, separated in time, is used to obtain weight matrices, which, when appropriately partitioned, constitute the familiar data matrices of multivariate statistical analysis. Standard multivariate statistical analysis of samples may be carried out in the tangent space to shape space at the Procrustes average shape. Linear discriminant function scores were used for assessing, approximately, the evolutionary relevance of shape change and change in size (based on standard distance measures) in the species from the upper and lower Oligocene and lower middle Miocene of Victoria, Australia. The findings indicate that a small, though consistent, change in non-uniform shape of the carapace has taken place. The analysis indicates that the evolution in shape of the carapace could have been caused by random genetic drift in moderately large populations, tempered by weak selection.     

A cranial morphometric study of deer (Mammalia, Cervidae) from Argentina using three-dimensional landmarks

This study examines morphological variation in the crania (n = 70) of eight cervid species from Argentina. Forty 3-dimensional landmarks were acquired on each adult cervid cranium. The data were analysed using Morphologika software. The co-ordinates were registered and scaled to remove size differences by Procrustes analysis, and then principal components analysis was applied to examine shape variation. Shape variation associated with each principal component can be visualised in the program. The first principal component correlates strongly with the centroid size of the crania and also with the body mass and height of each species. The larger species were distinguished by relatively longer snouts and relatively smaller brains. The smallerMazama andPudu species cluster closely on the first as well as the other principal components. Among the larger species, the twoHippocamelus species, which live at higher altitudes, were clearly distinguished from the lowland species,Ozotoceros bezoarticus andBlastocerus dichotomus, on the basis of cranial flexion and the orientation of the occipital region. FinallyO. bezoarticus andB. dichotomus were compared directly and small differences were noted in the orbital region. The shape data was used to produce a distance matrix and a phenogram, which we relate to some of currently accepted phylogenetic relationships of this group of cervids.     

Advances in Geometric Morphometrics

Geometric morphometrics is the statistical analysis of form based on Cartesian landmark coordinates. After separating shape from overall size, position, and orientation of the landmark configurations, the resulting Procrustes shape coordinates can be used for statistical analysis. Kendall shape space, the mathematical space induced by the shape coordinates, is a metric space that can be approximated locally by a Euclidean tangent space. Thus, notions of distance (similarity) between shapes or of the length and direction of developmental and evolutionary trajectories can be meaningfully assessed in this space. Results of statistical techniques that preserve these convenient properties—such as principal component analysis, multivariate regression, or partial least squares analysis—can be visualized as actual shapes or shape deformations. The Procrustes distance between a shape and its relabeled reflection is a measure of bilateral asymmetry. Shape space can be extended to form space by augmenting the shape coordinates with the natural logarithm of Centroid Size, a measure of size in geometric morphometrics that is uncorrelated with shape for small isotropic landmark variation. The thin-plate spline interpolation function is the standard tool to compute deformation grids and 3D visualizations. It is also central to the estimation of missing landmarks and to the semilandmark algorithm, which permits to include outlines and surfaces in geometric morphometric analysis. The powerful visualization tools of geometric morphometrics and the typically large amount of shape variables give rise to a specific exploratory style of analysis, allowing the identification and quantification of previously unknown shape features.     

Effects of the evolution of the Serra do Mar mountains on the shape of the geographically isolated populations of Aegla schmitti Hobbs III, 1979 (Decapoda: Anomura)

The Aeglidae family is a useful model for evaluating the shape variation to present small geographically isolated populations and present a carapace structure that allows precise placement of landmarks. We analysed variations of the size and shape of the carapace of seven populations of A. schmitti. We used 18 bidimensional anatomical landmarks on the carapace of adult males. Size variation between the populations was analysed through a univariate analysis of variance, while the shape variation was analysed through a multivariate analysis of variance using the configurations aligned by the generalized Procrustes analysis. The ordination of the populations was investigated through a canonical variables analysis. Individuals from the populations differed in the carapace size, in this case, we observed the largest individuals in the Iguaçu river basin, followed by Ribeira do Iguape and the Coastal basin. Carapace shape also differed between the populations of A. schmitti, with the exception of the Mato Grande and Piraquara rivers: populations from the Iguaçu River presented a larger and rounder carapace when compared with the slenderer carapace of the Coastal basin, and the Ribeira do Iguape animals presented intermediate shapes. Data from the present study show great congruency with the draining pattern and geological history of the region.     

Geographic variation for wing shape in Drosophila serrata

Geographic variation in wing shape in female Drosophila serrata was examined by characterizing isofemale strains from 19 localities collected along a transect on the eastern coast of Australia. Shape variation was analyzed by Procrustes superimposition of landmark data followed by canonical variate analysis. The first extracted canonical variate showed a nonlinear association with latitude and accounted for 43% of the variance. There was a sharp increase in this variate at low latitudes as well as a gradual increase at high latitudes. These shape changes were associated with two landmarks at the edge of the wing. There was also a linear change in wing aspect. The isofemale heritability for two measures of shape was around 30%. Allometric relationships were weak both between localities and among isofemale strains within localities. The possibility that wing shape parameters are under selection independent of wing size is discussed.     

Effect of different larval rearing temperatures on the productivity (Ro) and morphology of the malaria vector Anopheles superpictus Grassi (Diptera: Culicidae) using geometric morphometrics

Temperature affects both the biology and morphology of mosquito vectors. Geometric morphometrics is a useful new tool for capturing and analyzing differences in shape and size in many morphological parameters, including wings. We have used this technique for capturing the differences in the wings of the malaria vector Anopheles superpictus, using cohorts reared at six different constant temperatures (15°, 20°, 25°, 27°, 30°, and 35° C) and also searched for potential correlations with the life tables of the species. We studied wing shape in both male and female adults, using 22 landmarks on the wing in relation to ecological parameters, including the development rate. The ecological zero was calculated as 9.93° C and the thermal constant as 296.34 day‐degrees. The rearing temperature affects egg, larval, and pupal development and also the total time from egg to adult. As rearing temperatures increased, longevity decreased in both sexes. In An. superpictus, R_o value and productivity correlated with the statistically significant gradual deformations in the wing shape related to size in both sexes. These deformations directly linked to differences in immature rearing temperatures. Analysis using PCA and UPGMA phenograms showed that although wings of females became narrower dorsoventrally as the temperature increased, they became broader in males. Comparisons of the wing landmarks indicated the medial part of the wing was most affected by larval rearing temperatures, showing relatively more deformations. Algorithmic values of the life tables were determined in correlation with the results of geometric morphometrics. Comparisons of centroid sizes in the cohorts showed that overall wing size became smaller in both sexes in response to higher rearing temperatures.     

An Integrated Approach for Landmark-Based Resistant Shape Analysis in 3D

The study of shape changes in morphology has seen a significant renovation in the last 20 years, particularly as a consequence of the development of geometric morphometric methods based on Cartesian coordinates of points. In order to extract information about shape differences when Cartesian coordinates are used, it is necessary to establish a common reference frame or system for all specimens to be compared. Therefore, a central issue in coordinate-based methods is which criterion should be used to align these configurations of points, since shape differences highly depend on those alignments. This is usually accomplished by aligning the configurations in a way that the sum of squared distances between coordinates of homologous points (landmarks) is minimized: the least-squares superimposition method. However, it is widely recognized that this method has some limitations when shape differences are not homogeneous across landmarks. Here we present an integrated approach for the resistant shape comparison of 3D landmark sets. It includes a new ordinary resistant Procrustes superimposition and its corresponding generalized resistant Procrustes version. In addition, they are combined with existing resistant multivariate statistical techniques for depicting the results. We demonstrate, by using both simulated and real datasets, that resistant Procrustes better detects and measures localized shape variation whenever present in up to half but one of the landmarks. The resistant Procrustes results are highly concordant with a priori biological information, and might dramatically improve the quality of inferences on patterns of shape variation.     

Geometric estimates of heritability in biological shape

The recently developed geometric morphometrics methods represent an important contribution of statistics and geometry to the study of biological shapes. We propose simple protocols using shape distances that incorporate geometric techniques into linear quantitative genetic models that should provide insights into the contribution of genetics to shape variation in organisms. The geometric approaches use Procrustes distances in a curved shape space and distances in tangent spaces within and among families to estimate shape heritability. We illustrate the protocols with an example of wing shape variation in the honeybee, Apis mellifera. The heritability of overall shape variation was small, but some localized components depicting shape changes on distal wing regions showed medium to large heritabilities. The genetic variance-covariance matrix of the geometric shape variables was significantly correlated with the phenotypic shape variance-covariance matrix. A comparison of the results of geometric methods with the traditional multivariate analysis of interlandmark distances indicated that even with a larger dimensionality, the interlandmark distances were not as rich in shape information as the landmark coordinates. Quantitative genetics studies of shape should greatly benefit from the application of geometric methods.