RESUMEN
The contribution of non-additive genetic effects to the genetic architecture of fitness and to the evolutionary potential of populations has been a topic of theoretical and empirical interest for a long time. Yet, the empirical study of these effects in natural populations remains scarce, perhaps because measuring dominance and epistasis relies heavily on experimental line crosses. In this study, we explored the contribution of dominance and epistasis in natural alpine populations of Arabidopsis thaliana for 2 fitness traits, the dry biomass and the estimated number of siliques, measured in a greenhouse. We found that, on average, crosses between inbred lines of A. thaliana led to mid-parent heterosis for dry biomass but outbreeding depression for an estimated number of siliques. While heterosis for dry biomass was due to dominance, we found that outbreeding depression for an estimated number of siliques could be attributed to the breakdown of beneficial epistatic interactions. We simulated and discussed the implication of these results for the adaptive potential of the studied populations, as well as the use of line-cross analyses to detect non-additive genetic effects.
Asunto(s)
Arabidopsis , Epistasis Genética , Arabidopsis/genética , Vigor Híbrido , Cruzamientos Genéticos , Biomasa , Modelos GenéticosRESUMEN
Physiological status can influence social behavior, which in turn can affect physiology and health. Previously, we reported that tumor growth in Drosophila virgin females depends on the social context, but did not investigate the underlying physiological mechanisms. Here, we sought to characterize the signal perceived between tumorous flies, ultimately discovering that the tumor suppressive effect varies depending on reproductive status. Firstly, we show that the tumor suppressive effect is neither dependent on remnant pheromone-like products nor on the microbiota. Transcriptome analysis of the heads of these tumorous flies reveals social-dependent gene-expression changes related to nervous-system activity, suggesting that a cognitive-like relay might mediate the tumor suppressive effect. The transcriptome also reveals changes in the expression of genes related to mating behavior. Surprisingly, we observed that this social-dependent tumor-suppressive effect is lost in fertilized females. After mating, Drosophila females change their behavior-favoring offspring survival-in response to peptides transferred via the male ejaculate, a phenomenon called "male manipulation". Remarkably, the social-dependent tumor suppressive effect is restored in females mated by sex-peptide deficient males. Since male manipulation has likely been selected to favor male gene transmission, our findings indicate that this evolutionary trait impedes social-dependent tumor growth slowdown.
Asunto(s)
Proteínas de Drosophila , Drosophila melanogaster , Animales , Femenino , Masculino , Drosophila melanogaster/metabolismo , Proteínas de Drosophila/metabolismo , Conducta Sexual Animal/fisiología , Péptidos/metabolismo , ReproducciónRESUMEN
The evolution of gene expression is constrained by the topology of gene regulatory networks, as co-expressed genes are likely to have their expressions affected together by mutations. Conversely, co-expression can also be an advantage when genes are under joint selection. Here, we assessed theoretically whether correlated selection (selection for a combination of traits) was able to affect the pattern of correlated gene expressions and the underlying gene regulatory networks. We ran individual-based simulations, applying a stabilizing correlated fitness function to three genetic architectures: a quantitative genetics (multilinear) model featuring epistasis and pleiotropy, a quantitative genetics model where each genes has an independent mutational structure, and a gene regulatory network model, mimicking the mechanisms of gene expression regulation. Simulations showed that correlated mutational effects evolved in the three genetic architectures as a response to correlated selection, but the response in gene networks was specific. The intensity of gene co-expression was mostly explained by the regulatory distance between genes (largest correlations being associated to genes directly interacting with each other), and the sign of co-expression was associated with the nature of the regulation (transcription activation or inhibition). These results concur to the idea that gene network topologies could partly reflect past selection patterns on gene expression.
Asunto(s)
Redes Reguladoras de Genes , Modelos Genéticos , Regulación de la Expresión Génica , Fenotipo , Mutación , Selección GenéticaRESUMEN
Transposable elements (TEs) are self-reproducing selfish DNA sequences that can invade the genome of virtually all living species. Population genetics models have shown that TE copy numbers generally reach a limit, either because the transposition rate decreases with the number of copies (transposition regulation) or because TE copies are deleterious, and thus purged by natural selection. Yet, recent empirical discoveries suggest that TE regulation may mostly rely on piRNAs, which require a specific mutational event (the insertion of a TE copy in a piRNA cluster) to be activated - the so-called TE regulation "trap model". We derived new population genetics models accounting for this trap mechanism, and showed that the resulting equilibria differ substantially from previous expectations based on a transposition-selection equilibrium. We proposed three sub-models, depending on whether or not genomic TE copies and piRNA cluster TE copies are selectively neutral or deleterious, and we provide analytical expressions for maximum and equilibrium copy numbers, as well as cluster frequencies for all of them. In the full neutral model, the equilibrium is achieved when transposition is completely silenced, and this equilibrium does not depend on the transposition rate. When genomic TE copies are deleterious but not cluster TE copies, no long-term equilibrium is possible, and active TEs are eventually eliminated after an active incomplete invasion stage. When all TE copies are deleterious, a transposition-selection equilibrium exists, but the invasion dynamics is not monotonic, and the copy number peaks before decreasing. Mathematical predictions were in good agreement with numerical simulations, except when genetic drift and/or linkage disequilibrium dominates. Overall, the trap-model dynamics appeared to be substantially more stochastic and less repeatable than traditional regulation models.
Asunto(s)
Elementos Transponibles de ADN , ARN de Interacción con Piwi , Elementos Transponibles de ADN/genética , Genética de Población , Selección Genética , Mutación , Evolución MolecularRESUMEN
The domestication of plant species leads to repeatable morphological evolution, often referred to as the phenotypic domestication syndrome. Domestication is also associated with important genomic changes, such as the loss of genetic diversity compared with adequately large wild populations, and modifications of gene expression patterns. Here, we explored theoretically the effect of a domestication-like scenario on the evolution of gene regulatory networks. We ran population genetics simulations in which individuals were featured by their genotype (an interaction matrix encoding a gene regulatory network) and their gene expressions, representing the phenotypic level. Our domestication scenario included a population bottleneck and a selection switch mimicking human-mediated directional and canalizing selection, i.e., change in the optimal gene expression level and selection toward more stable expression across environments. We showed that domestication profoundly alters genetic architectures. Based on four examples of plant domestication scenarios, our simulations predict (1) a drop in neutral allelic diversity; (2) a change in gene expression variance that depends upon the domestication scenario; (3) transient maladaptive plasticity; (4) a deep rewiring of the gene regulatory networks, with a trend toward gain of regulatory interactions; and (5) a global increase in the genetic correlations among gene expressions, with a loss of modularity in the resulting coexpression patterns and in the underlying networks. We provide empirically testable predictions on the differences of genetic architectures between wild and domesticated forms. The characterization of such systematic evolutionary changes in the genetic architecture of traits contributes to define a molecular domestication syndrome.
Asunto(s)
Domesticación , Redes Reguladoras de Genes , Genética de Población , Genotipo , Humanos , Fenotipo , Selección GenéticaRESUMEN
In many studies, related individuals are phenotyped in order to infer how their genotype contributes to their phenotype, through the estimation of parameters such as breeding values or locus effects. When it is not possible to phenotype all the individuals, it is important to properly sample the population to improve the precision of the statistical analysis. This article studies how to optimize such sampling designs for pedigrees and association studies. Two sampling methods are developed, stratified sampling and D optimality. It is found that it is important to take account of mutation when sampling pedigrees with many generations: as the size of mutation effects increases, optimized designs sample more individuals in late generations. Optimized designs for association studies tend to improve the joint estimation of breeding values and locus effects, all the more as sample size is low and the genetic architecture of the trait is simple. When the trait is determined by few loci, they are reminiscent of classical experimental designs for regression models and tend to select homozygous individuals. When the trait is determined by many loci, locus effects may be difficult to estimate, even if an optimized design is used.
Asunto(s)
Modelos Genéticos , Sitios de Carácter Cuantitativo , Genotipo , Linaje , FenotipoRESUMEN
Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.
Asunto(s)
Adaptación Fisiológica , Variación Genética , Adaptación Fisiológica/genética , Evolución Biológica , Redes Reguladoras de Genes , Genotipo , FenotipoRESUMEN
Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.
Asunto(s)
Flujo Genético , Tasa de Mutación , Selección Genética , Zea mays/genética , Aptitud Genética , Polinización , Zea mays/fisiologíaRESUMEN
Glycolysis and fatty acid (FA) synthesis directs the production of energy-carrying molecules and building blocks necessary to support cell growth, although the absolute requirement of these metabolic pathways must be deeply investigated. Here, we used Drosophila genetics and focus on the TOR (Target of Rapamycin) signaling network that controls cell growth and homeostasis. In mammals, mTOR (mechanistic-TOR) is present in two distinct complexes, mTORC1 and mTORC2; the former directly responds to amino acids and energy levels, whereas the latter sustains insulin-like-peptide (Ilp) response. The TORC1 and Ilp signaling branches can be independently modulated in most Drosophila tissues. We show that TORC1 and Ilp-dependent overgrowth can operate independently in fat cells and that ubiquitous over-activation of TORC1 or Ilp signaling affects basal metabolism, supporting the use of Drosophila as a powerful model to study the link between growth and metabolism. We show that cell-autonomous restriction of glycolysis or FA synthesis in fat cells retrains overgrowth dependent on Ilp signaling but not TORC1 signaling. Additionally, the mutation of FASN (Fatty acid synthase) results in a drop in TORC1 but not Ilp signaling, whereas, at the cell-autonomous level, this mutation affects none of these signals in fat cells. These findings thus reveal differential metabolic sensitivity of TORC1- and Ilp-dependent growth and suggest that cell-autonomous metabolic defects might elicit local compensatory pathways. Conversely, enzyme knockdown in the whole organism results in animal death. Importantly, our study weakens the use of single inhibitors to fight mTOR-related diseases and strengthens the use of drug combination and selective tissue-targeting.
Asunto(s)
Proliferación Celular , Proteínas de Drosophila/metabolismo , Cuerpo Adiposo/metabolismo , Insulinas/metabolismo , Diana Mecanicista del Complejo 1 de la Rapamicina/metabolismo , Transducción de Señal , Animales , Metabolismo Basal , Drosophila melanogaster , Cuerpo Adiposo/citología , Cuerpo Adiposo/fisiología , Acido Graso Sintasa Tipo I/metabolismo , Ácidos Grasos/biosíntesis , GlucólisisRESUMEN
Anticipating the genetic and phenotypic changes induced by natural or artificial selection requires reliable estimates of trait evolvabilities (genetic variances and covariances). However, whether or not multivariate quantitative genetics models are able to predict precisely the evolution of traits of interest, especially fitness-related, life history traits, remains an open empirical question. Here, we assessed to what extent the response to bivariate artificial selection on both body size and maturity in the medaka Oryzias latipes, a model fish species, fits the theoretical predictions. Three lines (Large, Small, and Control lines) were differentially selected for body length at 75 days of age, conditional on maturity. As maturity and body size were phenotypically correlated, this selection procedure generated a bi-dimensional selection pattern on two life history traits. After removal of nonheritable trends and noise with a random effect ("animal") model, the observed selection response did not match the expected bidirectional response. For body size, Large and Control lines responded along selection gradients (larger body size and stasis, respectively), but, surprisingly, the Small did not evolve a smaller body length and remained identical to the Control line throughout the experiment. The magnitude of the empirical response was smaller than the theoretical prediction in both selected directions. For maturity, the response was opposite to the expectation (the Large line evolved late maturity compared to the Control line, while the Small line evolved early maturity, while the opposite pattern was predicted due to the strong positive genetic correlation between both traits). The mismatch between predicted and observed response was substantial and could not be explained by usual sources of uncertainties (including sampling effects, genetic drift, and error in G matrix estimates).
RESUMEN
Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.
RESUMEN
Anthropogenic perturbations such as harvesting often select against a large body size and are predicted to induce rapid evolution toward smaller body sizes and earlier maturation. However, body-size evolvability and, hence, adaptability to anthropogenic perturbations remain seldom evaluated in wild populations. Here, we use a laboratory experiment over 6 generations to measure the ability of wild-caught medaka fish (Oryzias latipes) to evolve in response to bidirectional size-dependent selection mimicking opposite harvest regimes. Specifically, we imposed selection against a small body size (Large line), against a large body size (Small line) or random selection (Control line), and measured correlated responses across multiple phenotypic, life-history, and endocrine traits. As expected, the Large line evolved faster somatic growth and delayed maturation, but also evolved smaller body sizes at hatch, with no change in average levels of pituitary gene expressions of luteinizing, follicle-stimulating, or growth hormones (GH). In contrast, the Small medaka line was unable to evolve smaller body sizes or earlier maturation, but evolved smaller body sizes at hatch and showed marginally significant signs of increased reproductive investment, including larger egg sizes and elevated pituitary GH production. Natural selection on medaka body size was too weak to significantly hinder the effect of artificial selection, indicating that the asymmetric body-size response to size-dependent selection reflected an asymmetry in body-size evolvability. Our results show that trait evolvability may be contingent upon the direction of selection and that a detailed knowledge of trait evolutionary potential is needed to forecast population response to anthropogenic change.
RESUMEN
Dissecting the genetic basis of intraspecific variations in life history traits is essential to understand their evolution, notably for potential biocontrol agents. Such variations are observed in the endoparasitoid Cotesia typhae (Hymenoptera: Braconidae), specialized on the pest Sesamia nonagrioides (Lepidoptera: Noctuidae). Previously, we identified two strains of C. typhae that differed significantly for life history traits on an allopatric host population. To investigate the genetic basis underlying these phenotypic differences, we used a quantitative trait locus (QTL) approach based on restriction site-associated DNA markers. The characteristic of C. typhae reproduction allowed us generating sisters sharing almost the same genetic content, named clonal sibship. Crosses between individuals from the two strains were performed to generate F2 and F8 recombinant CSS. The genotypes of 181 clonal sibships were determined as well as the phenotypes of the corresponding 4,000 females. Informative markers were then used to build a high-quality genetic map. These 465 markers spanned a total length of 1,300 cM and were organized in 10 linkage groups which corresponded to the number of C. typhae chromosomes. Three QTLs were detected for parasitism success and two for offspring number, while none were identified for sex ratio. The QTLs explained, respectively, 27.7% and 24.5% of the phenotypic variation observed. The gene content of the genomic intervals was investigated based on the genome of C. congregata and revealed 67 interesting candidates, as potentially involved in the studied traits, including components of the venom and of the symbiotic virus (bracovirus) shown to be necessary for parasitism success in related wasps.
Asunto(s)
Polydnaviridae , Avispas , Animales , Femenino , Ligamiento Genético , Fenotipo , Polydnaviridae/genética , Sitios de Carácter Cuantitativo/genética , Reproducción/genética , Avispas/genéticaRESUMEN
Microbial symbionts are ubiquitous associates of living organisms but their role in mediating reproductive isolation (RI) remains controversial. We addressed this knowledge gap by employing the Drosophila paulistorum-Wolbachia model system. Semispecies in the D. paulistorum species complex exhibit strong RI between each other and knockdown of obligate mutualistic Wolbachia bacteria in female D. paulistorum flies triggers loss of assortative mating behavior against males carrying incompatible Wolbachia strains. Here we set out to determine whether de novo RI can be introduced by Wolbachia-knockdown in D. paulistorum males. We show that Wolbachia-knockdown D. paulistorum males (i) are rejected as mates by wild type females, (ii) express altered sexual pheromone profiles, and (iii) are devoid of the endosymbiont in pheromone producing cells. Our findings suggest that changes in Wolbachia titer and tissue tropism can induce de novo premating isolation by directly or indirectly modulating sexual behavior of their native D. paulistorum hosts.
Asunto(s)
Drosophila/microbiología , Reproducción/fisiología , Simbiosis/genética , Animales , Conducta Animal , Evolución Biológica , Proteínas de Drosophila/metabolismo , Femenino , Masculino , Aislamiento Reproductivo , Atractivos Sexuales/metabolismo , Atractivos Sexuales/fisiología , Conducta Sexual Animal/fisiología , Especificidad de la Especie , Wolbachia/fisiologíaRESUMEN
Understanding the importance of nongenetic heredity in the evolutionary process is a major topic in modern evolutionary biology. We modified a classical gene-network model by allowing parental transmission of gene expression and studied its evolutionary properties through individual-based simulations. We identified ontogenetic time (i.e. the time gene networks have to stabilize before being submitted to natural selection) as a crucial factor in determining the evolutionary impact of this phenotypic inheritance. Indeed, fast-developing organisms display enhanced adaptation and greater robustness to mutations when evolving in presence of nongenetic inheritance (NGI). In contrast, in our model, long development reduces the influence of the inherited state of the gene network. NGI thus had a negligible effect on the evolution of gene networks when the speed at which transcription levels reach equilibrium is not constrained. Nevertheless, simulations show that intergenerational transmission of the gene-network state negatively affects the evolution of robustness to environmental disturbances for either fast- or slow-developing organisms. Therefore, these results suggest that the evolutionary consequences of NGI might not be sought only in the way species respond to selection, but also on the evolution of emergent properties (such as environmental and genetic canalization) in complex genetic architectures.
Asunto(s)
Evolución Molecular , Redes Reguladoras de Genes/genética , Modelos Biológicos , Animales , Simulación por Computador , Ambiente , Regulación de la Expresión Génica , PadresRESUMEN
In sexual species, mating success depends on the male's capacity to find sexual partners and on female receptivity to mating. Mating is under evolutionary constraints to prevent interspecific mating and to maximize the reproductive success of both sexes. In Drosophila melanogaster, female receptivity to mating is mainly controlled by Sex peptide (SP, i.e. Acp70A) produced by the male accessory glands with other proteins (Acps). The transfer of SP during copulation dramatically reduces female receptivity to mating and prevents remating with other males. To date, female postmating responses are well-known in D. melanogaster but have been barely investigated in closely-related species or strains exhibiting different mating systems (monoandrous versus polyandrous). Here, we describe the diversity of mating systems in two strains of D. melanogaster and the three species of the yakuba complex. Remating delay and sexual receptivity were measured in cross-experiments following SP orthologs or Acp injections within females. Interestingly, we discovered strong differences between the two strains of D. melanogaster as well as among the three species of the yakuba complex. These results suggest that reproductive behavior is under the control of complex sexual interactions between the sexes and evolves rapidly, even among closely-related species.
Asunto(s)
Proteínas de Drosophila/fisiología , Drosophila/fisiología , Péptidos/fisiología , Animales , Femenino , Masculino , Reproducción/fisiología , Conducta Sexual Animal/fisiología , Especificidad de la EspecieRESUMEN
Transposable elements (TEs) are repeated DNA sequences that can constitute a substantial part of genomes. Studying TEs' activity, interactions, and accumulation dynamics is thus of major interest to understand genome evolution. Here, we describe the transposition dynamics of cut-and-paste mariner elements during experimental (short- and longer-term) evolution in Drosophila melanogaster Flies with autonomous and nonautonomous mariner copies were introduced in populations containing no active mariner, and TE accumulation was tracked by quantitative PCR for up to 100 generations. Our results demonstrate that (i) active mariner elements are highly invasive and characterized by an elevated transposition rate, confirming their capacity to spread in populations, as predicted by the "selfish-DNA" mechanism; (ii) nonautonomous copies act as parasites of autonomous mariner elements by hijacking the transposition machinery produced by active mariner, which can be considered as a case of hyperparasitism; (iii) this behavior resulted in a failure of active copies to amplify which systematically drove the whole family to extinction in less than 100 generations. This study nicely illustrates how the presence of transposition-competitive variants can deeply impair TE dynamics and gives clues to the extraordinary diversity of TE evolutionary histories observed in genomes.
Asunto(s)
Elementos Transponibles de ADN , Evolución Molecular Dirigida , Drosophila melanogaster/genética , Genoma de los Insectos , Animales , Proteínas de Unión al ADN/genética , Femenino , Modelos Lineales , Masculino , Fenotipo , Filogenia , Reacción en Cadena de la Polimerasa , Especificidad de la Especie , Transposasas/genéticaRESUMEN
Understanding and predicting evolution is a central challenge in both population and quantitative genetics. The amount of genetic variance for quantitative traits available in a population conditions the particular way in which this population will (or will not) evolve under natural or artificial selection. Here, we explore the potential of gene-gene interactions (epistasis) to induce evolutionary plateaus at which evolutionary change virtually collapses for a number of generations, followed by the release of previously cryptic genetic variation. First, we demonstrate theoretically that a wide range of epistatic interactions has the potential to generate temporary decelerations in the course of response to selection. Second, we perform simulations to show that such microevolutionary plateaus may occur in selection responses under empirically based assumptions. Finally, we show that such events can be traced in artificial selection experiments, thus providing further empirical evidence for this phenomenon.
Asunto(s)
Evolución Biológica , Pollos/genética , Epistasis Genética , Ratones/genética , Selección Genética , Animales , Peso Corporal , Pollos/crecimiento & desarrollo , Variación Genética , Ratones/crecimiento & desarrollo , Modelos GenéticosRESUMEN
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).