Variable paralog expression underlies phenotype variation

Human faces are variable; we look different from one another. Craniofacial disorders further increase this variability. Here we used the zebrafish mef2ca mutant, which produces variable phenotypes, to understand craniofacial variation. Comparing different mef2ca alleles demonstrated that severity, measured by penetrance and expressivity, correlates with variation. Years of selective breeding for low and high penetrance produced strains that are either resilient, or sensitive, to the mef2ca mutation. Comparing these strains further demonstrates that severity correlates with variation. Gene expression studies indicated that selective breeding upregulated and downregulated mef2ca paralog expression in the low- and high-penetrance strains, respectively. We hypothesized that heritable paralog expression variation underlies mutant phenotype variation. In support, mutagenizing all mef2ca paralogs in the low-penetrance strain demonstrated modular buffering by paralogs. Specifically, some paralogs buffer severity while others buffer variability. We present a novel, mechanistic model for phenotypic variation where cryptic vestigial paralog expression modularly buffers development and contributes to evolution. These studies are a major step forward in understanding of the mechanisms of facial variation, including how some genetically resilient individuals can overcome a deleterious mutation.