Unifying framework explaining how parental regulatory divergence can drive gene expression in hybrids and allopolyploids

Hybridization and polyploidy are powerful evolutionary forces, inducing a range of phenotypic outcomes, including non-additive expression, subgenome dominance, deviations in genomic dosage, and transcriptome downsizing. The reasons for these patterns and whether they are universal adaptive responses to genome merging and doubling remain debated. To address this, we develop a thermodynamic model of gene expression based on transcription factor (TF)-promoter binding. Applied to hybridization between species with divergent gene expression levels, cell volumes, or euchromatic ratios, this model distinguishes the effects of hybridization from those of polyploidy. Our results align with empirical observations, suggesting that gene regulation patterns in hybrids and polyploids often stem from the constrained interplay between inherited diverged regulatory networks rather than from subsequent adaptive evolution. In addition, occurrence of certain phenotypic traits depend on specific assumptions about promoter-TF coevolution and their distribution within the hybrid's nucleoplasm, offering new research avenues to understand the underlying mechanisms. In summary, our model explains how the legacy of divergent species directly influences the phenotypic traits of hybrids and allopolyploids.