Testing the accuracy of 3D automatic landmarking via genome-wide association studies.

Various advances in 3D automatic phenotyping and landmark-based geometric morphometric methods have been made. While it is generally accepted that automatic landmarking compromises the capture of the biological variation, no studies have directly tested the actual impact of such landmarking approaches in analyses requiring a large number of specimens and for which the precision of phenotyping is crucial to extract an actual biological signal adequately. Here, we use a recently developed 3D atlas-based automatic landmarking method to test its accuracy in detecting QTLs associated with craniofacial development of the house mouse skull and lower jaws for a large number of specimens (circa 700) that were previously phenotyped via a semiautomatic landmarking method complemented with manual adjustment. We compare both landmarking methods with univariate and multivariate mapping of the skull and the lower jaws. We find that most significant SNPs and QTLs are not recovered based on the data derived from the automatic landmarking method. Our results thus confirm the notion that information is lost in the automated landmarking procedure although somewhat dependent on the analyzed structure. The automatic method seems to capture certain types of structures slightly better, such as lower jaws whose shape is almost entirely summarized by its outline and could be assimilated as a 2D flat object. By contrast, the more apparent 3D features exhibited by a structure such as the skull are not adequately captured by the automatic method. We conclude that using 3D atlas-based automatic landmarking methods requires careful consideration of the experimental question.

Wild versus lab house mice: Effects of age, diet, and genetics on molar geometry and topography.

Molar morphology is shaped by phylogenetic history and adaptive processes related to food processing. Topographic parameters of the occlusal surface, such as sharpness and relief, can be especially informative regarding diet preferences of a species. The occlusal surface can however be deeply modified by wear throughout an animal's life, potentially obliterating other signals. Age being difficult to assess in wild populations, especially small rodents, experimental studies of wear through age in laboratory populations may constitute a powerful way to assess its impact on molar geometry and topography, and to validate descriptors of molar morphology that could mitigate this issue. Molar morphology was therefore quantified using 3D geometric morphometrics and topographic estimates in four groups of house mice: wild-trapped mice, lab-bred offspring of these wild mice, typical laboratory mice, and their hybrids. Three descriptors of the molar morphology were considered: the surface of the whole molar row, the surface of the first upper molar, and a truncated template of the first upper molar mimicking advanced wear. Increasing wear with age was demonstrated in the different groups, with a more pronounced effect in the wild-trapped population. The geometry of the molar row is not only modified by wear, but also by the relative position of the late developing molars on the jaw due to loading during mastication. As a consequence, the alignment of the molars is modified in wild mice, showing a qualitative difference between wild animals and their lab-bred offspring. Results obtained from the lab should thus be transferred with caution to the interpretation of differences in wild populations. Topographic estimates computed for the first upper molar seems to provide more stable parameters than those based on the whole molar row, because issues related to non-planar occlusal surface along the molar row are discarded. The truncated template was proven efficient in discarding the wear effect to focus on genetic differences, allowing an efficient characterization of the hybridization signature between wild and lab mice. Dominance of the wild phenotype for the first molar shape supports that the lab strain evolved in a context of relaxation of the selective pressures related to nutrition.

Natural copy number variation of tandemly repeated regulatory SNORD RNAs leads to individual phenotypic differences in mice.

Genic copy number differences can have phenotypic consequences, but so far this has not been studied in detail in natural populations. Here, we analysed the natural variation of two families of tandemly repeated regulatory small nucleolar RNAs (SNORD115 and SNORD116) in the house mouse (Mus musculus). They are encoded within the Prader-Willi Syndrome gene region, known to be involved in behavioural, metabolic, and osteogenic functions in mammals. We determined that the copy numbers of these SNORD RNAs show substantial natural variation, both in wild-derived mice as well as in an inbred mouse strain (C57BL/6J). We show that copy number differences are subject to change across generations, making them highly variable and resulting in individual differences. In transcriptome data from brain samples, we found SNORD copy-number correlated regulation of possible target genes, including Htr2c, a predicted target gene of SNORD115, as well as Ankrd11, a predicted target gene of SNORD116. Ankrd11 is a chromatin regulator, which has previously been implicated in regulating the development of the skull. Based on morphometric shape analysis of the skulls of individual mice of the inbred strain, we show that shape measures correlate with SNORD116 copy numbers in the respective individuals. Our results suggest that the variable dosage of regulatory RNAs can lead to phenotypic variation between individuals that would typically have been ascribed to environmentally induced variation, while it is actually encoded in individual differences of copy numbers of regulatory molecules.

Comparative analysis of squamate brains unveils multi-level variation in cerebellar architecture associated with locomotor specialization.

Ecomorphological studies evaluating the impact of environmental and biological factors on the brain have so far focused on morphology or size measurements, and the ecological relevance of potential multi-level variations in brain architecture remains unclear in vertebrates. Here, we exploit the extraordinary ecomorphological diversity of squamates to assess brain phenotypic diversification with respect to locomotor specialization, by integrating single-cell distribution and transcriptomic data along with geometric morphometric, phylogenetic, and volumetric analysis of high-definition 3D models. We reveal significant changes in cerebellar shape and size as well as alternative spatial layouts of cortical neurons and dynamic gene expression that all correlate with locomotor behaviours. These findings show that locomotor mode is a strong predictor of cerebellar structure and pattern, suggesting that major behavioural transitions in squamates are evolutionarily correlated with mosaic brain changes. Furthermore, our study amplifies the concept of 'cerebrotype', initially proposed for vertebrate brain proportions, towards additional shape characters.

Bracketing phenogenotypic limits of mammalian hybridization.

An increasing number of mammalian species have been shown to have a history of hybridization and introgression based on genetic analyses. Only relatively few fossils, however, preserve genetic material, and morphology must be used to identify the species and determine whether morphologically intermediate fossils could represent hybrids. Because dental and cranial fossils are typically the key body parts studied in mammalian palaeontology, here we bracket the potential for phenotypically extreme hybridizations by examining uniquely preserved cranio-dental material of a captive hybrid between grey and ringed seals. We analysed how distinct these species are genetically and morphologically, how easy it is to identify the hybrids using morphology and whether comparable hybridizations happen in the wild. We show that the genetic distance between these species is more than twice the modern human-Neanderthal distance, but still within that of morphologically similar species pairs known to hybridize. By contrast, morphological and developmental analyses show grey and ringed seals to be highly disparate, and that the hybrid is a predictable intermediate. Genetic analyses of the parent populations reveal introgression in the wild, suggesting that grey-ringed seal hybridization is not limited to captivity. Taken together, we postulate that there is considerable potential for mammalian hybridization between phenotypically disparate taxa.

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.

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

This paper provides a step-by-step guide for the morphological analysis of corolla and the decomposition of corolla shape variation into its symmetric and asymmetric components. The shape and symmetric organisation of corolla are key traits in the developmental and evolutionary biology of flowering plants. The various spatial layout of petals can exhibit bilateral symmetry, rotational symmetry or more complex combination of symmetry types. Here, I describe a general landmark-based geometric morphometric framework for the full statistical shape analysis of corolla and exemplify its use with four fully worked out case studies including tissue treatment, imaging, landmark data collection, file formatting, and statistical analyses: (i) bilateral symmetry (Fedia graciliflora), (ii) two perpendicular axes of bilateral symmetry (Erysimum mediohispanicum), (iii) rotational symmetry (Vinca minor), and (iv) combined bilateral and rotational symmetry (Trillium undulatum). The necessary tools for such analyses are not implemented in standard morphometric software and they are therefore provided here as functions running in the R environment. Principal Component Analysis is used to separate symmetric and asymmetric components of variation, respectively, quantifying variation among and within individuals. For bilaterally symmetric flowers, only one component of left-right asymmetric variation is extracted, while flowers with more complex symmetric layout have components of asymmetric variation associated with each symmetry operator implied (e.g., left-right asymmetry and adaxial-abaxial asymmetry). Fundamental information on the genetic, developmental, and environmental determinants of shape variation can be inferred from this decomposition (e.g., directional asymmetry, fluctuating asymmetry) and further exploited to document patterns of canalization, developmental stability, developmental modularity and morphological integration. Even if symmetry and asymmetry are not the primary interest of a study on corolla shape variation, statistical and anatomical arguments support the use of the framework advocated. This didactic protocol will help both morphometricians and non-morphometricians to further understand the role of symmetry in the development, variation and adaptive evolution of flowers.

Differential tissue growth and cell adhesion alone drive early tooth morphogenesis: An ex vivo and in silico study.

From gastrulation to late organogenesis animal development involves many genetic and bio-mechanical interactions between epithelial and mesenchymal tissues. Ectodermal organs, such as hairs, feathers and teeth are well studied examples of organs whose development is based on epithelial-mesenchymal interactions. These develop from a similar primordium through an epithelial folding and its interaction with the mesenchyme. Despite extensive knowledge on the molecular pathways involved, little is known about the role of bio-mechanical processes in the morphogenesis of these organs. We propose a simple computational model for the biomechanics of one such organ, the tooth, and contrast its predictions against cell-tracking experiments, mechanical relaxation experiments and the observed tooth shape changes over developmental time. We found that two biomechanical processes, differential tissue growth and differential cell adhesion, were enough, in the model, for the development of the 3D morphology of the early tooth germ. This was largely determined by the length and direction of growth of the cervical loops, lateral folds of the enamel epithelium. The formation of these cervical loops was found to require accelerated epithelial growth relative to other tissues and their direction of growth depended on specific differential adhesion between the three tooth tissues. These two processes and geometrical constraints in early tooth bud also explained the shape asymmetry between the lateral cervical loops and those forming in the anterior and posterior of the tooth. By performing mechanical perturbations ex vivo and in silico we inferred the distribution and direction of tensile stresses in the mesenchyme that restricted cervical loop lateral growth and forced them to grow downwards. Overall our study suggests detailed quantitative explanations for how bio-mechanical processes lead to specific morphological 3D changes over developmental time.

The ecological origins of snakes as revealed by skull evolution.

The ecological origin of snakes remains amongst the most controversial topics in evolution, with three competing hypotheses: fossorial; marine; or terrestrial. Here we use a geometric morphometric approach integrating ecological, phylogenetic, paleontological, and developmental data for building models of skull shape and size evolution and developmental rate changes in squamates. Our large-scale data reveal that whereas the most recent common ancestor of crown snakes had a small skull with a shape undeniably adapted for fossoriality, all snakes plus their sister group derive from a surface-terrestrial form with non-fossorial behavior, thus redirecting the debate toward an underexplored evolutionary scenario. Our comprehensive heterochrony analyses further indicate that snakes later evolved novel craniofacial specializations through global acceleration of skull development. These results highlight the importance of the interplay between natural selection and developmental processes in snake origin and diversification, leading first to invasion of a new habitat and then to subsequent ecological radiations.

Geometric morphometrics reveals shifts in flower shape symmetry and size following gene knockdown of CYCLOIDEA and ANTHOCYANIDIN SYNTHASE.

BACKGROUND: While floral symmetry has traditionally been assessed qualitatively, recent advances in geometric morphometrics have opened up new avenues to specifically quantify flower shape and size using robust multivariate statistical methods. In this study, we examine, for the first time, the ability of geometric morphometrics to detect morphological differences in floral dorsoventral asymmetry following virus-induced gene silencing (VIGS). Using Fedia graciliflora Fisch. & Meyer (Valerianaceae) as a model, corolla shape of untreated flowers was compared using canonical variate analysis to knockdown phenotypes of CYCLOIDEA2A (FgCYC2A), ANTHOCYANIDIN SYNTHASE (FgANS), and empty vector controls. RESULTS: Untreated flowers and all VIGS treatments were morphologically distinct from each other, suggesting that VIGS may cause subtle shifts in floral shape. Knockdowns of FgCYC2A were the most dramatic, affecting the position of dorsal petals in relation to lateral petals, thereby resulting in more actinomorphic-like flowers. Additionally, FgANS knockdowns developed larger flowers with wider corolla tube openings. CONCLUSIONS: These results provide a method to quantify the role that specific genes play in the developmental pathway affecting the dorsoventral axis of symmetry in zygomorphic flowers. Additionally, they suggest that ANS may have an unintended effect on floral size and shape.

Development and evolution of segmentation assessed by geometric morphometrics: The centipede Strigamia maritima as a case study.

Using the centipede model species Strigamia maritima as a subject of study, we illustrate the potential of geometric morphometrics for investigating the development and evolution of segmentation, with a specific focus on post-embryonic segmental patterning. We show how these techniques can contribute detailed descriptive data for comparative purposes, but also precious information on some features of the developmental system that are considered relevant for the evolvability of a segmented body architecture, such as developmental stability and canalization. Morphometric analyses allow to separately investigate several sources of phenotypic variation along a segmented body axis, like constitutive and random segment heteronomy, both within and among individuals. Specifically, in S. maritima, the segmental pattern of ventral sclerite shapes mirrors that of their bilateral fluctuating asymmetry and among-individual variation in associating the most anterior and most posterior segments in diverging from the central ones. Also, among segments, there seems to be a correlation between fluctuating asymmetry and shape variation among individuals, suggesting that canalization and developmental stability are somehow associated. Overall, these associations might stem from a joint influence of the segmental position on the two processes of developmental buffering.

Modularity and developmental stability in segmented animals: variation in translational asymmetry in geophilomorph centipedes.

Does a modular body organization present a challenge for developmental control? We investigate the idea of a possible developmental cost of modularity by examining the relationship between modularity and developmental stability in a multi-segmented arthropod taxon: the geophilomorph centipedes. In a sample of eight species, we tested the correlation between developmental stability, estimated from measures of translational fluctuating asymmetry, and the number of trunk segments and some other morphological traits, both at the species and individual levels. We found sizeable differences in size and shape patterns of variation at the level of species. However, we did not find any clear evidence of correlation between fluctuating asymmetry and the number of trunk segments or the other morphological traits considered. Thus, our results provide no support to the idea of a possible trade-off between the cardinality of a modular system and the level of developmental precision in the phenotypic expression of its modules. The results of this exploratory study invite further investigations of patterns of translational fluctuating asymmetry in segmented animals and other modular organisms, as these have the potential to reveal features of developmental stability that cannot be captured by the study of bilateral asymmetry alone.

A new method to pair-match metacarpals using bilateral asymmetry and shape analysis.

Matching pairs of skeletal elements is of fundamental importance when re-associating skeletons from commingled settings. This study presents a new method that combines the use of traditional anthropometric size measurements and the analysis of shape, through geometric morphometrics, to form pairs from commingled settings. We tested the hypothesis that shape asymmetry among individuals is smaller than between individuals that share similar dimensions. The metacarpals were measured in their maximum length, and absolute means of asymmetry between right and left metacarpals were calculated. To analyze the shape characteristics, two-dimensional landmarks were recorded from digital images of 111 metacarpals. The specimens were analyzed using generalized Procrustes analysis and multivariate statistics. The accuracy of pair-matching bones using a combined method was of 100%, showing that the incorporation of geometric morphometrics techniques--in the assessment of shape in physical anthropology--is of particular use when comparing the shape differences/similarities between individuals.

Geometric morphometrics of corolla shape: dissecting components of symmetric and asymmetric variation in Erysimum mediohispanicum (Brassicaceae).

Symmetry is an important feature of floral structure, and floral symmetries are diverse and often complex. We use a new morphometric approach for analysing shapes with complex types of symmetry, which partitions shape variation into a component of symmetric variation and different components of asymmetry. This approach, based on the mathematical theory of symmetry groups, can be used for landmark configurations with any type of symmetry and is therefore promising as a general framework for morphometric analyses of floral symmetry and asymmetry. We demonstrate this approach by quantifying floral shape variation in a wild population of Erysimum mediohispanicum (Brassicaceae). Flowers of this species are disymmetric, so that the symmetry in the left-right and adaxial-abaxial directions can be considered separately and in combination. Both principal component analysis and Procrustes ANOVA indicate that symmetric variation accounts for most of the total variance and that adaxial-abaxial asymmetry is the dominant component of fluctuating asymmetry. Each component is associated with specific patterns of shape variation. These results illustrate the potential of the new method and suggest new areas for future research. The new morphometric approach is promising for further analyses of floral symmetry and asymmetry in evolutionary and developmental contexts.

Beyond bilateral symmetry: geometric morphometric methods for any type of symmetry.

BACKGROUND: Studies of symmetric structures have made important contributions to evolutionary biology, for example, by using fluctuating asymmetry as a measure of developmental instability or for investigating the mechanisms of morphological integration. Most analyses of symmetry and asymmetry have focused on organisms or parts with bilateral symmetry. This is not the only type of symmetry in biological shapes, however, because a multitude of other types of symmetry exists in plants and animals. For instance, some organisms have two axes of reflection symmetry (biradial symmetry; e.g. many algae, corals and flowers) or rotational symmetry (e.g. sea urchins and many flowers). So far, there is no general method for the shape analysis of these types of symmetry. RESULTS: We generalize the morphometric methods currently used for the shape analysis of bilaterally symmetric objects so that they can be used for analyzing any type of symmetry. Our framework uses a mathematical definition of symmetry based on the theory of symmetry groups. This approach can be used to divide shape variation into a component of symmetric variation among individuals and one or more components of asymmetry. We illustrate this approach with data from a colonial coral that has ambiguous symmetry and thus can be analyzed in multiple ways. Our results demonstrate that asymmetric variation predominates in this dataset and that its amount depends on the type of symmetry considered in the analysis. CONCLUSIONS: The framework for analyzing symmetry and asymmetry is suitable for studying structures with any type of symmetry in two or three dimensions. Studies of complex symmetries are promising for many contexts in evolutionary biology, such as fluctuating asymmetry, because these structures can potentially provide more information than structures with bilateral symmetry.