Genetic and environmental factors affecting cryptic variations in gene regulatory networks
© Iwasaki et al.; licensee BioMed Central Ltd. 2013
Received: 7 January 2013
Accepted: 16 April 2013
Published: 26 April 2013
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.
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.
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.
The speed and direction of phenotypic evolution depends on the quantity and quality of genetic variations [1, 2], and there have been a number of attempts to quantify such variations [2, 3]. Understanding how biological populations generate and maintain genetic variations that contribute to phenotypic variations is one of the most important issues in evolutionary biology. When phenotypic variations are caused by additive genetic variations, evolutionary responses to directional selection are expected to be readily halted because the genetic variation is exhausted in the process of selection [4, 5]. However, several experiments involving artificial selection have shown that responses to directional selection can continue for many generations [6, 7]. In addition, a population of organisms sometimes shows much larger phenotypic variance when it encounters novel environments compared with that when it lives in normal environments. Standing genetic variations that do not translate into phenotypic differences in the current genetic and environmental background but that can become visible in a different background are called cryptic genetic variations (CGVs). They are considered to contribute to evolution by generating phenotypic diversity in response to changes in the environment and genetic background [8–11].
Various mechanisms are involved in preventing genetic variations from manifesting as phenotypic variations. The ability to retain phenotypes and functions despite internal and external perturbations is called biological robustness . The molecular chaperone called HSP90 is a well-known example of such a mechanism [13, 14]. It neutralizes non-synonymous substitutions on DNA sequences by assisting the proper folding of polypeptide chains. Defects in HSP90 or environmental stresses beyond its capacity cause the emergence of hidden variations. A system that contributes to evolvability by hiding and releasing CGVs is referred to as an evolutionary capacitor [13, 15, 16]. In addition, many genes other than HSP90 are suggested to be involved in stabilizing developmental processes and regulating CGVs [14, 17]. Apart form the mechanisms involving a molecular chaperone, redundancy and modularity also produce robustness in organisms .
One study  proposed a model whereby the release of hidden genetic variation due to a change in the environment or genetic background can be explained by epistasis or genotype–environment interactions. In fact, epistatic interactions were found to be involved in a number of directional selection experiments . Recent studies have also revealed how genes interact with each other during early development (in multi-cellular organisms) and in response to environmental stimuli (in microorganisms) [21–23]. It has recently been suggested that a phenotype is not only the sum of gene effects but also the product of complex interactions among genes [24–26]. Such networks of gene interactions are termed gene regulatory networks (GRNs), and these are considered to play crucial roles in cell differentiation and specific biological functions by modulating the expression of different gene combinations as well as the extent, site, and time of gene expression [27, 28]. Epistasis and genotype–environment interactions are generic features of GRNs [29, 30], and the biological properties of networks of interacting genes are thus fundamental factors that produce CGVs [29–33].
To understand the nature and the evolution of CGVs, we need to explore the structure and mechanisms of GRNs that accumulate hidden variations and release these variations upon exposure to novel environmental stimuli. We also need to determine how GRNs themselves change during evolutionary processes and to explore the conditions under which evolved GRNs can produce CGVs.
CGVs should be accumulated through population genetic processes, such as mutations, genetic drift, and natural selection. However, most studies of the robustness of GRNs [34, 35] have not explicitly assumed gene frequency changes in populations, and in particular, no attempt has been made to examine how CGVs could be accumulated through GRNs that have evolved via population genetic processes even in studies that have discussed the evolution of robustness due to changing gene frequencies [36, 37]. To understand the role of GRNs in cryptic variation, we need to construct models that include GRNs that produce different phenotypes when exposed to different environmental stimuli and to examine the number of CGVs acquired through genetic drift and natural selection for each phenotype produced by these GRNs.
This study aimed to show how interactions between GRNs and the environment influence the accumulation and emergence of CGVs. For this purpose, we constructed an individual-based model comprising individuals with plastic GRNs, i.e., GRNs that produce different phenotypes in response to different environmental stimuli. We then determined which factors, e.g., strength of selection, the mode of mutations, the size and properties of GRNs, and the number of environments that the organisms encountered, affected the number of CGVs.
Outline of the model
The lifetime fitness of an individual was measured as the geometric mean of the fitness of the phenotypes for different environments, and the fitness for each environment was calculated on the basis of a fitness function for the phenotype in each environment. The individuals reproduced according to their lifetime fitness values. When an offspring was produced, mutations (cis- and trans-regulatory mutations) were expected to occur. Gene duplications and deletions were not assumed in this model.
A population of individuals was allowed to evolve over a number of generations until a mutation-drift balance was achieved. The genetic and phenotypic diversity and the network characteristics of these were then measured. In addition, phenotypic diversity was measured when the individuals were exposed to novel environmental signals, i.e., signals that differed from those found in the normal environments during the evolution of the population.
The individuals had V g e n o m e (20 by default) genes, which together constituted a GRN. Genes at the downstream terminals, which are known as differentiation gene batteries, have been suggested to represent genes expressed in a cell that fulfill a specific function [23, 41]. We assumed that four of the V g e n o m e genes were phenotypic genes that did not have any control over the other genes and that the rest were regulatory genes. Phenotypic genes assigned by us were twice the genes assigned in a previous study  to enable us to apply distinct pressures upon phenotypes in different environments. Each gene had a cis-regulatory region and a coding region. The coding region of each regulatory gene produced a transcription factor when it was expressed. Each cis-regulatory region was composed of L (20 by default) cis-elements that were potentially recognized by specific transcription factors or signals induced by environmental stimuli. Each cis-element and transcription factors had a recognition number (a c i s and a t r a n s , an integral number between 1 and 200). Each transcription factor bound to the cis-elements with the corresponding number and induced the expression of the gene. A default value of L and a range of recognition numbers were chosen so that the expected number of target genes for a regulatory gene (L×V g e n o m e /200) was 2.0 . Each cis-element and transcription factor also had an interaction coefficient (v c i s and v t r a n s , a real number between −1 and 1) that affected the strength of the transcriptional activation/repression. The absolute value of a coefficient was chosen from uniform random numbers between 0 and 1, and its sign was positive with the fixed probability so that the expected proportion of positive edges in a GRN became A. We refer to A as the positive edge bias and set its default value to 0.6 given the observation that 60% of the edges in the bacteria GRN are positive . In addition to the transcription factors, which were the products of the coding regions, we assumed that each environmental signal was another type of trans-element bound to the cis-elements, and we set the interaction coefficient to +1.0.
where β m a x and K m are constants that determine the maximum transcription rate and the threshold against regulatory input; both were set to 0.1 so that the normalized expression level ranged from 0.0 to 1.0 and was balanced with the degradation rate α (described below). We assumed that the genes without transcriptional regulation (x=0) were also expressed to some extent; therefore, the basal transcription rate β b a s a l was set to 0.001, which is consistent with the previous study by Tsuda and Kawata. However, this parameter may not have greatly affected the phenotypes observed in this study because we obtained quite similar results when β b a s a l =0.
where the time step dt was set to 1.0 in the numerical calculations. We regarded the expression of the gene to be constant when the variance of the expression level in 32 recent loops was <10−6. These values were determined for computational convenience. Individuals received input signals from environmental stimuli through a signaling pathway that functioned via trans-elements that activated the expression of genes with the corresponding cis-elements. In this model, each trans-element from an environmental stimulus also had a random trans-number (a t r a n s , an integral number between 1 and 200, the same as the transcription factors described above) and a constant interaction coefficient (v t r a n s , +1.0).
All the mutations in our simulation were point mutations on cis- or trans-elements, and they could be categorized into four types of changes: cis-number, cis-coefficient, trans-number, and trans-coefficient. Each type had the same mutation rate (μ), and the default value was set to 10−4 per individual per generation, which was in balance with the population size N (described below).
Phenotype and fitness
The level of gene expression reached a stable state after the GRN was activated by environmental signals. The phenotype of an individual was defined as the combination of the stable-state expression levels of four phenotypic genes. The phenotypic value of an individual was represented by a point in four-dimensional space with four axes, each of which denoted the expression level of a phenotypic gene. The individuals whose expression levels were very low (<0.01) or did not reach stable states after 5,000 loops of interaction were assumed to be dead. These values were determined for computational convenience.
where di,j is the Euclidean distance from the optimum defined by the environment j and the phenotypic value of the individual i, and σ is the parameter that determined the breadth of the foot of the fitness landscape. For simplicity of explanation, we used S=−l n2σ2 to represent the strength of selection and set it to 3.0 by default. When S was larger, the selection gradient became steeper and small changes in phenotypic values were subjected to stronger selective pressure, i.e., the mutations that produced phenotypic variations were supposed to be purged more quickly from the population.
Analysis of factors affecting the number of cryptic variations
Range and default parameters used in the simulations
Mutation rate [gene−1 generation−1]
V g e n o m e
Number of genes
Number of cis-regulatory elements
Positive edge bias
Number of environments
Strength of selection
The number of genes (V g e n o m e ), the number of cis-regulatory elements per gene (L), and positive edge ratio (A) were set as the initial parameter values given their effects on network properties such as order, weighted size, and density. The population size (N), mutation rate (μ), strength of stabilizing selection (S), and number of environments (H) were also set as the initial parameter values. The effects of each of these properties on the number of CGVs were examined by modifying them. The values for these parameters are shown in Table 1. Simulations were conducted 20 times for each parameter. The number of standing genetic variations (G), visible phenotypic variations (P n o r m a l ), and cryptic variations (P c r y p t i c ) was measured after the simulations (see the previous subsection for their definitions). We obtained regression coefficients from a GLM in which G, P n o r m a l , or P c r y p t i c was the response variable, and each of the parameters was an explanatory variable.
The structural properties of biological networks are considered to be responsible for biological processes such as environmental response and cell differentiation [44, 45]. The following network properties were measured using NetworkX 1.7, a Python language software package .
Gene duplication was not assumed in our model, and the number of genes that an individual had was determined by a parameter (V g e n o m e ). However, not all the genes were responsible for phenotypic expression in our model because how the genes interacted with each other depended on the randomly variable cis-regulatory elements of the genes. Thus, we defined core genes as the genes in a GRN that were connected to phenotypic genes and the network order (V c o r e ) as the number of core genes.
Network size and density
If the two genes showed a regulatory interaction, we regarded them as being neighbors. The number of neighbors of gene i was defined as degree k i . In our model, the expected density values were dependent on the number of genes and cis-regulatory elements of each gene.
Average clustering coefficient
The average clustering coefficient of a network tends to increase when the network has more locally dense compartments, i.e., clusters. Cancerous networks of oncogenes in humans are known to form a giant, highly clustered component in the whole signaling network .
Number of self-loops
When a gene product regulates its own expression, it is a form of autoregulation and is denoted by a loop in a network called a self-loop. It is known that negative feedback stabilizes its expression level and shortens the response time, while positive feedback delays the response and results in bistability . A mutagenesis experiment with yeast confirmed that genes with negative feedback loops are mutationally robust .
Positive edge ratio
Whether the expression of a gene was activated or repressed by the binding of a transcription factor was determined by the signs of their interaction coefficients. The expected value of the ratio of positive edges to negative edges was dependent on the bias parameter A: however, it randomly varied during the initial state and was changed by the evolution during the simulation. We measured the realized ratio of positive edges to negative edges at the end of the simulation.
Analysis of factors affecting network properties
Network properties evolved during the simulations. We examined the effects of the number of genes, number of cis-regulatory elements, mutation rate, population size, strength of stabilizing selection, and number of environments on the properties of the evolved networks by altering these parameter values. Simulations were conducted 20 times for each parameter, and the average values of the network properties were measured thereafter. A GLM was then created with each parameter as a response variable and the average network property as an explanatory variable.
Cryptic variations in GRNs
The effect of GRN properties on the accumulation of cryptic variations
Effects of population size, mutation rate, number of genes, number of cis-regulatory elements, and positive edge bias on cryptic variations
The effects of each parameter on population variation
Rank of P
P n o r m a l
P c r y p t i c
P n o r m a l
P c r y p t i c
V g e n o m e
P n o r m a l
V g e n o m e
P c r y p t i c
V g e n o m e
P n o r m a l
P c r y p t i c
P n o r m a l
P c r y p t i c
P n o r m a l
P c r y p t i c
P n o r m a l
P c r y p t i c
Effects of selective pressure and the number of within-generation environments on the accumulation of cryptic variations
Effects of altered parameters on network characteristics
The effects of varying the parameter values on the network properties of evolved populations were also analyzed Additional file 2: (Table S2). As expected, the mutation rate (μ) and population size (N) did not affect the network properties of the evolved populations. The network size significantly increased with increasing numbers of genes (V g e n o m e ) and cis-regulatory elements (L). The increasing number of genes (V g e n o m e ) caused an increase in network order and a decrease in density, whereas the increasing number of cis-regulatory elements (L) caused an increase in network density.
Factors affecting the number of cryptic variations
The accumulation and release of CGVs were predicted to be the outcome of epistasis and genotype–environment interactions by analytical population genetic modeling . Hermisson and Wagner (2004) assumed sexual reproduction and argued that the sex dependence of allelic effects is an important mechanism for conditional neutrality. In the context of a recent categorization , their model corresponds to “post-mutation evolvability under high recombination,” whereas ours corresponds to “pre-mutation evolvability under low recombination.” Here we showed that such epistatic behavior could be produced as a generic feature of GRNs even without sexual reproduction and without an explicit coefficient for hidden variation as expected from previous stochastic models [29, 34, 53–56]. The present study used individual-based models with explicit implementation of a population of GRNs that produce phenotypes responding to environmental stimuli, and our findings indicate that the cryptic phenotypic variations unveiled by a novel environmental stimulus can accumulate and be maintained in evolved populations. The signals induced by a novel environmental stimulus probably initiated different cascades of gene interactions, which led to novel phenotypic variations due to standing genetic variations that were invisible under normal circumstances.
In the present study, we used two methods to examine the effects of factors that affect the number of CGVs. First, because genetic variations and some network properties varied among the simulation replicates for the same parameter set values, we examined the effects of varying values of these properties. Second, we altered the parameter values and examined the effects of varying parameter values on the CGVs. In our simulations, the number of genes, V g e n o m e , was not assumed to be altered during the evolutionary process; therefore, the effect of V g e n o m e could be examined by altering the initial values of the parameters.
The number of CGVs accumulated in each of the populations increased with increases in the total numbers of genetic polymorphisms (Figure 4H) when the simulation replicates with the default parameters were analyzed. This is reasonable because the more polymorphic genotypes a population has, the more diverse phenotypes it is supposed to produce. This finding was also supported by the results of the simulations where only the population size and mutation rate were altered (Table 2). In addition, increasing the network order by altering the number of genes increased the number of total and cryptic genetic variations even under conditions where the mutation rate per individual (not per gene) was kept constant.
Under the condition of a fixed number of genes, the number of CGVs increased with increasing the weighted size of the GRNs (Figure 4C). The weighted size of a GRN was expected to increase when increasing the order, size, and ratio of the positive edges. In the second analysis, increasing the number of genes caused increases in the order and size of the GRNs followed by an increase in the number of CGVs. In addition, it was confirmed that increasing the ratio of positive edges by altering the positive edge bias caused increases in the weighted size and the number of CGVs in the GRNs. The genetic difference among networks with large weighted sizes may have been amplified through genetic interactions and resulted in larger differences in expression outcomes, which would have contributed to the emergence of CGVs.
These findings indicate that the number of genes could affect the number of CGVs in two ways. First, a population of GRNs with a large number of genes would be able to accumulate a larger number of genetic variations. Second, a higher network order would have a larger weighted size, which may contribute to the emergence of more phenotypic variations. Thus, we propose that the network order (i.e., the number of genes) is the most important factor in facilitating the accumulation and release of CGVs.
Effect of strength of selection and environmental variations
The accumulation of CGVs was affected by the strength of the selective pressure and the number of environments that the individuals experienced during their lifetime. Increasing the strength of selection and the number of environments decreased the genetic variation to a similar extent; however, the effects of these two environmental factors on the number of CGVs differed. This indicates that the quality of accumulated genetic variations may differ depending on the ancestral environments that the populations experienced through the generations.
We predicted that strong selection usually purges genetic variation within populations; therefore, increasing the strength of stabilizing selection decreases the number of genetic and phenotypic variations. The present findings showed that this prediction was only true for visible variations and not for cryptic variations. The number of CGVs increased despite the strong selective pressure. Not all mutations in the GRNs led to phenotypic differences in a normal environment. The mutations that caused phenotypically visible differences were purged by strong pressure, whereas those that did not affect the phenotype were free from selection and accumulated in populations through mutation-drift balance. Therefore, the genetic variations that did not contribute to the phenotypes in the normal environments were preserved in the populations regardless of the strength of the selective pressure.
The number of environments an individual was exposed to during a lifetime also decreased the accumulation of genetic variations. Here, the individuals were required to express different phenotypes in response to different environmental signals experienced during their lifetimes using appropriate regulatory pathways. When an individual experienced a larger number of environments, more of the GRN should be involved in producing phenotypic differences than that in normal environments. Therefore, exposure to a larger number of different environments would probably decrease the number of mutations accumulated through genetic drift. In other words, populations with high plasticity may have less evolvability than those with less plasticity, which would result in less phenotypic diversity in response to novel environments. An experimental study indicated that populations that evolved in rapidly changing environments may result in genotypes with high levels of pleiotropy and historical constraints, which in tern could impede the future evolution . This finding indicates that there may be intrinsic genetic costs associated with adaptation to variable environments [58–60]. These costs differ from other types of costs associated with plasticity because they cannot be directly detected by negative correlations between individual fitness and plasticity but can instead emerge as reduced genetic variability on an evolutionary timescale. However, some predictions [54–56] contradict this. One possible reason for this discrepancy is a difference in focus. Espinosa-Soto et al.  showed that perturbing the initial state of each GRN every generation shortened the time required to discover the genotype of an arbitrary novel phenotype. Fierst  measured the evolutionary variability of genetic systems toward novel phenotypic optima without changing the environmental cue. Draghi and Whitlock  emphasized the directionality of plasticity and evolvability. Here we investigated the quantitative tendencies of genetic and phenotypic variations and detected a negative effect of environmental heterogeneity on the number of variations present in a population.
On the other hand, environmental heterogeneity would be expected to have positive effects on the number of CGVs in the light of GRN evolution. The more heterogeneous the environment that individuals experienced during their lifetimes, the larger and denser the GRNs tended to become Additional file 2: (Table S2). This was probably because only small parts of the GRNs were used in low heterogeneity environments and then unused genes were then allowed to mutate and become disconnected from the GRNs without affecting their fitness. In other words, populations that evolve in static environments over long periods may lose the ability to sense and respond to the environment, and would therefore show less revealed phenotypic variations in novel environments. A previous theoretical study  showed that the evolution of complex GRNs could be remarkably promoted by the fixation of beneficial gene duplications under randomly fluctuating environmental conditions and that such GRNs tend to exhibit high mutational robustness and evolvability. Furthermore, they showed that large and complex GRNs could not be evolved in cyclically fluctuating environments. This indicates that predictable variable environments, including the heterogeneous environments that we used here, cannot promote the evolution of higher evolvability as well as of larger numbers of CGVs. Therefore, interactions with variable environments in which higher plasticity has evolved may promote the accumulation of CGVs by facilitating the evolution of larger GRNs in some cases. In turn, the expansion of GRNs could facilitate evolutionary adaptation to novel environments and niche construction. Thus, we propose that the interaction between GRNs and variable environments could be a sustainable source of evolvability.
Acquisition of novel traits promoted by cryptic genetic variations in GRNs
It is believed that CGVs contribute to evolutionary responses to environmental changes and to long-term selection by providing phenotypic variations in response to changes in the background [8, 9, 61]. Furthermore, some researchers argue that CGVs also promote the acquisition of novel traits because of their ability to accumulate and release multiple mutations in individuals and populations [62, 63].
Several studies support these arguments. For example, the black mutant strain of Manduca sexta, which was originally green, showed variations in thermosensitivity where some larvae turned green and the others stayed black in response to heat shock over a specific period . The authors of that study also artificially selected these variations and successfully established two lines: one for sensitivity and the other for insensitivity. CGVs that accumulate during strong stabilizing selection could contribute to phenotypic evolution in this manner.
In addition, a previous model study  predicted that GRNs are robust and evolvable enough to cover the broad genotypic space and to reach genotypes that produce novel phenotypes through recurrent neutral mutations. Our individual-based model revealed that CGVs could accumulate and produce phenotypic diversity to contribute to evolvability in the context of population genetic processes. It has been proposed that novel phenotypes do not necessarily require new genes. Rather, changes in the expression patterns of existing genes are important for phenotypic evolution [65–68]. Further simulations with properly defined novel phenotypes will help elucidate the evolution of novelty.
We constructed an individual-based model of GRNs that controlled gene expression in response to environmental stimuli. The model facilitated the analysis of network properties in the context of population genetics. It showed that populations of GRNs accumulate and release cryptic variations, the number of which varies depending on the properties of the GRNs and the environments to which they have been subjected across the generations. Our findings indicate that the expansion of GRNs and adaptation to novel environments are mutually facilitating, resulting in a sustainable sources of evolvability. This study thus provides important insight into the origins of biological diversity and complexity.
Cryptic genetic variation.
We are grateful to Takashi Makino and Yukari Ohno for their valuable comments on this study. We thank Satoshi Kato, Wakako Yamaguchi, Mitsuhiko Sato, and Satoshi Tamate for their assistance with the maintenance of the computer environment. WMI is supported by a Japan Society for the Promotion of Science (JSPS) Research Fellowship. This research was supported by a Grant-in-Aid for JSPS Fellows (10J56112) to WMI and partially by the Global COE Program “Center for ecosystem management adapting to global change” (J03) of the Ministry of Education, Culture, Sports, Science and Technology of Japan granted to MK.
- Lande R: Quantitative genetic analysis of multivariate evolution, applied to brain: body size allometry. Evolution. 1979, 33: 402-416. 10.2307/2407630.View Article
- Hansen TF, Houle D: Measuring and comparing evolvability and constraint in multivariate characters. J Evol Biol. 2008, 21 (5): 1201-1219. 10.1111/j.1420-9101.2008.01573.x.PubMedView Article
- Steppan S, Phillips P, Houle D: Comparative quantitative genetics: evolution of the G matrix. Trends Ecol Evol. 2002, 17 (7): 320-327. 10.1016/S0169-5347(02)02505-3.View Article
- Barton N, Keightley P: Understanding quantitative genetic variation. Nat Rev Genet. 2002, 3: 11-21.PubMedView Article
- Crow JF: Maintaining evolvability. J Genet. 2008, 87 (4): 349-353. 10.1007/s12041-008-0057-8.PubMedView Article
- Dudley JW: From means to QTL: The Illinois long-term selection experiment as a case study in quantitative genetics. Crop Science. 2007, 47: S20-S31.View Article
- Dunnington E, Siegel P: Long-term selection for 8-week body weight in chickens–direct and correlated responses. Theor Appl Genet. 1985, 71 (2): 305-313.PubMed
- Le Rouzic A, Carlborg O: Evolutionary potential of hidden genetic variation. Trends Ecol Evol. 2008, 23: 33-37. 10.1016/j.tree.2007.09.014.PubMedView Article
- Schlichting CD: Hidden reaction norms, cryptic genetic variation, and evolvability. Ann N Y Acad Sci. 2008, 1133: 187-203. 10.1196/annals.1438.010.PubMedView Article
- Lande R: Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation. J Evol Biol. 2009, 22 (7): 1435-1446. 10.1111/j.1420-9101.2009.01754.x.PubMedView Article
- Lauter N, Doebley J: Genetic variation for phenotypically invariant traits detected in teosinte: implications for the evolution of novel forms. Genet. 2002, 160: 333-342.
- Kitano H: Biological robustness. Nat Rev Genet. 2004, 5 (11): 826-837.PubMedView Article
- Rutherford SL, Lindquist S: Hsp90 as a capacitor for morphological evolution. Nature. 1998, 396 (6709): 336-342. 10.1038/24550.PubMedView Article
- Bergman A, Siegal ML: Evolutionary capacitance as a general feature of complex gene networks. Nature. 2003, 424 (6948): 549-552. 10.1038/nature01765.PubMedView Article
- Masel J, Siegal ML: Robustness: mechanisms and consequences. Trends Genet. 2009, 25 (9): 395-403. 10.1016/j.tig.2009.07.005.PubMedView Article
- Masel J, Trotter MV: Robustness and evolvability. Trends Genet. 2010, 26 (9): 406-14. 10.1016/j.tig.2010.06.002.PubMedView Article
- Takahashi KH: Multiple capacitors for natural genetic variation in Drosophila melanogaster. Mol Ecol. 2013, 22 (5): 1356-1365. 10.1111/mec.12091.PubMedView Article
- Wagner A: Distributed robustness versus redundancy as causes of mutational robustness. Bioessays. 2005, 27 (2): 176-188. 10.1002/bies.20170.PubMedView Article
- Hermisson J, Wagner GP: The population genetic theory of hidden variation and genetic robustness. Genetics. 2004, 168 (4): 2271-2284. 10.1534/genetics.104.029173.PubMedView Article
- Carlborg O, Jacobsson L, Ahgren P, Siegel P, Andersson L: Epistasis and the release of genetic variation during long-term selection. Nat Genet. 2006, 38 (4): 418-420. 10.1038/ng1761.PubMedView Article
- Farkas IJ, Wu C, Chennubhotla C, Bahar I, Oltvai ZN: Topological basis of signal integration in the transcriptional-regulatory network of the yeast, Saccharomyces cerevisiae. BMC Bioinformatics. 2006, 7: 478-10.1186/1471-2105-7-478.PubMedView Article
- Chen HC, Lee HC, Lin TY, Li WH, Chen BS: Quantitative characterization of the transcriptional regulatory network in the yeast cell cycle. Bioinformatics. 2004, 20 (12): 1914-1927. 10.1093/bioinformatics/bth178.PubMedView Article
- Davidson EH: The Regulatory Genome: Gene Regulatory Networks In Development And Evolution. 2006, Burlington, MA: Academic Press
- Erwin DH, Davidson EH: The evolution of hierarchical gene regulatory networks. Nat Rev Genet. 2009, 10 (2): 141-148.PubMedView Article
- Siegal ML, Promislow DEL, Bergman A: Functional and evolutionary inference in gene networks: does topology matter?. Genetica. 2007, 129: 83-103.PubMedView Article
- Wray GA: The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007, 8 (3): 206-216.PubMedView Article
- Kauffman SA: Metabolic stability and epigenesis in randomly constructed genetic nets. J Theor Biol. 1969, 22 (3): 437-467. 10.1016/0022-5193(69)90015-0.PubMedView Article
- Kauffman SA: The Origins of Order: Self-Organization and Selection in Evolution. 1993, New York: Oxford University Press
- Gjuvsland AB, Hayes BJ, Omholt SW, Carlborg O: Statistical epistasis is a generic feature of gene regulatory networks. Genetics. 2007, 175: 411-420.PubMedView Article
- Omholt S, Plahte E, Oyehaug L, Xiang K: Gene regulatory networks generating the phenomena of additivity, dominance and epistasis. Genetics. 2000, 155 (2): 969-980.PubMed
- Boldogkoi Z: Genetic background of phenotypic variation. Prog Nat Sci. 2007, 17 (10): 1119-1126.
- Weitz JS, Benfey PN, Wingreen NS: Evolution, interactions, and biological networks. PLoS Biol. 2007, 5: 10-12. 10.1371/journal.pbio.0050010.View Article
- Boldogkoi Z: Gene network polymorphism is the raw material of natural selection: the selfish gene network hypothesis. J Mol Evol. 2004, 59 (3): 340-357. 10.1007/s00239-004-2629-4.PubMedView Article
- Ciliberti S, Martin OC, Wagner A: Innovation and robustness in complex regulatory gene networks. Proc Natl Acad Sci U S A. 2007, 104 (34): 13591-13596. 10.1073/pnas.0705396104.PubMedView Article
- von Dassow G, Meir E, Munro EM, Odell GM: The segment polarity network is a robust developmental module. Nature. 2000, 406 (6792): 188-192. 10.1038/35018085.PubMedView Article
- Tsuda ME, Kawata M: Evolution of gene regulatory networks by fluctuating selection and intrinsic constraints. PLoS Comput Biol. 2010, 6 (8): e1000873-10.1371/journal.pcbi.1000873.PubMedView Article
- Ciliberti S, Martin OC, Wagner A: Robustness can evolve gradually in complex regulatory gene networks with varying topology. PLoS Comput Biol. 2007, 3 (2): e15-10.1371/journal.pcbi.0030015.PubMedView Article
- Lynch M: The origins of genome architecture. Sunderland. 2007, MA: Sinauer Associates
- Alon U: An introduction to systems biology: design principles of biological circuits. 2007, Boca Raton, FL: Chapman & Hall/CRC
- Andrecut M, Huang S, Kauffman SA: Heuristic approach to sparse approximation of gene regulatory networks. J Comput Biol. 2008, 15 (9): 1173-1186. 10.1089/cmb.2008.0087.PubMedView Article
- Britten RJ, Davidson EH: Gene regulation for higher cells: a theory. Science. 1969, 165 (3891): 349-357. 10.1126/science.165.3891.349.PubMedView Article
- Rabani M, Levin JZ, Fan L, Adiconis X, Raychowdhury R, Garber M, Gnirke A, Nusbaum C, Hacohen N, Friedman N, Amit I, Regev A: Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat Biotechnol. 2011, 29 (5): 436-442. 10.1038/nbt.1861.PubMedView Article
- R Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2012, Austria: Vienna, [http://www.R-project.org/]
- Davidson EH, Levine MS: Properties of developmental gene regulatory networks. Proc Natl Acad Sci U S A. 2008, 105 (51): 20063-20066. 10.1073/pnas.0806007105.PubMedView Article
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298 (5594): 824-827. 10.1126/science.298.5594.824.PubMedView Article
- Hagberg AA, Schult DA, Swart PJ: Exploring network structure, dynamics, and function using NetworkX. Proceedings of the 7th Python in Science Conference (SciPy2008). 2008, Pasadena, CA USA: SciPy.org, 11-15.
- Cui Q, Ma Y, Jaramillo M, Bari H, Awan A, Yang S, Zhang S, Liu L, Lu M, O’Connor-McCourt M, Purisima EO, Wang E: A map of human cancer signaling. Mol Syst Biol. 2007, 3: 152-PubMedView Article
- Denby CM, Im JH, Yu RC, Pesce CG, Brem RB: Negative feedback confers mutational robustness in yeast transcription factor regulation. Proc Natl Acad Sci U S A. 2012, 109 (10): 3874-3878. 10.1073/pnas.1116360109.PubMedView Article
- Newman M: Assortative mixing in networks. Phys Rev Lett. 2002, 89 (20): 208701-PubMedView Article
- Maslov S, Sneppen K: Specificity and stability in topology of protein networks. Science. 2002, 296 (5569): 910-913. 10.1126/science.1065103.PubMedView Article
- Payne JL, Dodds PS, Eppstein MJ: Information cascades on degree-correlated random networks. Phys Rev E. 2009, 80 (2): 026125-View Article
- Guo W, Rao MB: On control of the false discovery rate under no assumption of dependency. J Stat Plann Inference. 2008, 138 (10): 3176-3188. 10.1016/j.jspi.2008.01.003.View Article
- Hodgins-Davis A, Townsend JP: Evolving gene expression: from G to E to GxE. Trends Ecol Evol. 2009, 24 (12): 649-658. 10.1016/j.tree.2009.06.011.PubMedView Article
- Espinosa-Soto C, Martin OC, Wagner A: Phenotypic plasticity can facilitate adaptive evolution in gene regulatory circuits. BMC Evol Biol. 2011, 11: 5-10.1186/1471-2148-11-5.PubMedView Article
- Fierst JL: A history of phenotypic plasticity accelerates adaptation to a new environment. J Evol Biol. 2011, 24 (9): 1992-2001. 10.1111/j.1420-9101.2011.02333.x.PubMedView Article
- Draghi JA, Whitlock MC: Phenotypic plasticity facilitates mutational variance, genetic variance, and evolvability along the major axis of environmental variation. Evolution. 2012, 66 (9): 2891-2902. 10.1111/j.1558-5646.2012.01649.x.PubMedView Article
- Collins S, de Meaux J: Adaptation to different rates of environmental change in Chlamydomonas. Evolution. 2009, 63 (11): 2952-2965. 10.1111/j.1558-5646.2009.00770.x.PubMedView Article
- DeWitt T, Sih A, Wilson D: Costs and limits of phenotypic plasticity. Trends Ecol Evol. 1998, 13 (2): 77-81. 10.1016/S0169-5347(97)01274-3.PubMedView Article
- van Kleunen M, Fischer M: Constraints on the evolution of adaptive phenotypic plasticity in plants. New Phytol. 2005, 166: 49-60. 10.1111/j.1469-8137.2004.01296.x.PubMedView Article
- Auld JR, Agrawal AA, Relyea RA: Re-evaluating the costs and limits of adaptive phenotypic plasticity. Proc Biol Sci. 2010, 277 (1681): 503-511. 10.1098/rspb.2009.1355.PubMedView Article
- Wagner A: Robustness and evolvability: a paradox resolved. Proc Biol Sci. 2008, 275 (1630): 91-100. 10.1098/rspb.2007.1137.PubMedView Article
- West-Eberhard MJ: Developmental plasticity and the origin of species differences. Proc Natl Acad Sci U S A. 2005, 102 (Suppl 1): 6543-6549.PubMedView Article
- Pigliucci M, Murren CJ: Genetic assimilation and a possible evolutionary paradox: can macroevolution sometimes be so fast as to pass us by?. Evolution. 2003, 57 (7): 1455-1464.PubMedView Article
- Suzuki Y, Nijhout HF: Evolution of a polyphenism by genetic accommodation. Science. 2006, 311 (5761): 650-652. 10.1126/science.1118888.PubMedView Article
- Shubin N, Tabin C, Carroll S: Deep homology and the origins of evolutionary novelty. Nature. 2009, 457 (7231): 818-23. 10.1038/nature07891.PubMedView Article
- Prud’homme B, Gompel N, Carroll SB: Emerging principles of regulatory evolution. Proc Natl Acad Sci U S A. 2007, 104 (Suppl 1): 8605-8612.PubMedView Article
- Carroll SB, Grenier JK, Weatherbee SD: From DNA to diversity: molecular genetics and the evolution of animal design. Oxford. 2004, UK: Blackwell
- Shubin N, Tabin C, Carroll S: Fossils, genes and the evolution of animal limbs. Nature. 1997, 388 (6643): 639-648. 10.1038/41710.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.