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.     

Biometrics,biomathematics and the morphometric synthesis

At the core of contemporarymorphometrics—the quantitative study of biological shape variation—is a synthesis of two originally divergent methodological styles. One contributory tradition is the multivariate analysis of covariance matrices originally developed as biometrics and now dominant across a broad expanse of applied statistics. This approach, couched solely in the linear geometry of covariance structures, ignores biomathematical aspects of the original measurements. The other tributary emphasizes the direct visualization of changes in biological form. However, making objective the biological meaning of the features seen in those diagrams was always problematical; also, the representation of variation, as distinct from pairwise difference, proved infeasible. To combine these two variants of biomathematical modeling into a valid praxis for quantitative studies of biological shape was a goal earnestly sought though most of this century. That goal was finally achieved in the 1980s when techniques from mathematical statistics, multivariate biometrics, non-Euclidean geometry and computer graphics were combined in a coherent new system of tools for the complete regionalized quantitative analysis oflandmark points together with the biomedical images in which they are seen. In this morphometric synthesis, correspondence of landmarks (biologically labeled geometric points, like “bridge of the nose”) across specimens is taken as a biomathematical primitive. The shapes of configurations of landmarks are defined as equivalence classes with respect to the Euclidean similarity group and then represented as single points in David Kendall'sshape space, a Riemannian manifold with Procrustes distance as metric. All conventional multivariate strategies carry over to the study of shape variation and covariation when shapes are interpreted in the tangent space to the shape manifold at an average shape. For biomathematical interpretation of such analyses, one needs a basis for the tangent space compatible with the reality of local biotheoretical processes and explanations at many different geometric scales, and one needs graphics for visualizing average shape differences and other statistical contrasts there. Both of these needs are managed by thethin-plate spline, a deformation function that has an unusually helpful linear algebra. The spline also links the biometrics of landmarks to deformation analysis of the images from which the landmarks originally arose. This article reviews the history and principal tools of this synthesis in their biomathematical and biometrical context and demonstrates their usefulness in a study of focal neuroanatomical anomalies in schizophrenia.     

Comparative analyses of <Emphasis Type="Italic">Astacus leptodactylus</Emphasis> morphological characteristics from Croatia and Armenia

Narrow-clawed crayfish, Astacus leptodactylus is a native European freshwater crayfish species, also distributed in Croatian freshwater systems belonging to the Black Sea drainage. Its taxonomical status is still in the process of change and discussion, and the data on morphological, molecular, ecological and zoogeographical characteristics of this species are scarce. Therefore, comparative analyses of morphological characteristics were applied with the aim to contribute to the knowledge on the morphometrical and meristical characteristics of A. leptodactylus. Recent research proved that measurements of a large number of morphometrical characteristics, in combination with multivariate statistical analysis, could provide a good instrument for identification and differentiation between populations. In this research altogether 143 animals were analysed (121 from two Croatian populations and 22 from Armenia). 22 morphometrical characteristics and 4 meristical characteristics, per crayfish, were measured. It was found that males and females differ between populations in measured meristical and morphometrical characteristics. None of recorded meristical characteristics proved itself to be reliable characteristic for distinguishing populations. From measured morphometrical characteristics the most discriminant characteristics for separating males from different populations were those describing carapace shape and for females those that describe shape of the claws.     

Heterochrony and geometric morphometrics: a comparison of cranial growth in Pan paniscus versus Pan troglodytes

Heterochrony, the classic framework in which to study ontogeny and phylogeny, in essence relies on a univariate concept of shape. Though principal component (PC) plots of multivariate shape data seem to resemble classical bivariate allometric plots, the language of heterochrony cannot be translated directly into general multivariate methodology. We simulate idealized multivariate ontogenetic trajectories and explore their appearance in PC plots of shape space and size-shape space. Only if the trajectories of two related species lie along exactly the same path in shape space can the classic terminology of heterochrony apply and pure dissociation of size change against shape change be detected. Regional heterochrony--the variation of apparent heterochrony by region--implies a dissociation of local growth fields and cannot be identified in an overall PC analysis. We exemplify a geometric morphometric approach to these issues using adult and subadult crania of 48 Pan paniscus and 47 Pan troglodytes specimens. On each specimen, we digitized 47 landmarks and 144 semilandmarks on facial curves and the external neurocranial surface. We reject the hypothesis of global heterochrony in the cranium of Pan as well as regional heterochrony for the lower face, the upper face, and the neurocranium.     

Getting into Shape: An Empirical Comparison of Traditional Truss-Based Morphometric Methods with a Newer Geometric Method Applied to New World Cichlids

Body shape is a difficult, but important, trait to quantify. Researchers have traditionally used multivariate analysis of several linear measures ('trusses') across the body form to quantify shape. Newer geometric morphometric methods claim to better estimate shape because they analyze the geometry among the locations of all landmarks simultaneously rather than the linear distances between pairs of landmarks. We tested this claim by comparing the results of several traditional morphometric analyses against a newer geometric analysis involving thin-plate splines (TPS), all applied to a common data set of morphologically variable new world cichlids Amphilophus citrinellus and A. zaliosus. The TPS method yielded slightly stronger evidence of morphological differences among forms, although traditional methods also distinguished the two species. Perhaps our most important result was the idiosyncratic interpretation of shape variation among the traditional truss-based methods, whereas the generation of deformation grids using the TPS approach yielded clear and visually interpretable figures. Our results indicate that geometric morphometrics can be a more effective way to analyze and interpret body form, but also that traditional methods can be relied upon to provide statistical evidence of shape differences, although not necessarily accurate information about the nature of variation in shape.     

Ontogenetic changes in the carapace shape of the non-marine ostracod Eucypris virens (Jurine)

Developmental changes in carapace form (size+shape) during ontogeny have been explored in Eucypris virens (Crustacea, Ostracoda) using elliptic Fourier analysis. Clones from different geographic localities raised under controlled constant conditions (temperature and photoperiod) were used to characterize developmental pathways in the species. A larger data set including field populations and laboratory populations cultured under a range of environmental conditions were used to infer influence of environmental factors on carapace shape changes during ontogeny. Size changes between consecutive juvenile stages support empirical laws describing the doubling of ostracod volume at each moult. Ontogenetic changes point out the remarkable influence of environmental conditions on carapace shape.     

Testing the success of palaeontological methods in the delimitation of clam shrimp (Crustacea,Branchiopoda) on extant species

Fossil spinicaudatan taxonomy heavily relies on carapace features (size, shape, ornamentation) and palaeontologists have greatly refined methods to study and describe carapace variability. Whether carapace features alone are sufficient for distinguishing between species of a single genus has remained untested. In our study, we tested common palaeontological methods on 481 individuals of the extant Australian genus Ozestheria that have been previously assigned to ten species based on genetic analysis. All species are morphologically distinct based on geometric morphometrics (p ≤ 0.001), but they occupy overlapping regions in Ozestheria morphospace. Linear discriminant analysis of Fourier shape coefficients reaches a mean model performance of 93.8% correctly classified individuals over all possible 45 pairwise species comparisons. This can be further increased by combining the size and shape datasets. Nine of the ten examined species are clearly sexually dimorphic but male and female morphologies strongly overlap within species with little influence on model performance. Ornamentation is commonly species-diagnostic; seven ornamentation types are distinguished of which six are species-specific while one is shared by four species. A transformation of main ornamental features (e.g. from punctate to smooth) can occur among closely related species suggesting short evolutionary timescales. Our overall results support the taxonomic value of carapace features, which should also receive greater attention in the taxonomy of extant species. The extensive variation in carapace shape and ornamentation is noteworthy and several species would probably have been assigned to different genera or families if these had been fossils, bearing implications for the systematics of fossil Spinicaudata.     

DETECTING SHAPE VARIATION IN OAK LEAF MORPHOLOGY: A COMPARISON OF ROTATIONAL-FIT METHODS

Rotational-fit methods were used to examine shape variation in oak leaf morphology by reference to a set of 14 landmarks for each leaf. Within-tree, between-tree, and between-species variations were examined. Generalized least-squares and generalized resistant-fit analyses revealed patterns of landmark variation that could be related to leaf architecture. The two species (Quercus palustris and Q. velutina) did not illustrate the same patterns of within-tree variability nor did they reveal similar degrees of between-tree variability. Incorporating an affine transformation algorithm resulted in little insight in some comparisons but suggested strong uniform shape change in other comparisons (especially between species). Resistant-fit mean square statistics were equally variable 1) in upper vs. lower crown samples within trees, 2) between trees within each species, and 3) between species. In addition, resistant-fit mean square statistics were found to be a poor measure of similarity, whether derived by comparison to a single reference object or by way of pairwise comparisons. Fundamental species differences in leaf shape are suggested by relationships among particular sets of landmarks, although overall shape differences cannot be explained fully by these analyses.     

Character coupling for taxa discrimination: a critical appraisal of quadratic assignment procedures (QAP)1

Taxa can be characterized by character coupling represented in similarity matrices. The customary methods of testing equality of variance-covariance matrices are based upon the multinormality assumption which is, however, frequently unacceptable in reality. Quadratic assignment procedures (QAP) have proved to be an alternative. They represent a type of computer-based test and utilize a random-permutation strategy to discover significant pattern correspondences between matrices. A comparison of the applicability of both testing methods requires an example with underlying multinormality. The samples of two species of land snails (Pulmonata, Helicidae), i. e. Arianta arbustorum (n = 104) and Arianta chamaeleon (n = 36), fulfil this requirement. Four parameters of shape and two parameters of spiral change were determined in each shell. The data serve as the basis for similarity matrices (variance-covariance, product-moment and rank order correlations). The inspection of methods reveals that QAP are suitable for correlation matrices, but can be applied for variance-covariance matrices with limitations only. Nevertheless, they are recommended procedures in taxonomy and evolutionary biology. Straightforward application, independence from distributional assumptions, and the possibility to test hypotheses of character coupling are advantageous features. The snail species are significantly discriminated by character coupling. Also, their parameters of shape and spiral change are morphologically integrated in a different way.     

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.     

Clinal variation in carapace shape in the South American freshwater crab,Aegla uruguayana (Anomura: Aeglidae)

South America has been influenced by different geoclimatic events ever since its separation from Africa. The inland water fauna has evolved in response to the changing landscape. Currently, there are indications of variations in populations, occurring to different degrees that would indicate a clinal pattern in morphology. Among South America's fauna, the freshwater anomuran, Aegla, is an enigmatic group as a result of its endemicity and is composed of only one genus. Of all the species in this family, Aegla uruguayana has the broadest distribution. Its native habitats have been influenced by several marine transgressions during the Miocene–Quaternary Periods; thus, it is likely that their current distribution has been more recent. Its habitat spreads across a number of isolated basins and sub‐basins that display distinct degrees of isolation/connection, making clinal variation patterns in the morphology of this species possible. The present study aimed to evaluate the pattern of carapace shape variation in A. uruguayana and how it relates to the isolation and/or connection of populations from different basins and sub‐basins, allowing the determination of any extant clinal patterns. The specimens studied belong to 25 separate populations, representing all areas in which the species currently exists. A total of 523 crabs were analyzed. We identified 13 landmarks and four semi‐landmarks in the carapace. The aeglids were divided into seven size intervals to avoid an allometry effect. In each size category, shape relationships analyzed by principal component analysis suggest a geographical pattern corresponding to the distribution of the populations studied. An evaluation of covariation between body shape and geographical coordinates reveals a strong pattern and shows that population distribution had a significant effect on species morphology. Additionally, according to covariance analysis, the variation in shape was not associated with the environmental variables studied. We observed a clinal pattern throughout the species distribution, which could be attributed to genetic drift. It is possible that this process is being amplified by the geographical isolation of the basins, differences in environmental characteristics, and low dispersal ability. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 914–930.     

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.     

The pH-related morphological variations of two acidophilic species of Desmidiales (Viridiplantae) isolated from a lowland peat bog,Czech Republic

Morphological variation related to pH was investigated in two acidophilic desmid species (Euastrum binale var. gutwinskii and Staurastrum hirsutum) utilizing geometric morphometric methods. Clones isolated from acidic habitats were cultured using a range of pH values from 3.5 to 6.5. The plasticity of ensuing populations was quantified and illustrated by the general Procrustes superimposition of landmarks placed along the outline of cells and subsequent statistical analyses of shape data. In both species, there was a significant effect of pH on the morphology of cells. In Staurastrum hirsutum, the pH-related morphological change was accompanied by a decrease in the size of cells cultured at a higher pH. However, in Euastrum binale, cell size did not differ in relation to pH, but cell shape was characterized by a deepening of the incisions between cell lobes at higher pH. In both species, cell complexity based on surface-to-volume ratio was positively correlated with increasing pH. We conclude that by manipulating their surface-to-volume ratios, these desmid species can respond to pH variations in their environment.     

First global approach: morphological and biological variability in a genetically homogeneous population of the European pilchard, Sardina pilchardus (Walbaum, 1792) in the North Atlantic coast

The European pilchard Sardina pilchardus represents the most commercially relevant fisheries resource in many countries bordering north Atlantic coasts and the Mediterranean Sea, being especially significant along the coast of Morocco. The continuous exploitation of this pelagic species for several decades places Morocco as the leader in sardine production. However, the conditions of exploitation of this resource underwent a great change during the last recent years. In order to identify the populations of the European pilchard sardine (Sardina pilchardus, Walbaum, 1792) in the Atlantic coast of Morocco and Spain, we have combined the truss network data to conduct multivariate analysis with biologic parameters and genetic analysis based on Microsatellite and mitochondrial control region data. Sardine morphometrics data truss variables from 10 samples spanning the Atlantic coast of Morocco were analysed by multivariate analysis. Thirteen morphometric measurements and some biological parameters such as the sex and the age of fishes were made for each individual. Discriminant analysis on size-corrected truss variables and cluster analysis of mean fishes shape using landmark data indicate, that the shape of north Moroccan sardines is distinct from the shape of sardines from south Morocco. However the analysis of the mitochondrial region and four microsatellites loci (Sp2, Sp7, Sp8 and SpI5) demonstrated an homogeneity population in the Moroccan Atlantic coast, with a low but significant genetic differentiation, which followed an isolation-by-distance pattern according to Mantel test.     

Linear Discrimination,Ordination, and the Visualization of Selection Gradients in Modern Morphometrics

Linear discriminant analysis (LDA) is a multivariate classification technique frequently applied to morphometric data in various biomedical disciplines. Canonical variate analysis (CVA), the generalization of LDA for multiple groups, is often used in the exploratory style of an ordination technique (a low-dimensional representation of the data). In the rare case when all groups have the same covariance matrix, maximum likelihood classification can be based on these linear functions. Both LDA and CVA require full-rank covariance matrices, which is usually not the case in modern morphometrics. When the number of variables is close to the number of individuals, groups appear separated in a CVA plot even if they are samples from the same population. Hence, reliable classification and assessment of group separation require many more organisms than variables. A simple alternative to CVA is the projection of the data onto the principal components of the group averages (between-group PCA). In contrast to CVA, these axes are orthogonal and can be computed even when the data are not of full rank, such as for Procrustes shape coordinates arising in samples of any size, and when covariance matrices are heterogeneous. In evolutionary quantitative genetics, the selection gradient is identical to the coefficient vector of a linear discriminant function between the populations before vs. after selection. When the measured variables are Procrustes shape coordinates, discriminant functions and selection gradients are vectors in shape space and can be visualized as shape deformations. Except for applications in quantitative genetics and in classification, however, discriminant functions typically offer no interpretation as biological factors.     

Testing the cranial evolutionary allometric ‘rule’ in Galliformes

Recent comparative studies have indicated the existence of a common cranial evolutionary allometric (CREA) pattern in mammals and birds, in which smaller species have relatively smaller faces and bigger braincases than larger species. In these studies, cranial allometry was tested using a multivariate regression between shape (described using landmarks coordinates) and size (i.e. centroid size), after accounting for phylogenetic relatedness. Alternatively, cranial allometry can be determined by comparing the sizes of two anatomical parts using a bivariate regression analysis. In this analysis, a slope higher or lower than one indicates the existence of positive or negative allometry, respectively. Thus, in those species that support the CREA ‘rule’, positive allometry is expected for the association between face size and braincase size, which would indicate that larger species have disproportionally larger faces. In this study, I applied these two approaches to explore cranial allometry in 83 Galliformes (Aves, Galloanserae), ranging in mean body weight from 30 g to 2.5 kg. The multivariate regression between shape and centroid size revealed the existence of a significant allometric pattern resembling CREA, whereas the second analysis revealed a negative allometry for beak size and braincase size (i.e. contrary to the CREA ‘rule’, larger galliform species have disproportionally shorter beaks than smaller galliform species). This study suggests that the presence of CREA may be overestimated when using cranium size as the standard measurement.     

Morphological stasis and phylogenetic relationships in Tadpole shrimps, Triops (Crustacea: Notostraca)

The tadpole shrimp (Triops) is a well-known 'living fossil' whose fundamental morphology has been unchanged for over 170 million years. Thus, tadpole shrimps are suitable subjects for the study of morphological stasis. We were able to obtain samples of three species of Triops (T. granarius, T longicaudatus and T. cancriformis) from four regions in Japan. Taxonomic species were identified by diagnostic morphology. We inferred phylogenetic relationships between individual samples using mitochondrial 16S rRNA. Carapace shapes were compared among populations using shape coordinate methods. The phylogeny inferred from mtDNA shows that T. granarius is phylogenetically more similar to T. longicaudatus than to T. cancriformis. mtDNA sequences did not differ among the populations of T granarius. However, there were two distinct phylogenetic species within T. longicaudatus. In spite of the similarity in fundamental morphological characteristics among Triops species, mtDNA sequences of Triops showed marked differences among the four phylogenetic species. Among the populations of T. granarius , the carapace shape of the Fukuoka population was significantly different from those of other populations of T. granarius. The carapace shape of the Kagawa population of T. longicaudatus was more similar to those of the Shizuoka and Kagawa populations of T. granarius. The shape of the carapace of T. cancriformis was significantly different from those of T. granarius and T longicaudatus. Thus, taxonomic species, phylogenetic species and populations with similar carapace morphology did not correspond with each other. The present results indicate that most of the morphological change did not occur at the time of speciation (lineage separation) and that morphological stases are important evolutionary patterns, but they are not species-level properties.     

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.     

Geometric morphometrics of the forewing shape and size discriminate Plebeia species (Hymenoptera: Apidae) nesting in different substrates

Historically, studies evaluating morphological diversity in stingless bees (Hymenoptera, Apidae: Meliponini) by geometric morphometrics have been used to successfully discriminate taxa and/or populations. Moreover, the use of geometric morphometrics to evaluate phylogenetic morphological variation among stingless bee species has received less attention. Here, we used geometric morphometrics to assess taxonomic discrimination and putative phylogenetic signals for six diapausing stingless bee species (Plebeia) occurring in southern Brazil. In all, 12 landmarks were captured from forewings of P. droryana, P. saiqui, P. emerina, P. remota, P. nigriceps and P. wittmanni. Our data show that the centroid size of the forewings reliably discriminated, for example, between P. droryana and P. emerina from P. saiqui. Moreover, this trait does not have a significant phylogenetic signal. In turn, we found that the overall accuracy in discriminating between the six Plebeia species according to forewing shape was 84%, while the confusion matrix achieved 71%. Interestingly, our discriminant analysis separated Plebeia species nesting in tree cavities from those nesting under granitic rocks. The latter group has second cubital (landmarks = 5, 6, 7), first medial (landmarks = 2, 3, 8) and first submarginal cells (landmarks = 3, 4, 9, 10) that are larger than those of species nesting in trees. The forewing shape showed a strong phylogenetic signal, therefore suggesting that its variation may be due to an evolutionary history shared between Plebeia species studied here rather than to environmental features. This work sheds light on the value of forewing size and shape attributes in discriminating Plebeia species within same genus. We suggest that landmarks separating different taxonomic groups could be incorporated into dichotomous keys to help in identifying clades of complex resolution.