Genotype by Environment Interactions in Gene Regulation Underlie the Response to Soil Drying in the Model Grass Brachypodium distachyon.

Gene expression is a quantitative trait under the control of genetic and environmental factors and their interaction, so-called genotype and environment (G × E). Understanding the mechanisms driving G × E is fundamental for ensuring stable crop performance across environments and for predicting the response of natural populations to climate change. Gene expression is regulated through complex molecular networks, yet the interactions between genotype and environment in gene regulation are rarely considered, particularly at the genome scale. Current frameworks and experimental designs often lack power to explicitly test network rewiring or to systematically compare regulatory networks. Here, we leverage a highly replicated RNA-sequencing dataset to model genome-scale gene expression variation between two natural accessions of the model grass Brachypodium distachyon and their response to soil drying. We first identified genotypic, environmental, and G × E effects on physiological, metabolic, and gene expression traits. We identify patterns of conservation-or variation-in gene coexpression networks and link these coexpression features to physiological traits. We further develop predictions of gene-gene interactions using causal inference and screen for interactions specific to-or with higher affinity in-a single genotype, treatment, or their interaction, G × E. Our analyses identify variation in candidate gene regulatory networks that may shape the evolution of environmental response in B. distachyon. We highlight the environmentally dependent regulatory control of several metabolic traits shown previously to play a role in drought acclimation. The framework presented here provides a scalable approach for more complex comparisons, particularly with the growing availability of large datasets from technologies such as single-cell transcriptomics.