A system genetics analysis uncovers the regulatory variants controlling drought response in wheat

Plants activate a variable response to drought stress by modulating transcription of key genes. However, our knowledge of genetic variations governing gene expression in response to drought stress remains limited in natural germplasm. Here, we performed a comprehensive analysis of the transcriptional variability of 200 wheat accessions in response to drought stress by using a systems genetics approach integrating pan-transcriptome, co-expression networks, transcriptome-wide association study (TWAS), and expression quantitative trait loci (eQTLs) mapping. We identified 1621 genes and eight co-expression modules significantly correlated with wheat drought tolerance. We also defined 620 664 and 654 798 independent eQTLs associated with the expression of 17 429 and 18 080 eGenes under normal and drought stress conditions. Focusing on dynamic regulatory variants, we further identified 572 eQTL hotspots and constructed transcription factors governed drought-responsive network by the XGBoost model. Subsequently, by combining with genome-wide association study (GWAS), we uncovered a 369-bp insertion variant in the TaKCS3 promoter containing multiple cis-regulatory elements recognized by eQTL hotspot-associated transcription factors that enhance its transcription. Further functional analysis indicated that elevating TaKCS3 expression affects cuticular wax composition to reduce water loss during drought stress, and thereby increase drought tolerance. This study sheds light on the genome-wide genetic variants that influence dynamic transcriptional changes during drought stress and provides a valuable resource for the mining of drought-tolerant genes in the future.