The population genetic theory of hidden variation and genetic robustness

Second premolar agenesis is associated with mandibular form: a geometric morphometric analysis of mandibular cross-sections

The aim of this study was to compare mandibular form (i.e., size and shape) between patients with agenesis of the lower second premolar (P2) and a control group with no agenesis. Three hypotheses were tested: (H1) agenesis causes a change in mandibular morphology because of inadequate alveolar ridge development in the area of the missing tooth (mandibular plasticity); (H2) agenesis is caused by spatial limitations within the mandible (dental plasticity); and (H3) common genetic/epigenetic factors cause agenesis and affect mandibular form (pleiotropy). A geometric morphometric analysis was applied to cross-sectional images of computed tomography (CT) scans of three matched groups (n=50 each): (1) regularly erupted P2; (2) agenesis of P2 and the primary second molar in situ; and (3) agenesis of P2 and the primary second molar missing for >3 months. Cross-sections of the three areas of interest (first premolar, P2, first molar) were digitized with 23 landmarks and superimposed by a generalized Procrustes analysis. On average, the mandibular cross-sections were narrower and shorter in patients with P2 agenesis compared with that in the control group. Both agenesis groups featured a pronounced submandibular fossa. These differences extended at least one tooth beyond the agenesis-affected region. Taken together with the large interindividual variation that resulted in massively overlapping group distributions, these findings support genetic and/or epigenetic pleiotropy (H3) as the most likely origin of the observed covariation between mandibular form and odontogenesis. Clinically, reduced dimensions and greater variability of mandibular form, as well as a pronounced submandibular fossa, should be expected during the treatment planning of patients with P2 agenesis.

Modifiers of the Genotype-Phenotype Map: Hsp90 and Beyond

Disruption of certain genes alters the heritable phenotypic variation among individuals. Research on the chaperone Hsp90 has played a central role in determining the genetic basis of this phenomenon, which may be important to evolution and disease. Key studies have shown that Hsp90 perturbation modifies the effects of many genetic variants throughout the genome. These modifications collectively transform the genotype–phenotype map, often resulting in a net increase or decrease in heritable phenotypic variation. Here, we summarize some of the foundational work on Hsp90 that led to these insights, discuss a framework for interpreting this research that is centered upon the standard genetics concept of epistasis, and propose major questions that future studies in this area should address.

Why and how genetic canalization evolves in gene regulatory networks

Background

Genetic canalization reflects the capacity of an organism’s phenotype to remain unchanged in spite of mutations. As selection on genetic canalization is weak and indirect, whether or not genetic canalization can reasonably evolve in complex genetic architectures is still an open question. In this paper, we use a quantitative model of gene regulatory network to describe the conditions in which substantial canalization is expected to emerge in a stable environment.

Results

Through an individual-based simulation framework, we confirmed that most parameters associated with the network topology (complexity and size of the network) have less influence than mutational parameters (rate and size of mutations) on the evolution of genetic canalization. We also established that selecting for extreme phenotypic optima (nil or full gene expression) leads to much higher canalization levels than selecting for intermediate expression levels. Overall, constrained networks evolve less canalization than networks in which some genes could evolve freely (i.e. without direct stabilizing selection pressure on gene expression).

Conclusions

Taken together, these results lead us to propose a two-fold mechanism involved in the evolution of genetic canalization in gene regulatory networks: the shrinkage of mutational target (useless genes are virtually removed from the network) and redundancy in gene regulation (so that some regulatory factors can be lost without affecting gene expression).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0801-2) contains supplementary material, which is available to authorized users.

Selection Transforms the Landscape of Genetic Variation Interacting with Hsp90

The protein-folding chaperone Hsp90 has been proposed to buffer the phenotypic effects of mutations. The potential for Hsp90 and other putative buffers to increase robustness to mutation has had major impact on disease models, quantitative genetics, and evolutionary theory. But Hsp90 sometimes contradicts expectations for a buffer by potentiating rapid phenotypic changes that would otherwise not occur. Here, we quantify Hsp90’s ability to buffer or potentiate (i.e., diminish or enhance) the effects of genetic variation on single-cell morphological features in budding yeast. We corroborate reports that Hsp90 tends to buffer the effects of standing genetic variation in natural populations. However, we demonstrate that Hsp90 tends to have the opposite effect on genetic variation that has experienced reduced selection pressure. Specifically, Hsp90 tends to enhance, rather than diminish, the effects of spontaneous mutations and recombinations. This result implies that Hsp90 does not make phenotypes more robust to the effects of genetic perturbation. Instead, natural selection preferentially allows buffered alleles to persist and thereby creates the false impression that Hsp90 confers greater robustness.

Most biologists appreciate that natural selection filters new mutations (e.g., by eliminating deleterious ones), such that genetic variation in nature is biased. The idea that selection also skews the types of genetic interactions that exist in nature is less appreciated. For example, studies spanning diverse species have shown that the protein Hsp90, which helps other proteins to fold properly, tends to diminish the observable effects of genetic variation. This observation has led to the assumption that Hsp90 also buffers the effects of new mutations. This untested assumption has served as a rationale for cancer-treatment strategies and shaped our understanding of variation in complex traits. We measured the effects of new mutations on the shapes and sizes of individual yeast cells and found that Hsp90 does not tend to buffer these effects. Instead, Hsp90 interacts with new mutations in diverse ways, sometimes buffering, but more often enhancing mutational effects on cell shape and size. We conclude that selection preferentially allows buffered mutations to persist in natural populations. This result alters common perceptions about why cryptic (i.e., buffered) genetic variation exists and casts doubt on cancer-treatment strategies aiming to target presumed buffers of mutational effects.

Genetic structure of phenotypic robustness in the collaborative cross mouse diallel panel

Developmental stability and canalization describe the ability of developmental systems to minimize phenotypic variation in the face of stochastic micro-environmental effects, genetic variation and environmental influences. Canalization is the ability to minimize the effects of genetic or environmental effects, whereas developmental stability is the ability to minimize the effects of micro-environmental effects within individuals. Despite much attention, the mechanisms that underlie these two components of phenotypic robustness remain unknown. We investigated the genetic structure of phenotypic robustness in the collaborative cross (CC) mouse reference population. We analysed the magnitude of fluctuating asymmetry (FA) and among-individual variation of cranial shape in reciprocal crosses among the eight parental strains, using geometric morphometrics and a diallel analysis based on a Bayesian approach. Significant differences among genotypes were found for both measures, although they were poorly correlated at the level of individuals. An overall positive effect of inbreeding was found for both components of variation. The strain CAST/EiJ exerted a positive additive effect on FA and, to a lesser extent, among-individual variance. Sex- and other strain-specific effects were not significant. Neither FA nor among-individual variation was associated with phenotypic extremeness. Our results support the existence of genetic variation for both developmental stability and canalization. This finding is important because robustness is a key feature of developmental systems. Our finding that robustness is not related to phenotypic extremeness is consistent with theoretical work that suggests that its relationship to stabilizing selection is not straightforward.

Cryptic Genetic Variation in Evolutionary Developmental Genetics

Evolutionary developmental genetics has traditionally been conducted by two groups: Molecular evolutionists who emphasize divergence between species or higher taxa, and quantitative geneticists who study variation within species. Neither approach really comes to grips with the complexities of evolutionary transitions, particularly in light of the realization from genome-wide association studies that most complex traits fit an infinitesimal architecture, being influenced by thousands of loci. This paper discusses robustness, plasticity and lability, phenomena that we argue potentiate major evolutionary changes and provide a bridge between the conceptual treatments of macro- and micro-evolution. We offer cryptic genetic variation and conditional neutrality as mechanisms by which standing genetic variation can lead to developmental system drift and, sheltered within canalized processes, may facilitate developmental transitions and the evolution of novelty. Synthesis of the two dominant perspectives will require recognition that adaptation, divergence, drift and stability all depend on similar underlying quantitative genetic processes-processes that cannot be fully observed in continuously varying visible traits.

Diverse genetic architectures lead to the same cryptic phenotype in a yeast cross

Cryptic genetic variants that do not typically influence traits can interact epistatically with each other and mutations to cause unexpected phenotypes. To improve understanding of the genetic architectures and molecular mechanisms that underlie these interactions, we comprehensively dissected the genetic bases of 17 independent instances of the same cryptic colony phenotype in a yeast cross. In eight cases, the phenotype resulted from a genetic interaction between a de novo mutation and one or more cryptic variants. The number and identities of detected cryptic variants depended on the mutated gene. In the nine remaining cases, the phenotype arose without a de novo mutation due to two different classes of higher-order genetic interactions that only involve cryptic variants. Our results may be relevant to other species and disease, as most of the mutations and cryptic variants identified in our study reside in components of a partially conserved and oncogenic signalling pathway.

Cryptic genetic variants may not individually show discernible phenotypic effects, but collectively, these polymorphisms can lead to unexpected, genetically complex traits that might be relevant to evolution and disease. Here, the authors use large yeast populations to comprehensively dissect the genetic bases of 17 independent occurrences of a phenotype that arises due to combinations of epistatically interacting cryptic variants.

Looking Beyond the Genes: The Interplay Between Signaling Pathways and Mechanics in the Shaping and Diversification of Epithelial Tissues

The core of Evo-Devo lies in the intuition that the way tissues grow during embryonic development, the way they sustain their structure and function throughout lifetime, and the way they evolve are closely linked. Epithelial tissues are ubiquitous in metazoans, covering the gut and internal branched organs, as well as the skin and its derivatives (ie, teeth). Here, we discuss in vitro, in vivo, and in silico studies on epithelial tissues to illustrate the conserved, dynamical, and complex aspects of their development. We then explore the implications of the dynamical and nonlinear nature of development on the evolution of their size and shape at the phenotypic and genetic levels. In rare cases, when the interplay between signaling and mechanics is well understood at the cell level, it is becoming clear that the structure of development leads to covariation of characters, an integration which in turn provides some predictable structure to evolutionary changes. We suggest that such nonlinear systems are prone to genetic drift, cryptic genetic variation, and context-dependent mutational effects. We argue that experimental and theoretical studies at the cell level are critical to our understanding of the phenotypic and genetic evolution of epithelial tissues, including carcinomas.

Mitonuclear Epistasis for Development Time and Its Modification by Diet in Drosophila

Mitochondrial (mtDNA) and nuclear genes have to operate in a coordinated manner to maintain organismal function, and the regulation of this homeostasis presents a substantial source of potential epistatic (G × G) interactions. How these interactions shape the fitness landscape is poorly understood. Here we developed a novel mitonuclear epistasis model, using selected strains of the Drosophila Genetic Reference Panel (DGRP) and mitochondrial genomes from within Drosophila melanogaster and D. simulans to test the hypothesis that mtDNA × nDNA interactions influence fitness. In total we built 72 genotypes (12 nuclear backgrounds × 6 mtDNA haplotypes, with 3 from each species) to dissect the relationship between genotype and phenotype. Each genotype was assayed on four food environments. We found considerable variation in several phenotypes, including development time and egg-to-adult viability, and this variation was partitioned into genetic (G), environmental (E), and higher-order (G × G, G × E, and G × G × E) components. Food type had a significant impact on development time and also modified mitonuclear epistases, evidencing a broad spectrum of G × G × E across these genotypes. Nuclear background effects were substantial, followed by mtDNA effects and their G × G interaction. The species of mtDNA haplotype had negligible effects on phenotypic variation and there was no evidence that mtDNA variation has different effects on male and female fitness traits. Our results demonstrate that mitonuclear epistases are context dependent, suggesting the selective pressure acting on mitonuclear genotypes may vary with food environment in a genotype-specific manner.

Genetic assimilation: a review of its potential proximate causes and evolutionary consequences

BACKGROUND

Most, if not all, organisms possess the ability to alter their phenotype in direct response to changes in their environment, a phenomenon known as phenotypic plasticity. Selection can break this environmental sensitivity, however, and cause a formerly environmentally induced trait to evolve to become fixed through a process called genetic assimilation. Essentially, genetic assimilation can be viewed as the evolution of environmental robustness in what was formerly an environmentally sensitive trait. Because genetic assimilation has long been suggested to play a key role in the origins of phenotypic novelty and possibly even new species, identifying and characterizing the proximate mechanisms that underlie genetic assimilation may advance our basic understanding of how novel traits and species evolve.

SCOPE

This review begins by discussing how the evolution of phenotypic plasticity, followed by genetic assimilation, might promote the origins of new traits and possibly fuel speciation and adaptive radiation. The evidence implicating genetic assimilation in evolutionary innovation and diversification is then briefly considered. Next, the potential causes of phenotypic plasticity generally and genetic assimilation specifically are examined at the genetic, molecular and physiological levels and approaches that can improve our understanding of these mechanisms are described. The review concludes by outlining major challenges for future work.

CONCLUSIONS

Identifying and characterizing the proximate mechanisms involved in phenotypic plasticity and genetic assimilation promises to help advance our basic understanding of evolutionary innovation and diversification.

Multi-locus Genotypes Underlying Temperature Sensitivity in a Mutationally Induced Trait

Determining how genetic variation alters the expression of heritable phenotypes across conditions is important for agriculture, evolution, and medicine. Central to this problem is the concept of genotype-by-environment interaction (or 'GxE'), which occurs when segregating genetic variation causes individuals to show different phenotypic responses to the environment. While many studies have sought to identify individual loci that contribute to GxE, obtaining a deeper understanding of this phenomenon may require defining how sets of loci collectively alter the relationship between genotype, environment, and phenotype. Here, we identify combinations of alleles at seven loci that control how a mutationally induced colony phenotype is expressed across a range of temperatures (21, 30, and 37 °C) in a panel of yeast recombinants. We show that five predominant multi-locus genotypes involving the detected loci result in trait expression with varying degrees of temperature sensitivity. By comparing these genotypes and their patterns of trait expression across temperatures, we demonstrate that the involved alleles contribute to temperature sensitivity in different ways. While alleles of the transcription factor MSS11 specify the potential temperatures at which the trait can occur, alleles at the other loci modify temperature sensitivity within the range established by MSS11 in a genetic background- and/or temperature-dependent manner. Our results not only represent one of the first characterizations of GxE at the resolution of multi-locus genotypes, but also provide an example of the different roles that genetic variants can play in altering trait expression across conditions.

Novel genetic capacitors and potentiators for the natural genetic variation of sensory bristles and their trait specificity in Drosophila melanogaster

Cryptic genetic variation (CGV) is defined as the genetic variation that has little effect on phenotypic variation under a normal condition, but contributes to heritable variation under environmental or genetic perturbations. Genetic buffering systems that suppress the expression of CGV and store it in a population are called genetic capacitors, and the opposite systems are called genetic potentiators. One of the best-known candidates for a genetic capacitor and potentiator is the molecular chaperone protein, HSP90, and one of its characteristics is that it affects the genetic variation in various morphological traits. However, it remains unclear whether the wide-ranging effects of HSP90 on a broad range of traits are a general feature of genetic capacitors and potentiators. In the current study, I searched for novel genetic capacitors and potentiators for quantitative bristle traits of Drosophila melanogaster and then investigated the trait specificity of their genetic buffering effect. Three bristle traits of D. melanogaster were used as the target traits, and the genomic regions with genetic buffering effects were screened using the 61 genomic deficiencies examined previously for genetic buffering effects in wing shape. As a result, four and six deficiencies with significant effects on increasing and decreasing the broad-sense heritability of the bristle traits were identified, respectively. Of the 18 deficiencies with significant effects detected in the current study and/or by the previous study, 14 showed trait-specific effects, and four affected the genetic buffering of both bristle traits and wing shape. This suggests that most genetic capacitors and potentiators exert trait-specific effects, but that general capacitors and potentiators with effects on multiple traits also exist.

Transcriptional Derepression Uncovers Cryptic Higher-Order Genetic Interactions

Disruption of certain genes can reveal cryptic genetic variants that do not typically show phenotypic effects. Because this phenomenon, which is referred to as ‘phenotypic capacitance’, is a potential source of trait variation and disease risk, it is important to understand how it arises at the genetic and molecular levels. Here, we use a cryptic colony morphology trait that segregates in a yeast cross to explore the mechanisms underlying phenotypic capacitance. We find that the colony trait is expressed when a mutation in IRA2, a negative regulator of the Ras pathway, co-occurs with specific combinations of cryptic variants in six genes. Four of these genes encode transcription factors that act downstream of the Ras pathway, indicating that the phenotype involves genetically complex changes in the transcriptional regulation of Ras targets. We provide evidence that the IRA2 mutation reveals the phenotypic effects of the cryptic variants by disrupting the transcriptional silencing of one or more genes that contribute to the trait. Supporting this role for the IRA2 mutation, deletion of SFL1, a repressor that acts downstream of the Ras pathway, also reveals the phenotype, largely due to the same cryptic variants that were detected in the IRA2 mutant cross. Our results illustrate how higher-order genetic interactions among mutations and cryptic variants can result in phenotypic capacitance in specific genetic backgrounds, and suggests these interactions might reflect genetically complex changes in gene expression that are usually suppressed by negative regulation.

Some genetic polymorphisms have phenotypic effects that are masked under most conditions, but can be revealed by mutations or environmental change. The genetic and molecular mechanisms that suppress and uncover these cryptic genetic variants are important to understand. Here, we show that a single mutation in a yeast cross causes a major phenotypic change through its genetic interactions with two specific combinations of cryptic variants in six genes. This result suggests that in some cases cryptic variants themselves play roles in revealing their own phenotypic effects through their genetic interactions with each other and the mutations that reveal them. We also demonstrate that most of the genes harboring cryptic variation in our system are transcription factors, a finding that supports an important role for perturbation of gene regulatory networks in the uncovering of cryptic variation. As a final part of our study, we interrogate how a mutation exposes combinations of cryptic variants and obtain evidence that it does so by disrupting the silencing of one or more genes that must be expressed for the cryptic variants to exert their effects. To prove this point, we delete the transcriptional repressor that mediates this silencing and demonstrate that this deletion reveals a similar set of cryptic variants to the ones that were discovered in the initial mutant background. These findings advance our understanding of the genetic and molecular mechanisms that reveal cryptic variation.

Mapping quantitative trait loci in selected breeding populations: A segregation distortion approach

Quantitative trait locus (QTL) mapping is often conducted in line-crossing experiments where a sample of individuals is randomly selected from a pool of all potential progeny. QTLs detected from such an experiment are important for us to understand the genetic mechanisms governing a complex trait, but may not be directly relevant to plant breeding if they are not detected from the breeding population where selection is targeting for. QTLs segregating in one population may not necessarily segregate in another population. To facilitate marker-assisted selection, QTLs must be detected from the very population which the selection is targeting. However, selected breeding populations often have depleted genetic variation with small population sizes, resulting in low power in detecting useful QTLs. On the other hand, if selection is effective, loci controlling the selected trait will deviate from the expected Mendelian segregation ratio. In this study, we proposed to detect QTLs in selected breeding populations via the detection of marker segregation distortion in either a single population or multiple populations using the same selection scheme. Simulation studies showed that QTL can be detected in strong selected populations with selected population sizes as small as 25 plants. We applied the new method to detect QTLs in two breeding populations of rice selected for high grain yield. Seven QTLs were identified, four of which have been validated in advanced generations in a follow-up study. Cloned genes in the vicinity of the four QTLs were also reported in the literatures. This mapping-by-selection approach provides a new avenue for breeders to improve breeding progress. The new method can be applied to breeding programs not only in rice but also in other agricultural species including crops, trees and animals.

The origins of mutational robustness

Biological systems are resistant to genetic changes; a property known as mutational robustness, the origin of which remains an open question. In recent years, researchers have explored emergent properties of biological systems and mechanisms of genetic redundancy to reveal how mutational robustness emerges and persists. Several mechanisms have been proposed to explain the origin of mutational robustness, including molecular chaperones and gene duplication. The latter has received much attention, but its role in robustness remains controversial. Here, I examine recent findings linking genetic redundancy through gene duplication and mutational robustness. Experimental evolution and genome resequencing have made it possible to test the role of gene duplication in tolerating mutations at both the coding and regulatory levels. This evidence as well as previous findings on regulatory reprogramming of duplicates support the role of gene duplication in the origin of robustness.

Predicting human genetic interactions from cancer genome evolution

Synthetic Lethal (SL) genetic interactions play a key role in various types of biological research, ranging from understanding genotype-phenotype relationships to identifying drug-targets against cancer. Despite recent advances in empirical measuring SL interactions in human cells, the human genetic interaction map is far from complete. Here, we present a novel approach to predict this map by exploiting patterns in cancer genome evolution. First, we show that empirically determined SL interactions are reflected in various gene presence, absence, and duplication patterns in hundreds of cancer genomes. The most evident pattern that we discovered is that when one member of an SL interaction gene pair is lost, the other gene tends not to be lost, i.e. the absence of co-loss. This observation is in line with expectation, because the loss of an SL interacting pair will be lethal to the cancer cell. SL interactions are also reflected in gene expression profiles, such as an under representation of cases where the genes in an SL pair are both under expressed, and an over representation of cases where one gene of an SL pair is under expressed, while the other one is over expressed. We integrated the various previously unknown cancer genome patterns and the gene expression patterns into a computational model to identify SL pairs. This simple, genome-wide model achieves a high prediction power (AUC = 0.75) for known genetic interactions. It allows us to present for the first time a comprehensive genome-wide list of SL interactions with a high estimated prediction precision, covering up to 591,000 gene pairs. This unique list can potentially be used in various application areas ranging from biotechnology to medical genetics.

The devil is in the details: genetic variation in introduced populations and its contributions to invasion

The influence of genetic variation on invasion success has captivated researchers since the start of the field of invasion genetics 50 years ago. We review the history of work on this question and conclude that genetic variation-as surveyed with molecular markers-appears to shape invasion rarely. Instead, there is a significant disconnect between marker assays and ecologically relevant genetic variation in introductions. We argue that the potential for adaptation to facilitate invasion will be shaped by the details of genotypes affecting phenotypes, and we highlight three areas in which we see opportunities to make powerful new insights. (i) The genetic architecture of adaptive variation. Traits shaped by large-effect alleles may be strongly impacted by founder events yet more likely to respond to selection when genetic drift is strong. Large-effect loci may be especially relevant for traits involved in biotic interactions. (ii) Cryptic genetic variation exposed during invasion. Introductions have strong potential to uncover masked variation due to alterations in genetic and ecological environments. (iii) Genetic interactions during admixture of multiple source populations. As divergence among sources increases, positive followed by increasingly negative effects of admixture should be expected. Although generally hypothesized to be beneficial during invasion, admixture is most often reported among sources of intermediate divergence, supporting the possibility that incompatibilities among divergent source populations might be limiting their introgression. Finally, we note that these details of invasion genetics can be coupled with comparative demographic analyses to link genetic changes to the evolution of invasiveness itself.

The unforeseen challenge: from genotype-to-phenotype in cell populations

Biological cells present a paradox, in that they show simultaneous stability and flexibility, allowing them to adapt to new environments and to evolve over time. The emergence of stable cell states depends on genotype-to-phenotype associations, which essentially reflect the organization of gene regulatory modes. The view taken here is that cell-state organization is a dynamical process in which the molecular disorder manifests itself in a macroscopic order. The genome does not determine the ordered cell state; rather, it participates in this process by providing a set of constraints on the spectrum of regulatory modes, analogous to boundary conditions in physical dynamical systems. We have developed an experimental framework, in which cell populations are exposed to unforeseen challenges; novel perturbations they had not encountered before along their evolutionary history. This approach allows an unbiased view of cell dynamics, uncovering the potential of cells to evolve and develop adapted stable states. In the last decade, our experiments have revealed a coherent set of observations within this framework, painting a picture of the living cell that in many ways is not aligned with the conventional one. Of particular importance here, is our finding that adaptation of cell-state organization is essentially an efficient exploratory dynamical process rather than one founded on random mutations. Based on our framework, a set of concepts underlying cell-state organization-exploration evolving by global, non-specific, dynamics of gene activity-is presented here. These concepts have significant consequences for our understanding of the emergence and stabilization of a cell phenotype in diverse biological contexts. Their implications are discussed for three major areas of biological inquiry: evolution, cell differentiation and cancer. There is currently no unified theoretical framework encompassing the emergence of order, a stable state, in the living cell. Hopefully, the integrated picture described here will provide a modest contribution towards a physics theory of the cell.

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.

A simulation study of gene-by-environment interactions in GWAS implies ample hidden effects

The switch to a modern lifestyle in recent decades has coincided with a rapid increase in prevalence of obesity and other diseases. These shifts in prevalence could be explained by the release of genetic susceptibility for disease in the form of gene-by-environment (GxE) interactions. Yet, the detection of interaction effects requires large sample sizes, little replication has been reported, and a few studies have demonstrated environmental effects only after summing the risk of GWAS alleles into genetic risk scores (GRSxE). We performed extensive simulations of a quantitative trait controlled by 2500 causal variants to inspect the feasibility to detect gene-by-environment interactions in the context of GWAS. The simulated individuals were assigned either to an ancestral or a modern setting that alters the phenotype by increasing the effect size by 1.05–2-fold at a varying fraction of perturbed SNPs (from 1 to 20%). We report two main results. First, for a wide range of realistic scenarios, highly significant GRSxE is detected despite the absence of individual genotype GxE evidence at the contributing loci. Second, an increase in phenotypic variance after environmental perturbation reduces the power to discover susceptibility variants by GWAS in mixed cohorts with individuals from both ancestral and modern environments. We conclude that a pervasive presence of gene-by-environment effects can remain hidden even though it contributes to the genetic architecture of complex traits.

Estimating directional epistasis

Epistasis, i.e., the fact that gene effects depend on the genetic background, is a direct consequence of the complexity of genetic architectures. Despite this, most of the models used in evolutionary and quantitative genetics pay scant attention to genetic interactions. For instance, the traditional decomposition of genetic effects models epistasis as noise around the evolutionarily-relevant additive effects. Such an approach is only valid if it is assumed that there is no general pattern among interactions—a highly speculative scenario. Systematic interactions generate directional epistasis, which has major evolutionary consequences. In spite of its importance, directional epistasis is rarely measured or reported by quantitative geneticists, not only because its relevance is generally ignored, but also due to the lack of simple, operational, and accessible methods for its estimation. This paper describes conceptual and statistical tools that can be used to estimate directional epistasis from various kinds of data, including QTL mapping results, phenotype measurements in mutants, and artificial selection responses. As an illustration, I measured directional epistasis from a real-life example. I then discuss the interpretation of the estimates, showing how they can be used to draw meaningful biological inferences.

Mutation of the RDR1 gene caused genome-wide changes in gene expression, regional variation in small RNA clusters and localized alteration in DNA methylation in rice

BACKGROUND

Endogenous small (sm) RNAs (primarily si- and miRNAs) are important trans/cis-acting regulators involved in diverse cellular functions. In plants, the RNA-dependent RNA polymerases (RDRs) are essential for smRNA biogenesis. It has been established that RDR2 is involved in the 24 nt siRNA-dependent RNA-directed DNA methylation (RdDM) pathway. Recent studies have suggested that RDR1 is involved in a second RdDM pathway that relies mostly on 21 nt smRNAs and functions to silence a subset of genomic loci that are usually refractory to the normal RdDM pathway in Arabidopsis. Whether and to what extent the homologs of RDR1 may have similar functions in other plants remained unknown.

RESULTS

We characterized a loss-of-function mutant (Osrdr1) of the OsRDR1 gene in rice (Oryza sativa L.) derived from a retrotransposon Tos17 insertion. Microarray analysis identified 1,175 differentially expressed genes (5.2% of all expressed genes in the shoot-tip tissue of rice) between Osrdr1 and WT, of which 896 and 279 genes were up- and down-regulated, respectively, in Osrdr1. smRNA sequencing revealed regional alterations in smRNA clusters across the rice genome. Some of the regions with altered smRNA clusters were associated with changes in DNA methylation. In addition, altered expression of several miRNAs was detected in Osrdr1, and at least some of which were associated with altered expression of predicted miRNA target genes. Despite these changes, no phenotypic difference was identified in Osrdr1 relative to WT under normal condition; however, ephemeral phenotypic fluctuations occurred under some abiotic stress conditions.

CONCLUSIONS

Our results showed that OsRDR1 plays a role in regulating a substantial number of endogenous genes with diverse functions in rice through smRNA-mediated pathways involving DNA methylation, and which participates in abiotic stress response.

The genotype-phenotype map of yeast complex traits: basic parameters and the role of natural selection

Most phenotypic traits are controlled by many genes, but a global picture of the genotype-phenotype map (GPM) is lacking. For example, in no species do we know generally how many genes affect a trait and how large these effects are. It is also unclear to what extent GPMs are shaped by natural selection. Here we address these fundamental questions using the reverse genetic data of 220 morphological traits in 4,718 budding yeast strains, each of which lacks a nonessential gene. We show that 1) the proportion of genes affecting a trait varies from <1% to >30%, averaging 6%, 2) most traits are impacted by many more small-effect genes than large-effect genes, and 3) the mean effect of all nonessential genes on a trait decreases precipitously as the estimated importance of the trait to fitness increases. An analysis of 3,116 yeast gene expression traits in 754 gene-deletion strains reveals a similar pattern. These findings illustrate the power of genome-wide reverse genetics in genotype-phenotype mapping, uncover an enormous range of genetic complexity of phenotypic traits, and suggest that the GPM of cellular organisms has been shaped by natural selection for mutational robustness.

Transposable elements in cancer as a by-product of stress-induced evolvability

Transposable elements (TEs) are ubiquitous in eukaryotic genomes. Barbara McClintock's famous notion of TEs acting as controlling elements modifying the genetic response of an organism upon exposure to stressful environments has since been solidly supported in a series of model organisms. This requires the TE activity response to possess an element of specificity and be targeted toward certain parts of the genome. We propose that a similar TE response is present in human cells, and that this stress response may drive the onset of human cancers. As such, TE-driven cancers may be viewed as an evolutionary by-product of organisms' abilities to genetically adapt to environmental stress.

Epistasis and natural selection shape the mutational architecture of complex traits

The evolutionary trajectories of complex traits are constrained by levels of genetic variation as well as genetic correlations among traits. As the ultimate source of all genetic variation is mutation, the distribution of mutations entering populations profoundly affects standing variation and genetic correlations. Here we use an individual-based simulation model to investigate how natural selection and gene interactions (that is, epistasis) shape the evolution of mutational processes affecting complex traits. We find that the presence of epistasis allows natural selection to mould the distribution of mutations, such that mutational effects align with the selection surface. Consequently, novel mutations tend to be more compatible with the current forces of selection acting on the population. These results suggest that in many cases mutational effects should be seen as an outcome of natural selection rather than as an unbiased source of genetic variation that is independent of other evolutionary processes.

Tissue culture-induced heritable genomic variation in rice, and their phenotypic implications

BACKGROUND

Somaclonal variation generally occurs in plants regenerated from tissue culture. However, fundamental issues regarding molecular characteristics, mutation rates and mutation spectra of plant somatic variation as well as their phenotypic relevance have been addressed only recently. Moreover, these studies have reported highly discrepant results in different plant species and even in the same plant genotype.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated heritable genomic variation induced by tissue culture in rice by whole genome re-sequencing of an extensively selfed somaclonal line (TC-reg-2008) and its wild type (WT) donor (cv. Hitomebore). We computed the overall mutation rate, single nucleotide polymorphisms (SNPs), small scale insertions/deletions (Indels) and mobilization of transposable elements (TEs). We assessed chromosomal distribution of the various types of genomic variations, tested correlations between SNPs and Indels, and examined concomitancy between TE activity and its cytosine methylation states. We also performed gene ontology (GO) analysis of genes containing nonsynonymous mutations and large-effect mutations, and assayed effects of the genomic variations on phenotypes under both normal growing condition and several abiotic stresses. We found that heritable somaclonal genomic variation occurred extensively in rice. The genomic variations distributed non-randomly across each of the 12 rice chromosomes, and affected a large number of functional genes. The phenotypic penetrance of the genomic variations was condition-dependent.

CONCLUSIONS/SIGNIFICANCE

Tissue culture is a potent means to generate heritable genetic variations in rice, which bear distinct difference at least in space (chromosomal distribution) from those occurred under natural settings. Our findings have provided new information regarding the mutation rate and spectrum as well as chromosomal distribution pattern of somaclonal variation in rice. Our data also suggest that rice possesses a strong capacity to canalize genetic variations under normal growing conditions to maintain phenotypic robustness, which however can be released by certain abiotic stresses to generate variable phenotypes.

Genotypic-specific variance in Caenorhabditis elegans lifetime fecundity

Organisms live in heterogeneous environments, so strategies that maximze fitness in such environments will evolve. Variation in traits is important because it is the raw material on which natural selection acts during evolution. Phenotypic variation is usually thought to be due to genetic variation and/or environmentally induced effects. Therefore, genetically identical individuals in a constant environment should have invariant traits. Clearly, genetically identical individuals do differ phenotypically, usually thought to be due to stochastic processes. It is now becoming clear, especially from studies of unicellular species, that phenotypic variance among genetically identical individuals in a constant environment can be genetically controlled and that therefore, in principle, this can be subject to selection. However, there has been little investigation of these phenomena in multicellular species. Here, we have studied the mean lifetime fecundity (thus a trait likely to be relevant to reproductive success), and variance in lifetime fecundity, in recently-wild isolates of the model nematode Caenorhabditis elegans. We found that these genotypes differed in their variance in lifetime fecundity: some had high variance in fecundity, others very low variance. We find that this variance in lifetime fecundity was negatively related to the mean lifetime fecundity of the lines, and that the variance of the lines was positively correlated between environments. We suggest that the variance in lifetime fecundity may be a bet-hedging strategy used by this species.

Cryptic genetic variation: evolution's hidden substrate

Cryptic genetic variation (CGV) is invisible under normal conditions, but it can fuel evolution when circumstances change. In theory, CGV can represent a massive cache of adaptive potential or a pool of deleterious alleles that are in need of constant suppression. CGV emerges from both neutral and selective processes, and it may inform about how human populations respond to change. CGV facilitates adaptation in experimental settings, but does it have an important role in the real world? Here, we review the empirical support for widespread CGV in natural populations, including its potential role in emerging human diseases and the growing evidence of its contribution to evolution.

Phenotypic plasticity and epigenetic marking: an assessment of evidence for genetic accommodation

The relationship between genotype (which is inherited) and phenotype (the target of selection) is mediated by environmental inputs on gene expression, trait development, and phenotypic integration. Phenotypic plasticity or epigenetic modification might influence evolution in two general ways: (1) by stimulating evolutionary responses to environmental change via population persistence or by revealing cryptic genetic variation to selection, and (2) through the process of genetic accommodation, whereby natural selection acts to improve the form, regulation, and phenotypic integration of novel phenotypic variants. We provide an overview of models and mechanisms for how such evolutionary influences may be manifested both for plasticity and epigenetic marking. We point to promising avenues of research, identifying systems that can best be used to address the role of plasticity in evolution, as well as the need to apply our expanding knowledge of genetic and epigenetic mechanisms to our understanding of how genetic accommodation occurs in nature. Our review of a wide variety of studies finds widespread evidence for evolution by genetic accommodation.

Nonlinear effects of temperature on body form and developmental canalization in the threespine stickleback

Theoretical models predict that nonlinear environmental effects on the phenotype also affect developmental canalization, which in turn can influence the tempo and course of organismal evolution. Here, we used an oceanic population of threespine stickleback (Gasterosteus aculeatus) to investigate temperature-induced phenotypic plasticity of body size and shape using a paternal half-sibling, split-clutch experimental design and rearing offspring under three different temperature regimes (13, 17 and 21 °C). Body size and shape of 466 stickleback individuals were assessed by a set of 53 landmarks and analysed using geometric morphometric methods. At approximately 100 days, individuals differed significantly in both size and shape across the temperature groups. However, the temperature-induced differences between 13 and 17 °C (mainly comprising relative head and eye size) deviated considerably from those between 17 and 21 °C (involving the relative size of the ectocoracoid, the operculum and the ventral process of the pelvic girdle). Body size was largest at 17 °C. For both size and shape, phenotypic variance was significantly smaller at 17 °C than at 13 and 21 °C, indicating that development is most stable at the intermediate temperature matching the conditions encountered in the wild. Higher additive genetic variance at 13 and 21 °C indicates that the plastic response to temperature had a heritable basis. Understanding nonlinear effects of temperature on development and the underlying genetics are important for modelling evolution and for predicting outcomes of global warming, which can lead not only to shifts in average morphology but also to destabilization of development.

Rapid evolution of quantitative traits: theoretical perspectives

An increasing number of studies demonstrate phenotypic and genetic changes in natural populations that are subject to climate change, and there is hope that some of these changes will contribute to avoiding species extinctions ('evolutionary rescue'). Here, we review theoretical models of rapid evolution in quantitative traits that can shed light on the potential for adaptation to a changing climate. Our focus is on quantitative-genetic models with selection for a moving phenotypic optimum. We point out that there is no one-to-one relationship between the rate of adaptation and population survival, because the former depends on relative fitness and the latter on absolute fitness. Nevertheless, previous estimates that sustainable rates of genetically based change usually do not exceed 0.1 haldanes (i.e., phenotypic standard deviations per generation) are probably correct. Survival can be greatly facilitated by phenotypic plasticity, and heritable variation in plasticity can further speed up genetic evolution. Multivariate selection and genetic correlations are frequently assumed to constrain adaptation, but this is not necessarily the case and depends on the geometric relationship between the fitness landscape and the structure of genetic variation. Similar conclusions hold for adaptation to shifting spatial gradients. Recent models of adaptation in multispecies communities indicate that the potential for rapid evolution is strongly influenced by interspecific competition.

Pleiotropic effects of a mitochondrial-nuclear incompatibility depend upon the accelerating effect of temperature in Drosophila

Interactions between mitochondrial and nuclear gene products that underlie eukaryotic energy metabolism can cause the fitness effects of mutations in one genome to be conditional on variation in the other genome. In ectotherms, the effects of these interactions are likely to depend upon the thermal environment, because increasing temperature accelerates molecular rates. We find that temperature strongly modifies the pleiotropic phenotypic effects of an incompatible interaction between a Drosophila melanogaster polymorphism in the nuclear-encoded, mitochondrial tyrosyl-transfer (t)RNA synthetase and a D. simulans polymorphism in the mitochondrially encoded tRNA(Tyr). The incompatible mitochondrial-nuclear genotype extends development time, decreases larval survivorship, and reduces pupation height, indicative of decreased energetic performance. These deleterious effects are ameliorated when larvae develop at 16° and exacerbated at warmer temperatures, leading to complete sterility in both sexes at 28°. The incompatible genotype has a normal metabolic rate at 16° but a significantly elevated rate at 25°, consistent with the hypothesis that inefficient energy metabolism extends development in this genotype at warmer temperatures. Furthermore, the incompatibility decreases metabolic plasticity of larvae developed at 16°, indicating that cooler development temperatures do not completely mitigate the deleterious effects of this genetic interaction. Our results suggest that the epistatic fitness effects of metabolic mutations may generally be conditional on the thermal environment. The expression of epistatic interactions in some environments, but not others, weakens the efficacy of selection in removing deleterious epistatic variants from populations and may promote the accumulation of incompatibilities whose fitness effects will depend upon the environment in which hybrids occur.

Histone variant HTZ1 shows extensive epistasis with, but does not increase robustness to, new mutations

Biological systems produce phenotypes that appear to be robust to perturbation by mutations and environmental variation. Prior studies identified genes that, when impaired, reveal previously cryptic genetic variation. This result is typically interpreted as evidence that the disrupted gene normally increases robustness to mutations, as such robustness would allow cryptic variants to accumulate. However, revelation of cryptic genetic variation is not necessarily evidence that a mutationally robust state has been made less robust. Demonstrating a difference in robustness requires comparing the ability of each state (with the gene perturbed or intact) to suppress the effects of new mutations. Previous studies used strains in which the existing genetic variation had been filtered by selection. Here, we use mutation accumulation (MA) lines that have experienced minimal selection, to test the ability of histone H2A.Z (HTZ1) to increase robustness to mutations in the yeast Saccharomyces cerevisiae. HTZ1, a regulator of chromatin structure and gene expression, represents a class of genes implicated in mutational robustness. It had previously been shown to increase robustness of yeast cell morphology to fluctuations in the external or internal microenvironment. We measured morphological variation within and among 79 MA lines with and without HTZ1. Analysis of within-line variation confirms that HTZ1 increases microenvironmental robustness. Analysis of between-line variation shows the morphological effects of eliminating HTZ1 to be highly dependent on the line, which implies that HTZ1 interacts with mutations that have accumulated in the lines. However, lines without HTZ1 are, as a group, not more phenotypically diverse than lines with HTZ1 present. The presence of HTZ1, therefore, does not confer greater robustness to mutations than its absence. Our results provide experimental evidence that revelation of cryptic genetic variation cannot be assumed to be caused by loss of robustness, and therefore force reevaluation of prior claims based on that assumption.

Genetic and environmental factors affecting cryptic variations in gene regulatory networks

BACKGROUND

Cryptic genetic variation (CGV) is considered to facilitate phenotypic evolution by producing visible variations in response to changes in the internal and/or external environment. Several mechanisms enabling the accumulation and release of CGVs have been proposed. In this study, we focused on gene regulatory networks (GRNs) as an important mechanism for producing CGVs, and examined how interactions between GRNs and the environment influence the number of CGVs by using individual-based simulations.

RESULTS

Populations of GRNs were allowed to evolve under various stabilizing selections, and we then measured the number of genetic and phenotypic variations that had arisen. Our results showed that CGVs were not depleted irrespective of the strength of the stabilizing selection for each phenotype, whereas the visible fraction of genetic variation in a population decreased with increasing strength of selection. On the other hand, increasing the number of different environments that individuals encountered within their lifetime (i.e., entailing plastic responses to multiple environments) suppressed the accumulation of CGVs, whereas the GRNs with more genes and interactions were favored in such heterogeneous environments.

CONCLUSIONS

Given the findings that the number of CGVs in a population was largely determined by the size (order) of GRNs, we propose that expansion of GRNs and adaptation to novel environments are mutually facilitating and sustainable sources of evolvability and hence the origins of biological diversity and complexity.

Asymmetric flies: the control of developmental noise in Drosophila

What are the sources of phenotypic variation and which factors shape this variation are fundamental questions of developmental and evolutionary biology. Despite this simple formulation and intense research, controversy remains. Three points are particularly discussed: (1) whether adaptive developmental mechanisms buffering variation exist at all; (2) if yes, do they involve specific genes and processes, i.e., different from those involved in the development of the traits that are buffered?; and (3) whether different mechanisms specifically buffer the various sources of variation, i.e., genetic, environmental and stochastic, or whether a generalist process buffers them all at once. We advocate that experimental work integrating different levels of analysis will improve our understanding of the origin of phenotypic variation and thus help answering these contentious questions. In this paper, we first survey the current views on these issues, highlighting potential sources of controversy. We then focus on the stochastic part of phenotypic variation, as measured by fluctuating asymmetry, and on current knowledge about the genetic basis of developmental stability. We report our recent discovery that an individual gene, Cyclin G, plays a central role-adaptive or not-in developmental stability in Drosophila. ( 1) We discuss the implications of this discovery on the regulation of organ size and shape, and finally point out open questions.

Directional selection causes decanalization in a group I ribozyme

A canalized genotype is robust to environmental or genetic perturbations. Canalization is expected to result from stabilizing selection on a well-adapted phenotype. Decanalization, the loss of robustness, might follow periods of directional selection toward a new optimum. The evolutionary forces causing decanalization are still unknown, in part because it is difficult to determine the fitness effects of mutations in populations of organisms with complex genotypes and phenotypes. Here, we report direct experimental measurements of robustness in a system with a simple genotype and phenotype, the catalytic activity of an RNA enzyme. We find that the robustness of a population of RNA enzymes decreases during a period of directional selection in the laboratory. The decrease in robustness is primarily caused by the selective sweep of a genotype that is decanalized relative to the wild-type, both in terms of mutational robustness and environmental robustness (thermodynamic stability). Our results experimentally demonstrate that directional selection can cause decanalization on short time scales, and demonstrate co-evolution of mutational and environmental robustness.

The relationship between phenotypic and environmental variation: do physiological responses reduce interindividual differences?

What is the effect of a variable environment on phenotypic variation? Does the physiological response to a new environment increase or decrease the differences among individuals? We provide a speculative hypothesis suggesting that the induction of a physiological response to environmental change minimizes phenotypic differences among individuals in outbred genetically variable populations. Although this suggestion runs counter to the general idea that environmental variation induces phenotypic variation, we provide evidence that this is not always the case. One explanation for this counterintuitive hypothesis is that in a variable environment, the physiological mechanism that maintains homeostasis changes the concentrations of active transcription factors (TFs). This change in TFs reduces the effectiveness of nucleotide polymorphisms in TF binding sites and thus reduces the variation among individuals in mRNA expression and in the phenotypes affected by these mRNA transcripts. Thus, there are fewer differences among individuals in a variable environment compared with the variation observed in a constant environment. Our conjecture is that the physiological mechanisms that maintain homeostasis in response to environmental variation canalize phenotypic variation. If our hypothesis is correct, then the physiological canalization of gene expression in a variable environment hides genetic variation and thereby reduces the evolutionary costs of polymorphism. This hypothesis provides a new perspective on the mechanisms by which high levels of genetic variation can persist in real-world populations.

Genetic divergence and the genetic architecture of complex traits in chromosome substitution strains of mice

Background

The genetic architecture of complex traits strongly influences the consequences of inherited mutations, genetic engineering, environmental and genetic perturbations, and natural and artificial selection. But because most studies are under-powered, the picture of complex traits is often incomplete. Chromosome substitution strains (CSSs) are a unique paradigm for these genome surveys because they enable statistically independent, powerful tests for the phenotypic effects of each chromosome on a uniform inbred genetic background. A previous CSS survey in mice and rats revealed many complex trait genes (QTLs), large phenotypic effects, extensive epistasis, as well as systems properties such as strongly directional phenotypic changes and genetically-determined limits on the range of phenotypic variation. However, the unusually close genetic relation between the CSS progenitor strains in that study raised questions about the impact of genetic divergence: would greater divergence between progenitor strains, with the corresponding changes in gene regulation and protein function, lead to significantly more distinctive phenotypic features, or alternatively would epistasis and systems constraints, which are pervasive in CSSs, limit the range of phenotypic variation regardless of the extent of DNA sequence variation?

Results

We analyzed results for an extensive survey of traits in two new panels of CSSs where the donor strains were derived from inbred strains with more distant origins and discovered a strong similarity in genetic and systems properties among the three CSS panels, regardless of divergence time.

Conclusion

Our results argue that DNA sequence differences between host and donor strains did not substantially affect the architecture of complex traits, and suggest instead that strong epistasis buffered the phenotypic effects of genetic divergence, thereby constraining the range of phenotypic variation.