  Biological systems must be robust for stable function against perturbations,
but robustness alone is not sufficient. The ability to switch between
appropriate states (phenotypes) in response to different conditions is
essential for biological functions. How are robustness and plasticity
simultaneously acquired through evolution? We examine the evolution of
genotypes that realize plastic switching between two phenotypes upon external
inputs, as well as stationary expressions of phenotypes. We introduce a
statistical physics model consisting of spins, with active and regulatory
sites, which are distinct from each other. We represent the phenotype and
genotype as spin configurations and the spin-spin interactions, respectively.
The fitness for selection is given so that it takes a higher value as more of
the active sites take two requested spin configurations depending on the
regulation. We numerically evolve the interaction matrix by changing them with
mutations and selection of those with higher fitness. Our numerical simulations
show that characteristic genotypes evolve slightly above the transition
temperature between replica symmetric and replica symmetry breaking phase.
These genotypes shape two spin configurations separately depending on the
regulation, where the two phenotypes are dominantly represented by the
genotypes' first and second eigenmodes, and smooth switching of two phenotypes
are achieved by following a one-dimensional path connecting the two phenotypes.
Upon changes in regulations, spin configurations are attracted to this path,
which allows for robust and plastic switching between the two phenotypes. The
statistical-physics analysis show that the free energy landscape has a valley
along the switching path. Our finding indicates that the compatibility of the
robustness and plasticity is acquired by the evolution of the
low-dimensionality in the phenotype space.
