Hybrid incompatibility can result from gene misregulation produced by divergence in trans-acting regulatory factors and their cis-regulatory targets. However, change in trans-acting factors may be constrained by pleiotropy, which would in turn limit the evolution of incompatibility. We employed a mechanistically explicit bioenergetic model of gene expression wherein parameter combinations (number of transcription factor molecules, energetic properties of binding to the regulatory site, and genomic background size) determine the shape of the genotype-phenotype (G-P) map, and interacting allelic variants of mutable cis and trans sites determine the phenotype along that map. Misregulation occurs when the phenotype differs from its optimal value. We simulated a pleiotropic regulatory pathway involving a positively selected and a conserved trait regulated by a shared transcription factor (TF), with two populations evolving in parallel. Pleiotropic constraints shifted evolution in the positively selected trait to its cis-regulatory locus. We nevertheless found that the TF genotypes often evolved, accompanied by compensatory evolution in the conserved trait, and both traits contributed to hybrid misregulation. Compensatory evolution resulted in "developmental system drift," whereby the regulatory basis of the conserved phenotype changed although the phenotype itself did not. Pleiotropic constraints became stronger and in some cases prohibitive when the bioenergetic properties of the molecular interaction produced a G-P map that was too steep. Likewise, compensatory evolution slowed and hybrid misregulation was not evident when the G-P map was too shallow. A broad pleiotropic "sweet spot" nevertheless existed where evolutionary constraints were moderate to weak, permitting substantial hybrid misregulation in both traits. None of these pleiotropic constraints manifested when the TF contained nonrecombining domains independently regulating the respective traits.A DAPTIVE changes in phenotype often occur through changes in gene regulation (Wray 2007;Carroll 2008). Regulatory networks map an organism's genotype to its phenotype through developmental and physiological processes (Wilkins 2002). Such networks consist of interacting loci that can respond to selection by changing the expression levels of individual genes in space and time. In addition to gene interaction (epistasis) (Phillips 2008), gene networks are also characterized by pleiotropy, where single genetic loci have manifold effects (Gibson 1996;Paaby and Rockman 2013).Gene interactions are pervasive in the evolution of hybrid incompatibility, an important form of reproductive isolation between species (Coyne and Orr 2004). The leading model for the evolution of hybrid incompatibility, the BatesonDobzhansky-Muller (BDM) model, requires interactions between at least two genetic loci (Bateson 1909;Dobzhansky 1937;Muller 1942). Interpopulation divergence in regulatory interactions may be expected to result in hybrid incompatibility due to misregulation of th...