The dominance of gene expression controlled by trans-eQTL hotspots contributes to phenotypic heterosis in maize

Abstract

Heterosis, or hybrid vigor, is a key phenomenon in genetics research and agricultural production, and has been primarily attributed to non-additive genetic effects such as dominance -- a prevailing consensus shaped by decades of empirical research and theoretical debate. Although dominance may arguably arise from distal modifiers, their selective advantage is debated due to presumably small individual effects. To address this long-standing question, particularly how genetic dominance manifests at the transcriptomic level and contributes to phenotypic heterosis, we integrated transcriptomic and phenotypic data from a large population of maize hybrids and their inbred parents. We found that [~] 30% of the expressed seedling genes in a significant proportion of hybrids exhibited expression patterns deviating from the average of the two parents, indicative of non-additivity. Further analysis suggests that while hybrid gene expression per se is primarily regulated by cis-eQTLs, expression dominance (or non-additivity) is disproportionately controlled by trans-eQTLs. These trans-eQTLs cluster into hotspots that regulate the non-additivity of hundreds of target genes, mostly within co-expression networks, and are notably enriched for transcription factors (TFs). Focusing on one such hotspot, we functionally validated a classical maize gene ZmR1, a basic helix-loop-helix (bHLH) TF associated with multiple seedling trait heterosis, as a candidate regulator of expression dominance across hundreds of genes. Overexpression of ZmR1 enhances expression dominance of downstream genes and increases phenotypic heterosis in both seedling and adult traits. Further experiments confirmed its direct regulatory role in modulating genes involved in anthocyanin biosynthesis and lignin metabolism, driving transcriptome-level dominance. These results provide empirical support for the modifier hypothesis under an omnigenic model, suggesting that heterosis arises not from the modification of a single genes inheritance but through the coordinated regulation of hundreds of phenotype-associated genes, thereby helping to reconcile the long-standing debate over the genetic basis of dominance in heterosis.