Transcriptomic and Metabolomic Analysis of Seedling-Stage Soybean Responses to PEG-Simulated Drought Stress

Soybean is an important crop grown worldwide, and drought stress seriously affects the yield and quality of soybean. Therefore, it is necessary to elucidate the molecular mechanisms underlying soybean resistance to drought stress. In this study, RNA-seq technology and ultra-performance liquid chromatography–tandem mass spectrometry were used to analyze the transcriptome and metabolome changes in soybean leaves at the seedling stage under drought stress. The results showed that there were 4790 and 3483 DEGs (differentially expressed genes) and 156 and 124 DAMs (differentially expressed metabolites), respectively, in the HN65CK vs. HN65S0 and HN44CK vs. HN44S0 comparison groups. Comprehensive analysis of transcriptomic and metabolomic data reveals metabolic regulation of seedling soybean in response to drought stress. Some candidate genes such as LOC100802571, LOC100814585, LOC100777350 and LOC100787920, LOC100800547, and LOC100785313 showed different expression trends between the two cultivars, which may cause differences in drought resistance. Secondly, a large number of flavonoids were identified, and the expression of Monohydroxy-trimethoxyflavone-O-(6″-malonyl)glucoside was upregulated between the two varieties. Finally, several key candidate genes and metabolites involved in isoflavone biosynthesis and the TCA cycle were identified, suggesting that these metabolic pathways play important roles in soybean response to drought. Our study deepens the understanding of soybean drought resistance mechanisms and provides references for soybean drought resistance breeding.