Hsp90 is one of the most conserved and abundant molecular chaperones and is an essential component of the protective stress response; however, its roles in abiotic stress responses in soybean (Glycine max) remain obscure. Here, 12 GmHsp90 genes from soybean were identified and found to be expressed and to function differentially under abiotic stresses. The 12 GmHsp90 genes were isolated and named GmHsp90A1–GmHsp90A6, GmHsp90B1, GmHsp90B2, GmHsp90C1.1, GmHsp90C1.2, GmHsp90C2.1 and GmHsp90C2.2 based on their characteristics and high homology to other Hsp90s according to a new nomenclature system. Quantitative real-time PCR expression data revealed that all the genes exhibited higher transcript levels in leaves and could be strongly induced under heat, osmotic and salt stress but not cold stress. Overexpression of five typical genes (GmHsp90A2, GmHsp90A4, GmHsp90B1, GmHsp90C1.1 and GmHsp90C2.1) in Arabidopsis thaliana provided useful evidences that GmHsp90 genes can decrease damage of abiotic stresses. In addition, an abnormal accumulation of proline was detected in some transgenic Arabidopsis plants suggested overexpressing GmHsp90s may affect the synthesis and response system of proline. Our work represents a systematic determination of soybean genes encoding Hsp90s, and provides useful evidence that GmHsp90 genes function differently in response to abiotic stresses and may affect the synthesis and response system of proline.

Background: The WRKY proteins are a superfamily of transcription factors and members play essential roles in the modulation of diverse physiological processes, such as growth, development, senescence and response to biotic and abiotic stresses. However, the biological roles of the majority of the WRKY family members remains poorly understood in soybean relative to the research progress in model plants. Results: In this study, we identified and characterized GmWRKY40, which is a group IIc WRKY gene. Transient expression analysis revealed that the GmWRKY40 protein is located in the nucleus of plant cells. Expression of GmWRKY40 was strongly induced in soybean following infection with Phytophthora sojae, or treatment with methyl jasmonate, ethylene, salicylic acid, and abscisic acid. Furthermore, soybean hairy roots silencing GmWRKY40 enhanced susceptibility to P. sojae infection compared with empty vector transgenic roots. Moreover, suppression of GmWRKY40 decreased the accumulation of reactive oxygen species (ROS) and modified the expression of several oxidation-related genes. Yeast two-hybrid experiment combined with RNA-seq analysis showed that GmWRKY40 interacted with 8 JAZ proteins with or without the WRKY domain or zinc-finger domain of GmWRKY40, suggesting there were different interaction patterns among these interacted proteins. Conclusions: Collectively, these results suggests that GmWRKY40 functions as a positive regulator in soybean plants response to P. sojae through modulating hydrogen peroxide accumulation and JA signaling pathway.

Plant height (PH) is an important agronomic trait and is closely related to yield in soybean [Glycine max (L.) Merr.]. Previous studies have identified many QTLs for PH. Due to the complex genetic background of PH in soybean, there are few reports on its fine mapping. In this study, we used a mapping population derived from a cross between a chromosome segment substitution line CSSL3228 (donor N24852 (G. Soja), a receptor NN1138–2 (G. max)) and NN1138–2 to fine map a wild soybean allele of greater PH by QTL-seq and linkage mapping. We identified a QTL for PH in a 1.73 Mb region on soybean chromosome 13 through QTL-seq, which was confirmed by SSR marker-based classical QTL mapping in the mapping population. The linkage analysis showed that the QTLs of PH were located between the SSR markers BARCSOYSSR_13_1417 and BARCSOYSSR_13_1421 on chromosome 13, and the physical distance was 69.3 kb. RT-PCR and sequence analysis of possible candidate genes showed that Glyma.13 g249400 revealed significantly higher expression in higher PH genotypes, and the gene existed 6 differences in the amino acids encoding between the two parents. Data presented here provide support for Glyma.13 g249400 as a possible candidate genes for higher PH in wild soybean line N24852.

Seedling drought stress is one of the most important constraints affecting soybean yield and quality. To unravel the molecular mechanisms under soybean drought tolerance, we conducted comprehensive comparative transcriptome analyses of drought-tolerant genotype Jindou 21 (JD) and drought-sensitive genotype Tianlong No.1 (N1) seedlings that had been exposed to drought treatment. A total of 6038 and 4112 differentially expressed genes (DEGs) were identified in drought-tolerant JD and drought-sensitive N1, respectively. Subsequent KEGG pathway analyses showed that numerous DEGs in JD are predominately involved in signal transduction pathways, including plant hormone signaling pathway, calcium signaling pathway, and MAPK signaling pathway. Interestingly, JA and BR plant hormone signal transduction pathways were found specifically participating in drought-tolerant JD. Meanwhile, the differentially expressed CPKs, CIPKs, MAPKs, and MAP3Ks of calcium and MAPK signaling pathway were only identified in JD. The number of DEGs involved in transcription factors (TFs) is larger in JD than that of in N1. Moreover, some differently expressed transcriptional factor genes were only identified in drought-tolerant JD, including FAR1, RAV, LSD1, EIL, and HB-PHD. In addition, this study suggested that JD could respond to drought stress by regulating the cell wall remodeling and stress-related protein genes such as EXPs, CALSs, CBPs, BBXs, and RD22s. JD is more drought tolerant than N1 owing to more DEGs being involved in multiple signal transduction pathways (JA, BR, calcium, MAPK signaling pathway), stress-related TFs, and proteins. The above valuable genes and pathways will deepen the understanding of the molecular mechanisms under drought stress in soybean.

Soybean is one of the most important economic and oil crops across the world. Phytophthora root rot (PRR), caused by Phytophthora sojae (P. sojae), is a major disease in most soybean-growing regions worldwide. Here, we investigated metabolic changes in hypocotyls of two soybean lines, Nannong 10-1 (resistant line, R) and 06-070583 (susceptible line, S), at two time points (12 and 36 hpi) after P. sojae infection and metabolic differences between the R line and the S line. In total, 90 differentially accumulated metabolites (DAMs) were identified after P. sojae infection; the levels of 50 metabolites differed between the R line and the S line. There are 28 DAMs that not only differentially accumulated between the R line and the S line but also differentially accumulated after P. sojae infection. Based on the changes of these DAMs in response to P. sojae infection in different lines and at different timepoints, and the differences in the contents of these DAMs between the R line and the S line, we speculated that DAMs, including sugars (monosaccharides and oligosaccharides), organic acids (oxalic acid, cumic acid), amino acid derivatives, and other secondary metabolites (mannitol, octanal, hypoxanthine, and daidzein etc.) may participate in the metabolic-level defense response of soybean to P. sojae. In this study, an integrated pathway-level analysis of transcriptomics (obtained by RNA-Seq) and metabolomics data illustrated the poor connections and interdependencies between the metabolic and transcriptional responses of soybean to P. sojae infection. This work emphasizes the value of metabolomic studies of plant–pathogen interactions and paves the way for future research of critical metabolic determinants of the soybean-P. sojae interaction.