Center for Genomics and Systems Biology, New York University, New York, United States
Simon Niels Groen
Center for Genomics and Systems Biology, New York University, New York, United States; Department of Nematology and Department of Botany & Plant Sciences, University of California, Riverside, Riverside, United States; Center for Plant Cell Biology, Institute for Integrative Genome Biology, University of California, Riverside, Riverside, United States
Maricris L Zaidem
Center for Genomics and Systems Biology, New York University, New York, United States; Department of Biology, University of Oxford, Oxford, United Kingdom
Andres Godwin C Sajise
International Rice Research Institute, Los Baños, Philippines
Irina Calic
Department of Biological Sciences, Fordham University, Bronx, United States; Inari Agriculture Nv, Gent, Belgium
Mignon Natividad
International Rice Research Institute, Los Baños, Philippines
Kenneth McNally
International Rice Research Institute, Los Baños, Philippines
Georgina V Vergara
International Rice Research Institute, Los Baños, Philippines; Institute of Crop Science, University of the Philippines, Los Baños, Philippines
Rahul Satija
Center for Genomics and Systems Biology, New York University, New York, United States; New York Genome Center, New York, United States
Steven J Franks
Department of Biological Sciences, Fordham University, Bronx, United States
Rakesh K Singh
International Rice Research Institute, Los Baños, Philippines; International Center for Biosaline Agriculture, Dubai, United Arab Emirates
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 Oryza sativa ssp. indica (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 trans-eQTLs rather than cis-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 cis-trans reinforcement is more common than cis-trans 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 cis- and trans-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.