Systems genomics of salinity stress response in rice

Populations can adapt to stressful environments through changes in gene expression. However, the fitness effect of gene expression in mediating stress response and adaptation remains largely unexplored. Here, we use an integrative field dataset obtained from 780 plants of <i>Oryza sativa</i> ssp. <i>indica</i> (rice) grown in a field experiment under normal or moderate salt stress conditions to examine selection and evolution of gene expression variation under salinity stress conditions. We find that salinity stress induces increased selective pressure on gene expression. Further, we show that <i>trans</i>-eQTLs rather than <i>cis</i>-eQTLs are primarily associated with rice's gene expression under salinity stress, potentially via a few master-regulators. Importantly, and contrary to the expectations, we find that <i>cis-trans</i> reinforcement is more common than <i>cis-trans</i> compensation which may be reflective of rice diversification subsequent to domestication. We further identify genetic fixation as the likely mechanism underlying this compensation/reinforcement. Additionally, we show that <i>cis</i>- and <i>trans</i>-eQTLs are under balancing and purifying selection, respectively, giving us insights into the evolutionary dynamics of gene expression variation. By examining genomic, transcriptomic, and phenotypic variation across a rice population, we gain insights into the molecular and genetic landscape underlying adaptive salinity stress responses, which is relevant for other crops and other stresses.