Transcriptomic Analysis on Developing Seed Uncovers Candidate Genes Associated with Seed Storage Protein in Soybean

Soybean [Glycine max (L.) Merr.] is a globally significant crop that provides essential meal protein and vegetable oil for human consumption. The protein content in soybean seeds is a critical factor that affects nutrition regarding human dietary needs as well as livestock feed. Therefore, identifying the key genes that affect the soybean seed protein content is one of the major goals in soybean research. To identify candidate genes and related pathways involved in soybean seed storage protein during seed development, an RNA-seq analysis was conducted in two soybean varieties that differ in protein content. A series of pathways related to seed protein metabolism, including “Photosynthesis”, “TCA cycle”, and “Starch and sucrose metabolism” pathways, were identified through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Seven candidate genes exhibiting two different gene regulation patterns were identified, six of which are directly related to the seed storage protein pathway, and one of which is related to the carbon binding pathway. An integrated analysis of transcriptomic and candidate gene expression trend suggested that 40 days after flowering (DAF) might be a crucial period for seed protein accumulation in soybean. Through a Weighted Gene Co-expression Network Analysis (WGCNA), two modules and two novel hub genes were found, which may be highly correlated with seed protein development. These findings might be valuable for a complete understanding of the genetic basis of seed protein content and lay a theoretical foundation for future gene functional identification and breeding efforts in soybean.