The developmental-genetics of canalization

Dosage-dependent effects of FGFR2W290R mutation on craniofacial shape and cellular dynamics of the basicranial synchondroses

Craniosynostosis is a common yet complex birth defect, characterized by premature fusion of the cranial sutures that can be syndromic or nonsyndromic. With over 180 syndromic associations, reaching genetic diagnoses and understanding variations in underlying cellular mechanisms remains a challenge. Variants of FGFR2 are highly associated with craniosynostosis and warrant further investigation. Using the missense mutation FGFR2W290R , an effective mouse model of Crouzon syndrome, craniofacial features were analyzed using geometric morphometrics across developmental time (E10.5-adulthood, n = 665 total). Given the interrelationship between the cranial vault and basicranium in craniosynostosis patients, the basicranium and synchondroses were analyzed in perinates. Embryonic time points showed minimal significant shape differences. However, hetero- and homozygous mutant perinates and adults showed significant differences in shape and size of the cranial vault, face, and basicranium, which were associated with cranial doming and shortening of the basicranium and skull. Although there were also significant shape and size differences associated with the basicranial bones and clear reductions in basicranial ossification in cleared whole-mount samples, there were no significant alterations in chondrocyte cell shape, size, or orientation along the spheno-occipital synchondrosis. Finally, shape differences in the cranial vault and basicranium were interrelated at perinatal stages. These results point toward the possibility that facial shape phenotypes in craniosynostosis may result in part from pleiotropic effects of the causative mutations rather than only from the secondary consequences of the sutural defects, indicating a novel direction of research that may shed light on the etiology of the broad changes in craniofacial morphology observed in craniosynostosis syndromes.

A meta-analysis of Boolean network models reveals design principles of gene regulatory networks

Gene regulatory networks (GRNs) play a central role in cellular decision-making. Understanding their structure and how it impacts their dynamics constitutes thus a fundamental biological question. GRNs are frequently modeled as Boolean networks, which are intuitive, simple to describe, and can yield qualitative results even when data are sparse. We assembled the largest repository of expert-curated Boolean GRN models. A meta-analysis of this diverse set of models reveals several design principles. GRNs exhibit more canalization, redundancy, and stable dynamics than expected. Moreover, they are enriched for certain recurring network motifs. This raises the important question why evolution favors these design mechanisms.

The behavioral origins of phylogenic responses and ontogenic habits

An examination of innate behavior and its possible origins suggests parallels with the formation of habitual behavior. Inflexible but adaptive responses-innate reflexive behavior, Pavlovian conditioned responses, and operant habits-may have evolved from variable behavior in phylogeny and ontogeny. This form of "plasticity-first" scientific narrative was unpopular post-Darwin but has recently gained credibility in evolutionary biology. The present article seeks to identify originating events and contingencies contributing to such inflexible but adaptive behavior at both phylogenic and ontogenic levels of selection. In ontogeny, the development of inflexible performance (i.e., habit) from variable operant behavior is reminiscent of the genetic accommodation of initially variable phylogenic traits. The effects characteristic of habit (e.g., unresponsiveness to reinforcer devaluation) are explicable as the result of a conflict between behaviors at distinct levels of selection. The present interpretation validates the practice of seeking hard analogies between evolutionary biology and operant behavior. Finding such parallels implies the validity of a claim that organismal behavior, both innate and learned, is a product of selection by consequences. A complete and coherent account of organismal behavior may ultimately focus on functional selective histories in much the same way evolutionary biology does with its subject matter.

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.

mtDNA Single-Nucleotide Variants Associated with Type 2 Diabetes

Type 2 diabetes (T2D) is a chronic systemic disease with a complex etiology, characterized by insulin resistance and mitochondrial dysfunction in various cell tissues. To explore this relationship, we conducted a secondary analysis of complete mtDNA sequences from 1261 T2D patients and 1105 control individuals. Our findings revealed significant associations between certain single-nucleotide polymorphisms (SNPs) and T2D. Notably, the variants m.1438A>G (rs2001030) (controls: 32 [27.6%], T2D: 84 [72.4%]; OR: 2.46; 95%CI: 1.64-3.78; p < 0.001), m.14766C>T (rs193302980) (controls: 498 [36.9%], T2D: 853 [63.1%]; OR: 2.57, 95%CI: 2.18-3.04, p < 0.001), and m.16519T>C (rs3937033) (controls: 363 [43.4%], T2D: 474 [56.6%]; OR: 1.24, 95%CI: 1.05-1.47, p = 0.012) were significantly associated with the likelihood of developing diabetes. The variant m.16189T>C (rs28693675), which has been previously documented in several studies across diverse populations, showed no association with T2D in our analysis (controls: 148 [13.39] T2D: 171 [13.56%]; OR: 1.03; 95%CI: 0.815-1.31; p = 0.83). These results provide evidence suggesting a link between specific mtDNA polymorphisms and T2D, possibly related to association rules, topological patterns, and three-dimensional conformations associated with regions where changes occur, rather than specific point mutations in the sequence.

Ethanol Drives Evolution of Hsp90-Dependent Robustness by Redundancy in Yeast Domestication

Protein folding promotes and constrains adaptive evolution. We uncover this surprising duality in the role the protein-folding chaperone Hsp90 plays in mediating the interplay between proteome and the genome which acts to maintain the integrity of yeast metabolism in the face of proteotoxic stressors in anthropic niches. Of great industrial relevance, ethanol concentrations generated by fermentation in the making of beer and bread disrupt critical Hsp90-dependent nodes of metabolism and exert strong selective pressure for increased copy number of key genes encoding components of these nodes, yielding the classical genetic signatures of beer and bread domestication. This work establishes a mechanism of adaptive canalization in an ecology of major economic significance and highlights Hsp90-contingent variation as an important source of phantom heritability in complex traits.

Zebrafish crocc2 mutants exhibit divergent craniofacial shape, misregulated variability, and aberrant cartilage morphogenesis

BACKGROUND

Phenotypic variation is of paramount importance in development, evolution, and human health; however, the molecular mechanisms that influence organ shape and shape variability are not well understood. During craniofacial development, the behavior of skeletal precursors is regulated by both biochemical and environmental inputs, and the primary cilia play critical roles in transducing both types of signals. Here, we examine a gene that encodes a key constituent of the ciliary rootlets, crocc2, and its role in cartilage morphogenesis in larval zebrafish.

RESULTS

Geometric morphometric analysis of crocc2 mutants revealed altered craniofacial shapes and expanded variation. At the cellular level, we observed altered chondrocyte shapes and planar cell polarity across multiple stages in crocc2 mutants. Notably, cellular defects were specific to areas that experience direct mechanical input. Cartilage cell number, apoptosis, and bone patterning were not affected in crocc2 mutants.

CONCLUSIONS

Whereas "regulatory" genes are widely implicated in patterning the craniofacial skeleton, genes that encode "structural" aspects of the cell are increasingly implicated in shaping the face. Our results add crocc2 to this list, and demonstrate that it affects craniofacial geometry and canalizes phenotypic variation. We propose that it does so via mechanosensing, possibly through the ciliary rootlet. If true, this would implicate a new organelle in skeletal development and evolution.

Robust morphogenesis by chaotic dynamics

This research illustrates that complex dynamics of gene products enable the creation of any prescribed cellular differentiation patterns. These complex dynamics can take the form of chaotic, stochastic, or noisy chaotic dynamics. Based on this outcome and previous research, it is established that a generic open chemical reactor can generate an exceptionally large number of different cellular patterns. The mechanism of pattern generation is robust under perturbations and it is based on a combination of Turing's machines, Turing instability and L. Wolpert's gradients. These results can help us to explain the formidable adaptive capacities of biochemical systems.

Canalization of Phenotypes-When the Transcriptome is Constantly but Weakly Perturbed

Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

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.

Developmental disruption and restoration of brain synaptome architecture in the murine Pax6 neurodevelopmental disease model

Neurodevelopmental disorders of genetic origin delay the acquisition of normal abilities and cause disabling phenotypes. Nevertheless, spontaneous attenuation and even complete amelioration of symptoms in early childhood and adolescence can occur in many disorders, suggesting that brain circuits possess an intrinsic capacity to overcome the deficits arising from some germline mutations. We examined the molecular composition of almost a trillion excitatory synapses on a brain-wide scale between birth and adulthood in mice carrying a mutation in the homeobox transcription factor Pax6, a neurodevelopmental disorder model. Pax6 haploinsufficiency had no impact on total synapse number at any age. By contrast, the molecular composition of excitatory synapses, the postnatal expansion of synapse diversity and the acquisition of normal synaptome architecture were delayed in all brain regions, interfering with networks and electrophysiological simulations of cognitive functions. Specific excitatory synapse types and subtypes were affected in two key developmental age-windows. These phenotypes were reversed within 2-3 weeks of onset, restoring synapse diversity and synaptome architecture to the normal developmental trajectory. Synapse subtypes with rapid protein turnover mediated the synaptome remodeling. This brain-wide capacity for remodeling of synapse molecular composition to recover and maintain the developmental trajectory of synaptome architecture may help confer resilience to neurodevelopmental genetic disorders.

Morphological integration and evolutionary potential of the primate shoulder: Variation among taxa and implications for genetic covariances with the basicranium, pelvis, and arm

Within the primate order, the morphology of the shoulder girdle is immensely variable and has been shown to reflect the functional demands of the upper limb. The observed morphological variation among extant primate taxa consequently has been hypothesized to be driven by selection for different functional demands. Evolutionary analyses of the shoulder girdle often assess this anatomical region, and its traits, individually, therefore implicitly assuming independent evolution of the shoulder girdle. However, the primate shoulder girdle has developmental and functional covariances with the basicranium and pelvic girdle that have been shown to potentially influence its evolution. It is unknown whether these relationships are similar or even present across primate taxa, and how they may affect morphological variation among primates. This study evaluates the strength of covariance and evolutionary potential across four anatomical regions: shoulder girdle, basicranium, pelvis, and distal humerus. Measures of morphological integration and evolutionary potential (conditioned covariance and evolutionary flexibility) are assessed across eight anthropoid primate taxa. Results demonstrate a consistent pattern of morphological constraint within paired anatomical regions across primates. Differences in evolutionary flexibility are observed among primate genera, with humans having the highest evolutionary potential overall. This pattern does not follow functional differences, but rather a separation between monkeys and apes. Therefore, evolutionary hypotheses of primate shoulder girdle morphological variation that evaluate functional demands alone may not account for the effect of these relationships. Collectively, our findings suggest differences in genetic covariance among anatomical regions may have contributed to the observable morphological variation among taxa.

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 ).

Metabolic and immunological responses of Drosophila melanogaster to dietary restriction and bacterial infection differ substantially between genotypes in a population

To respond to changing environmental conditions, a population may either shift toward better‐adapted genotypes or adapt on an individual level. The present work aimed to quantify the relevance of these two processes by comparing the responses of defined Drosophila melanogaster populations to different stressors. To do this, we infected two homogeneous populations (isofemale lines), which differ significantly in fitness, and a synthetic heterogeneous population with a specific pathogen and/or exposed them to food restriction. Pectobacterium carotovorum was used to infect Drosophila larvae either fed standard or protein‐restricted diet. In particular, the two homogeneous groups, which diverged in their fitness, showed considerable differences in all parameters assessed (survivorship, protein and lipid contents, phenol‐oxidase (PO) activity, and antibacterial rate). Under fully nutritious conditions, larvae of the homogeneous population with low fitness exhibited lower survivorship and protein levels, as well as higher PO activity and antibacterial rate compared with the fitter population. A protein‐restricted diet and bacterial infection provoked a decrease in survivorship, and antibacterial rate in most populations. Bacterial infection elicited an opposite response in protein and lipid content in both isofemale lines tested. Interestingly, the heterogeneous population showed a complex response pattern. The response of the heterogeneous population followed the fit genotype in terms of survival and antibacterial activity but followed the unfit genotype in terms of PO activity. In conclusion, our results show that defined genotypes exhibit highly divergent responses to varying stressors that are difficult to predict. Furthermore, the responses of heterogeneous populations do not follow a fixed pattern showing a very high degree of plasticity and differences between different genotypes.

Canalization of the Polygenic Risk for Common Diseases and Traits in the UK Biobank Cohort

Since organisms develop and thrive in the face of constant perturbations due to environmental and genetic variation, species may evolve resilient genetic architectures. We sought evidence for this process, known as canalization, through a comparison of the prevalence of phenotypes as a function of the polygenic score (PGS) across environments in the UK Biobank cohort study. Contrasting seven diseases and three categorical phenotypes with respect to 151 exposures in 408,925 people, the deviation between the prevalence-risk curves was observed to increase monotonically with the PGS percentile in one-fifth of the comparisons, suggesting extensive PGS-by-Environment (PGS×E) interaction. After adjustment for the dependency of allelic effect sizes on increased prevalence in the perturbing environment, cases where polygenic influences are greater or lesser than expected are seen to be particularly pervasive for educational attainment, obesity, and metabolic condition type-2 diabetes. Inflammatory bowel disease analysis shows fewer interactions but confirms that smoking and some aspects of diet influence risk. Notably, body mass index has more evidence for decanalization (increased genetic influence at the extremes of polygenic risk), whereas the waist-to-hip ratio shows canalization, reflecting different evolutionary pressures on the architectures of these weight-related traits. An additional 10 % of comparisons showed evidence for an additive shift of prevalence independent of PGS between exposures. These results provide the first widespread evidence for canalization protecting against disease in humans and have implications for personalized medicine as well as understanding the evolution of complex traits. The findings can be explored through an R shiny app at https://canalization-gibsonlab.shinyapps.io/rshiny/.

The nature of coxofemoral joint pathology across family Canidae

We evaluated coxofemoral joints from museum specimens of: Vulpes lagopus; Vulpes vulpes; Vulpes velox; Nyctereutes procyonoides; Urocyon cinereoargenteus; Aenocyon [Canis] dirus; Canis latrans; Canis lupus lupus; Canis lupus familiaris; C. l. familiaris × latrans; and Canis dingo. Acetabular components included: fossa; articular surface; medial and lateral articular margins; and periarticular surfaces. Acetabular components variably revealed: osteophyte-like features; varying appearance of articular margin rims (especially contour changes); rough bone surfaces (especially fossa and articular surface); and surface wear. Proximal femoral components included: articular surface; articular margin; periarticular surfaces; and joint capsule attachment. Femoral components variably revealed: rough bone surface; bone loss; articular margin osteophyte-like features; caudal post-developmental mineralized prominence; and enthesophytes along the joint capsule attachment. Non-metric multidimensional scaling was used to analyze right-left asymmetric relationships between observed traits, across taxa. Significantly different acetabular trait asymmetry involved only C. latrans-C. l. familiaris; V. vulpes-N. procyonoides, and U. cinereoargenteus-N. procyonoides. There were no significant lateralized differences in proximal femoral traits involving modern canids, ancient and modern C. l. familiaris, or modern vulpines. Thus, the observations were strongly bilateral. We hypothesized high similarity of traits across taxa. The data confirm the hypothesis and strongly suggest broad and deep morphological and mechanistic conservation that almost certainly pre-existed (at least) all modern canids. Further zoological studies are needed to evaluate phylogenic implications in greater detail.

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.

Unexpected variation of human molar size patterns

A tenet of mammalian, including primate dental evolution, is the Inhibitory Cascade Model, where first molar (M1) size predicts in a linear cline the size and onset time of the second (M2) and third (M3) molars: a larger M1 portends a progressively smaller and later-developing M2 and M3. In contemporary modern Homo sapiens, later-developing M3s are less likely to erupt properly. The Inhibitory Cascade Model is also used to predict molar sizes of extinct taxa, including fossil Homo. The extent to which Inhibitory Cascade Model predictions hold in contemporary H. sapiens molars is unclear, including whether this tenet informs about molar initiation, development, and eruption. We tested these questions here. In our radiographic sample of 323 oral quadrants and molar rows from contemporary humans based on mesiodistal crown lengths, we observed the distribution of molar proportions with a central tendency around parity (M1 = M2 = M3) that parsed into 13 distinct molar size ratio patterns. These patterns presented at different frequencies (e.g., M1 > M2 > M3 in about one-third of cases) that reflected whether the molar row was located in the maxilla or mandible and included both linear (e.g., M1 < M2 < M3) and nonlinear molar size ratio progressions (e.g., M1 > M2 < M3). Up to four patterns were found in the same subject's mouth. Lastly, M1 size alone does not predict M3 size, developmental timing, or eruption; rather, M2 size is integral to predicting M3 size. Our study indicates that human molar size is genetically 'softwired' and sensitive to factors local to the human upper jaw vs. lower jaw. The lack of a single stereotypical molar size ratio for contemporary H. sapiens suggests that predictions of fossil H. sapiens molar sizes using the Inhibitory Cascade Model must be made with caution.

The cichlid oral and pharyngeal jaws are evolutionarily and genetically coupled

Evolutionary constraints may significantly bias phenotypic change, while "breaking" from such constraints can lead to expanded ecological opportunity. Ray-finned fishes have broken functional constraints by developing two jaws (oral-pharyngeal), decoupling prey capture (oral jaw) from processing (pharyngeal jaw). It is hypothesized that the oral and pharyngeal jaws represent independent evolutionary modules and this facilitated diversification in feeding architectures. Here we test this hypothesis in African cichlids. Contrary to our expectation, we find integration between jaws at multiple evolutionary levels. Next, we document integration at the genetic level, and identify a candidate gene, smad7, within a pleiotropic locus for oral and pharyngeal jaw shape that exhibits correlated expression between the two tissues. Collectively, our data show that African cichlid evolutionary success has occurred within the context of a coupled jaw system, an attribute that may be driving adaptive evolution in this iconic group by facilitating rapid shifts between foraging habitats, providing an advantage in a stochastic environment such as the East African Rift-Valley.

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.

Organ geometry channels reproductive cell fate in the Arabidopsis ovule primordium

In multicellular organisms, sexual reproduction requires the separation of the germline from the soma. In flowering plants, the female germline precursor differentiates as a single spore mother cell (SMC) as the ovule primordium forms. Here, we explored how organ growth contributes to SMC differentiation. We generated 92 annotated 3D images at cellular resolution in Arabidopsis. We identified the spatio-temporal pattern of cell division that acts in a domain-specific manner as the primordium forms. Tissue growth models uncovered plausible morphogenetic principles involving a spatially confined growth signal, differential mechanical properties, and cell growth anisotropy. Our analysis revealed that SMC characteristics first arise in more than one cell but SMC fate becomes progressively restricted to a single cell during organ growth. Altered primordium geometry coincided with a delay in the fate restriction process in katanin mutants. Altogether, our study suggests that tissue geometry channels reproductive cell fate in the Arabidopsis ovule primordium.

Growing points in attachment disorganization: looking back to advance forward

In this special issue paper we reflect on the next generation of attachment research with a focus on disorganization, a central but still poorly understood topic in this area. We suggest that progress will be facilitated by a return to attachment theory's evolutionary roots, and to the emphasis on biological function that inspired Bowlby's original thinking. Increased interdisciplinary cross-fertilization and collaborations would enable novel and generative research on some of the long-standing questions surrounding attachment disorganization. Accordingly, we present an agenda for future research that encompasses contributions of modern ethology and neurobiology, novel hypotheses based on the concept of adaptive decanalization, connections with neurodevelopmental vulnerability and risk for mental disorders such as schizophrenia, and the possibility of sex differences in the behavioral manifestations of attachment disorganization. We believe that these avenues of theory and research offer exciting potential for innovative work in attachment disorganization in the years ahead.

Exploring the Causality Between Hypothyroidism and Non-alcoholic Fatty Liver: A Mendelian Randomization Study

The etiology of non-alcoholic fatty liver disease (NAFLD) involves complex interaction of genetic and environmental factors. A large number of observational studies have shown that hypothyroidism contributes to a high risk of NAFLD. However, the exact causality is still unknown. Due to the progress of genome-wide association study (GWAS) and the discovery of Mendelian randomization (MR), it is possible to explore the causality between the two diseases. In this study, in order to research into the influence of intermediate phenotypes on outcome, nine independent genetic variants of hypothyroidism obtained from the GWAS were used as instrumental variables (IVs) to perform MR analysis on NAFLD. Since there was no heterogeneity between IVs (P = 0.70), a fixed-effects model was used. The correlation between hypothyroidism and NAFLD was evaluated by using inverse-variance weighted (IVW) method and weighted median method. Then the sensitivity test was analyzed. The results showed that there was a high OR (1.7578; 95%CI 1.1897–2.5970; P = 0.0046) and a low intercept (−0.095; P = 0.431). None of the genetic variants drove the overall result (P < 0.01). Simply, we proved for the first time that the risk of NAFLD increases significantly on patients with hypothyroidism. Furthermore, we explained possible causes of NAFLD caused by hypothyroidism.

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.

Weak genetic signal for phenotypic integration implicates developmental processes as major regulators of trait covariation

Phenotypic integration is an important metric that describes the degree of covariation among traits in a population, and is hypothesized to arise due to selection for shared functional processes. Our ability to identify the genetic and/or developmental underpinnings of integration is marred by temporally overlapping cell-, tissue- and structure-level processes that serve to continually 'overwrite' the structure of covariation among traits through ontogeny. Here, we examine whether traits that are integrated at the phenotypic level also exhibit a shared genetic basis (e.g. pleiotropy). We micro-CT scanned two hard tissue traits, and two soft tissue traits (mandible, pectoral girdle, atrium and ventricle, respectively) from an F5 hybrid population of Lake Malawi cichlids, and used geometric morphometrics to extract 3D shape information from each trait. Given the large degree of asymmetric variation that may reflect developmental instability, we separated symmetric from asymmetric components of shape variation. We then performed quantitative trait loci (QTL) analysis to determine the degree of genetic overlap between shapes. While we found ubiquitous associations among traits at the phenotypic level, except for a handful of notable exceptions, our QTL analysis revealed few overlapping genetic regions. Taken together, this indicates developmental interactions can play a large role in determining the degree of phenotypic integration among traits, and likely obfuscate the genotype to phenotype map, limiting our ability to gain a comprehensive picture of the genetic contributors responsible for phenotypic divergence.

Developing a theoretical evolutionary framework to solve the mystery of parturition initiation

Eutherian mammals have characteristic lengths of gestation that are key for reproductive success, but relatively little is known about the processes that determine the timing of parturition, the process of birth, and how they are coordinated with fetal developmental programs. This issue remains one of biology's great unsolved mysteries and has significant clinical relevance because preterm birth is the leading cause of infant and under 5 year old child mortality worldwide. Here, we consider the evolutionary influences and potential signaling mechanisms that maintain or end pregnancy in eutherian mammals and use this knowledge to formulate general theoretical evolutionary models. These models can be tested through evolutionary species comparisons, studies of experimental manipulation of gestation period and birth timing, and human clinical studies. Understanding how gestation time and parturition are determined will shed light on this fundamental biological process and improve human health through the development of therapies to prevent preterm birth.

CLAVATA Signaling Ensures Reproductive Development in Plants across Thermal Environments

The ability to thrive in diverse environments requires that species maintain development and reproduction despite dynamic conditions. Many developmental processes are stabilized through robust signaling pathways that cooperatively ensure proper development.¹ During reproduction, plants like Arabidopsis thaliana continuously generate flowers on growing indeterminate inflorescences.² Flower primordia initiation and outgrowth depends on the hormone auxin and is robust across diverse environments.^3-6 Here, we show that reproductive development under different thermal conditions requires the integration of multiple pathways regulating auxin-dependent flower production. In colder/ambient temperatures, the receptor complex CLAVATA2/CORYNE (CLV2/CRN) is necessary for continuous flower outgrowth during inflorescence development. CLV2/CRN signaling is independent of CLAVATA1 (CLV1)-related receptor signaling but involves the CLAVATA3 INSENSITIVE RECEPTOR KINASE (CIK) family co-receptors, with higher order cik mutant combinations phenocopying clv2/crn flower outgrowth defects. Developing crn inflorescences display reduced auxin signaling, and restoration of auxin biosynthesis is sufficient to restore flower outgrowth in colder and ambient temperatures. In contrast, at higher temperatures, both clv2/crn signaling and heat-induced auxin biosynthesis via YUCCA family genes are synergistically required to maintain flower development. Our work reveals a novel mechanism integrating peptide hormone and auxin signaling in the regulation of flower development across diverse thermal environments.

Phenotypes to remember: Evolutionary developmental memory capacity and robustness

There is increased awareness of the possibility of developmental memories resulting from evolutionary learning. Genetic regulatory and neural networks can be modelled by analogous formalism raising the important question of productive analogies in principles, processes and performance. We investigate the formation and persistence of various developmental memories of past phenotypes asking how the number of remembered past phenotypes scales with network size, to what extent memories stored form by Hebbian-like rules, and how robust these developmental "devo-engrams" are against networks perturbations (graceful degradation). The analogy between neural and genetic regulatory networks is not superficial in that it allows knowledge transfer between fields that used to be developed separately from each other. Known examples of spectacular phenotypic radiations could partly be accounted for in such terms.

Schizophrenia: Developmental Variability Interacts with Risk Factors to Cause the Disorder: Nonspecific Variability-Enhancing Factors Combine with Specific Risk Factors to Cause Schizophrenia

A new etiological model is proposed for schizophrenia that combines variability-enhancing nonspecific factors acting during development with more specific risk factors. This model is better suited than the current etiological models of schizophrenia, based on the risk factors paradigm, for predicting and/or explaining several important findings about schizophrenia: high co-morbidity rates, low specificity of many risk factors, and persistence in the population of the associated genetic polymorphisms. Compared with similar models, e.g., de-canalization, common psychopathology factor, sexual-selection, or differential sensitivity to the environment, this proposal is more general and integrative. Recently developed research methods have proven the existence of genetic and environmental factors that enhance developmental variability. Applying such methods to newly collected or already available data can allow for testing the hypotheses upon which this model is built. If validated, this model may change the understanding of the etiology of schizophrenia, the research models, and preventionbrk paradigms.

Canalization and Robustness in Human Genetics and Disease

Canalization refers to the evolution of populations such that the number of individuals who deviate from the optimum trait, or experience disease, is minimized. In the presence of rapid cultural, environmental, or genetic change, the reverse process of decanalization may contribute to observed increases in disease prevalence. This review starts by defining relevant concepts, drawing distinctions between the canalization of populations and robustness of individuals. It then considers evidence pertaining to three continuous traits and six domains of disease. In each case, existing genetic evidence for genotype-by-environment interactions is insufficient to support a strong inference of decanalization, but we argue that the advent of genome-wide polygenic risk assessment now makes an empirical evaluation of the role of canalization in preventing disease possible. Finally, the contributions of both rare and common variants to congenital abnormality and adult onset disease are considered in light of a new kerplunk model of genetic effects.

A Registration and Deep Learning Approach to Automated Landmark Detection for Geometric Morphometrics

Geometric morphometrics is the statistical analysis of landmark-based shape variation and its covariation with other variables. Over the past two decades, the gold standard of landmark data acquisition has been manual detection by a single observer. This approach has proven accurate and reliable in small-scale investigations. However, big data initiatives are increasingly common in biology and morphometrics. This requires fast, automated, and standardized data collection. We combine techniques from image registration, geometric morphometrics, and deep learning to automate and optimize anatomical landmark detection. We test our method on high-resolution, micro-computed tomography images of adult mouse skulls. To ensure generalizability, we use a morphologically diverse sample and implement fundamentally different deformable registration algorithms. Compared to landmarks derived from conventional image registration workflows, our optimized landmark data show up to a 39.1% reduction in average coordinate error and a 36.7% reduction in total distribution error. In addition, our landmark optimization produces estimates of the sample mean shape and variance-covariance structure that are statistically indistinguishable from expert manual estimates. For biological imaging datasets and morphometric research questions, our approach can eliminate the time and subjectivity of manual landmark detection whilst retaining the biological integrity of these expert annotations.

Embryonic geometry underlies phenotypic variation in decanalized conditions

During development, many mutations cause increased variation in phenotypic outcomes, a phenomenon termed decanalization. Phenotypic discordance is often observed in the absence of genetic and environmental variations, but the mechanisms underlying such inter-individual phenotypic discordance remain elusive. Here, using the anterior-posterior (AP) patterning of the Drosophila embryo, we identified embryonic geometry as a key factor predetermining patterning outcomes under decanalizing mutations. With the wild-type AP patterning network, we found that AP patterning is robust to variations in embryonic geometry; segmentation gene expression remains reproducible even when the embryo aspect ratio is artificially reduced by more than twofold. In contrast, embryonic geometry is highly predictive of individual patterning defects under decanalized conditions of either increased bicoid (bcd) dosage or bcd knockout. We showed that the phenotypic discordance can be traced back to variations in the gap gene expression, which is rendered sensitive to the geometry of the embryo under mutations.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Hsp90 and complex birth defects: A plausible mechanism for the interaction of genes and environment

How genes and environment interact to cause birth defects is not well understood, but key to developing new strategies to modify risk. The threshold model has been proposed to represent this complex interaction. This model stipulates that while environmental exposure or genetic mutation alone may not result in a defect, factors in combination increase phenotypic variability resulting in more individuals crossing the disease threshold where birth defects manifest. Many environmental factors that contribute to birth defects induce widespread cellular stress and misfolding of proteins. Yet, the impact of the stress response on the threshold model is not typically considered in discephering the etiology of birth defects. This mini-review will explore a potential mechanism for gene-environment interactions co-opted from studies of evolution. This model stipulates that heat shock proteins that mediate the stress response induced by environmental factors can influence the number of individuals that cross disease thresholds resulting in increased incidence of birth defects. Studies in the field of evolutionary biology have demonstrated that heat shock proteins and Hsp90 in particular provide a link between environmental stress, genotype and phenotype. Hsp90 is a highly expressed molecular chaperone that assists a wide variety of protein clients with folding and conformational changes needed for proper function. Hsp90 also chaperones client proteins with potentially deleterious amino acid changes to suppress variation caused by genetic mutations. However, upon exposure to stress, Hsp90 abandons its normal physiological clients and is diverted to assist with the misfolded protein response. This can impact the activity of signaling pathways that involve Hsp90 clients as well as unmask suppressed protein variation, essentially creating complex traits in a single step. In this capacity Hsp90 acts as an evolutionary capacitor allowing stored variation to accumulate and then become expressed in times of stress. This mechanism provides a substrate which natural selection can act upon at the population level allowing survival of the species with selective pressure. However, at the level of the individual, this mechanism can result in simultaneous expression of deleterious variants as well as reduced activity of a variety of Hsp90 chaperoned pathways, potentially shifting phenotypic variability over the disease threshold resulting in birth defects.

Tempos and modes of collectivity in the history of life

Collective integration and processing of information have increased through the history of life, through both the formation of aggregates in which the entities may have very different properties and which jointly coarse-grained environmental variables (ranging from widely varying metabolism in microbial consortia to the ecological diversity of species on reefs) and through collectives of similar entities (such as cells within an organism or social groups). Such increases have been implicated in significant transitions in the history of life, including aspects of the origin of life, the generation of pangenomes among microbes and microbial communities such as stromatolites, multicellularity and social insects. This contribution provides a preliminary overview of the dominant modes of collective information processing in the history of life, their phylogenetic distribution and extent of convergence, and the effects of new modes for integrating and acting upon information on the tempo of evolutionary change.

Testing the evolutionary constraints of metamorphosis: The ontogeny of head shape in Triturus newts

In vertebrates with complex, biphasic, life cycles, larvae have a distinct morphology and ecological preferences compared to metamorphosed juveniles and adults. In amphibians, abrupt and rapid metamorphic changes transform aquatic larvae to terrestrial juveniles. The main aim of this study is to test whether, relative to larval stages, metamorphosis (1) resets the pattern of variation between ontogenetic stages and species, (2) constrains intraspecific morphological variability, and (3) similar to the "hour-glass" model reduces morphological disparity. We explore postembryonic ontogenetic trajectories of head shape (from hatching to completed metamorphosis) of two well-defined, morphologically distinct Triturus newts species and their F1 hybrids. Variation in head shape is quantified and compared on two levels: dynamic (across ontogenetic stages) and static (at a particular stage). Our results show that the ontogenetic trajectories diverge early during development and continue to diverge throughout larval stages and metamorphosis. The high within-group variance and the largest disparity level (between-group variance) characterize the metamorphosed stage. Hence, our results indicate that metamorphosis does not canalize head shape variation generated during larval development and that metamorphosed phenotype is not more constrained relative to larval ones. Therefore, metamorphosis cannot be regarded as a developmental constraint, at least not for salamander head shape.

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.

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.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.

Decanalizing thinking on genetic canalization

The concept of genetic canalization has had an abiding influence on views of complex-trait evolution. A genetically canalized system has evolved to become less sensitive to the effects of mutation. When a gene product that supports canalization is compromised, the phenotypic impacts of a mutation should be more pronounced. This expected increase in mutational effects not only has important consequences for evolution, but has also motivated strategies to treat disease. However, recent studies demonstrate that, when putative agents of genetic canalization are impaired, systems do not behave as expected. Here, we review the evidence that is used to infer whether particular gene products are agents of genetic canalization. Then we explain how such inferences often succumb to a converse error. We go on to show that several candidate agents of genetic canalization increase the phenotypic impacts of some mutations while decreasing the phenotypic impacts of others. These observations suggest that whether a gene product acts as a 'buffer' (lessening mutational effects) or a 'potentiator' (increasing mutational effects) is not a fixed property of the gene product but instead differs for the different mutations with which it interacts. To investigate features of genetic interactions that might predispose them toward buffering versus potentiation, we explore simulated gene-regulatory networks. Similarly to putative agents of genetic canalization, the gene products in simulated networks also modify the phenotypic effects of mutations in other genes without a strong overall tendency towards lessening or increasing these effects. In sum, these observations call into question whether complex traits have evolved to become less sensitive (i.e., are canalized) to genetic change, and the degree to which trends exist that predict how one genetic change might alter another's impact. We conclude by discussing approaches to address these and other open questions that are brought into focus by re-thinking genetic canalization.