Developmental nonlinearity drives phenotypic robustness

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

Deciphering the origin of developmental stability: The role of intracellular expression variability in evolutionary conservation

Progress in evolutionary developmental biology (evo-devo) has deepened our understanding of how intrinsic properties of embryogenesis, along with natural selection and population genetics, shape phenotypic diversity. A focal point of recent empirical and theoretical research is the idea that highly developmentally stable phenotypes are more conserved in evolution. Previously, we demonstrated that in Japanese medaka (Oryzias latipes), embryonic stages and genes with high stability, estimated through whole-embryo RNA-seq, are highly conserved in subsequent generations. However, the precise origin of the stability of gene expression levels evaluated at the whole-embryo level remained unclear. Such stability could be attributed to two distinct sources: stable intracellular expression levels or spatially stable expression patterns. Here we demonstrate that stability observed in whole-embryo RNA-seq can be attributed to stability at the cellular level (low variability in gene expression at the cellular levels). We quantified the intercellular variations in expression levels and spatial gene expression patterns for seven key genes involved in patterning dorsoventral and rostrocaudal regions during early development in medaka. We evaluated intracellular variability by counting transcripts and found its significant correlation with variation observed in whole-embryo RNA-seq data. Conversely, variation in spatial gene expression patterns, assessed through intraindividual left-right asymmetry, showed no correlation. Given the previously reported correlation between stability and conservation of expression levels throughout embryogenesis, our findings suggest a potential general trend: the stability or instability of developmental systems-and the consequent evolutionary diversity-may be primarily anchored in intrinsic fundamental elements such as the variability of intracellular states.

Evo-Devo Mechanobiology: The Missing Link

While the modern framework of evolutionary development (evo-devo) has been decidedly genetic, historic analyses have also considered the importance of mechanics in the evolution of form. With the aid of recent technological advancements in both quantifying and perturbing changes in the molecular and mechanical effectors of organismal shape, how molecular and genetic cues regulate the biophysical aspects of morphogenesis is becoming increasingly well studied. As a result, this is an opportune time to consider how the tissue-scale mechanics that underlie morphogenesis are acted upon through evolution to establish morphological diversity. Such a focus will enable a field of evo-devo mechanobiology that will serve to better elucidate the opaque relations between genes and forms by articulating intermediary physical mechanisms. Here, we review how the evolution of shape is measured and related to genetics, how recent strides have been made in the dissection of developmental tissue mechanics, and how we expect these areas to coalesce in evo-devo studies in the future.

Developing the genotype-to-phenotype relationship in evolutionary theory: A primer of developmental features

For decades, there have been repeated calls for more integration across evolutionary and developmental biology. However, critiques in the literature and recent funding initiatives suggest this integration remains incomplete. We suggest one way forward is to consider how we elaborate the most basic concept of development, the relationship between genotype and phenotype, in traditional models of evolutionary processes. For some questions, when more complex features of development are accounted for, predictions of evolutionary processes shift. We present a primer on concepts of development to clarify confusion in the literature and fuel new questions and approaches. The basic features of development involve expanding a base model of genotype-to-phenotype to include the genome, space, and time. A layer of complexity is added by incorporating developmental systems, including signal-response systems and networks of interactions. The developmental emergence of function, which captures developmental feedbacks and phenotypic performance, offers further model elaborations that explicitly link fitness with developmental systems. Finally, developmental features such as plasticity and developmental niche construction conceptualize the link between a developing phenotype and the external environment, allowing for a fuller inclusion of ecology in evolutionary models. Incorporating aspects of developmental complexity into evolutionary models also accommodates a more pluralistic focus on the causal importance of developmental systems, individual organisms, or agents in generating evolutionary patterns. Thus, by laying out existing concepts of development, and considering how they are used across different fields, we can gain clarity in existing debates around the extended evolutionary synthesis and pursue new directions in evolutionary developmental biology. Finally, we consider how nesting developmental features in traditional models of evolution can highlight areas of evolutionary biology that need more theoretical attention.

Toward a universal measure of robustness across model organs and systems

The development of an individual must be capable of resisting the harmful effects of internal and external perturbations. This capacity, called robustness, can make the difference between normal variation and disease. Some systems and organs are more resilient in their capacity to correct the effects of internal disturbances such as mutations. Similarly, organs and organisms differ in their capacity to be resilient against external disturbances, such as changes in temperature. Furthermore, all developmental systems must be somewhat flexible to permit evolutionary change, and understanding robustness requires a comparative framework. Over the last decades, most research on developmental robustness has been focusing on specific model systems and organs. Hence, we lack tools that would allow cross-species and cross-organ comparisons. Here, we emphasize the need for a uniform framework to experimentally test and quantify robustness across study systems and suggest that the analysis of fluctuating asymmetry might be a powerful proxy to do so. Such a comparative framework will ultimately help to resolve why and how organs of the same and different species differ in their sensitivity to internal (e.g., mutations) and external (e.g., temperature) perturbations and at what level of biological organization buffering capacities exist and therefore create robustness of the developmental system.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

Complexities of recapitulating polygenic effects in natural populations: replication of genetic effects on wing shape in artificially selected and wild-caught populations of Drosophila melanogaster

Identifying the genetic architecture of complex traits is important to many geneticists, including those interested in human disease, plant and animal breeding, and evolutionary genetics. Advances in sequencing technology and statistical methods for genome-wide association studies have allowed for the identification of more variants with smaller effect sizes, however, many of these identified polymorphisms fail to be replicated in subsequent studies. In addition to sampling variation, this failure to replicate reflects the complexities introduced by factors including environmental variation, genetic background, and differences in allele frequencies among populations. Using Drosophila melanogaster wing shape, we ask if we can replicate allelic effects of polymorphisms first identified in a genome-wide association studies in three genes: dachsous, extra-macrochaete, and neuralized, using artificial selection in the lab, and bulk segregant mapping in natural populations. We demonstrate that multivariate wing shape changes associated with these genes are aligned with major axes of phenotypic and genetic variation in natural populations. Following seven generations of artificial selection along the dachsous shape change vector, we observe genetic differentiation of variants in dachsous and genomic regions containing other genes in the hippo signaling pathway. This suggests a shared direction of effects within a developmental network. We also performed artificial selection with the extra-macrochaete shape change vector, which is not a part of the hippo signaling network, but showed a largely shared direction of effects. The response to selection along the emc vector was similar to that of dachsous, suggesting that the available genetic diversity of a population, summarized by the genetic (co)variance matrix (G), influenced alleles captured by selection. Despite the success with artificial selection, bulk segregant analysis using natural populations did not detect these same variants, likely due to the contribution of environmental variation and low minor allele frequencies, coupled with small effect sizes of the contributing variants.

Evolutionary mechanisms modulating the mammalian skull development

Mammals possess impressive craniofacial variation that mirrors their adaptation to diverse ecological niches, feeding behaviour, physiology and overall lifestyle. The spectrum of craniofacial geometries is established mainly during embryonic development. The formation of the head represents a sequence of events regulated on genomic, molecular, cellular and tissue level, with each step taking place under tight spatio-temporal control. Even minor variations in timing, position or concentration of the molecular drivers and the resulting events can affect the final shape, size and position of the skeletal elements and the geometry of the head. Our knowledge of craniofacial development increased substantially in the last decades, mainly due to research using conventional vertebrate model organisms. However, how developmental differences in head formation arise specifically within mammals remains largely unexplored. This review highlights three evolutionary mechanisms acknowledged to modify ontogenesis: heterochrony, heterotopy and heterometry. We present recent research that links changes in developmental timing, spatial organization or gene expression levels to the acquisition of species-specific skull morphologies. We highlight how these evolutionary modifications occur on the level of the genes, molecules and cellular processes, and alter conserved developmental programmes to generate a broad spectrum of skull shapes characteristic of the class Mammalia. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Fgf8 regulates first pharyngeal arch segmentation through pouch-cleft interactions

Introduction: The pharyngeal arches are transient developmental structures that, in vertebrates, give rise to tissues of the head and neck. A critical process underlying the specification of distinct arch derivatives is segmentation of the arches along the anterior-posterior axis. Formation of ectodermal-endodermal interfaces is a key mediator of this process, and although it is essential, mechanisms regulating the establishment of these interfaces vary between pouches and between taxa. Methods: Here, we focus on the patterning and morphogenesis of epithelia associated with the first pharyngeal arch, the first pharyngeal pouch (pp1) and the first pharyngeal cleft (pc1), and the role of Fgf8 dosage in these processes in the mouse model system. Results: We find that severe reductions of Fgf8 levels disrupt both pp1 and pc1 development. Notably, out-pocketing of pp1 is largely robust to Fgf8 reductions, however, pp1 extension along the proximal-distal axis fails when Fgf8 is low. Our data indicate that Fgf8 is required for specification of regional identity in both pp1 and pc1, for localized changes in cell polarity, and for elongation and extension of both pp1 and pc1. Discussion: Based on Fgf8-mediated changes in tissue relationships between pp1 and pc1, we hypothesize that extension of pp1 requires physical interaction with pc1. Overall, our data indicate a critical role for the lateral surface ectoderm in segmentation of the first pharyngeal arch that has previously been under-appreciated.

Precise modulation of transcription factor levels identifies features underlying dosage sensitivity

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

Morphology, evolution, and the whole organism imperative: Why evolutionary questions need multi-trait evolutionary quantitative genetics

Since Washburn's New Physical Anthropology, researchers have sought to understand the complexities of morphological evolution among anatomical regions in human and non-human primates. Researchers continue, however, to preferentially use comparative and functional approaches to examine complex traits, but these methods cannot address questions about evolutionary process and often conflate function with fitness. Moreover, researchers also tend to examine anatomical elements in isolation, which implicitly assumes independent evolution among different body regions. In this paper, we argue that questions asked in primate evolution are best examined using multiple anatomical regions subjected to model-bound methods built from an understanding of evolutionary quantitative genetics. A nascent but expanding number of studies over the last two decades use this approach, examining morphological integration, evolvability, and selection modeling. To help readers learn how to use these methods, we review fundamentals of evolutionary processes within a quantitative genetic framework, explore the importance of neutral evolutionary theory, and explain the basics of evolutionary quantitative genetics, namely the calculation of evolutionary potential for multiple traits in response to selection. Leveraging these methods, we demonstrate their use to understand non-independence in possible evolutionary responses across the limbs, limb girdles, and basicranium of humans. Our results show that model-bound quantitative genetic methods can reveal unexpected genetic covariances among traits that create a novel but measurable understanding of evolutionary complexity among multiple traits. We advocate for evolutionary quantitative genetic methods to be a standard whenever appropriate to keep studies of primate morphological evolution relevant for the next seventy years and beyond.

Fgf8 regulates first pharyngeal arch segmentation through pouch-cleft interactions

The pharyngeal arches are transient developmental structures that, in vertebrates, give rise to tissues of the head and neck. A critical process underlying the specification of distinct arch derivatives is segmentation of the arches along the anterior-posterior axis. Out-pocketing of the pharyngeal endoderm between the arches is a key mediator of this process, and although it is essential, mechanisms regulating out-pocketing vary between pouches and between taxa. Here, we focus on the patterning and morphogenesis of epithelia associated with the first pharyngeal arch, the first pharyngeal pouch (pp1) and the first pharyngeal cleft (pc1), and the role of Fgf8 dosage in these processes. We find that severe reductions of Fgf8 levels disrupt both pp1 and pc1 development. Notably, out-pocketing of pp1 is largely robust to Fgf8 reductions, however, pp1 extension along the proximal-distal axis fails when Fgf8 is low. Our data indicate that extension of pp1 requires physical interaction with pc1, and that multiple aspects of pc1 morphogenesis require Fgf8 . In particular, Fgf8 is required for specification of regional identity in both pp1 and pc1, for localized changes in cell polarity, and for elongation and extension of both pp1 and pc1. Overall, our data indicate a critical role for the lateral surface ectoderm in segmentation of the first pharyngeal arch that has previously been under-appreciated.

Stability in gene expression and body-plan development leads to evolutionary conservation

BACKGROUND

Phenotypic evolution is mainly explained by selection for phenotypic variation arising from factors including mutation and environmental noise. Recent theoretical and experimental studies have suggested that phenotypes with greater developmental stability tend to have a constant phenotype and gene expression level within a particular genetic and environmental condition, and this positively correlates with stronger evolutionary conservation, even after the accumulation of genetic changes. This could reflect a novel mechanism that contributes to evolutionary conservation; however, it remains unclear whether developmental stability is the cause, or whether at least it contributes to their evolutionary conservation. Here, using Japanese medaka lines, we tested experimentally whether developmental stages and gene expression levels with greater stability led to their evolutionary conservation.

RESULTS

We first measured the stability of each gene expression level and developmental stage (defined here as the whole embryonic transcriptome) in the inbred F0 medaka population. We then measured their evolutionary conservation in the F3 generation by crossing the F0 line with the distantly related Japanese medaka line (Teradomori), followed by two rounds of intra-generational crossings. The results indicated that the genes and developmental stages that had smaller variations in the F0 generation showed lower diversity in the hybrid F3 generation, which implies a causal relationship between stability and evolutionary conservation.

CONCLUSIONS

These findings suggest that the stability in phenotypes, including the developmental stages and gene expression levels, leads to their evolutionary conservation; this most likely occurs due to their low potential to generate phenotypic variation. In addition, since the highly stable developmental stages match with the body-plan-establishment stage, it also implies that the developmental stability potentially contributed to the strict conservation of animal body plan.

Covariation of brain and skull shapes as a model to understand the role of crosstalk in development and evolution

Covariation among discrete phenotypes can arise due to selection for shared functions, and/or shared genetic and developmental underpinnings. The consequences of such phenotypic integration are far-reaching and can act to either facilitate or limit morphological variation. The vertebrate brain is known to act as an "organizer" of craniofacial development, secreting morphogens that can affect the shape of the growing neurocranium, consistent with roles for pleiotropy in brain-neurocranium covariation. Here, we test this hypothesis in cichlid fishes by first examining the degree of shape integration between the brain and the neurocranium using three-dimensional geometric morphometrics in an F₅ hybrid population, and then genetically mapping trait covariation using quantitative trait loci (QTL) analysis. We observe shape associations between the brain and the neurocranium, a pattern that holds even when we assess associations between the brain and constituent parts of the neurocranium: the rostrum and braincase. We also recover robust genetic signals for both hard- and soft-tissue traits and identify a genomic region where QTL for the brain and braincase overlap, implicating a role for pleiotropy in patterning trait covariation. Fine mapping of the overlapping genomic region identifies a candidate gene, notch1a, which is known to be involved in patterning skeletal and neural tissues during development. Taken together, these data offer a genetic hypothesis for brain-neurocranium covariation, as well as a potential mechanism by which behavioral shifts may simultaneously drive rapid change in neuroanatomy and craniofacial morphology.

Pleiotropic functions of chordin gene causing drastic morphological changes in ornamental goldfish

Breeders and fanciers have established many peculiar morphological phenotypes in ornamental goldfish. Among them, the twin-tail and dorsal-finless phenotypes have particularly intrigued early and recent researchers, as equivalent morphologies are extremely rare in nature. These two mutated phenotypes appeared almost simultaneously within a short time frame and were fixed in several strains. However, little is known about how these two different mutations could have co-occurred during such a short time period. Here, we demonstrate that the chordin gene, a key factor in dorsal-ventral patterning, is responsible not only for the twin-tail phenotype but also for the dorsal-finless phenotype. Our F2 backcrossing and functional analyses revealed that the penetrance/expressivity of the dorsal-finless phenotype can be suppressed by the wild-type allele of chdS. Based on these findings, we propose that chdS^wt may have masked the expression of the dorsal-finless phenotype, acting as a capacitor buffering gene to allow accumulation of genetic mutations. Once this gene lost its original function in the twin-tail goldfish lineages, the dorsal-finless phenotype could be highly expressed. Thus, this study experimentally demonstrates that the rapid genetic fixation of morphological mutations during a short domestication time period may be related to the robustness of embryonic developmental mechanisms.

Turing's turtles all the way down: A conserved role of EDAR in the carapacial ridge suggests a deep homology of prepatterns across ectodermal appendages

The scutes of the turtle shell are epidermal shields that begin their formation during the early stages of shell development. Like other skin appendages, turtle scutes are hypothesized to be patterned by reaction-diffusion systems. We have previously established ex vivo and in silico systems to study these mechanisms experimentally and have further shown that mathematical models can explain the dynamics of the induction of turtle scute primordia and the generation of final scute architecture. Using these foundations, we expand our current knowledge and test the roles of ectodysplasin and activin signaling in the development of turtle scutes. We find that these molecules play important roles in the prepatterning of scute primordia along the carapacial ridge and show that blocking Edar signaling may lead to a complete loss of marginal scute primordia. We show that it is possible to reproduce these observations using simple mathematical modeling, thereby suggesting a stabilizing role for ectodysplasin within the reaction-diffusion mechanisms. Finally, we argue that our findings further entrench turtle scutes within a class of developmental systems composed of hierarchically nested reaction-diffusion mechanisms, which is conserved across ectodermal organs.

Functional constraints channel mandible shape ontogenies in rodents

In mammals, postnatal growth plays an essential role in the acquisition of the adult shape. During this period, the mandible undergoes many changing functional constraints, leading to spatialization of bone formation and remodelling to accommodate various dietary and behavioural changes. The interactions between the bone, muscles and teeth drive this developmental plasticity, which, in turn, could lead to convergences in the developmental processes constraining the directionality of ontogenies, their evolution and thus the adult shape variation. To test the importance of the interactions between tissues in shaping the ontogenetic trajectories, we compared the mandible shape at five postnatal stages on three rodents: the house mouse, the Mongolian gerbil and the golden hamster, using geometric morphometrics. After an early shape differentiation, by both longer gestation and allometric scaling in gerbils or early divergence of postnatal ontogeny in hamsters in comparison with the mouse, the ontogenetic trajectories appear more similar around weaning. The changes in muscle load associated with new food processing and new behaviours at weaning seem to impose similar physical constraints on the mandible, driving the convergences of the ontogeny at that stage despite an early anatomical differentiation. Nonetheless, mice present a rather different timing compared with gerbils or hamsters.

Breaking constraints: The development and evolution of extreme fin morphology in the Bramidae

The developmental process establishes the foundation upon which natural selection may act. In that same sense, it is inundated with numerous constraints that work to limit the directions in which a phenotype may respond to selective pressures. Extreme phenotypes have been used in the past to identify tradeoffs and constraints and may aid in recognizing how alterations to the Baupläne can influence the trajectories of lineages. The Bramidae, a family of Scombriformes consisting of 20 extant species, are unique in that five species greatly deviate from the stout, ovaloid bodies that typify the bramids. The Ptericlinae, or fanfishes, are instead characterized by relatively elongated body plans and extreme modifications to their medial fins. Here, we explore the development of Bramidae morphologies and examine them through a phylogenetic lens to investigate the concepts of developmental and evolutionary constraints. Contrary to our predictions that the fanfishes had been constrained by inherited properties of an ancestral state, we find that the fanfishes exhibit both increased rates of trait evolution and differ substantially from the other bramids in their developmental trajectories. Conversely, the remaining bramid genera differ little, both among one another and in comparison, to the sister family Caristiidae. In all, our data suggest that the fanfishes have broken constraints, thereby allowing them to mitigate trade‐offs on distinctive aspects of morphology.

Nucleolin loss of function leads to aberrant Fibroblast Growth Factor signaling and craniofacial anomalies

Ribosomal RNA (rRNA) transcription and ribosome biogenesis are global processes required for growth and proliferation of all cells, yet perturbation of these processes in vertebrates leads to tissue-specific defects termed ribosomopathies. Mutations in rRNA transcription and processing proteins often lead to craniofacial anomalies; however, the cellular and molecular reasons for these defects are poorly understood. Therefore, we examined the function of the most abundant nucleolar phosphoprotein, Nucleolin (Ncl), in vertebrate development. ncl mutant (ncl-/-) zebrafish present with craniofacial anomalies such as mandibulofacial hypoplasia. We observed that ncl-/- mutants exhibited decreased rRNA synthesis and p53-dependent apoptosis, consistent with a role in ribosome biogenesis. However, we found that Nucleolin also performs functions not associated with ribosome biogenesis. We discovered that the half-life of fgf8a mRNA was reduced in ncl-/- mutants, which perturbed Fgf signaling, resulting in misregulated Sox9a-mediated chondrogenesis and Runx2-mediated osteogenesis. Consistent with this model, exogenous FGF8 treatment significantly rescued the cranioskeletal phenotype in ncl-/- zebrafish, suggesting that Nucleolin regulates osteochondroprogenitor differentiation. Our work has therefore uncovered tissue-specific functions for Nucleolin in rRNA transcription and post-transcriptional regulation of growth factor signaling during embryonic craniofacial development.

MusMorph, a database of standardized mouse morphology data for morphometric meta-analyses

Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N = 10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase ( www.facebase.org , https://doi.org/10.25550/3-HXMC ) and GitHub ( https://github.com/jaydevine/MusMorph ).

Extend mixed models to multilayer neural networks for genomic prediction including intermediate omics data

With the growing amount and diversity of intermediate omics data complementary to genomics (e.g. DNA methylation, gene expression, and protein abundance), there is a need to develop methods to incorporate intermediate omics data into conventional genomic evaluation. The omics data help decode the multiple layers of regulation from genotypes to phenotypes, thus forms a connected multilayer network naturally. We developed a new method named NN-MM to model the multiple layers of regulation from genotypes to intermediate omics features, then to phenotypes, by extending conventional linear mixed models ("MM") to multilayer artificial neural networks ("NN"). NN-MM incorporates intermediate omics features by adding middle layers between genotypes and phenotypes. Linear mixed models (e.g. pedigree-based BLUP, GBLUP, Bayesian Alphabet, single-step GBLUP, or single-step Bayesian Alphabet) can be used to sample marker effects or genetic values on intermediate omics features, and activation functions in neural networks are used to capture the nonlinear relationships between intermediate omics features and phenotypes. NN-MM had significantly better prediction performance than the recently proposed single-step approach for genomic prediction with intermediate omics data. Compared to the single-step approach, NN-MM can handle various patterns of missing omics measures and allows nonlinear relationships between intermediate omics features and phenotypes. NN-MM has been implemented in an open-source package called "JWAS".

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Proteostasis deregulation as a driver of C9ORF72 pathogenesis

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are two related neurodegenerative disorders that display overlapping features. The hexanucleotide repeat expansion GGGGCC (G₄ C₂ ) in C9ORF72 gene has been causally linked to both ALS and FTD emergence, thus opening a novel potential therapeutic target for disease intervention. The main driver of C9ORF72 pathology is the disruption of distinct cellular processes involved in the function of the proteostasis network. Here we discuss main findings relating to the induction of neurodegeneration by C9ORF72 mutation and proteostasis deregulation, highlighting the role of the endoplasmic reticulum stress, nuclear transport, and autophagy in the disease process. We further discuss possible points of intervention to target proteostasis mediators to treat C9ORF72-linked ALS/FTD.

Commonalities and evolutionary divergences of mandible shape ontogenies in rodents

In mammals, significant changes take place during postnatal growth, linked to changes in diet (from sucking to gnawing). During this period, mandible development is highly interconnected with muscle growth and the epigenetic interactions between muscle and bone control the spatialization of bone formation and remodelling in response to biomechanical strain. This mechanism contributes to postnatal developmental plasticity and may have influenced the course of evolutionary divergences between species and clades. We sought to model postnatal changes at a macroevolutionary scale by analysing ontogenetic trajectories of mandible shape across 16 species belonging mainly to two suborders of Rodents, Myomorpha and Hystricomorpha, which differ in muscle attachments, tooth growth and life-history traits. Myomorpha species present a much stronger magnitude of changes over a shorter growth period. Among Hystricomorpha, part of the observed adult shape is set up prenatally, and most postnatal trajectories are genus-specific, which agrees with nonlinear developmental trajectories over longer gestational periods. Beside divergence at large scale, we find some collinearities between evolutionary and developmental trajectories. A common developmental trend was also observed, leading to enlargement of the masseter fossa during postnatal growth. The tooth growth, especially hypselodonty, seems to be a major driver of divergences of postnatal trajectories. These muscle- and tooth-related effects on postnatal trajectories suggest opportunities for developmental plasticity in the evolution of the mandible shape, opportunities that may have differed across Rodent clades.

Connexin 43 contributes to phenotypic robustness of the mouse skull

BACKGROUND

We compared skull shape and variation among genetically modified mice that exhibit different levels of connexin43 (Cx43) channel function, to determine whether Cx43 contributes to craniofacial phenotypic robustness. Specifically, we used two heterozygous mutant mouse models (G60S/+ and I130T/+) that, when compared to their wildtype counterparts, have an ~80% and ~50% reduction in Cx43 function, respectively.

RESULTS

Both mutant strains showed significant differences in skull shape compared to wildtype littermates and while these differences were more severe in the G60S/+ mouse, shape differences were localized to similar regions of the skull in both mutants. However, increased skull shape variation was observed in G60S/+ mutants only. Additionally, covariation of skull structures was disrupted in the G60S/+ mutants only, indicating that while a 50% reduction in Cx43 function is sufficient to cause a shift in mean skull shape, the threshold for Cx43 function for disrupting craniofacial phenotypic robustness is lower.

CONCLUSIONS

Collectively, our results indicate Cx43 can contribute to phenotypic robustness of the skull through a nonlinear relationship between Cx43 gap junctional function and phenotypic outcomes.

Ontogenetic and in silico models of spatial-packing in the hypermuscular mouse skull

Networks linking single genes to multiple phenotypic outcomes can be founded on local anatomical interactions as well as on systemic factors like biochemical products. Here we explore the effects of such interactions by investigating the competing spatial demands of brain and masticatory muscle growth within the hypermuscular myostatin-deficient mouse model and in computational simulations. Mice that lacked both copies of the myostatin gene (-/-) and display gross hypermuscularity, and control mice that had both copies of the myostatin gene (+/+) were sampled at 1, 7, 14 and 28 postnatal days. A total of 48 mice were imaged with standard as well as contrast-enhanced microCT. Size metrics and landmark configurations were collected from the image data and were analysed alongside in silico models of tissue expansion. Findings revealed that: masseter muscle volume was smaller in -/- mice at day 1 but became, and remained thereafter, larger by 7 days; -/- endocranial volumes begin and remained smaller; -/- enlargement of the masticatory muscles was associated with caudolateral displacement of the calvarium, lateral displacement of the zygomatic arches, and slight dorsal deflection of the face and basicranium. Simulations revealed basicranial retroflexion (flattening) and dorsal deflection of the face associated with muscle expansion and abrogative covariations of basicranial flexion and ventral facial deflection associated with endocranial expansion. Our findings support the spatial-packing theory and highlight the importance of understanding the harmony of competing spatial demands that can shape and maintain mammalian skull architecture during ontogeny.

A role for the fusogen eff-1 in epidermal stem cell number robustness in Caenorhabditis elegans

Developmental patterning in Caenorhabditis elegans is known to proceed in a highly stereotypical manner, which raises the question of how developmental robustness is achieved despite the inevitable stochastic noise. We focus here on a population of epidermal cells, the seam cells, which show stem cell-like behaviour and divide symmetrically and asymmetrically over post-embryonic development to generate epidermal and neuronal tissues. We have conducted a mutagenesis screen to identify mutants that introduce phenotypic variability in the normally invariant seam cell population. We report here that a null mutation in the fusogen eff-1 increases seam cell number variability. Using time-lapse microscopy and single molecule fluorescence hybridisation, we find that seam cell division and differentiation patterns are mostly unperturbed in eff-1 mutants, indicating that cell fusion is uncoupled from the cell differentiation programme. Nevertheless, seam cell losses due to the inappropriate differentiation of both daughter cells following division, as well as seam cell gains through symmetric divisions towards the seam cell fate were observed at low frequency. We show that these stochastic errors likely arise through accumulation of defects interrupting the continuity of the seam and changing seam cell shape, highlighting the role of tissue homeostasis in suppressing phenotypic variability during development.

Developmental plasticity reveals hidden fish phenotypes and enables morphospace diversification

The establishment of a given phenotype is only one expression from a range of hidden developmental possibilities. Developmental plasticity at hidden reaction norms might elicit phenotypic diversification under new developmental environments. Current discussion benefits from empirical analyses that integrate multiple environmental stimuli to evaluate how plastic responses may shape phenotypic variation. We raised Megaleporinus macrocephalus fish in different environmental settings to address contributions of developmental plasticity for emergence of new phenotypes and subsequent morphospace diversification. Plastic morphotypes were evaluated at two complementary scales, the M. macrocephalus morphospace and the higher taxonomic level of Anostomidae family. Morphospace analyses demonstrated that developmental plasticity quickly releases distinct head morphotypes that were hidden in the parental monomorphic population. Plastic morphotypes occupied discrete and previously unfilled morphospace regions, a result obtained from comparisons with a control population and in analyses including several Anostomidae species. Plastic responses involved adjustments in shape and relative position of head bonesets, and fish raised under specific environmental combinations rescued phenotypic patterns described for different genera. Therefore, developmental plasticity possibly contributes to adaptive radiation in Anostomidae. Results illustrate how plastic responses enable morphospace diversification and contribute to evolution.

A model of developmental canalization, applied to human cranial form

Developmental mechanisms that canalize or compensate perturbations of organismal development (targeted or compensatory growth) are widely considered a prerequisite of individual health and the evolution of complex life, but little is known about the nature of these mechanisms. It is even unclear if and how a “target trajectory” of individual development is encoded in the organism’s genetic-developmental system or, instead, emerges as an epiphenomenon. Here we develop a statistical model of developmental canalization based on an extended autoregressive model. We show that under certain assumptions the strength of canalization and the amount of canalized variance in a population can be estimated, or at least approximated, from longitudinal phenotypic measurements, even if the target trajectories are unobserved. We extend this model to multivariate measures and discuss reifications of the ensuing parameter matrix. We apply these approaches to longitudinal geometric morphometric data on human postnatal craniofacial size and shape as well as to the size of the frontal sinuses. Craniofacial size showed strong developmental canalization during the first 5 years of life, leading to a 50% reduction of cross-sectional size variance, followed by a continual increase in variance during puberty. Frontal sinus size, by contrast, did not show any signs of canalization. Total variance of craniofacial shape decreased slightly until about 5 years of age and increased thereafter. However, different features of craniofacial shape showed very different developmental dynamics. Whereas the relative dimensions of the nasopharynx showed strong canalization and a reduction of variance throughout postnatal development, facial orientation continually increased in variance. Some of the signals of canalization may owe to independent variation in developmental timing of cranial components, but our results indicate evolved, partly mechanically induced mechanisms of canalization that ensure properly sized upper airways and facial dimensions.

Developmental mechanisms that canalize or compensate perturbations of organismal development are a prerequisite of individual health and the evolution of complex life. However, surprisingly little is known about these mechanisms, partly because the “target trajectories” of individual development cannot be directly observed. Here we develop a statistical model of developmental canalization that allows one to estimate the strength of canalization and the amount of canalized variance in a population even if the target trajectories are unobserved. We applied these approaches to data on human postnatal craniofacial growth. Whereas overall craniofacial size was strongly canalized during the first 5 years of age, frontal sinus size did not show any signs of canalization. The relative dimensions of the nasopharynx showed strong canalization and a reduction of variance throughout postnatal development, whereas other shape features, such as facial orientation, continually increased in variance. Our results indicate evolved, partly mechanically induced mechanisms of canalization that ensure properly sized upper airways and facial dimensions.

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

Effects of Reduced Connexin43 Function on Mandibular Morphology and Osteogenesis in Mutant Mouse Models of Oculodentodigital Dysplasia

Mutations in the gene encoding the gap-junctional protein connexin43 (Cx43) are the cause of the human disease oculodentodigital dysplasia (ODDD). The mandible is often affected in this disease, with clinical reports describing both mandibular overgrowth and conversely, retrognathia. These seemingly opposing observations underscore our relative lack of understanding of how ODDD affects mandibular morphology. Using two mutant mouse models that mimic the ODDD phenotype (I130T/+ and G60S/+), we sought to uncover how altered Cx43 function may affect mandibular development. Specifically, mandibles of newborn mice were imaged using micro-CT, to enable statistical comparisons of shape. Tissue-level comparisons of key regions of the mandible were conducted using histomorphology, and we quantified the mRNA expression of several cartilage and bone cell differentiation markers. Both G60S/+ and I130T/+ mutant mice had altered mandibular morphology compared to their wildtype counterparts, and the morphological effects were similarly localized for both mutants. Specifically, the biggest phenotypic differences in mutant mice were focused in regions exposed to mechanical forces, such as alveolar bone, muscular attachment sites, and articular surfaces. Histological analyses revealed differences in ossification of the intramembranous bone of the mandibles of both mutant mice compared to their wildtype littermates. However, chondrocyte organization within the secondary cartilages of the mandible was unaffected in the mutant mice. Overall, our results suggest that the morphological differences seen in G60S/+ and I130T/+ mouse mandibles are due to delayed ossification and suggest that mechanical forces may exacerbate the effects of ODDD on the skeleton.

The development, patterning and evolution of neural crest cell differentiation into cartilage and bone

Neural crest cells are a vertebrate-specific migratory, multipotent cell population that give rise to a diverse array of cells and tissues during development. Cranial neural crest cells, in particular, generate cartilage, bone, tendons and connective tissue in the head and face as well as neurons, glia and melanocytes. In this review, we focus on the chondrogenic and osteogenic potential of cranial neural crest cells and discuss the roles of Sox9, Runx2 and Msx1/2 transcription factors and WNT, FGF and TGFβ signaling pathways in regulating neural crest cell differentiation into cartilage and bone. We also describe cranioskeletal defects and disorders arising from gain or loss-of-function of genes that are required for patterning and differentiation of cranial neural crest cells. Finally, we discuss the evolution of skeletogenic potential in neural crest cells and their function as a conduit for intraspecies and interspecies variation, and the evolution of craniofacial novelties.

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.

Selective breeding modifies mef2ca mutant incomplete penetrance by tuning the opposing Notch pathway

Deleterious genetic mutations allow developmental biologists to understand how genes control development. However, not all loss of function genetic mutants develop phenotypic changes. Many deleterious mutations only produce a phenotype in a subset of mutant individuals, a phenomenon known as incomplete penetrance. Incomplete penetrance can confound analyses of gene function and our understanding of this widespread phenomenon remains inadequate. To better understand what controls penetrance, we capitalized on the zebrafish mef2ca mutant which produces craniofacial phenotypes with variable penetrance. Starting with a characterized mef2ca loss of function mutant allele, we used classical selective breeding methods to generate zebrafish strains in which mutant-associated phenotypes consistently appear with low or high penetrance. Strikingly, our selective breeding for low penetrance converted the mef2ca mutant allele behavior from homozygous lethal to homozygous viable. Meanwhile, selective breeding for high penetrance converted the mef2ca mutant allele from fully recessive to partially dominant. Comparing the selectively-bred low- and high-penetrance strains revealed that the strains initially respond similarly to the mutation, but then gene expression differences between strains emerge during development. Thus, altered temporal genetic circuitry can manifest through selective pressure to modify mutant penetrance. Specifically, we demonstrate differences in Notch signaling between strains, and further show that experimental manipulation of the Notch pathway phenocopies penetrance changes occurring through selective breeding. This study provides evidence that penetrance is inherited as a liability-threshold trait. Our finding that vertebrate animals can overcome a deleterious mutation by tuning genetic circuitry complements other reported mechanisms of overcoming deleterious mutations such as transcriptional adaptation of compensatory genes, alternative mRNA splicing, and maternal deposition of wild-type transcripts, which are not observed in our system. The selective breeding approach and the resultant genetic circuitry change we uncovered advances and expands our current understanding of genetic and developmental resilience.

Some deleterious gene mutations only affect a subset of genetically mutant animals. This widespread phenomenon, known as mutant incomplete penetrance, complicates discovery of causative gene mutations in both model organisms and human disease. This study utilized the zebrafish mef2ca transcription factor mutant that produces craniofacial skeleton defects with incomplete penetrance. Selectively breeding zebrafish families for low- or high-penetrance mutants for many generations created different zebrafish strains with consistently low or high penetrance. Comparing these strains allowed us to gain insight into the mechanisms that control penetrance. Specifically, genes under the control of mef2ca are initially similarly expressed between the two strains, but differences between strains emerge during development. We found that genetic manipulation of these downstream genes mimics the effects of our selective breeding. Thus, selective breeding for penetrance can change the genetic circuitry downstream of the mutated gene. We propose that small differences in gene circuitry between individuals is one mechanism underlying susceptibility or resilience to genetic mutations.

How the diversity of the faces arises

BACKGROUND

The evolution of the face is crucial for each species to adapt to different diets, environments, and in some species, to promote social interaction. The diversity in the shapes of the face results from divergence in the process of facial development that begins during early embryonic development.

HIGHLIGHTS

Here we review the recent advancements in the understanding of the genetic, epigenetic, molecular, and cellular basis of facial diversity. We also review the robustness of facial development and how it relates to the evolution of the face. Finally, we discuss the current challenges in achieving a deeper understanding of facial diversity.

CONCLUSION

We have gained much knowledge with respect to cis-regulatory elements, gene expression, cellular behavior, and the physical forces in facial development in the past two decades. Significant interdisciplinary work is needed to integrate these varied pieces of information into a complete picture of how the diversity of faces arises.

Developmental plasticity associated with early structural integration and evolutionary patterns: Examples of developmental bias and developmental facilitation in the skeletal system

The relation of developmental plasticity to evolutionary diversification is a key component of evolutionary theory involving developmental bias, but the basis of the relationship varies among traits and among taxa. Here I review some scenarios of how structural integration during early organogenesis could influence this relationship. When condensations are highly integrated and dependent on each other during early organogenesis, both plasticity and evolution are restricted, for example size proportions in molar tooth rows and phalanges within a digit. When similar condensations develop and remain separate (in tracheal cartilages and feather buds), they show high levels of variation and diversity in number but not in shape and size, at least at early stages. When non-similar structures form separately and then integrate while still undergoing patterning, high levels of plasticity (in number, size, shape; in rib uncinate processes) or new dimensions of ecologically-significant variation (cusp offset, in mammal teeth) are seen. Although each of these structural integration scenarios is unique, the modulation of evolvability is detectable and informative. Parsing the influence of structural integration at these developmental levels, rather than later-stage structural correlations or only through genetic covariation, may be necessary to advance understanding of evolvability of the phenotype.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse

BACKGROUND

Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression.

RESULTS

Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape.

CONCLUSIONS

Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

Epistasis regulates the developmental stability of the mouse craniofacial shape

Fluctuating asymmetry is a classic concept linked to organismal development. It has traditionally been used as a measure of developmental instability, which is the inability of an organism to buffer environmental fluctuations during development. Developmental stability has a genetic component that influences the final phenotype of the organism and can lead to congenital disorders. According to alternative hypotheses, this genetic component might be either the result of additive genetic effects or a by-product of developmental gene networks. Here we present a genome-wide association study of the genetic architecture of fluctuating asymmetry of the skull shape in mice. Geometric morphometric methods were applied to quantify fluctuating asymmetry: we estimated fluctuating asymmetry as Mahalanobis distances to the mean asymmetry, correcting first for genetic directional asymmetry. We applied the marginal epistasis test to study epistasis among genomic regions. Results showed no evidence of additive effects but several interacting regions significantly associated with fluctuating asymmetry. Among the candidate genes overlapping these interacting regions we found an over-representation of genes involved in craniofacial development. A gene network is likely to be associated with skull developmental stability, and genes originally described as buffering genes (e.g., Hspa2) might occupy central positions within these networks, where regulatory elements may also play an important role. Our results constitute an important step in the exploration of the molecular roots of developmental stability and the first empirical evidence about its genetic architecture.

Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Mouse Skull Mean Shape and Shape Robustness Rely on Different Genetic Architectures and Different Loci

The genetic architecture of skull shape has been extensively studied in mice and the results suggest a highly polygenic and additive basis. In contrast few studies have explored the genetic basis of the skull variability. Canalization and developmental stability are the two components of phenotypic robustness. They have been proposed to be emergent properties of the genetic networks underlying the development of the trait itself, but this hypothesis has been rarely tested empirically. Here we use outbred mice to investigate the genetic architecture of canalization of the skull shape by implementing a genome-wide marginal epistatic test on 3D geometric morphometric data. The same data set had been used previously to explore the genetic architecture of the skull mean shape and its developmental stability. Here, we address two questions: (1) Are changes in mean shape and changes in shape variance associated with the same genomic regions? and (2) Do canalization and developmental stability rely on the same loci and genetic architecture and do they involve the same patterns of shape variation? We found that unlike skull mean shape, among-individual shape variance and fluctuating asymmetry (FA) show a total lack of additive effects. They are both associated with complex networks of epistatic interactions involving many genes (protein-coding and regulatory elements). Remarkably, none of the genomic loci affecting mean shape contribute these networks despite their enrichment for genes involved in craniofacial variation and diseases. We also found that the patterns of shape FA and individual variation are largely similar and rely on similar multilocus epistatic genetic networks, suggesting that the processes channeling variation within and among individuals are largely common. However, the loci involved in these two networks are completely different. This in turn underlines the difference in the origin of the variation at these two levels, and points at buffering processes that may be specific to each level.

Developmental processes regulate craniofacial variation in disease and evolution

Variation in development mediates phenotypic differences observed in evolution and disease. Although the mechanisms underlying phenotypic variation are still largely unknown, recent research suggests that variation in developmental processes may play a key role. Developmental processes mediate genotype-phenotype relationships and consequently play an important role regulating phenotypes. In this review, we provide an example of how shared and interacting developmental processes may explain convergence of phenotypes in spliceosomopathies and ribosomopathies. These data also suggest a shared pathway to disease treatment. We then discuss three major mechanisms that contribute to variation in developmental processes: genetic background (gene-gene interactions), gene-environment interactions, and developmental stochasticity. Finally, we comment on evolutionary alterations to developmental processes, and the evolution of disease buffering mechanisms.

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.