A Step-by-Step Guide for Geometric Morphometrics of Floral Symmetry

Characterizing symmetry transitions in systems with dynamic morphology

The accurate quantification of symmetry is a key goal in biological inquiries because symmetry can affect biological performance and can reveal insights into development and evolutionary history. Recently, we proposed a versatile measure of symmetry, transformation information (TI), which provides an entropy-based measure of deviations from exact symmetry with respect to a parameterized family of transformations. Here we develop this measure further to quantify approximate symmetries and maximal symmetries represented by critical points in TI as a function of a transformation parameter. This framework allows us to characterize the evolution of symmetry by tracking qualitative changes with respect to these critical points. We apply TI to increasingly complex settings, from mathematically tractable probability distributions to differential equation models with emergent behaviors that are inspired by developmental biology and formulated in both static and growing domains. Our analysis of the qualitative changes in symmetry properties indicates a potential pathway toward a general mathematical framework for characterizing symmetry transitions akin to bifurcation theory for dynamical systems. The results reveal deep connections between observed symmetry transitions, subtle changes in morphology, and the underlying mechanisms that govern the dynamics of the system.

From Auto-Kirigami to Drumheads: Suspended, Self-Tearing, and Strained Graphene Nanostructures Formed by Nanoindentation

Nanoindentation of substrate-supported graphene can produce auto-kirigami (AK) structures: spontaneously folded and extended self-tearing nanoribbons up to several micrometers in length. However, the mechanisms governing their formation and yield are poorly understood. Here, we study graphene AK through statistical analysis of high-throughput experiments involving hundreds-fold arrays of indents on highly uniform regions of exfoliated monolayer and bilayer graphene, with no applied oscillation (in contrast with prior work). Post-mortem atomic force microscopy analysis reveals a baseline AK formation rate of 13-61% for monolayers and 0-17% for bilayers depending on inter-indent pitch. Force-distance curves of each type of nanostructure showed no appreciable differences. Moreover, graphene can remain intact after indentation, permitting formation of unbroken graphene suspended over or conformed within indents. Inter-indent pitch affects the absolute and relative formation rates of these nanostructures, attributed to indentation-induced tensile graphene strain. This advances the understanding of mechanisms for controlled formation of nanostructures, including twisted bilayers of graphene and other van der Waals materials.

Deep learning to capture leaf shape in plant images: Validation by geometric morphometrics

Plant leaves play a pivotal role in automated species identification using deep learning (DL). However, achieving reproducible capture of leaf variation remains challenging due to the inherent "black box" problem of DL models. To evaluate the effectiveness of DL in capturing leaf shape, we used geometric morphometrics (GM), an emerging component of eXplainable Artificial Intelligence (XAI) toolkits. We photographed Ranunculus auricomus leaves directly in situ and after herbarization. From these corresponding leaf images, we automatically extracted DL features using a neural network and digitized leaf shapes using GM. The association between the extracted DL features and GM shapes was then evaluated using dimension reduction and covariation models. DL features facilitated the clustering of leaf images by source populations in both in situ and herbarized leaf image datasets, and certain DL features were significantly associated with biological leaf shape variation as inferred by GM. DL features also enabled leaf classification into morpho-phylogenomic groups within the intricate R. auricomus species complex. We demonstrated that simple in situ leaf imaging and DL reproducibly captured leaf shape variation at the population level, while combining this approach with GM provided key insights into the shape information extracted from images by computer vision, a necessary prerequisite for reliable automated plant phenotyping.

Quantitative Analysis of Fish Morphology Through Landmark and Outline-based Geometric Morphometrics with Free Software

Morphology underpins key biological and evolutionary processes that remain elusive. This is in part due to the limitations in robustly and quantitatively analyzing shapes within and between groups in an unbiased and high-throughput manner. Geometric morphometrics (GM) has emerged as a widely employed technique for studying shape variation in biology and evolution. This study presents a comprehensive workflow for conducting geometric morphometric analysis of fish morphology. The step-by-step manual provides detailed instructions for using popular free software, such as the TPS series, MorphoJ, ImageJ, and R, to carry out generalized Procrustes analysis (GPA), principal component analysis (PCA), discriminant function analysis (DFA), canonical variate analysis (CVA), mean shape analysis, and thin plate spline analysis (TPS). The Momocs package in R is specifically utilized for in-depth analysis of fish outlines. In addition, selected functions from the dplyr package are used to assist in the analysis. The full process of fish outline analysis is covered, including extracting outline coordinates, converting and scaling data, defining landmarks, creating data objects, analyzing outline differences, and visualizing results. In conclusion, the current protocol compiles a detailed method for evaluating fish shape variation based on landmarks and outlines. As the field of GM continues to evolve and related software develops rapidly, the limitations associated with morphological analysis of fish are expected to decrease. Interoperable data formats and analytical methods may facilitate the sharing of morphological data and help resolve related scientific problems. The convenience of this protocol allows for fast and effective morphological analysis. Furthermore, this detailed protocol could be adapted to assess image-based differences across a broader range of species or to analyze morphological data of the same species from different origins. Key features • This protocol provides a comprehensive set of commonly used GM-analyzing methods and visualizing skills plus supporting information to help assess the appropriate analysis method • By incorporating both landmarks and outlines, this protocol facilitates a thorough analysis of two-dimensional shape variation in fish, covering a wide range of morphological features • The simplified workflow and detailed procedures make it accessible for non-experienced users to successfully complete the analysis while also providing valuable insights for experienced users Graphical overview Workflow for conducting geometric morphometrics analysis on fish. The steps include image acquisition as data sources, digitization of fish morphology using landmark-based methods, analysis of shape variation characteristics, and visualization of the results in relation to biological interpretation. Largemouth bass (Micropterus salmoides) is used as an example in the schematic representation.

Osteon shape variation in the femoral diaphysis: A geometric-morphometric approach on human cortical bone microstructure in an elderly sample

Geometric morphometrics (GMM) have been applied to understand morphological variation in biological structures. However, research studying cortical bone through geometric histomorphometrics (GHMM) is scarce. This research aims to develop a landmark-based GHMM protocol to depict osteonal shape variation in the femoral diaphysis, exploring the role of age and biomechanics in bone microstructure. Proximal, midshaft, and distal anatomical segments from the femoral diaphysis of six individuals were assessed, with 864 secondary intact osteons from eight periosteal sampling areas being manually landmarked. Observer error was tested using Procrustes ANOVA. Average osteonal shape and anatomical segment-specific variation were explored using principal component analysis. Osteon shape differences between segments were examined using canonical variate analysis (CVA). Sex differences were assessed through Procrustes ANOVA and discriminant function analysis (DFA). The impact of osteonal size on osteonal shape was investigated. High repeatability and reproducibility in osteon shape landmarking were reported. The average osteon shape captured was an elliptical structure, with PC1 reflecting more circular osteons. Significant differences in osteon shape were observed between proximal and distal segments according to CVA. Osteon shape differed between males and females, with DFA showing 52% cross-validation accuracies. No effect of size on shape was reported. Osteonal shape variation observed in this study might be explained by the elderly nature of the sample as well as biomechanical and physiological mechanisms playing different roles along the femoral diaphysis. Although a larger sample is needed to corroborate these findings, this study contributes to the best of our knowledge on human microanatomy, proposing a novel GHMM approach.

Molecular and morphological evidence of hybridization between two dimorphic sympatric species of Fuchsia (Onagraceae)

Hybridization is commonly reported in angiosperms, generally based on morphology, and in few cases confirmed by molecular markers. Fuchsia has a long tradition of ornamental cultivars with different hybrids produced by artificial crosses, so natural hybridization between sympatric Fuchsia species could be common. Natural hybridization between F. microphylla and F. thymifolia was tested using six newly developed microsatellites for F. microphylla in addition to other molecular markers with codominant and maternal inheritance. Geometric morphometrics of leaves and floral structures were also used to identify putative hybrids. Hybrids showed a different degree of genetic admixture between both parental species. Chloroplast DNA (cpDNA) sequences indicated that hybridization occurs in both directions, in fact, some of the hybrids showed new haplotypes for cpDNA and ITS (internal transcriber spacer of nuclear ribosomal RNA genes) sequences. The morphology of hybrid individuals varied between the two parental species, but they could be better identified by their leaves and floral tubes. Our study is the first to confirm the hybridization in natural populations of Fuchsia species and suggests that hybridization has probably occurred repeatedly throughout the entire distribution of the species. Phylogeographic analysis of both species will be essential to understanding the impact of hybridization throughout their complete distribution.

A geometric morphometric analysis of variation in Australian frog ilia and taxonomic interpretations

Anurans including frogs and toads exhibit an ilium that is often regarded as taxonomically diagnostic. The ilium, one of the three paired bones that make up the pelvic girdle, has been important in the fossil record for identifying anuran morphotypes. Osteological collections for Australian frogs are rare in herpetological museums, and skeletonizing whole-bodied specimens requires destroying soft tissue morphology which is valuable to anuran specialists working on living species. Computed tomography scans provide the opportunity to study anuran osteology without the loss of soft tissues. Our study, based on microcomputed tomography scans of extant Australian frogs from the public repository Morphosource and from museum collections focuses on the morphological differences between Australian frogs at the familial and generic levels using geometric morphometrics to compare the diagnostic shape of the ilium. Principal component analysis (PCA) and canonical variate analysis (CVA) were conducted to assess differences in the ilium between supraspecific groups of Australian frogs. The canonical variates analysis accurately predicted group membership (i.e., the correct family) with up to 76.2% success for cross-validated predictions and 100% of original group predictions. While the sample was limited to familial and generic level analyses, our research shows that ilial morphology in Australian frogs is taxonomically informative. This research provides a guide for identifying Australian anurans, including fossils, as well as new information relevant to considerations about their phylogenetic relationships, and the potential use of the fossil record to enhance efforts to conserve threatened living frog species.

Geometric Morphometrics sheds light on the systematics affinities of two enigmatic dwarf Neotropical sedges (Carex, Cyperaceae)

Geometric morphometrics (GM) is a powerful analytical tool that enables complete quantification of shapes. Its use in Botany has a great potential for complementing plant evolutionary and ecological studies. Taxonomic delimitation in Carex has been complicated due to reduction of characters and frequent homoplasy. This problem is more marked in cases where the species exhibit dwarfism. South America is the continent with the least understood Carex flora. The systematic relationships of some bizarre-looking groups were not unraveled until molecular phylogenetic studies resolved their relationships. In particular, there are two species only known from their type material whose affinities remain uncertain: Carexherteri and C.hypsipedos. These two taxa are acaulescent plants that respectively grow in the Uruguayan pampa and Peruvian high-altitude meadows. Recently, both species were ascribed to the Carexphalaroides group (subgen. Psyllophorae, sect. Junciformes) due to superficial morphological similarities, such as the androgynous peduncled spikes. However, their character combination is also coincident for its circumscription to sect. Abditispicae species. Nevertheless, in the absence of confirmation from molecular analyses, their placement must be considered preliminary until additional data can be provided. In this work we employ for the first time geometric morphometrics (GM) tools to assess the systematic affinities of two taxonomically problematic sedge species based on fruit shape. We compared utricle morphology of C.herteri and C.hypsipedos with that of C.phalaroides group and species in sect. Abditispicae. To this end we used GM and traditional morphometric approaches. Utricle shape variation along with other morphological features support the exclusion of these two species from the C.phalaroides gr. and, at the same time, show clear affinities of C.herteri to sect. Abditispicae. Carexhypsipedos remains as an incertae sedis species. Our work shows the potential utility of GM for the exploration of systematic affinities in sedges and in other graminoids.

3D pollination biology using micro-computed tomography and geometric morphometrics in Theobroma cacao

Premise

Imaging technologies that capture three-dimensional (3D) variation in floral morphology at micro- and nano-resolutions are increasingly accessible. In herkogamous flowers, such as those of Theobroma cacao, structural barriers between anthers and stigmas represent bottlenecks that restrict pollinator size and access to reproductive organs. To study the unresolved pollination biology of cacao, we present a novel application of micro-computed tomography (micro-CT) using floral dimensions to quantify pollinator functional size limits.

Methods

We generated micro-CT data sets from field-collected flowers and museum specimens of potential pollinators. To compare floral variation, we used 3D Slicer to place landmarks on the surface models and performed a geometric morphometric (GMM) analysis using geomorph R. We identified the petal side door (an opening between the petal hoods and filament) as the main bottleneck for pollinator access. We compared its mean dimensions with proposed pollinators to identify viable candidates.

Results

We identified three levels of likelihood for putative pollinators based on the number of morphological (body) dimensions that fit through the petal side door. We also found floral reward microstructures whose presence and location were previously unclear.

Discussion

Using micro-CT and GMM to study the 3D pollination biology of cacao provides new evidence for predicting unknown pollinators. Incorporating geometry and floral rewards will strengthen plant-pollinator trait matching models for cacao and other species.

A novel method for assessment of human midpalatal sutures using CBCT-based geometric morphometrics and complexity scores

INTRODUCTION

Management of dentofacial deficiencies requires knowledge about sutural morphology and complexity. The present study assesses midpalatal sutural morphology based on human cone-beam computed tomography (CBCT) using geometric morphometrics (GMM) and complexity scores. The study is the first to apply a sutural complexity score to human CBCT datasets and demonstrates the potential such a score has to improve objectiveness and comparability when analysing the midpalatal suture.

MATERIALS AND METHODS

CBCTs of various age and sex groups were analysed retrospectively (n = 48). For the geometric morphometric analysis, landmark acquisition and generalised Procrustes superimposition were combined with principal component analysis to detect variability in sutural shape patterns. For complexity analysis, a windowed short-time Fourier transform with a power spectrum density (PSD) calculation was applied to resampled superimposed semi-landmarks.

RESULTS

According to the GMM, younger patients exhibited comparable sutural patterns. With increasing age, the shape variation increased among the samples. The principal components did not sufficiently capture complexity patterns, so an additional methodology was applied to assess characteristics such as sutural interdigitation. According to the complexity analysis, the average PSD complexity score was 1.465 (standard deviation = 0.010). Suture complexity increased with patient age (p < 0.0001), but was not influenced by sex (p = 0.588). The intra-class correlation coefficient exceeded 0.9, indicating intra-rater reliability.

CONCLUSION

Our study demonstrated that GMM applied to human CBCTs can reveal shape variations and allow the comparison of sutural morphologies across samples. We demonstrate that complexity scores can be applied to study human sutures captured in CBCTs and complement GMM for a comprehensive sutural analysis.

Secular Trends in the Size and Shape of the Scapula among the Portuguese between the 19th and the 21st Centuries

Potential secular changes in the human scapula are fundamentally unbeknownst, with most of the preceding anatomical studies focusing on long-term changes in the long bones and the skull. As such, the cardinal purpose of this study pertains to the evaluation of secular trends on the shape and size of the scapula in a time period spanning from the 19th to the early 21st centuries. The study sample included 211 individuals (100 males and 111 females) from the Coimbra Identified Skeletal Collection and the 21st Century Identified Skeletal Collection. The size and shape of the scapula were evaluated using geometric morphometrics. Results show secular changes over a relatively short period of time in both the shape and size of the scapula in Portuguese nationals. Shape changes were observed in both sexes but expressed minimally, while a significant negative trend in the size of the scapula was detected in males. Scapular size decrement in males conceivably echoes general trends of the overall anatomy towards a narrower body associated with higher standards of living that include enhanced nutrition and universal healthcare, among other factors.

Patterned proliferation orients tissue-wide stress to control root vascular symmetry in Arabidopsis

Symmetric tissue alignment is pivotal to the functions of plant vascular tissue, such as long-distance molecular transport and lateral organ formation. During the vascular development of the Arabidopsis roots, cytokinins initially determine cell-type boundaries among vascular stem cells and subsequently promote cell proliferation to establish vascular tissue symmetry. Although it is unknown whether and how the symmetry of initially defined boundaries is progressively refined under tissue growth in plants, such boundary shapes in animal tissues are regulated by cell fluidity, e.g., cell migration and intercalation, lacking in plant tissues. Here, we uncover that cell proliferation during vascular development produces anisotropic compressive stress, smoothing, and symmetrizing cell arrangement of the vascular-cell-type boundary. Mechanistically, the GATA transcription factor HANABA-TARANU cooperates with the type-B Arabidopsis response regulators to form an incoherent feedforward loop in cytokinin signaling. The incoherent feedforward loop fine-tunes the position and frequency of vascular cell proliferation, which in turn restricts the source of mechanical stress to the position distal and symmetric to the boundary. By combinatorial analyses of mechanical simulations and laser cell ablation, we show that the spatially constrained environment of vascular tissue efficiently entrains the stress orientation among the cells to produce a tissue-wide stress field. Together, our data indicate that the localized proliferation regulated by the cytokinin signaling circuit is decoded into a globally oriented mechanical stress to shape the vascular tissue symmetry, representing a reasonable mechanism controlling the boundary alignment and symmetry in tissue lacking cell fluidity.

On-Demand Dynamic Chirality Selection in Flower Corolla-like Micropillar Arrays

Chiral morphology has been intensively studied in various fields including biology, organic chemistry, pharmaceuticals, and optics. On-demand and dynamic chiral inversion not only cannot be realized in most intrinsically chiral materials but also has mostly been limited to chemical or light-induced methods. Herein, we report reversible real-time magneto-mechanical chiral inversion of a three-dimensional (3D) micropillar array between achiral, clockwise, and counterclockwise chiral arrangements. Inspired by the flower corolla, achiral arrays of five and six radially arranged semicylindrical micropillars were employed as model systems to investigate the dynamic symmetry properties of arrays consisting of odd and even numbers of micropillars, respectively. Each micropillar underwent twisting actuation with a different twisting angle depending on the angle with the magnetic field direction and magnetic flux density, thereby collectively changing the chirality from the achiral to chiral state. Importantly, the morphological handedness of the micropillars was inverted within a few seconds by manipulating the direction of the magnetic field. A chiral morphology consisting of magnetically twisted micropillars was shape-fixed by the introduction of a polymeric binder. This binder could be simply washed off to return the shape-fixed twisted micropillars to their initial straight state. Magnetically programmable and reproducible 3D flower corolla-like micropillar arrays are expected to expand the potential of shape-reconfigurable devices that require real-time chiral manipulation in ambient environments.

Shape asymmetry - what's new?

Studies of shape asymmetry have become increasingly abundant as the methods of geometric morphometrics have gained widespread use. Most of these studies have focussed on fluctuating asymmetry and have largely obtained similar results as more traditional analyses of asymmetry in distance measurements, but several notable differences have also emerged. A key difference is that shape analyses provide information on the patterns, not just the amount of variation, and therefore tend to be more sensitive. Such analyses have shown that apparently symmetric structures in animals consistently show directional asymmetry for shape, but not for size. Furthermore, the long-standing prediction that phenotypic plasticity in response to environmental heterogeneity can contribute to fluctuating asymmetry has been confirmed for the first time for the shape of flower parts (but not for size). Finally, shape analyses in structures with complex symmetry, such as many flowers, can distinguish multiple types of directional asymmetry, generated by distinct direction-giving factors, which combine to the single component observable in bilaterally symmetric structures. While analyses of shape asymmetry are broadly compatible with traditional analyses of asymmetry, they incorporate more detailed morphological information, particularly for structures with complex symmetry, and therefore can reveal subtle biological effects that would otherwise not be apparent. This makes them a promising tool for a wide range of studies in the basic and applied life sciences.

Directional asymmetry and direction-giving factors: Lessons from flowers with complex symmetry

Directional asymmetry is a systematic difference between the left and right sides for structures with bilateral symmetry or a systematic differentiation among repeated parts for complex symmetry. This study explores factors that produce directional asymmetry in the flower of Iris pumila, a structure with complex symmetry that makes it possible to investigate multiple such factors simultaneously. The shapes and sizes of three types of floral organs, the falls, standards, and style branches, were quantified using the methods of geometric morphometrics. For each flower, this study recorded the compass orientations of floral organs as well as their anatomical orientations relative to the two spathes subtending each flower. To characterize directional asymmetry at the whole-flower level, differences in the average sizes and shapes according to compass orientation and relative orientation were computed, and the left-right asymmetry was also evaluated for each individual organ. No size or shape differences within flowers were found in relation to anatomical position; this may relate to the terminal position of flowers in Iris pumila, suggesting that there may be no adaxial-abaxial polarity, which is very prominent in many other taxa. There was clear directional asymmetry of shape in relation to compass orientation, presumably driven by a consistent environmental gradient such as solar irradiance. There was also clear directional asymmetry between left and right halves of every floral organ, most likely related to the arrangement of organs in the bud. These findings indicate that different factors are acting to produce directional asymmetry at different levels. In conventional analyses not recording flower orientations, these effects would be impossible to disentangle from each other and would probably be included as part of fluctuating asymmetry.

Palatal morphology predicts the paleobiology of early salamanders

Ecological preferences and life history strategies have enormous impacts on the evolution and phenotypic diversity of salamanders, but the yet established reliable ecological indicators from bony skeletons hinder investigations into the paleobiology of early salamanders. Here, we statistically demonstrate by using time-calibrated cladograms and geometric morphometric analysis on 71 specimens in 36 species, that both the shape of the palate and many non-shape covariates particularly associated with vomerine teeth are ecologically informative in early stem- and basal crown-group salamanders. Disparity patterns within the morphospace of the palate in ecological preferences, life history strategies, and taxonomic affiliations were analyzed in detail, and evolutionary rates and ancestral states of the palate were reconstructed. Our results show that the palate is heavily impacted by convergence constrained by feeding mechanisms and also exhibits clear stepwise evolutionary patterns with alternative phenotypic configurations to cope with similar functional demand. Salamanders are diversified ecologically before the Middle Jurassic and achieved all their present ecological preferences in the Early Cretaceous. Our results reveal that the last common ancestor of all salamanders share with other modern amphibians a unified biphasic ecological preference, and metamorphosis is significant in the expansion of ecomorphospace of the palate in early salamanders.

Phylogenetic, ecological and intraindividual variability patterns in grass phytolith shape

BACKGROUND AND AIMS

Grass silica short cell (GSSC) phytoliths appear to be the most reliable source of fossil evidence for tracking the evolutionary history and paleoecology of grasses. In recent years, modern techniques that quantitatively assess phytolith shape variation have widened opportunities for the classification of grass fossil phytoliths. However, phylogenetic, ecological and intraindividual variability patterns in phytolith shape remain largely unexplored.

METHODS

The full range of intraindividual phytolith shape variation [3650 two-dimensional (2-D) outlines] from 73 extant grass species, 48 genera, 18 tribes and eight subfamilies (particularly Pooideae) was analysed using geometric morphometric analysis based on semi-landmarks spanning phytolith outlines.

KEY RESULTS

The 2-D phytolith shape is mainly driven by deep-time diversification of grass subfamilies. There is distinct phytolith shape variation in early-diverging lineages of Pooideae (Meliceae, Stipeae). The amount of intraindividual variation in phytolith shape varies among species, resulting in a remarkable pattern across grass phylogeny.

CONCLUSIONS

The phylogenetic pattern in phytolith shape was successfully revealed by applying geometric morphometrics to 2-D phytolith shape outlines, strengthening the potential of phytoliths to track the evolutionary history and paleoecology of grasses. Geometric morphometrics of 2-D phytolith shape is an excellent tool for analysis requiring large numbers of phytolith outlines, making it useful for quantitative palaeoecological reconstruction.

Geometric Morphometric Analysis of the Humerus in New and Old World Vultures

The vulture guild is composed of two distinct groups, Old and New World, that provide a unique insight into how morphology varies among convergent species. All vultures are considered to be large birds of prey that utilize a style of flight called thermal soaring to search and feed primarily on carrion. Even though this flight style is exhibited among all 23 species, slight variations in their skeletal morphology may relate to their differences in ecology. We hypothesized that vulture humeral morphology varies in relation to these organisms' habitat, average body mass, courtship displays, and migratory behavior. To address this hypothesis, we used three-dimensional geometric morphometrics to measure the overall shape differences of vulture humeri. Humeral morphology was found to vary most by habitat association and migratory tendency. The humeri of vultures that inhabit forested areas exhibit features that suggest increased flapping flight compared to those in open and mountainous regions. Migratory species were found to possess more robust features near the glenohumeral joint. We found these (and other features) have some utility for predicting ecology and behavior, but we suggest that further investigation into skeletal and muscular wing elements may reveal greater understanding of the habits of extinct vulture species. This article is protected by copyright. All rights reserved.

Patterns of leaf morphological variation in Quercus frainetto Ten. growing on different soil types in Serbia

Leaf morphology is at a certain level defined by the ways in which plants adapt to different habitats, especially in large trees. In this study, morphological variations in leaf size and shape of the Hungarian oak (Quercus frainetto Ten.) growing on different soil types (lithic leptosol, vertisol, cambisol) were investigated in the central part of Serbia (Sumadija). The information on soil type was obtained using a digitalized soil map of the Republic of Serbia, while leaf traits were characterized by geometric morphometric methods. Landmark analysis and leaf measurements showed significant differences among the analyzed groups, with individuals growing on nutrient-poor, shallow soils having smaller leaves with greater lobation. The observed differences suggest that the levels of soil productivity influence variations in leaf patterns. More studies on a larger sample size and along a broader spatial scale are needed to fully understand the differences in the patterns of leaf morphological variation in Q. frainetto.

Identification of approximate symmetries in biological development

Virtually all forms of life, from single-cell eukaryotes to complex, highly differentiated multicellular organisms, exhibit a property referred to as symmetry. However, precise measures of symmetry are often difficult to formulate and apply in a meaningful way to biological systems, where symmetries and asymmetries can be dynamic and transient, or be visually apparent but not reliably quantifiable using standard measures from mathematics and physics. Here, we present and illustrate a novel measure that draws on concepts from information theory to quantify the degree of symmetry, enabling the identification of approximate symmetries that may be present in a pattern or a biological image. We apply the measure to rotation, reflection and translation symmetries in patterns produced by a Turing model, as well as natural objects (algae, flowers and leaves). This method of symmetry quantification is unbiased and rigorous, and requires minimal manual processing compared to alternative measures. The proposed method is therefore a useful tool for comparison and identification of symmetries in biological systems, with potential future applications to symmetries that arise during development, as observed in vivo or as produced by mathematical models. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Characteristics of the Donkey's Dorsal Profile in Relation to Its Functional Body Condition Assessment

Simple Summary

Even though animal posture is known to reflect an emotional state, the presence of chronic postures associated with poor welfare has already been investigated in horses. Measuring, quantifying, and comparing postures requires an application of an objective tool, such as geometric morphometrics. This study aimed to use geometric morphometrics to characterize the dorsal profile of donkeys in relation to their individual features. Forty donkeys were photographed and characterized using the body condition score (BCS), fatty neck score (FNS), dental condition score (DCS), sex, and breed. Then, photographs were analyzed using geometric morphometrics and the differences in dorsal profile between the examined criteria were tested. Within the entire set of donkey dorsal profiles, deformation related to BCS and FNS was observed. FNS measurement seems to have the strongest influence on a donkey’s dorsal profile among the examined criteria. Concluding, the donkeys’ body condition affects their dorsal profile, and both FNS and BCS measurements should be considered when the donkeys’ dorsal profiles are investigated. However, to evaluate the link between the dorsal profile and the welfare state of donkeys, more studies are required.

Abstract

As the breeding of donkeys has increased due to different types of use, welfare evaluation importance increases. This equid’s welfare state has been described using body condition indicators and the geometric morphometrics method. However, the dorsal profile has not yet been assessed in donkeys. In this study, the body condition score (BCS), fatty neck score (FNS), dental condition score (DCS), sex, and breed were used as criteria of dorsal profile deformations. Photographs of 40 donkeys were analyzed using geometric morphometrics. Within the entire set of dorsal profiles, the variance of the first three principal components (PCs) was PC1 = 37.41%, PC2 = 23.43%, and PC3 = 13.34%. The dorsal profiles displayed deformation as an effect of FNS and BCS on size (FNS p = 0.012; BCS p = 0.024) and shape (FNS p < 0.0001; BCS p < 0.0001), rather than as an effect of DCS (p < 0.0001), sex (p = 0.0264), and breed (p < 0.0001) only on shape. The highest distances among the categories (Mahalanobis distances: MD ≥ 13.26; Procrustes distances: PD ≥ 0.044) were noted for FNS. The lowest distances were noted between jennets and males (MD = 4.58; PD = 0.012) and between BCS 1 and BCS 2 (MD = 4.70; PD = 0.018). Donkeys’ body condition affects their dorsal profile and both FNS and BCS measurements should be considered when a donkey’s dorsal profile is investigated.

Ancient Starch Remains Reveal the Vegetal Diet of the Neolithic Late Dawenkou Culture in Jiangsu, East China

The Liangwangcheng site, located in Pizhou County, Xuzhou City, northern Jiangsu Province, is one of the most important Neolithic Dawenkou Culture archeological sites in the Haidai area of China’s eastern seaboard. In recent years, archaeobotanical studies in the Haidai area, mainly focusing on Shandong Province, have yielded fruitful results, while relatively few such studies have been undertaken in northern Jiangsu Province. Here, we report the results of dental residue analysis conducted on 31 individual human skulls unearthed from the Late Dawenkou Culture Liangwangcheng site. The starch granules extracted from these residue samples indicate that foxtail and broomcorn millet, rice, roots and tubers, and legumes comprised the vegetal diet of Liangwangcheng’s occupants. Evidence suggests that mixed rice–millet agriculture played a definite role, with the coexistence of gathering as an economic element. According to archaeobotanical evidence from surrounding cotemporaneous sites, the Late Neolithic human groups that lived in the lower Huang-Huai River drainage shared similar subsistence patterns. Our results provide new evidence for a more comprehensive understanding of plant resource utilization and agricultural development in northern Jiangsu during the Dawenkou period.

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones

PREMISE

The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data.

METHODS

We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape.

RESULTS

Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations.

CONCLUSIONS

The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

Three-dimensional quantification of twisting in the Arabidopsis petiole

Organisms have a variety of three-dimensional (3D) structures that change over time. These changes include twisting, which is 3D deformation that cannot happen in two dimensions. Twisting is linked to important adaptive functions of organs, such as adjusting the orientation of leaves and flowers in plants to align with environmental stimuli (e.g. light, gravity). Despite its importance, the underlying mechanism for twisting remains to be determined, partly because there is no rigorous method for quantifying the twisting of plant organs. Conventional studies have relied on approximate measurements of the twisting angle in 2D, with arbitrary choices of observation angle. Here, we present the first rigorous quantification of the 3D twisting angles of Arabidopsis petioles based on light sheet microscopy. Mathematical separation of bending and twisting with strict definition of petiole cross-sections were implemented; differences in the spatial distribution of bending and twisting were detected via the quantification of angles along the petiole. Based on the measured values, we discuss that minute degrees of differential growth can result in pronounced twisting in petioles.

Male sterility significantly elevates shape variation and fluctuating asymmetry of zygomorphic corolla in gynodioecious Glechoma hederacea (Lamiaceae)

Female flowers of gynodioecious plants usually have smaller corollas than bisexual flowers. This difference is explained as a developmental consequence of stamen abortion and as a result of stronger selection for larger corolla size in hermaphroditic flowers that solely ensure male function within populations. This study evaluated whether the size difference of zygomorphic corollas in a widely distributed gynodioecious herb Glechoma hederacea is accompanied by variation in shape and bilateral fluctuating asymmetry of sexually differentiated flowers. Geometric morphometric analyses of bilateral symmetry in the shapes of corolla lower lips showed that male-sterile flowers were significantly more plastic and asymmetric, implying that they may be subjected to weaker stabilizing selection for corolla shape in comparison to hermaphrodites. These results illustrated that sexual differentiation may be an important factor contributing to bilateral fluctuating asymmetry in the shape of zygomorphic flowers.

Morphometric Analysis of Coptotermes spp. Soldier Caste (Blattodea: Rhinotermitidae) in Indonesia and Evidence of Coptotermes gestroi Extreme Head-Capsule Shapes

Linear and geometric morphometrics approaches were conducted to analyze the head capsule (HC) shape of collected soldier caste specimens of Coptotermes from various locations in Indonesia. The soldiers' morphology was observed and measured. The results of the principal component analysis of the group of all species showed two important groups of variables, i.e., the body size and setae characteristics of the pronotum and head. The multicollinearity of the morphometric variables showed the importance of body measurements as well as important alternative characteristics such as the pronotum setae (PrS) and HC setae. Four trends of HC shape were observed across the species. Interestingly, three extreme shapes were depicted by geometric morphometrics of the C. gestroi HC. The phylogenetic tree inferred from 12S and 16S mitochondrial gene fragments showed high confidence for C. gestroi populations. The lateral expansion of the posterior part of the HC across the species was in accordance with the increasing of the number of hairlike setae on the pronotum and HC. These differences among species might be associated with mandible-force-related defensive labor and sensitivity to environmental stressors.

Inter- and intraspecific variation in grass phytolith shape and size: a geometric morphometrics perspective

BACKGROUND AND AIMS

The relative contributions of inter- and intraspecific variation to phytolith shape and size have only been investigated in a limited number of studies. However, a detailed understanding of phytolith variation patterns among populations or even within a single plant specimen is of key importance for the correct taxonomic identification of grass taxa in fossil samples and for the reconstruction of vegetation and environmental conditions in the past. In this study, we used geometric morphometric analysis for the quantification of different sources of phytolith shape and size variation.

METHODS

We used landmark-based geometric morphometric methods for the analysis of phytolith shapes in two extant grass species (Brachypodium pinnatum and B. sylvaticum). For each species, 1200 phytoliths were analysed from 12 leaves originating from six plants growing in three populations. Phytolith shape and size data were subjected to multivariate Procrustes analysis of variance (ANOVA), multivariate regression, principal component analysis and linear discriminant analysis.

KEY RESULTS

Interspecific variation largely outweighed intraspecific variation with respect to phytolith shape. Individual phytolith shapes were classified with 83 % accuracy into their respective species. Conversely, variation in phytolith shapes within species but among populations, possibly related to environmental heterogeneity, was comparatively low.

CONCLUSIONS

Our results imply that phytolith shape relatively closely corresponds to the taxonomic identity of closely related grass species. Moreover, our methodological approach, applied here in phytolith analysis for the first time, enabled the quantification and separation of variation that is not related to species discrimination. Our findings strengthen the role of grass phytoliths in the reconstruction of past vegetation dynamics.

Insect morphometry is reproducible under average investigation standards

Morphometric research is being applied to a growing number and variety of organisms. Discoveries achieved via morphometric approaches are often considered highly transferable, in contrast to the tacit and idiosyncratic interpretation of discrete character states. The reliability of morphometric workflows in insect systematics has never been a subject of focused research, but such studies are sorely needed. In this paper, we assess the reproducibility of morphometric studies of ants where the mode of data collection is a shared routine.We compared datasets generated by eleven independent gaugers, that is, collaborators, who measured 21 continuous morphometric traits on the same pool of individuals according to the same protocol. The gaugers possessed a wide range of morphometric skills, had varying expertise among insect groups, and differed in their facility with measuring equipment. We used intraclass correlation coefficients (ICC) to calculate repeatability and reproducibility values (i.e., intra- and intergauger agreements), and we performed a multivariate permutational multivariate analysis of variance (PERMANOVA) using the Morosita index of dissimilarity with 9,999 iterations.The calculated average measure of intraclass correlation coefficients of different gaugers ranged from R = 0.784 to R = 0.9897 and a significant correlation was found between the repeatability and the morphometric skills of gaugers (p = 0.016). There was no significant association with the magnification of the equipment in the case of these rather small ants. The intergauger agreement, that is the reproducibility, varied between R = 0.872 and R = 0.471 (mean R = 0.690), but all gaugers arrived at the same two-species conclusion. A PERMANOVA test revealed no significant gauger effect on species identity (R 2 = 0.69, p = 0.58).Our findings show that morphometric studies are reproducible when observers follow the standard protocol; hence, morphometric findings are widely transferable and will remain a valuable data source for alpha taxonomy.

Linking Genes to Shape in Plants Using Morphometrics

A transition from qualitative to quantitative descriptors of morphology has been facilitated through the growing field of morphometrics, representing the conversion of shapes and patterns into numbers. The analysis of plant form at the macromorphological scale using morphometric approaches quantifies what is commonly referred to as a phenotype. Quantitative phenotypic analysis of individuals with contrasting genotypes in turn provides a means to establish links between genes and shapes. The path from a gene to a morphological phenotype is, however, not direct, with instructive information progressing both across multiple scales of biological complexity and through nonintuitive feedback, such as mechanical signals. In this review, we explore morphometric approaches used to perform whole-plant phenotyping and quantitative approaches in capture processes in the mesoscales, which bridge the gaps between genes and shapes in plants. Quantitative frameworks involving both the computational simulation and the discretization of data into networks provide a putative path to predicting emergent shape from underlying genetic programs.

Symmetric and Asymmetric Components of Shape Variation in the Diatom Genus Frustulia (Bacillariophyta)

Irregularities in cell division can produce asymmetry in symmetric structures, such as outlines of diatom cells, which can reflect genetic, environmental, or random variability in developmental processes. This study examined 12 phylogenetic lineages of the diatom genus Frustulia using landmark-based geometric morphometrics to assess the variation between cell segments separated by apical and transapical axes. Although asymmetric variation within cells differed in some lineages, these irregularities most likely did not reflect the evolutionary history of the lineages. The intraclonal phenotypic plasticity of diatom frustules was induced rather by nongenetic factors, i.e., inherited valve abnormalities, constraints of siliceous cell walls, and random developmental instability during morphogenesis. The positive correlations between the symmetric and asymmetric components of shape variation indicated that the morphogenesis of diatom cell walls affected irregularities within cells and variability among the symmetrized cells to a similar extent.

Within-individual phenotypic plasticity in flowers fosters pollination niche shift

Phenotypic plasticity, the ability of a genotype of producing different phenotypes when exposed to different environments, may impact ecological interactions. We study here how within-individual plasticity in Moricandia arvensis flowers modifies its pollination niche. During spring, this plant produces large, cross-shaped, UV-reflecting lilac flowers attracting mostly long-tongued large bees. However, unlike most co-occurring species, M. arvensis keeps flowering during the hot, dry summer due to its plasticity in key vegetative traits. Changes in temperature and photoperiod in summer trigger changes in gene expression and the production of small, rounded, UV-absorbing white flowers that attract a different assemblage of generalist pollinators. This shift in pollination niche potentially allows successful reproduction in harsh conditions, facilitating M. arvensis to face anthropogenic perturbations and climate change.

A Binary-Based Matrix Model for Malus Corolla Symmetry and Its Variational Significance

Floral symmetry (corolla symmetry) has important biological significance in plant genetics and evolution. However, it is often multi-dimensional and difficult to quantify. Here, we constructed a multi-dimensional data matrix [X Y Z] by extracting three qualitative variables with binary properties (X: corolla regularity of interval and coplanarity; Y: petal regularity of shape and size; Z: petal local regularity of curling and wrinkle) from different dimensions of petals (overall to individual, and then to the local): all petals (corolla), individual petals, and local areas of petals. To quantitatively express the degree of Malus corolla symmetry, these variables were then combined with weight assignments (X: 2² > Y: 2¹ > Z: 2⁰) based on their contributions to the corolla symmetry and the algorithm rule of converting binary to decimal values, which facilitated the unification of qualitative and quantitative analyses. Our results revealed significant reduction in degrees of Malus corolla symmetry along the direction of local to overall. Species showed higher degree of corolla symmetry than cultivars; however, taxa with stronger corolla symmetry might not necessarily be species. These findings provide new insights into the circumscription of Malus controversial species. The matrix model should be reference for future evaluation of angiosperm flower symmetry (lack of corolla fusion).

Capturing variation in floral shape: a virtual3D based morphospace for Pelargonium

Background

Variation in floral shapes has long fascinated biologists and its modelling enables testing of evolutionary hypotheses. Recent comparative studies that explore floral shape have largely ignored 3D floral shape. We propose quantifying floral shape by using geometric morphometrics on a virtual3D model reconstructed from 2D photographical data and demonstrate its performance in capturing shape variation.

Methods

This approach offers unique benefits to complement established imaging techniques (i) by enabling adequate coverage of the potential morphospace of large and diverse flowering-plant clades; (ii) by circumventing asynchronicity in anthesis of different floral parts; and (iii) by incorporating variation in copy number of floral organs within structures. We demonstrate our approach by analysing 90 florally-diverse species of the Southern African genus Pelargonium (Geraniaceae). We quantify Pelargonium floral shapes using 117 landmarks and show similarities in reconstructed morphospaces for nectar tube, corolla (2D datasets), and a combined virtual3D dataset.

Results

Our results indicate that Pelargonium species differ in floral shape, which can also vary extensively within a species. PCA results of the reconstructed virtual3D floral models are highly congruent with the separate 2D morphospaces, indicating it is an accurate, virtual, representation of floral shape. Through our approach, we find that adding the third dimension to the data is crucial to accurately interpret the manner of, as well as levels of, shape variation in flowers.

Gynodioecy in the common spindle tree (Euonymus europaeus L.) involves differences in the asymmetry of corolla shapes between sexually differentiated flowers

Gynodioecy is typically associated with a smaller perianth size in purely pistillate flowers than in hermaphrodite flowers. However, it is unclear whether this size differentiation is associated with any differences in flower shape between the two sexual groups. A geometric morphometric analysis of the symmetry of tetrameric corolla shapes was used in the study of Euonymus europaeus L., Darwin’s classical system of floral sexual differentiation. I investigated whether there are any shape differences between the female and bisexual flowers, with respect to both purely symmetric variation involving coordinated shape changes of the four petals and asymmetry among petals within flowers. The corolla shapes of the female and bisexual flowers and the variability among flowers within each sexual group were very similar in the purely symmetric components of shape variation. However, the female flowers were significantly more asymmetric with respect to both the lateral and transversal asymmetry of their corolla shapes. This is the first study to apply geometric morphometrics in the analysis of morphological patterns in a sexually differentiated gynodioecious plant system. The results showed that subtle shape differences in corolla asymmetry differ between the sexual groups and indicate diverging developmental or selection signals between the sexes.

Morphological Analysis of Size and Shape (MASS): An integrative software program for morphometric analyses of leaves

PREMISE

Morphometric analysis is a common approach for comparing and categorizing botanical samples; however, completing a suite of analyses using existing tools may require a multi-stage, multi-program process. To facilitate streamlined analysis within a single program, Morphological Analysis of Size and Shape (MASS) for leaves was developed. Its utility is demonstrated using exemplar leaf samples from Acer saccharum, Malus domestica, and Lithospermum.

METHODS

Exemplar samples were obtained from across a single tree (Acer saccharum), three trees in the same species (Malus domestica), and online, digitized herbarium specimens (Lithospermum). MASS was used to complete simple geometric measurements of samples, such as length and area, as well as geometric morphological analyses including elliptical Fourier and Procrustes analyses. Principal component analysis (PCA) of data was also completed within the same program.

RESULTS

MASS is capable of making desired measurements and analyzing traditional morphometric data as well as landmark and outline data.

DISCUSSION

Using MASS, differences were observed among leaves of the three studied taxa, but only in Malus domestica were differences statistically significant or correlated with other morphological features. In the future, MASS could be applied for analysis of other two-dimensional organs and structures. MASS is available for download at https://github.com/gillianlynnryan/MASS.

Geometric morphometrics of nested symmetries unravels hierarchical inter- and intra-individual variation in biological shapes

Symmetry is a pervasive feature of organismal shape and the focus of a large body of research in Biology. Here, we consider complex patterns of symmetry where a phenotype exhibits a hierarchically structured combination of symmetries. We extend the Procrustes ANOVA for the analysis of nested symmetries and the decomposition of the overall morphological variation into components of symmetry (among-individual variation) and asymmetry (directional and fluctuating asymmetry). We illustrate its use with the Aristotle's lantern, the masticatory apparatus of 'regular' sea urchins, a complex organ displaying bilateral symmetry nested within five-fold rotational symmetry. Our results highlight the importance of characterising the full symmetry of a structure with nested symmetries. Higher order rotational symmetry appears strongly constrained and developmentally stable compared to lower level bilateral symmetry. This contrast between higher and lower levels of asymmetry is discussed in relation to the spatial pattern of the lantern morphogenesis. This extended framework is applicable to any biological object exhibiting nested symmetries, regardless of their type (e.g., bilateral, rotational, translational). Such cases are extremely widespread in animals and plants, from arthropod segmentation to angiosperm inflorescence and corolla shape. The method therefore widens the research scope on developmental instability, canalization, developmental modularity and morphological integration.

Geometric morphometrics of nested symmetries unravels hierarchical inter- and intra-individual variation in biological shapes

Symmetry is a pervasive feature of organismal shape and the focus of a large body of research in Biology. Here, we consider complex patterns of symmetry where a phenotype exhibits a hierarchically structured combination of symmetries. We extend the Procrustes ANOVA for the analysis of nested symmetries and the decomposition of the overall morphological variation into components of symmetry (among-individual variation) and asymmetry (directional and fluctuating asymmetry). We illustrate its use with the Aristotle’s lantern, the masticatory apparatus of ‘regular’ sea urchins, a complex organ displaying bilateral symmetry nested within five-fold rotational symmetry. Our results highlight the importance of characterising the full symmetry of a structure with nested symmetries. Higher order rotational symmetry appears strongly constrained and developmentally stable compared to lower level bilateral symmetry. This contrast between higher and lower levels of asymmetry is discussed in relation to the spatial pattern of the lantern morphogenesis. This extended framework is applicable to any biological object exhibiting nested symmetries, regardless of their type (e.g., bilateral, rotational, translational). Such cases are extremely widespread in animals and plants, from arthropod segmentation to angiosperm inflorescence and corolla shape. The method therefore widens the research scope on developmental instability, canalization, developmental modularity and morphological integration.