Developmental constraint through negative pleiotropy in the zygomatic arch

Variable Craniofacial Shape and Development among Multiple Cave-Adapted Populations of Astyanax mexicanus

Astyanax mexicanus is a freshwater fish species with blind cave morphs and sighted surface morphs. Like other troglodytic species, independently evolved cave-dwelling A. mexicanus populations share several stereotypic phenotypes, including the expansion of certain sensory systems, as well as the loss of eyes and pigmentation. Here, we assess the extent to which there is also parallelism in craniofacial development across cave populations. Since multiple forces may be acting upon variation in the A. mexicanus system, including phylogenetic history, selection, and developmental constraint, several outcomes are possible. For example, eye regression may have triggered a conserved series of compensatory developmental events, in which case we would expect to observe highly similar craniofacial phenotypes across cave populations. Selection for cave-specific foraging may also lead to the evolution of a conserved craniofacial phenotype, especially in regions of the head directly associated with feeding. Alternatively, in the absence of a common axis of selection or strong developmental constraints, craniofacial shape may evolve under neutral processes such as gene flow, drift, and bottlenecking, in which case patterns of variation should reflect the evolutionary history of A. mexicanus. Our results found that cave-adapted populations do share certain anatomical features; however, they generally did not support the hypothesis of a conserved craniofacial phenotype across caves, as nearly every pairwise comparison was statistically significant, with greater effect sizes noted between more distantly related cave populations with little gene flow. A similar pattern was observed for developmental trajectories. We also found that morphological disparity was lower among all three cave populations versus surface fish, suggesting eye loss is not associated with increased variation, which would be consistent with a release of developmental constraint. Instead, this pattern reflects the relatively low genetic diversity within cave populations. Finally, magnitudes of craniofacial integration were found to be similar among all groups, meaning that coordinated development among anatomical units is robust to eye loss in A. mexicanus. We conclude that, in contrast to many conserved phenotypes across cave populations, global craniofacial shape is more variable, and patterns of shape variation are more in line with population structure than developmental architecture or selection.

Flexible conservatism in the skull modularity of convergently evolved myrmecophagous placental mammals

Background

The skull of placental mammals constitutes one of the best studied systems for phenotypic modularity. Several studies have found strong evidence for the conserved presence of two- and six-module architectures, while the strength of trait correlations (integration) has been associated with major developmental processes such as somatic growth, muscle-bone interactions, and tooth eruption. Among placentals, ant- and termite-eating (myrmecophagy) represents an exemplar case of dietary convergence, accompanied by the selection of several cranial morphofunctional traits such as rostrum elongation, tooth loss, and mastication loss. Despite such drastic functional modifications, the covariance patterns of the skull of convergently evolved myrmecophagous placentals are yet to be studied in order to assess the potential consequences of this dietary shift on cranial modularity.

Results

Here, we performed a landmark-based morphometric analysis of cranial covariance patterns in 13 species of myrmecophagous placentals. Our analyses reveal that most myrmecophagous species present skulls divided into six to seven modules (depending on the confirmatory method used), with architectures similar to those of non-myrmecophagous placentals (therian six modules). Within-module integration is also similar to what was previously described for other placentals, suggesting that most covariance-generating processes are conserved across the clade. Nevertheless, we show that extreme rostrum elongation and tooth loss in myrmecophagid anteaters have resulted in a shift in intermodule correlations in the proximal region of the rostrum. Namely, the naso-frontal and maxillo-palatine regions are strongly correlated with the oro-nasal module, suggesting an integrated rostrum conserved from pre-natal developmental processes. In contrast, the similarly toothless pangolins show a weaker correlation between the anterior rostral modules, resembling the pattern of toothed placentals.

Conclusions

These results reveal that despite some integration shifts related to extreme functional and morphological features of myrmecophagous skulls, cranial modular architectures have conserved the typical mammalian scheme.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02030-9.

The genetic basis of neurocranial size and shape across varied lab mouse populations

Brain and skull tissues interact through molecular signalling and mechanical forces during head development, leading to a strong correlation between the neurocranium and the external brain surface. Therefore, when brain tissue is unavailable, neurocranial endocasts are often used to approximate brain size and shape. Evolutionary changes in brain morphology may have resulted in secondary changes to neurocranial morphology, but the developmental and genetic processes underlying this relationship are not well understood. Using automated phenotyping methods, we quantified the genetic basis of endocast variation across large genetically varied populations of laboratory mice in two ways: (1) to determine the contributions of various genetic factors to neurocranial form and (2) to help clarify whether a neurocranial variation is based on genetic variation that primarily impacts bone development or on genetic variation that primarily impacts brain development, leading to secondary changes in bone morphology. Our results indicate that endocast size is highly heritable and is primarily determined by additive genetic factors. In addition, a non-additive inbreeding effect led to founder strains with lower neurocranial size, but relatively large brains compared to skull size; suggesting stronger canalization of brain size and/or a general allometric effect. Within an outbred sample of mice, we identified a locus on mouse chromosome 1 that is significantly associated with variation in several positively correlated endocast size measures. Because the protein-coding genes at this locus have been previously associated with brain development and not with bone development, we propose that genetic variation at this locus leads primarily to variation in brain volume that secondarily leads to changes in neurocranial globularity. We identify a strain-specific missense mutation within Akt3 that is a strong causal candidate for this genetic effect. Whilst it is not appropriate to generalize our hypothesis for this single locus to all other loci that also contribute to the complex trait of neurocranial skull morphology, our results further reveal the genetic basis of neurocranial variation and highlight the importance of the mechanical influence of brain growth in determining skull morphology.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Integration and the Developmental Genetics of Allometry

Allometry refers to the ways in which organismal shape is associated with size. It is a special case of integration, or the tendency for traits to covary, in that variation in size is ubiquitous and evolutionarily important. Allometric variation is so commonly observed that it is routinely removed from morphometric analyses or invoked as an explanation for evolutionary change. In this case, familiarity is mistaken for understanding because rarely do we know the mechanisms by which shape correlates with size or understand their significance. As with other forms of integration, allometric variation is generated by variation in developmental processes that affect multiple traits, resulting in patterns of covariation. Given this perspective, we can dissect the genetic and developmental determinants of allometric variation. Our work on the developmental and genetic basis for allometric variation in craniofacial shape in mice and humans has revealed that allometric variation is highly polygenic. Different measures of size are associated with distinct but overlapping patterns of allometric variation. These patterns converge in part on a common genetic basis. Finally, environmental modulation of size often generates variation along allometric trajectories, but the timing of genetic and environmental perturbations can produce deviations from allometric patterns when traits are differentially sensitive over developmental time. These results question the validity of viewing allometry as a singular phenomenon distinct from morphological integration more generally.

The effect of automated landmark identification on morphometric analyses

Morphometric analysis of anatomical landmarks allows researchers to identify specific morphological differences between natural populations or experimental groups, but manually identifying landmarks is time-consuming. We compare manually and automatically generated adult mouse skull landmarks and subsequent morphometric analyses to elucidate how switching from manual to automated landmarking will impact morphometric analysis results for large mouse (Mus musculus) samples (n = 1205) that represent a wide range of 'normal' phenotypic variation (62 genotypes). Other studies have suggested that the use of automated landmarking methods is feasible, but this study is the first to compare the utility of current automated approaches to manual landmarking for a large dataset that allows the quantification of intra- and inter-strain variation. With this unique sample, we investigated how switching to a non-linear image registration-based automated landmarking method impacts estimated differences in genotype mean shape and shape variance-covariance structure. In addition, we tested whether an initial registration of specimen images to genotype-specific averages improves automatic landmark identification accuracy. Our results indicated that automated landmark placement was significantly different than manual landmark placement but that estimated skull shape covariation was correlated across methods. The addition of a preliminary genotype-specific registration step as part of a two-level procedure did not substantially improve on the accuracy of one-level automatic landmark placement. The landmarks with the lowest automatic landmark accuracy are found in locations with poor image registration alignment. The most serious outliers within morphometric analysis of automated landmarks displayed instances of stochastic image registration error that are likely representative of errors common when applying image registration methods to micro-computed tomography datasets that were initially collected with manual landmarking in mind. Additional efforts during specimen preparation and image acquisition can help reduce the number of registration errors and improve registration results. A reduction in skull shape variance estimates were noted for automated landmarking methods compared with manual landmarking. This partially reflects an underestimation of more extreme genotype shapes and loss of biological signal, but largely represents the fact that automated methods do not suffer from intra-observer landmarking error. For appropriate samples and research questions, our image registration-based automated landmarking method can eliminate the time required for manual landmarking and have a similar power to identify shape differences between inbred mouse genotypes.