Geometric Morphometrics on Gene Expression Patterns Within Phenotypes: A Case Example on Limb Development

High-Throughput 3D Phenotyping of Plant Shoot Apical Meristems From Tissue-Resolution Data

Confocal imaging is a well-established method for investigating plant phenotypes on the tissue and organ level. However, many differences are difficult to assess by visual inspection and researchers rely extensively on ad hoc manual quantification techniques and qualitative assessment. Here we present a method for quantitatively phenotyping large samples of plant tissue morphologies using triangulated isosurfaces. We successfully demonstrate the applicability of the approach using confocal imaging of aerial organs in Arabidopsis thaliana. Automatic identification of flower primordia using the surface curvature as an indication of outgrowth allows for high-throughput quantification of divergence angles and further analysis of individual flowers. We demonstrate the throughput of our method by quantifying geometric features of 1065 flower primordia from 172 plants, comparing auxin transport mutants to wild type. Additionally, we find that a paraboloid provides a simple geometric parameterisation of the shoot inflorescence domain with few parameters. We utilise parameterisation methods to provide a computational comparison of the shoot apex defined by a fluorescent reporter of the central zone marker gene CLAVATA3 with the apex defined by the paraboloid. Finally, we analyse the impact of mutations which alter mechanical properties on inflorescence dome curvature and compare the results with auxin transport mutants. Our results suggest that region-specific expression domains of genes regulating cell wall biosynthesis and local auxin transport can be important in maintaining the wildtype tissue shape. Altogether, our results indicate a general approach to parameterise and quantify plant development in 3D, which is applicable also in cases where data resolution is limited, and cell segmentation not possible. This enables researchers to address fundamental questions of plant development by quantitative phenotyping with high throughput, consistency and reproducibility.

Regulatory encoding of quantitative variation in spatial activity of a Drosophila enhancer

Developmental enhancers control the expression of genes prefiguring morphological patterns. The activity of an enhancer varies among cells of a tissue, but collectively, expression levels in individual cells constitute a spatial pattern of gene expression. How the spatial and quantitative regulatory information is encoded in an enhancer sequence is elusive. To link spatial pattern and activity levels of an enhancer, we used systematic mutations of the yellow spot enhancer, active in developing Drosophila wings, and tested their effect in a reporter assay. Moreover, we developed an analytic framework based on the comprehensive quantification of spatial reporter activity. We show that the quantitative enhancer activity results from densely packed regulatory information along the sequence, and that a complex interplay between activators and multiple tiers of repressors carves the spatial pattern. Our results shed light on how an enhancer reads and integrates trans-regulatory landscape information to encode a spatial quantitative pattern.

Species specificity and sexual dimorphism in tooth shape among the three sympatric haplochromine species in Lake Kivu cichlids

Tooth shape is used to differentiate between morphologically similar species of vertebrates, including fish. This study aimed to quantify tooth shape of three sympatric species: Haplochromis kamiranzovu, H. insidiae, and H. astatodon endemic to Lake Kivu, whose existing identification criteria are currently only qualitative. A quantitative tooth shape analysis was performed based on digitized tooth outline data with a subsequent elliptic Fourier analysis to test for differences among the three species. We looked at crown shape and size differences within H. kamiranzovu and H. insidiae at geographical, habitat, and gender levels. No comparison at habitat level was done for H. astatodon because it is found only in littoral zone. The analysis revealed significant tooth shape differences among the three species. Haplochromis astatodon had a significantly longer major cusp height and a longer and larger minor cusp than that of H. insidiae. It had also a longer major cusp height and a longer and larger minor cusp than that of H. kamiranzovu. Tooth shape differences of H. kamiranzovu and H. insidiae species were not significantly different between littoral and pelagic fish (p > .05) while differences were significant between southern and northern Lake Kivu populations (p < .05). Tooth sizes in H. kamiranzovu and H. insidiae were significantly different, both in height and width as well as in their ratios, and this was true at sex and geographic levels (p < .05), but not at habitat level (p > .05). Tooth shape was also significantly different with sharp teeth for males compared with females of southern populations versus northern ones. These shape- and size-related differences between sexes suggest differences in the foraging strategies toward available food resources in the lake habitat. Further research should explain the genetic basis of the observed pattern.

Quantification of gene expression patterns to reveal the origins of abnormal morphogenesis

The earliest developmental origins of dysmorphologies are poorly understood in many congenital diseases. They often remain elusive because the first signs of genetic misregulation may initiate as subtle changes in gene expression, which are hard to detect and can be obscured later in development by secondary effects. Here, we develop a method to trace back the origins of phenotypic abnormalities by accurately quantifying the 3D spatial distribution of gene expression domains in developing organs. By applying Geometric Morphometrics to 3D gene expression data obtained by Optical Projection Tomography, we determined that our approach is sensitive enough to find regulatory abnormalities that have never been detected previously. We identified subtle but significant differences in the gene expression of a downstream target of a Fgfr2 mutation associated with Apert syndrome, demonstrating that these mouse models can further our understanding of limb defects in the human condition. Our method can be applied to different organ systems and models to investigate the etiology of malformations.

Our development in the womb is complex. Genes need to switch on and off in a precise order, controlling the activity of millions of cells as they work together to form different tissues. For everything to happen smoothly, cells must use instructions provided by each gene exactly at the correct moment and in the correct place. In this biological assembly line, the slightest change can lead to a defect.

Certain genetic mutations can change when and where cells use particular genes, and this can cause errors in development. These kinds of mutations are a common cause of birth defects, but we cannot always pinpoint how they begin. For example, a single mutation in a gene called FGFR2 causes malformations in the head, the heart and the limbs in a rare disease called Apert syndrome.

The first signs that development has gone wrong can be subtle changes in the use of certain genes, impossible to detect with standard methods. As development continues, other processes can mask the impact of problems with certain genes. Ultimately, changes alter the shape of the developing embryo.

Genetically engineered mouse models can mimic the gene defects that cause disease in humans. But current methods are not sensitive enough to detect the very first signs of defects. Now, Martínez-Abadías et al. developed a new method to detect these subtle changes and reveal the precise moment when development starts to go wrong.

In mice, a specific mutation in the FGFR2 gene affects the activity of a series of other genes. To track the levels of one of these genes, Martínez-Abadías et al. marked mouse embryos using a chemical label. Scanning the embryos then revealed the pattern of the cells using the gene during the earliest stages of development. In mice carrying a mutation in the FGFR2 gene, subtle changes in gene expression began just a few hours after their limbs start to develop. But it took another half a day to see the effects of these changes on the shape and size of the growing limbs. This approach revealed changes in gene expression before any problems with development were visible by eye.

Tracking subtle changes in the way cells use genes could allow us to detect the origins of embryo malformations before they appear, pointing at the best moment to start a treatment. With further development, the model could extend to other genes, proteins, animal models and diseases.

Quantitative Morphological Variation in the Developing Drosophila Wing

Quantitative genetic variation in morphology is pervasive in all species and is the basis for the evolution of differences among species. The measurement of morphological form in adults is now beginning to be combined with comparable measurements of form during development. Here we compare the shape of the developing wing to its adult form in a holometabolous insect, Drosophila melanogaster. We used protein expression patterns to measure shape in the developing precursors of the final adult wing. Three developmental stages were studied: late larval third instar, post-pupariation and in the adult fly. We studied wild-type animals in addition to mutants of two genes (shf and ds) that have known effects on adult wing shape and size. Despite experimental noise related to the difficulty of comparing developing structures, we found consistent differences in wing shape and size at each developmental stage between genotypes. Quantitative comparisons of variation arising at different developmental stages with the variation in the final structure enable us to determine when variation arises, and to generate hypotheses about the causes of that variation. In addition we provide linear rules allowing us to link wing morphology in the larva, with wing morphology in the pupa. Our approach provides a framework to analyze quantitative morphological variation in the developing fly wing. This framework should help to characterize the natural variation of the larval and pupal wing shape, and to measure the contribution of the processes occurring during these developmental stages to the natural variation in adult wing morphology.

The developmental-genetics of canalization

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Elliptical Fourier analysis: fundamentals, applications, and value for forensic anthropology

The numerical description of skeletal morphology enables forensic anthropologists to conduct objective, reproducible, and structured tests, with the added capability of verifying morphoscopic-based analyses. One technique that permits comprehensive quantification of outline shape is elliptical Fourier analysis. This curve fitting technique allows a form's outline to be approximated via the sum of multiple sine and cosine waves, permitting the profile perimeter of an object to be described in a dense (continuous) manner at a user-defined level of precision. A large amount of shape information (the entire perimeter) can thereby be collected in contrast to other methods relying on sparsely located landmarks where information falling in between the landmarks fails to be acquired. First published in 1982, elliptical Fourier analysis employment in forensic anthropology from 2000 onwards reflects a slow uptake despite large computing power that makes its calculations easy to conduct. Without hurdles arising from calculation speed or quantity, the slow uptake may partly reside with the underlying mathematics that on first glance is extensive and potentially intimidating. In this paper, we aim to bridge this gap by pictorially illustrating how elliptical Fourier harmonics work in a simple step-by-step visual fashion to facilitate universal understanding and as geared towards increased use in forensic anthropology. We additionally provide a short review of the method's utility for osteology, a summary of past uses in forensic anthropology, and software options for calculations that largely save the user the trouble of coding customized routines.