Evolution of dominance in gene expression pattern associated with phenotypic robustness

BACKGROUND: Mendelian inheritance is a fundamental law of genetics. When we consider two genomes in a diploid cell, a heterozygote’s phenotype is dominated by a particular homozygote according to the law of dominance. Classical Mendelian dominance is concerned with which proteins are dominant, and is usually based on simple genotype–phenotype relationship in which one gene regulates one phenotype. However, in reality, some interactions between genes can exist, resulting in deviations from Mendelian dominance. Whether and how Mendelian dominance is generalized to the phenotypes of gene expression determined by gene regulatory networks (GRNs) remains elusive. RESULTS: Here, by using the numerical evolution of diploid GRNs, we discuss whether the dominance of phenotype evolves beyond the classical Mendelian case of one-to-one genotype–phenotype relationship. We examine whether complex genotype–phenotype relationship can achieve Mendelian dominance at the expression level by a pair of haplotypes through the evolution of the GRN with interacting genes. This dominance is defined via a pair of haplotypes that differ from each other but have a common phenotype given by the expression of target genes. We numerically evolve the GRN model for a diploid case, in which two GRN matrices are added to give gene expression dynamics and simulate evolution with meiosis and recombination. Our results reveal that group Mendelian dominance evolves even under complex genotype–phenotype relationship. Calculating the degree of dominance shows that it increases through the evolution, correlating closely with the decrease in phenotypic fluctuations and the increase in robustness to initial noise. We also demonstrate that the dominance of gene expression patterns evolves concurrently. This evolution of group Mendelian dominance and pattern dominance is associated with phenotypic robustness against meiosis-induced genome mixing, whereas sexual recombination arising from the mixing of genomes from the parents further enhances dominance and robustness. Due to this dominance, the robustness to genetic differences increases, while optimal fitness is sustained to a significant difference between the two genomes. CONCLUSION: Group Mendelian dominance and gene-expression pattern dominance are achieved associated with the increase in phenotypic robustness to noise. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12862-021-01841-6.