Transcriptomic and Metabolomic Analyses Reveal the Key Genes Related to Shade Tolerance in Soybean

Soybean (Glycine max) is an important crop, rich in proteins, vegetable oils and several other phytochemicals, which is often affected by light during growth. However, the specific regulatory mechanisms of leaf development under shade conditions have yet to be understood. In this study, the transcriptome and metabolome sequencing of leaves from the shade-tolerant soybean ‘Nanxiadou 25′ under natural light (ND1) and 50% shade rate (SHND1) were carried out, respectively. A total of 265 differentially expressed genes (DEGs) were identified, including 144 down-regulated and 121 up-regulated genes. Meanwhile, KEGG enrichment analysis of DEGs was performed and 22 DEGs were significantly enriched in the top five pathways, including histidine metabolism, riboflavin metabolism, vitamin B6 metabolism, glycerolipid metabolism and cutin, suberine and wax biosynthesis. Among all the enrichment pathways, the most DEGs were enriched in plant hormone signaling pathways with 19 DEGs being enriched. Transcription factors were screened out and 34 differentially expressed TFs (DETFs) were identified. Weighted gene co-expression network analysis (WGCNA) was performed and identified 10 core hub genes. Combined analysis of transcriptome and metabolome screened out 36 DEGs, and 12 potential candidate genes were screened out and validated by quantitative real-time polymerase chain reaction (qRT-PCR) assay, which may be related to the mechanism of shade tolerance in soybean, such as ATP phosphoribosyl transferase (ATP-PRT2), phosphocholine phosphatase (PEPC), AUXIN-RESPONSIVE PROTEIN (IAA17), PURPLE ACID PHOSPHATASE (PAP), etc. Our results provide new knowledge for the identification and function of candidate genes regulating soybean shade tolerance and provide valuable resources for the genetic dissection of soybean shade tolerance molecular breeding.