Multi-omics atlas of combinatorial abiotic stress responses in wheat

Field-grown crops rarely experience growth conditions in which yield can be maximized. Environmental stresses occur in combination, with advancements in crop tolerance further complicated by its polygenic nature. Strategic targeting of causal genes is required to meet future crop production needs. Here, we employed a systems biology approach in wheat (Triticum aestivum L.) to investigate physio-metabolic adjustments and transcriptome reprogramming involved in acclimations to heat, drought, salinity and all combinations therein. A significant shift in magnitude and complexity of plant response was evident across stress scenarios based on the agronomic losses, increased proline concentrations and 8.7-fold increase in unique differentially expressed transcripts (DETs) observed under the triple stress condition. Transcriptome data from all stress treatments were assembled into an online, open access eFP browser for visualizing gene expression during abiotic stress. Weighted gene co-expression network analysis revealed 152 hub genes of which 32% contained the ethylene-responsive element binding factor-associated amphiphilic repression (EAR) transcriptional repression motif. Cross-referencing against the 31 DETs common to all stress treatments isolated TaWRKY33 as a leading candidate for greater plant tolerance to combinatorial stresses. Integration of our findings with available literature on gene functional characterization allowed us to further suggest flexible gene combinations for future adaptive gene stacking in wheat. Our approach demonstrates the strength of robust multi-omics-based data resources for gene discovery in complex environmental conditions. Accessibility of such datasets will promote cross-validation of candidate genes across studies and aid in accelerating causal gene validation for crop resiliency.