Gene network simulations provide testable predictions for the molecular domestication syndrome

Abstract The domestication of plant and animal species lead 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 and modifications of gene expression patterns. Here, we explored theoretically the effect of domestication at the genomic level by characterizing the impact of a domestication-like scenario on 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 (change in the optimal gene expression level) mimicking canalizing selection, i.e. evolution towards more stable expression to parallel enhanced environmental stability in man-made habitat. We showed that domestication profoundly alters genetic architectures. Based on the well-documented example of the maize (Zea mays ssp. mays) domestication, our simulations predicted (i) a drop in neutral allelic diversity, (ii) a change in gene expression variance that depended upon the domestication scenario, (iii) transient maladaptive plasticity, (iv) a deep rewiring of the gene regulatory networks, with a trend towards gain of regulatory interactions between genes, and (v) 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. Extending the range of parameters, we provide empirically testable predictions on the differences of genetic architectures between wild and domesticated and forms. The characterization of such systematic evolutionary changes in the genetic architecture of traits contributes to define a molecular domestication syndrome.