Canalization: a molecular genetic perspective

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.

Robust linear DNA degradation supports replication-initiation-defective mutants in Escherichia coli

RecBCD helicase/nuclease supports replication fork progress via recombinational repair or linear DNA degradation, explaining recBC mutant synthetic lethality with replication elongation defects. Since replication initiation defects leave chromosomes without replication forks, these should be insensitive to the recBCD status. Surprisingly, we found that both Escherichia coli dnaA46(Ts) and dnaC2(Ts) initiation mutants at semi-permissive temperatures are also recBC-colethal. Interestingly, dnaA46 recBC lethality suppressors suggest underinitiation as the problem, while dnaC2 recBC suppressors signal overintiation. Using genetic and physical approaches, we studied the dnaA46 recBC synthetic lethality, for the possibility that RecBCD participates in replication initiation. Overproduced DnaA46 mutant protein interferes with growth of dnaA+ cells, while the residual viability of the dnaA46 recBC mutant depends on the auxiliary replicative helicase Rep, suggesting replication fork inhibition by the DnaA46 mutant protein. The dnaA46 mutant depends on linear DNA degradation by RecBCD, rather than on recombinational repair. At the same time, the dnaA46 defect also interacts with Holliday junction-moving defects, suggesting reversal of inhibited forks. However, in contrast to all known recBC-colethals, which fragment their chromosomes, the dnaA46 recBC mutant develops no chromosome fragmentation, indicating that its inhibited replication forks are stable. Physical measurements confirm replication inhibition in the dnaA46 mutant shifted to semi-permissive temperatures, both at the level of elongation and initiation, while RecBCD gradually restores elongation and then initiation. We propose that RecBCD-catalyzed resetting of inhibited replication forks allows replication to displace the "sticky" DnaA46(Ts) protein from the chromosomal DNA, mustering enough DnaA for new initiations.

Leveraging genetic discoveries for sleep to determine causal relationships with common complex traits

Sleep occurs universally and is a biological necessity for human functioning. The consequences of diminished sleep quality impact physical and physiological systems such as neurological, cardiovascular, and metabolic processes. In fact, people impacted by common complex diseases experience a wide range of sleep disturbances. It is challenging to uncover the underlying molecular mechanisms responsible for decreased sleep quality in many disease systems owing to the lack of suitable sleep biomarkers. However, the discovery of a genetic component to sleep patterns has opened a new opportunity to examine and understand the involvement of sleep in many disease states. It is now possible to use major genomic resources and technologies to uncover genetic contributions to many common diseases. Large scale prospective studies such as the genome wide association studies (GWAS) have successfully revealed many robust genetic signals associated with sleep-related traits. With the discovery of these genetic variants, a major objective of the community has been to investigate whether sleep-related traits are associated with disease pathogenesis and other health complications. Mendelian Randomization (MR) represents an analytical method that leverages genetic loci as proxy indicators to establish causal effect between sleep traits and disease outcomes. Given such variants are randomly inherited at birth, confounding bias is eliminated with MR analysis, thus demonstrating evidence of causal relationships that can be used for drug development and to prioritize clinical trials. In this review, we outline the results of MR analyses performed to date on sleep traits in relation to a multitude of common complex diseases.

Morphological canalization, integration, and plasticity in response to population density in Abutilon theophrasti: Influences of soil conditions and growth stages

Phenotypic integration and developmental canalization have been hypothesized to constrain the degree of phenotypic plasticity, but little evidence exists, probably due to the lack of studies on the relationships among the three processes, especially for plants under different environments. We conducted a field experiment by subjecting plants of Abutilon theophrasti to three densities, under infertile and fertile soil conditions, and analyzing correlations among canalization, integration, and plasticity in a variety of measured morphological traits after 50 and 70 days, to investigate the relationships among the three variables in response to density and how these responses vary with soil conditions and growth stages. Results showed trait canalization decreased and phenotypic integration and the degree of plasticity (absolute plasticity) in traits increased with density. Phenotypic integration often positively correlated with absolute plasticity, whereas correlations between trait canalization and plasticity were insignificant in most cases, with a few positive ones between canalization and absolute plasticity at low and medium densities. As plants grew, these correlations intensified in infertile soil and attenuated in fertile soil. Our findings suggested the complexity of the relationship between canalization and plasticity: Decreased canalization is more likely to facilitate active plastic responses under more favorable conditions, whereas increased level of integration should mainly be an outcome of plastic responses. Soil conditions and growth stage may affect responses of these correlations to density via modifying plant size, competition strength, and plastic responses in traits. We also predicted that decreased canalization can be advantageous or disadvantageous, and the lack of response to stress may demonstrate a stronger ability of adaptation than passive response, thus should be adaptive plasticity as active response.

Size changes in miR‑21 knockout mice: Geometric morphometrics on teeth, alveolar bone and mandible

MicroRNA‑21 (miR‑21) is a small non‑coding RNA that is differentially expressed during tooth development, particularly during amelogenesis. Although orthodontic tooth movement and the innate immune response are impaired, miR‑21 knockout mice demonstrate no obvious skeletal phenotype. However, the consequence of miR‑21 knockout on tooth phenotype and corresponding alveolar bone is unknown. The current study utilized landmark‑based geometric morphometrics to identify anatomical dissimilarities of the three lower and upper molars, and the corresponding alveolar bone, in miR‑21 knockout and wild‑type control mice. The anatomical structures were visualized by microcomputer tomography. A total of 36 and 38 landmarks were placed on mandibular and maxillary molars, respectively. For the alveolar bone, 16 landmarks were selected on both anatomical sites. General Procrustes analysis revealed significantly smaller molars and dimensions of the alveolar bone in the mandible of the miR‑21 knockout mice when compared with wild‑type controls (P=0.03 and P=0.04, respectively). The overall dimension of the mandible was reduced by the lack of miR‑21 (P=0.02). In the maxilla, the dimension of the alveolar bone was significant (P=0.02); however, this was not observed in the molars (P=0.36). Based on principal component analysis, no changes in shape for any of the anatomical sites were observed. Dental and skeletal jaw length were calculated and no prognathism was identified. However, the fluctuating asymmetry of the molars in the mandible and the maxilla was reduced in the miR‑21 knockout mice by 38 and 27%, respectively. Taken together, the results of the present study revealed that the molars in the mandible and the dimension of the respective alveolar bone were smaller in miR‑21 mice compared with wild‑type littermates, suggesting that miR‑21 influences tooth development.

Decreased sensitivity to the anorectic effects of leptin in mice that lack a Pomc-specific neural enhancer

Enhancer redundancy has been postulated to provide a buffer for gene expression against genetic and environmental perturbations. While work in Drosophila has identified functionally overlapping enhancers, work in mammalian models has been limited. Recently, we have identified two partially redundant enhancers, nPE1 and nPE2, that drive proopiomelanocortin gene expression in the hypothalamus. Here we demonstrate that deletion of nPE1 produces mild obesity while knockout of nPE2 has no discernible metabolic phenotypes. Additionally, we show that acute leptin administration has significant effects on nPE1 knockout mice, with food intake and body weight change significantly impacted by peripheral leptin treatment. nPE1 knockout mice became less responsive to leptin treatment over time as percent body weight change increased over 2 week exposure to peripheral leptin. Both Pomc and Agrp mRNA were not differentially affected by chronic leptin treatment however we did see a decrease in Pomc and Agrp mRNA in both nPE1 and nPE2 knockout calorie restricted mice as compared to calorie restricted PBS-treated WT mice. Collectively, these data suggest dynamic regulation of Pomc by nPE1 such that mice with nPE1 knockout become less responsive to the anorectic effects of leptin treatment over time. Our results also support our earlier findings in which nPE2 may only be critical in adult mice that lack nPE1, indicating that these neural enhancers work synergistically to influence metabolism.

Integrative omics analysis reveals relationships of genes with synthetic lethal interactions through a pan-cancer analysis

Synthetic lethality is thought to play an important role in anticancer therapies. Herein, to understand the potential distributions and relationships between synthetic lethal interactions between genes, especially for pairs deriving from different sources, we performed an integrative analysis of genes at multiple molecular levels. Based on inter-species phylogenetic conservation of synthetic lethal interactions, gene pairs from yeast and humans were analyzed; a total of 37,588 candidate gene pairs containing 7,816 genes were collected. Of these, 49.74% of genes had 2–10 interactions, 22.93% were involved in hallmarks of cancer, and 21.61% were identified as core essential genes. Many genes were shown to have important biological roles via functional enrichment analysis, and 65 were identified as potentially crucial in the pathophysiology of cancer. Gene pairs with dysregulated expression patterns had higher prognostic values. Further screening based on mutation and expression levels showed that remaining gene pairs were mainly derived from human predicted or validated pairs, while most predicted pairs from yeast were filtered from analysis. Genes with synthetic lethality were further analyzed with their interactive microRNAs (miRNAs) at the isomiR level which have been widely studied as negatively regulatory molecules. The miRNA–mRNA interaction network revealed that many synthetic lethal genes contributed to the cell cycle (seven of 12 genes), cancer pathways (five of 12 genes), oocyte meiosis, the p53 signaling pathway, and hallmarks of cancer. Our study contributes to the understanding of synthetic lethal interactions and promotes the application of genetic interactions in further cancer precision medicine.

Anatomical similarity between the Sost-knockout mouse and sclerosteosis in humans

Sclerosteosis, a rare autosomal recessive genetic disorder caused by a mutation of the Sost gene, manifests in the facial skeleton by gigantism, facial distortion, mandibular prognathism, cranial nerve palsy, and, in extreme cases, compression of the medulla oblongata. Mice lacking sclerostin reflect some symptoms of sclerosteosis, but this is the first report of the effect on the facial skeleton. We used geometric morphometrics (GMM) to analyze the deformations of the murine facial skeleton from the wild‐type to the Sost gene knockout. Landmark coordinates were obtained by surface reconstructions from micro‐computed tomography. Centroid size, principal component scores in shape space and form space, and asymmetry were computed by the standard GMM formulas, and dental and skeletal jaw lengths were examined as ratios. We show here that, compared to wild type controls, mice lacking Sost have larger centroid size (effect size, p‐value: 4.59, <.001), higher mean asymmetry (1.14, .065), dental and skeletal mandibular prognathism (1.36, .010 and 5.92, <.001), a smaller foramen magnum (−1.71, .015), and calvaria that are more highly curved (form space p = 4.09, .002; shape space p = 12.82, .002). These features of mice lacking sclerostin largely correspond to the changes of the facial skeleton observed in sclerosteosis. This alignment further supports claims that the Sost gene plays a fundamental role in bony facial development in rodents and humans alike.

Covariation of the endocranium and splanchnocranium during great ape ontogeny

That great ape endocranial shape development persists into adolescence indicates that the splanchnocranium succeeds brain growth in driving endocranial development. However, the extent of this splanchnocranial influence is unknown. We applied two-block partial least squares analyses of Procrustes shape variables on an ontogenetic series of great ape crania to explore the covariation of the endocranium (the internal braincase) and splanchnocranium (face, or viscerocranium). We hypothesized that a transition between brain growth and splanchnocranial development in the establishment of final endocranial form would be manifest as a change in the pattern of shape covariation between early and adolescent ontogeny. Our results revealed a strong pattern of covariation between endocranium and splanchnocranium, indicating that chimpanzees, gorillas, and orangutans share a common tempo and mode of morphological integration from the eruption of the deciduous dentition onwards to adulthood: a reflection of elongating endocranial shape and continuing splanchnocranial prognathism. Within this overarching pattern, we noted that species variation exists in magnitude and direction, and that the covariation between the splanchnocranium and endocranium is somewhat weaker in early infancy compared to successive age groups. When correcting our covariation analyses for allometry, we found that an ontogenetic signal remains, signifying that allometric variation alone is insufficient to account for all endocranial-splanchnocranial developmental integration. Finally, we assessed the influence of the cranial base, which acts as the interface between the face and endocranium, on the shape of the vault using thin-plate spline warping. We found that not all splanchnocranial shape changes during development are tightly integrated with endocranial shape. This suggests that while the developmental expansion of the brain is the main driver of endocranial shape during early ontogeny, endocranial development from infancy onwards is moulded by the splanchnocranium in conjunction with the neurocranium.

Ankyrin domain encoding genes from an ancient horizontal transfer are functionally integrated into Nasonia developmental gene regulatory networks

BACKGROUND

How regulatory networks incorporate additional components and how novel genes are functionally integrated into well-established developmental processes are two important and intertwined questions whose answers have major implications for understanding the evolution of development. We recently discovered a set of lineage-restricted genes with strong and specific expression patterns along the dorsal-ventral (DV) axis of the embryo of the wasp Nasonia that may serve as a powerful system for addressing these questions. We sought to both understand the evolutionary history of these genes and to determine their functions in the Nasonia DV patterning system.

RESULTS

We have found that the novel DV genes are part of a large family of rapidly duplicating and diverging ankyrin domain-encoding genes that originated most likely by horizontal transfer from a prokaryote in a common ancestor of the wasp superfamily Chalcidoidea. We tested the function of those ankyrin-encoding genes expressed along the DV axis and found that they participate in early embryonic DV patterning. We also developed a new wasp model system (Melittobia) and found that some functional integration of ankyrin genes have been preserved for over 90 million years.

CONCLUSIONS

Our results indicate that regulatory networks can incorporate novel genes that then become necessary for stable and repeatable outputs. Even a modest role in developmental networks may be enough to allow novel or duplicate genes to be maintained in the genome and become fully integrated network components.

Canalization and genetic assimilation: Reassessing the radicality of the Waddingtonian concept of inheritance of acquired characters

Genetic assimilation is often mixed up with the Baldwin effect. For Waddington, genetic assimilation was both a phenomenon and a specific mechanism of adaptive evolution which was grounded in the concept of canalization. This theoretical link between canalization and genetic assimilation, which was pivotal in Waddington's view, has been weakened since the early 1960s. The aim of the present article is to emphasize the specificity and to reassess the possible radicality of Waddington's proposal. What he claimed to have elaborated was an actual and genuine mechanism of inheritance of acquired characters that did not rely on soft Lamarckian inheritance. Consequently his "theory" of genetic assimilation, unlike the Baldwin effect, might not be as easily integrated in the framework of the Modern Synthesis.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Partially redundant enhancers cooperatively maintain Mammalian pomc expression above a critical functional threshold

Cell-specific expression of many genes is conveyed by multiple enhancers, with each individual enhancer controlling a particular expression domain. In contrast, multiple enhancers drive similar expression patterns of some genes involved in embryonic development, suggesting regulatory redundancy. Work in Drosophila has indicated that functionally overlapping enhancers canalize development by buffering gene expression against environmental and genetic disturbances. However, little is known about regulatory redundancy in vertebrates and in genes mainly expressed during adulthood. Here we study nPE1 and nPE2, two phylogenetically conserved mammalian enhancers that drive expression of the proopiomelanocortin gene (Pomc) to the same set of hypothalamic neurons. The simultaneous deletion of both enhancers abolished Pomc expression at all ages and induced a profound metabolic dysfunction including early-onset extreme obesity. Targeted inactivation of either nPE1 or nPE2 led to very low levels of Pomc expression during early embryonic development indicating that both enhancers function synergistically. In adult mice, however, Pomc expression is controlled additively by both enhancers, with nPE1 being responsible for ∼80% and nPE2 for ∼20% of Pomc transcription. Consequently, nPE1 knockout mice exhibit mild obesity whereas nPE2-deficient mice maintain a normal body weight. These results suggest that nPE2-driven Pomc expression is compensated by nPE1 at later stages of development, essentially rescuing the earlier phenotype of nPE2 deficiency. Together, these results reveal that cooperative interactions between the enhancers confer robustness of Pomc expression against gene regulatory disturbances and preclude deleterious metabolic phenotypes caused by Pomc deficiency in adulthood. Thus, our study demonstrates that enhancer redundancy can be used by genes that control adult physiology in mammals and underlines the potential significance of regulatory sequence mutations in common diseases.

On the Nature and Evolutionary Impact of Phenotypic Robustness Mechanisms

Biologists have long observed that physiological and developmental processes are insensitive, or robust, to many genetic and environmental perturbations. A complete understanding of the evolutionary causes and consequences of this robustness is lacking. Recent progress has been made in uncovering the regulatory mechanisms that underlie environmental robustness in particular. Less is known about robustness to the effects of mutations, and indeed the evolution of mutational robustness remains a controversial topic. The controversy has spread to related topics, in particular the evolutionary relevance of cryptic genetic variation. This review aims to synthesize current understanding of robustness mechanisms and to cut through the controversy by shedding light on what is and is not known about mutational robustness. Some studies have confused mutational robustness with non-additive interactions between mutations (epistasis). We conclude that a profitable way forward is to focus investigations (and rhetoric) less on mutational robustness and more on epistasis.

Conversion of the thymus into a bipotent lymphoid organ by replacement of FOXN1 with its paralog, FOXN4

The thymus is a lymphoid organ unique to vertebrates, and it provides a unique microenvironment that facilitates the differentiation of immature hematopoietic precursors into mature T cells. We subjected the evolutionary trajectory of the thymic microenvironment to experimental analysis. A hypothetical primordial form of the thymus was established in mice by replacing FOXN1, the vertebrate-specific master regulator of thymic epithelial cell function, with its metazoan ancestor, FOXN4, thereby resetting the regulatory and coding changes that have occurred since the divergence of these two paralogs. FOXN4 exhibited substantial thymopoietic activity. Unexpectedly, histological changes and a functional imbalance between the lymphopoietic cytokine IL7 and the T cell specification factor DLL4 within the reconstructed thymus resulted in coincident but spatially segregated T and B cell development. Our results identify an evolutionary mechanism underlying the conversion of a general lymphopoietic organ to a site of exclusive T cell generation.

Contribution of nonohnologous duplicated genes to high habitat variability in mammals

The mechanism by which genetic systems affect environmental adaptation is a focus of considerable attention in the fields of ecology, evolution, and conservation. However, the genomic characteristics that constrain adaptive evolution have remained unknown. A recent study showed that the proportion of duplicated genes in whole Drosophila genomes correlated with environmental variability within habitat, but it remains unclear whether the correlation is observed even in vertebrates whose genomes including a large number of duplicated genes generated by whole-genome duplication (WGD). Here, we focus on fully sequenced mammalian genomes that experienced WGD in early vertebrate lineages and show that the proportion of small-scale duplication (SSD) genes in the genome, but not that of WGD genes, is significantly correlated with habitat variability. Moreover, species with low habitat variability have a higher proportion of lost duplicated genes, particularly SSD genes, than those with high habitat variability. These results indicate that species that inhabit variable environments may maintain more SSD genes in their genomes and suggest that SSD genes are important for adapting to novel environments and surviving environmental changes. These insights may be applied to predicting invasive and endangered species.

Role of mef2ca in developmental buffering of the zebrafish larval hyoid dermal skeleton

Phenotypic robustness requires a process of developmental buffering that is largely not understood, but which can be disrupted by mutations. Here we show that in mef2ca(b1086) loss of function mutant embryos and early larvae, development of craniofacial hyoid bones, the opercle (Op) and branchiostegal ray (BR), becomes remarkably unstable; the large magnitude of the instability serves as a positive attribute to learn about features of this developmental buffering. The OpBR mutant phenotype variably includes bone expansion and fusion, Op duplication, and BR homeosis. Formation of a novel bone strut, or a bone bridge connecting the Op and BR together occurs frequently. We find no evidence that the phenotypic stability in the wild type is provided by redundancy between mef2ca and its co-ortholog mef2cb, or that it is related to the selector (homeotic) gene function of mef2ca. Changes in dorsal-ventral patterning of the hyoid arch also might not contribute to phenotypic instability in mutants. However, subsequent development of the bone lineage itself, including osteoblast differentiation and morphogenetic outgrowth, shows marked variation. Hence, steps along the developmental trajectory appear differentially sensitive to the loss of buffering, providing focus for the future study.

Lines of least resistance and genetic architecture of house mouse (Mus musculus) mandible shape

Evolution along "lines of least resistance (LLR)" is a well-established phenomenon, which is ultimately based on variational properties of the developmental system. However, it has remained unclear which aspects of development are responsible for observed distributions of variation. This article suggests that house mouse mandible shape may present an opportunity for future research in this field. The genetic architecture of mouse mandible shape has been investigated in the laboratory, and its evolution has been surveyed in studies on natural variation. Here, I ask whether evolutionary diversification of mandible shape follows "LLR" by comparing principal directions of genetic and evolutionary variation, and I assess the potential contribution of specific genomic regions to evolutionary divergence along these directions. The role of two aspects of development, QTL number (presumably reflecting mutational target size) and canalization, is comparatively assessed. Overall, our results demonstrate a strong tendency of genetic and evolutionary systems to vary along similar directions, suggesting that mandible shape evolves along LLR at the level of populations, subspecies and species. At the level of genetic factors, effects bearing similarity to directions of evolution are significantly overrepresented, providing support for the mutational target size hypothesis. A role for canalization is not supported by a clear correlation between directions of evolution and size of genetic effects; however, the evidence for canalization remains ambiguous. These results provide some insights into how the developmental system may shape the variational properties of genetic systems and thus influence the direction of evolution.

The evolution of intelligent developmental systems

This chapter aims to understand the relations between the evolution and development of complex cognitive functions by emphasizing the context of complex, changeable environments. What evolves and develops in such contexts cannot be achieved by linear deterministic processes based on stable "codes". Rather, what is needed, even in the molecular ensembles of single-cell organisms, are "intelligent" systems with nonlinear dynamic processing, sensitive to informational structures, not just elements, in environments. This is the view emerging in recent molecular biology. The research is also constructing a new "biologic" of both evolution and development, providing a clearer rationale for transitions into more complex forms, including epigenetic, physiological, nervous, cognitive, and human sociocognitive forms. This chapter explains how these transitions form a nested hierarchical system in which the dynamics within and between levels creates emergent abilities so often underestimated or even demeaned in previous accounts, especially regarding human cognition. The implications of the view for human development in modern societies are also briefly considered.

Genetic redundancy strengthens the circadian clock leading to a narrow entrainment range

Circadian clocks are internal timekeepers present in almost all organisms. Driven by a genetic network of highly conserved structure, they generate self-sustained oscillations that entrain to periodic external signals such as the 24 h light-dark cycle. Vertebrates possess multiple, functionally overlapping homologues of the core clock genes. Furthermore, vertebrate clocks entrain to a range of periods three times as narrow as that of other organisms. We asked whether genetic redundancies play a role in governing entrainment properties and analysed locomotor activity rhythms of genetically modified mice lacking one set of clock homologues. Exposing them to non-24 h light-dark cycles, we found that the mutant mice have a wider entrainment range than the wild types. Spectral analysis furthermore revealed nonlinear phenomena of periodically forced self-sustained oscillators for which the entrainment range relates inversely to oscillator amplitude. Using the forced oscillator model to explain the observed differences in entrainment range between mutant and wild-type mice, we sought to quantify the overall oscillator amplitude of their clocks from the activity rhythms and found that mutant mice have weaker circadian clocks than wild types. Our results suggest that genetic redundancy strengthens the circadian clock leading to a narrow entrainment range in vertebrates.

Duplication and retention biases of essential and non-essential genes revealed by systematic knockdown analyses

When a duplicate gene has no apparent loss-of-function phenotype, it is commonly considered that the phenotype has been masked as a result of functional redundancy with the remaining paralog. This is supported by indirect evidence showing that multi-copy genes show loss-of-function phenotypes less often than single-copy genes and by direct tests of phenotype masking using select gene sets. Here we take a systematic genome-wide RNA interference approach to assess phenotype masking in paralog pairs in the Caenorhabditis elegans genome. Remarkably, in contrast to expectations, we find that phenotype masking makes only a minor contribution to the low knockdown phenotype rate for duplicate genes. Instead, we find that non-essential genes are highly over-represented among duplicates, leading to a low observed loss-of-function phenotype rate. We further find that duplicate pairs derived from essential and non-essential genes have contrasting evolutionary dynamics: whereas non-essential genes are both more often successfully duplicated (fixed) and lost, essential genes are less often duplicated but upon successful duplication are maintained over longer periods. We expect the fundamental evolutionary duplication dynamics presented here to be broadly applicable.

Duplicate genes occur in all organisms. It has been found that mutations in duplicate genes cause defects much less often than when single copy genes are mutated. It is widely believed that this is due to functional redundancy—that is, the two genes can carry out similar functions so that the non-mutated duplicate gene can cover for or “mask” the phenotype of the mutation in the first duplicate. To determine whether this hypothesis is true, it is necessary to test systematically whether defects indeed occur in the organism when both duplicate genes are inhibited. We have for the first time carried out such an analysis in a multicellular organism, the nematode Caenorhabditis elegans. In contrast to expectations, we observed that when both copies of duplicate genes are inhibited deleterious effects are very rare. We show that this is because duplicate genes are much more often non-essential compared to genes where there is only a single copy. Non-essential genes are also lost from the genome much more often than essential genes. However, when essential genes are duplicated, they remain present in the genome over longer periods. Our results give a framework to explain the evolutionary dynamics of duplications in the genome.

Gene loss, thermogenesis, and the origin of birds

Compared to related taxa, birds have exceptionally enlarged and diversified skeletal muscles, features that are closely associated with skeletal diversification and are commonly explained by a diversity of avian ecological niches and locomotion types. The thermogenic muscle hypothesis (TMH) for the origin of birds proposes that such muscle hyperplasia and the associated skeletal innovations are instead the consequence of the avian clade originating from an ancestral population that underwent several successive episodes of loss of genes associated with thermogenesis, myogenesis, and skeletogenesis. Direct bird ancestors met this challenge with a combination of behavioral strategies (e.g., brooding of nestlings) and acquisition of a variety of adaptations for enhanced nonshivering thermogenesis in skeletal muscle. The latter include specific biochemical alterations promoting muscle heat generation and dramatic expansion of thigh and breast muscle mass. The TMH proposes that such muscle hyperplasia facilitated bipedality, freeing upper limbs for new functions (e.g., flight, swimming), and, by altering the mechanical environment of embryonic development, generated skeletal novelties, sometimes abruptly, that became distinctive features of the avian body plan.

Habitat variability correlates with duplicate content of Drosophila genomes

The factors limiting the habitat range of species are crucial in understanding their biodiversity and response to environmental change. Yet the genetic and genomic architectures that produce genetic variation to enable environmental adaptation have remained poorly understood. Here we show that the proportion of duplicated genes (P_D) in the whole genomes of fully sequenced Drosophila species is significantly correlated with environmental variability within the habitats measured by the climatic envelope and habitat diversity. Furthermore, species with a low P_D tend to lose the duplicated genes owing to their faster evolution. These results indicate that the rapid relaxation of functional constraints on duplicated genes resulted in a low P_D for species with lower habitat diversity, and suggest that the maintenance of duplicated genes gives organisms an ecological advantage during evolution. We therefore propose that the P_D in a genome is related to adaptation to environmental variation.

A symmetric dual feedback system provides a robust and entrainable oscillator

Many organisms have evolved molecular clocks to anticipate daily changes in their environment. The molecular mechanisms by which the circadian clock network produces sustained cycles have extensively been studied and transcriptional-translational feedback loops are common structures to many organisms. Although a simple or single feedback loop is sufficient for sustained oscillations, circadian clocks implement multiple, complicated feedback loops. In general, different types of feedback loops are suggested to affect the robustness and entrainment of circadian rhythms.

To reveal the mechanism by which such a complex feedback system evolves, we quantify the robustness and light entrainment of four competing models: the single, semi-dual, dual, and redundant feedback models. To extract the global properties of those models, all plausible kinetic parameter sets that generate circadian oscillations are searched to characterize their oscillatory features. To efficiently perform such analyses, we used the two-phase search (TPS) method as a fast and non-biased search method and quasi-multiparameter sensitivity (QMPS) as a fast and exact measure of robustness to uncertainty of all kinetic parameters.

So far the redundant feedback model has been regarded as the most robust oscillator, but our extensive analysis corrects or overcomes this hypothesis. The dual feedback model, which is employed in biology, provides the most robust oscillator to multiple parameter perturbations within a cell and most readily entrains to a wide range of light-dark cycles. The kinetic symmetry between the dual loops and their coupling via a protein complex are found critically responsible for robust and entrainable oscillations. We first demonstrate how the dual feedback architecture with kinetic symmetry evolves out of many competing feedback systems.

Use of genetic markers and gene-diet interactions for interrogating population-level causal influences of diet on health

Differences in diet appear to contribute substantially to the burden of disease in populations, and therefore changes in diet could lead to major improvements in public health. This is predicated on the reliable identification of causal effects of nutrition on health, and unfortunately nutritional epidemiology has deficiencies in terms of identifying these. This is reflected in the many cases where observational studies have suggested that a nutritional factor is protective against disease, and randomized controlled trials have failed to verify this. The use of genetic variants as proxy measures of nutritional exposure-an application of the Mendelian randomization principle-can contribute to strengthening causal inference in this field. Genetic variants are not subject to bias due to reverse causation (disease processes influencing exposure, rather than vice versa) or recall bias, and if obvious precautions are applied are not influenced by confounding or attenuation by errors. This is illustrated in the case of epidemiological studies of alcohol intake and various health outcomes, through the use of genetic variants related to alcohol metabolism (in ALDH2 and ADH1B). Examples from other areas of nutritional epidemiology and of the informative nature of gene-environment interactions interpreted within the Mendelian randomization framework are presented, and the potential limitations of the approach addressed.

Gene regulation: enhancers stepping out of the shadow

The expression of many animal genes has been shown to be controlled by two--rather than one--enhancers with similar regulatory content. Such enhancer redundancy ensures robustness of gene expression under adverse environmental or genetic conditions.

Mendelian Randomization for Strengthening Causal Inference in Observational Studies: Application to Gene × Environment Interactions

Identification of environmentally modifiable factors causally influencing disease risk is fundamental to public-health improvement strategies. Unfortunately, observational epidemiological studies are limited in their ability to reliably identify such causal associations, reflected in the many cases in which conventional epidemiological studies have apparently identified associations that randomized controlled trials have failed to verify. The use of genetic variants as proxy measures of exposure -an application of the Mendelian randomization principle-can contribute to strengthening causal inference. Genetic variants are not subject to bias due to reverse causation (disease processes influencing exposure, rather than vice versa) or recall bias, and if simple precautions are applied, they are not influenced by confounding or attenuation by errors. The principles of Mendelian randomization are illustrated with specific reference to studies of the effects of alcohol intake on various health-related outcomes through the utilization of genetic variants related to alcohol metabolism (in ALDH2 and ADH1B). Ways of incorporating Gene × Environment interactions into the Mendelian randomization framework are developed, and the strengths and limitations of the approach discussed.

Buffering by gene duplicates: an analysis of molecular correlates and evolutionary conservation

Background

One mechanism to account for robustness against gene knockouts or knockdowns is through buffering by gene duplicates, but the extent and general correlates of this process in organisms is still a matter of debate. To reveal general trends of this process, we provide a comprehensive comparison of gene essentiality, duplication and buffering by duplicates across seven bacteria (Mycoplasma genitalium, Bacillus subtilis, Helicobacter pylori, Haemophilus influenzae, Mycobacterium tuberculosis, Pseudomonas aeruginosa, Escherichia coli), and four eukaryotes (Saccharomyces cerevisiae (yeast), Caenorhabditis elegans (worm), Drosophila melanogaster (fly), Mus musculus (mouse)).

Results

In nine of the eleven organisms, duplicates significantly increase chances of survival upon gene deletion (P-value ≤ 0.05), but only by up to 13%. Given that duplicates make up to 80% of eukaryotic genomes, the small contribution is surprising and points to dominant roles of other buffering processes, such as alternative metabolic pathways. The buffering capacity of duplicates appears to be independent of the degree of gene essentiality and tends to be higher for genes with high expression levels. For example, buffering capacity increases to 23% amongst highly expressed genes in E. coli. Sequence similarity and the number of duplicates per gene are weak predictors of the duplicate's buffering capacity. In a case study we show that buffering gene duplicates in yeast and worm are somewhat more similar in their functions than non-buffering duplicates and have increased transcriptional and translational activity.

Conclusion

In sum, the extent of gene essentiality and buffering by duplicates is not conserved across organisms and does not correlate with the organisms' apparent complexity. This heterogeneity goes beyond what would be expected from differences in experimental approaches alone. Buffering by duplicates contributes to robustness in several organisms, but to a small extent – and the relatively large amount of buffering by duplicates observed in yeast and worm may be largely specific to these organisms. Thus, the only common factor of buffering by duplicates between different organisms may be the by-product of duplicate retention due to demands of high dosage.

Gene duplications, robustness and evolutionary innovations

Mutational robustness facilitates evolutionary innovations. Gene duplications are unique kinds of mutations, in that they generally increase such robustness. The frequent association of gene duplications in regulatory networks with evolutionary innovation is thus a special case of a general mechanism linking innovation to robustness. The potential power of this mechanism to promote evolutionary innovations on large time scales is illustrated here with several examples. These include the role of gene duplications in the vertebrate radiation, flowering plant evolution and heart development, which encompass some of the most striking innovations in the evolution of life.

On the origins of novelty in development and evolution

The origin of novel traits is what draws many to evolutionary biology, yet our understanding of the mechanisms that underlie the genesis of novelty remains limited. Here I review definitions of novelty including its relationship to homology. I then discuss how ontogenetic perspectives may allow us to move beyond current roadblocks in our understanding of the mechanics of innovation. Specifically, I explore the roles of canalization, plasticity and threshold responses during development in generating a reservoir of cryptic genetic variation free to drift and accumulate in natural populations. Environmental or genetic perturbations that exceed the buffering capacity of development can then release this variation, and, through evolution by genetic accommodation, result in rapid diversification, recurrence of lost phenotypes as well as the origins of novel features. I conclude that, in our quest to understand the nature of innovation, the nature of development deserves to take center stage.

Strengthening causal inference in cardiovascular epidemiology through Mendelian randomization

Observational studies have contributed in a major way to understanding modifiable determinants of cardiovascular disease risk, but several examples exist of factors that were identified in observational studies as potentially protecting against coronary heart disease, that in randomized controlled trials had no such effect. The likely reason for misleading findings from observational epidemiological studies is that associations are influenced by confounding, bias, and reverse causation--where disease influences a risk factor, rather than vice versa. Mendelian randomization utilizes genetic variants that serve as proxy measures for modifiable risk factors to allow estimation of the causal influence of the modifiable risk factor in question. We present examples of the use of the Mendelian randomization approach and discuss both the limitations and potentials of this strategy.

Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology

Observational epidemiological studies suffer from many potential biases, from confounding and from reverse causation, and this limits their ability to robustly identify causal associations. Several high-profile situations exist in which randomized controlled trials of precisely the same intervention that has been examined in observational studies have produced markedly different findings. In other observational sciences, the use of instrumental variable (IV) approaches has been one approach to strengthening causal inferences in non-experimental situations. The use of germline genetic variants that proxy for environmentally modifiable exposures as instruments for these exposures is one form of IV analysis that can be implemented within observational epidemiological studies. The method has been referred to as 'Mendelian randomization', and can be considered as analogous to randomized controlled trials. This paper outlines Mendelian randomization, draws parallels with IV methods, provides examples of implementation of the approach and discusses limitations of the approach and some methods for dealing with these.

Recombination and the evolution of mutational robustness

Mutational robustness is the degree to which a phenotype, such as fitness, is resistant to mutational perturbations. Since most of these perturbations will tend to reduce fitness, robustness provides an immediate benefit for the mutated individual. However, robust systems decay due to the accumulation of deleterious mutations that would otherwise have been cleared by selection. This decay has received very little theoretical attention. At equilibrium, a population or asexual lineage is expected to have a mutation load that is invariant with respect to the selection coefficient of deleterious alleles, so the benefit of robustness (at the level of the population or asexual lineage) is temporary. However, previous work has shown that robustness can be favoured when robustness loci segregate independently of the mutating loci they act upon. We examine a simple two-locus model that allows for intermediate rates of recombination and inbreeding to show that increasing the effective recombination rate allows for the evolution of greater mutational robustness.

The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration

Although it is well known that many mutations influence phenotypic variability as well as the mean, the underlying mechanisms for variability effects are very poorly understood. The brachymorph (bm) phenotype results from an autosomal recessive mutation in the phosphoadenosine-phosphosulfate synthetase 2 gene (Papps2). A major cranial manifestation is a dramatic reduction in the growth of the chondrocranium which results from undersulfation of glycosaminoglycans (GAGs) in the cartilage matrix. We found that this reduction in the growth of the chondrocranium is associated with an altered pattern of craniofacial shape variation, a significant increase in phenotypic variance and a dramatic increase in morphological integration for craniofacial shape. Both effects are largest in the basicranium. The altered variation pattern indicates that the mutation produces developmental influences on shape that are not present in the wildtype. As the mutation dramatically reduces sulfation of GAGs, we infer that this influence is variation among individuals in the degree of sulfation, or variable expressivity of the mutation. This variation may be because of genetic variation at other loci that influence sulfation, environmental effects, or intrinsic effects. We infer that chondrocranial development exhibits greater sensitivity to variation in the sulfation of chondroitin sulfate when the degree of sulfation is low. At normal levels, sulfation probably contributes minimally to phenotypic variation. This case illustrates canalization in a particular developmental-genetic context.

Genetic epidemiology and public health: hope, hype, and future prospects

Genetic epidemiology is a rapidly expanding research field, but the implications of findings from such studies for individual or population health are unclear. The use of molecular genetic screening currently has some legitimacy in certain monogenic conditions, but no established value with respect to common complex diseases. Personalised medical care based on molecular genetic testing is also as yet undeveloped for common diseases. Genetic epidemiology can contribute to establishing the causal nature of environmentally modifiable risk factors, through the application of mendelian randomisation approaches and thus contribute to appropriate preventive strategies. Technological and other advances will allow the potential of genetic epidemiology to be revealed over the next few years, and the establishment of large population-based resources for such studies (biobanks) should contribute to this endeavour.

The Evolutionary Genetics of Canalization

Evolutionary genetics has recently made enormous progress in understanding how genetic variation maps into phenotypic variation. However why some traits are phenotypically invariant despite apparent genetic and environmental changes has remained a major puzzle. In the 1940s, Conrad Hal Waddington coined the concept and term "canalization" to describe the robustness of phenotypes to perturbation; a similar concept was proposed by Waddington's contemporary Ivan Ivanovich Schmalhausen. This paper reviews what has been learned about canalization since Waddington. Canalization implies that a genotype's phenotype remains relatively invariant when individuals of a particular genotype are exposed to different environments (environmental canalization) or when individuals of the same single- or multilocus genotype differ in their genetic background (genetic canalization). Consequently, genetic canalization can be viewed as a particular kind of epistasis, and environmental canalization and phenotypic plasticity are two aspects of the same phenomenon. Canalization results in the accumulation of phenotypically cryptic genetic variation, which can be released after a "decanalizing" event. Thus, canalized genotypes maintain a cryptic potential for expressing particular phenotypes, which are only uncovered under particular decanalizing environmental or genetic conditions. Selection may then act on this newly released genetic variation. The accumulation of cryptic genetic variation by canalization may therefore increase evolvability at the population level by leading to phenotypic diversification under decanalizing conditions. On the other hand, under canalizing conditions, a major part of the segregating genetic variation may remain phenotypically cryptic; canalization may therefore, at least temporarily, constrain phenotypic evolution. Mechanistically, canalization can be understood in terms of transmission patterns, such as epistasis, pleiotropy, and genotype by environment interactions, and in terms of genetic redundancy, modularity, and emergent properties of gene networks and biochemical pathways. While different forms of selection can favor canalization, the requirements for its evolution are typically rather restrictive. Although there are several methods to detect canalization, there are still serious problems with unambiguously demonstrating canalization, particularly its adaptive value.

Distributed robustness versus redundancy as causes of mutational robustness

A biological system is robust to mutations if it continues to function after genetic changes in its parts. Such robustness is pervasive on different levels of biological organization, from macromolecules to genetic networks and whole organisms. I here ask which of two possible causes of such robustness are more important on a genome-wide scale, for systems whose parts are genes, such as metabolic and genetic networks. The first of the two causes is redundancy of a system's parts: A gene may be dispensable if the genome contains redundant, back-up copies of the gene. The second cause, distributed robustness, is more poorly understood. It emerges from the distributed nature of many biological systems, where many (and different) parts contribute to system functions. I will here discuss evidence suggesting that distributed robustness is equally or more important for mutational robustness than gene redundancy. This evidence comes from the functional divergence of redundant genes, as well as from large-scale gene deletion studies. I also ask whether one can quantify the extent to which redundancy or distributed robustness contribute to mutational robustness.

Genetic analysis of an overlapping functional requirement for L1- and NCAM-type proteins during sensory axon guidance in Drosophila

L1- and NCAM-type cell adhesion molecules represent distinct protein families that function as specific receptors for different axon guidance cues. However, both L1 and NCAM proteins promote axonal growth by inducing neuronal tyrosine kinase activity and are coexpressed in subsets of axon tracts in arthropods and vertebrates. We have studied the functional requirements for the Drosophila L1- and NCAM-type proteins, Neuroglian (Nrg) and Fasciclin II (FasII), during postembryonic sensory axon guidance. The rescue of the Neuroglian loss-of-function (LOF) phenotype by transgenically expressed L1- and NCAM-type proteins demonstrates a functional interchangeability between these proteins in Drosophila photoreceptor pioneer axons, where both proteins are normally coexpressed. In contrast, the ectopic expression of Fasciclin II in mechanosensory neurons causes a strong enhancement of the axonal misguidance phenotype. Moreover, our findings demonstrate that this functionally redundant specificity to mediate axon guidance has been conserved in their vertebrate homologs, L1-CAM and NCAM.

Development of foregut and proboscis in the buccinid neogastropodNassarius mendicus: Evolutionary opportunity exploited by a developmental module

This article extends previous descriptions of foregut development and metamorphosis in neogastropods by providing data on the buccinid Nassarius mendicus, a species with a feeding larva. Histological sections showed that, like many other gastropods, the postmetamorphic buccal cavity and radular sac of N. mendicus differentiate during the larval stage from a ventral outpocketing of the distal larval esophagus. However, in N. mendicus the outpocketing also gives rise to the entire anterior esophagus and valve of Leiblein, suggesting that both these structures may be evolutionary derivatives of the gastropod buccal cavity. Scanning electron microscopy and three-dimensional reconstructions of section profiles revealed that the distal larval esophagus and larval mouth are completely destroyed at metamorphosis. The postmetamorphic mouth is formed as a new orifice. Furthermore, epithelia covering the proboscis and proboscis sac arise from preexisting epidermal epithelium of the larval head, an interpretation that contradicts an earlier suggestion on the origin of these epithelial elements in neogastropods with a feeding larval stage. These results, when compared to foregut development in other gastropods, lead me to propose that the gastropod buccal cavity and buccal mass is a developmental module. Canalized development of this module may have been important to the "evolvability" of the complex gastropod foregut, because it allowed a silent developmental novelty to arise (secondary formation of the postmetamorphic mouth) without disrupting development of the whole module. Nevertheless, this silent novelty might have subsequently facilitated dramatic evolutionary change by allowing the elaborate foregut structure of predatory, postmetamorphic neogastropods to arise in late stage larvae without compromising larval feeding.

Perspective: Evolution and detection of genetic robustness

Robustness is the invariance of phenotypes in the face of perturbation. The robustness of phenotypes appears at various levels of biological organization, including gene expression, protein folding, metabolic flux, physiological homeostasis, development, and even organismal fitness. The mechanisms underlying robustness are diverse, ranging from thermodynamic stability at the RNA and protein level to behavior at the organismal level. Phenotypes can be robust either against heritable perturbations (e.g., mutations) or nonheritable perturbations (e.g., the weather). Here we primarily focus on the first kind of robustness--genetic robustness--and survey three growing avenues of research: (1) measuring genetic robustness in nature and in the laboratory; (2) understanding the evolution of genetic robustness: and (3) exploring the implications of genetic robustness for future evolution.

Under cover: causes, effects and implications of Hsp90-mediated genetic capacitance

The environmentally responsive molecular chaperone Hsp90 assists the maturation of many key regulatory proteins. An unexpected consequence of this essential biochemical function is that genetic variation can accumulate in genomes and can remain phenotypically silent until Hsp90 function is challenged. Notably, this variation can be revealed by modest environmental change, establishing an environmentally responsive exposure mechanism. The existence of diverse cryptic polymorphisms with a plausible exposure mechanism in evolutionarily distant lineages has implications for the pace and nature of evolutionary change. Chaperone-mediated storage and release of genetic variation is undoubtedly rooted in protein-folding phenomena. As we discuss, proper protein folding crucially affects the trajectory from genotype to phenotype. Indeed, the impact of protein quality-control mechanisms and other fundamental cellular processes on evolution has heretofore been overlooked. A true understanding of evolutionary processes will require an integration of current evolutionary paradigms with the many new insights accruing in protein science.

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.

Epistasis in polygenic traits and the evolution of genetic architecture under stabilizing selection

We consider the effects of epistasis in a polygenic trait in the balance of mutation and stabilizing selection. The main issues are the genetic variation maintained in equilibrium and the evolution of the mutational effect distribution. The model assumes symmetric mutation and a continuum of alleles at all loci. Epistasis is modeled proportional to pairwise products of the single-locus effects. A general analytical formalism is developed. Assuming linkage equilibrium, we derive results for the equilibrium mutation load and the genetic and mutational variance in the house of cards and the Gaussian approximation. The additive genetic variation maintained in mutation-selection balance is reduced by any pattern of the epistatic interactions. The mutational variance, in contrast, is often increased. Large differences in mutational effects among loci emerge, and a negative correlation among (standard mean) locus mutation effects and mutation rates is predicted. Contrary to the common view since Waddington, we find that stabilizing selection in general does not lead to canalization of the trait. We propose that canalization as a target of selection instead occurs at the genic level. Here, primarily genes with a high mutation rate are buffered, often at the cost of decanalization of other genes. An intuitive interpretation of this view is given in the discussion.

From epigenesis to epigenetics: the case of C. H. Waddington

One continuous thread in this volume is the name of Conrad H. Waddington (1905-1975), the developmental biologist known as the inventor of the term epigenetics. After some biographical notes on his life, this article explores the meaning of the Waddingtonian equation and the context wherein it was developed. This equation holds that epigenesis + genetics = epigenetics, and refers in retrospect to the debate on epigenesis versus preformationism in neoclassical embryology. Whereas Waddington actualized this debate by linking epigenesis to developmental biology and preformation to genetics, thereby stressing the importance of genetic action in causal embryology, today's epigenetics more and more offers the possibility to enfeeble biological thinking in terms of genes only, as it expands the gene-centric view in biology by introducing a flexible and pragmatically oriented hierarchy of crucial genomic contexts that go beyond the organism.