GWAS and WGCNA Analysis Uncover Candidate Genes Associated with Oil Content in Soybean.

Soybean vegetable oil is an important source of the human diet. However, the analysis of the genetic mechanism leading to changes in soybean oil content is still incomplete. In this study, a total of 227 soybean materials were applied and analyzed by a genome-wide association study (GWAS). There are 44 quantitative trait nucleotides (QTNs) that were identified as associated with oil content. A total of six, four, and 34 significant QTN loci were identified in Xiangyang, Hulan, and Acheng, respectively. Of those, 26 QTNs overlapped with or were near the known oil content quantitative trait locus (QTL), and 18 new QTNs related to oil content were identified. A total of 594 genes were located near the peak single nucleotide polymorphism (SNP) from three tested environments. These candidate genes exhibited significant enrichment in tropane, piperidine, and pyridine alkaloid biosynthesiss (ko00960), ABC transporters (ko02010), photosynthesis-antenna proteins (ko00196), and betalain biosynthesis (ko00965). Combined with the GWAS and weighted gene co-expression network analysis (WGCNA), four candidate genes (Glyma.18G300100, Glyma.11G221100, Glyma.13G343300, and Glyma.02G166100) that may regulate oil content were identified. In addition, Glyma.18G300100 was divided into two main haplotypes in the studied accessions. The oil content of haplotype 1 is significantly lower than that of haplotype 2. Our research findings provide a theoretical basis for improving the regulatory mechanism of soybean oil content.

Integrating Genome-Wide Association Study, Transcriptome and Metabolome Reveal Novel QTL and Candidate Genes That Control Protein Content in Soybean.

Protein content (PC) is crucial to the nutritional quality of soybean [Glycine max (L.) Merrill]. In this study, a total of 266 accessions were used to perform a genome-wide association study (GWAS) in three tested environments. A total of 23,131 high-quality SNP markers (MAF ≥ 0.02, missing data ≤ 10%) were identified. A total of 40 association signals were significantly associated with PC. Among them, five novel quantitative trait nucleotides (QTNs) were discovered, and another 32 QTNs were found to be overlapping with the genomic regions of known quantitative trait loci (QTL) related to soybean PC. Combined with GWAS, metabolome and transcriptome sequencing, 59 differentially expressed genes (DEGs) that might control the change in protein content were identified. Meantime, four commonly upregulated differentially abundant metabolites (DAMs) and 29 commonly downregulated DAMs were found. Remarkably, the soybean gene Glyma.08G136900, which is homologous with Arabidopsis hydroxyproline-rich glycoproteins (HRGPs), may play an important role in improving the PC. Additionally, Glyma.08G136900 was divided into two main haplotype in the tested accessions. The PC of haplotype 1 was significantly lower than that of haplotype 2. The results of this study provided insights into the genetic mechanisms regulating protein content in soybean.

Multi-omics analysis reveals novel loci and a candidate regulatory gene of unsaturated fatty acids in soybean (Glycine max (L.) Merr).

BACKGROUND: Soybean is a major oil crop; the nutritional components of soybean oil are mainly controlled by unsaturated fatty acids (FA). Unsaturated FAs mainly include oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). The genetic architecture of unsaturated FAs in soybean seeds has not been fully elucidated, although many independent studies have been conducted. A 3 V multi-locus random single nucleotide polymorphism (SNP)-effect mixed linear model (3VmrMLM) was established to identify quantitative trait loci (QTLs) and QTL-by-environment interactions (QEIs) for complex traits. RESULTS: In this study, 194 soybean accessions with 36,981 SNPs were calculated using the 3VmrMLM model. As a result, 94 quantitative trait nucleotides (QTNs) and 19 QEIs were detected using single-environment (QTN) and multi-environment (QEI) methods. Three significant QEIs, namely rs4633292, rs39216169, and rs14264702, overlapped with a significant single-environment QTN. CONCLUSIONS: For QTNs and QEIs, further haplotype analysis of candidate genes revealed that the Glyma.03G040400 and Glyma.17G236700 genes were beneficial haplotypes that may be associated with unsaturated FAs. This result provides ideas for the identification of soybean lipid-related genes and provides insights for breeding high oil soybean varieties in the future.

Genome-wide identification and analysis of glyceraldehyde-3-phosphate dehydrogenase family reveals the role of GmGAPDH14 to improve salt tolerance in soybean (Glycine max L.).

Introduction: Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential key enzyme in the glycolytic pathway and plays an important role in stress responses. Although GAPDH family genes have been found in different plant species, the determination of their gene family analysis and their functional roles in soybean are still unknown. Methods: In this study, gene sequence and expression data were obtained using online tools, and systematic evolution, expression profile analysis, and qRT-PCR analysis were conducted. Results and Discussion: Here a total of 16 GmGAPDH genes were identified on nine chromosomes, which were classified into three clusters. Additionally, all GmGAPDH genes harbor two highly conserved domains, including Gp_dh_N (PF00044) and Gp_dh_C (PF02800). The qRTPCR analysis also showed that most GmGAPDH genes significantly responded to multiple abiotic stresses, including NaHCO3, polyethylene glycol, cold, and salt. Among them, GmGAPDH14 was extraordinarily induced by salt stress. The GmGAPDH14 gene was cloned and overexpressed through soybean hair roots. The overexpressed transgenic soybean plants of the GmGAPDH14 gene have also shown better growth than that of control plants. Moreover, the overexpressed transgenic plants of GmGAPDH14 gene had higher activities of superoxide dismutase but lower malonaldehyde (MDA) content than those of control plants under salt stress. Meanwhile, a total of four haplotypes were found for the GmGAPDH14 gene, and haplotypes 2, 3, and 4 were beneficial for the tolerance of soybean to salt stress. These results suggest that the GmGAPDH14 gene might be involved in the process of soybean tolerance to salt stress. The results of this study will be valuable in understanding the role of GAPDH genes in the abiotic stress response of soybean.

Combined analysis of the metabolome and transcriptome provides insight into seed oil accumulation in soybean.

BACKGROUND: Soybean (Glycine max (L.) Merr) is an important source of human food, animal feed, and bio-energy. Although the genetic network of lipid metabolism is clear in Arabidopsis, the understanding of lipid metabolism in soybean is limited. RESULTS: In this study, 30 soybean varieties were subjected to transcriptome and metabolome analysis. In total, 98 lipid-related metabolites were identified, including glycerophospholipid, alpha-linolenic acid, linoleic acid, glycolysis, pyruvate, and the sphingolipid pathway. Of these, glycerophospholipid pathway metabolites accounted for the majority of total lipids. Combining the transcriptomic and metabolomic analyses, we found that 33 lipid-related metabolites and 83 lipid-related genes, 14 lipid-related metabolites and 17 lipid-related genes, and 12 lipid-related metabolites and 25 lipid-related genes were significantly correlated in FHO (five high-oil varieties) vs. FLO (five low-oil varieties), THO (10 high-oil varieties) vs. TLO (10 low-oil varieties), and HO (15 high-oil varieties) vs. LO (15 low-oil varieties), respectively. CONCLUSIONS: The GmGAPDH and GmGPAT genes were significantly correlated with lipid metabolism genes, and the result revealed the regulatory relationship between glycolysis and oil synthesis. These results improve our understanding of the regulatory mechanism of soybean seed oil improvement.

Genome-Wide Identification of Candidate Genes Underlying Soluble Sugar Content in Vegetable Soybean (Glycine max L.) via Association and Expression Analysis.

Soluble sugar is a major indicator of the intrinsic quality of vegetable soybean [Glycine max (L.) Merr. ]. The improvement of soluble sugar content in soybean is very important due to its healthcare functions for humans. The genetic mechanism of soluble sugar in soybean is unclear. In this study, 278 diverse soybean accessions were utilized to identify the quantitative trait nucleotides (QTNs) for total soluble sugar content in soybean seeds based on a genome-wide association study (GWAS). A total of 25,921 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected for GWAS. Totally, thirteen QTNs associated with total soluble sugar content were identified, which were distributed on ten chromosomes. One hundred and fifteen genes near the 200-kb flanking region of these identified QTNs were considered candidate genes associated with total soluble sugar content in soybean seed. Gene-based association analysis and haplotype analysis were utilized to further identify the effect of candidate genes on total soluble sugar content. Totally, 84 SNPs from seventeen genes across four chromosomes were significantly associated with the total soluble sugar content. Among them, three SNPs from Glyma.02G292900 were identified at two locations, and other eighty-one SNPs from sixteen genes were detected at three locations. Furthermore, expression level analysis of candidate genes revealed that Glyma.02G293200 and Glyma.02G294900 were significantly positively associated with soluble sugar content and Glyma.02G294000 was significantly negatively associated with soluble sugar content. Six genes (i.e., Glyma.02G292600, Glyma.02G292700, Glyma.02G294000, Glyma.02G294300, Glyma.02G294400, and Glyma.15G264200) identified by GWAS were also detected by the analysis of differential expression genes based on soybean germplasms with higher and lower soluble sugar content. Among them, Glyma.02G294000 is the only gene that was identified by gene-based association analysis with total soluble sugar content and was considered an important candidate gene for soluble sugar content. These candidate genes and beneficial alleles would be useful for improving the soluble sugar content of soybean.

Membrane Lipids' Metabolism and Transcriptional Regulation in Maize Roots Under Cold Stress.

Low temperature is one of the major abiotic stresses that restrict the growth and development of maize seedlings. Membrane lipid metabolism and remodeling are key strategies for plants to cope with temperature stresses. In this study, an integrated lipidomic and transcriptomic analysis was performed to explore the metabolic changes of membrane lipids in the roots of maize seedlings under cold stress (5°C). The results revealed that major extraplastidic phospholipids [phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidic acid (PA), and phosphatidylinositol (PI)] were dominant membrane lipids in maize root tissues, accounting for more than 70% of the total lipids. In the transcriptome data of maize roots under cold stress, a total of 189 lipid-related differentially expressed genes (DEGs) were annotated and classified into various lipid metabolism pathways, and most of the DEGs were enriched in the "Eukaryotic phospholipid synthesis" (12%), "Fatty acid elongation" (12%), and "Phospholipid signaling" (13%) pathways. Under low temperature stress, the molar percentage of the most abundant phospholipid PC decreased around 10%. The significantly up-regulated expression of genes encoding phospholipase [phospholipase D (PLD)] and phosphatase PAP/LPP genes implied that PC turnover was triggered by cold stress mainly via the PLD pathway. Consequently, as the central product of PC turnover, the level of PA increased drastically (63.2%) compared with the control. The gene-metabolite network and co-expression network were constructed with the prominent lipid-related DEGs to illustrate the modular regulation of metabolic changes of membrane lipids. This study will help to explicate membrane lipid remodeling and the molecular regulation mechanism in field crops encountering low temperature stress.

Identification of Fatty Acid Desaturases in Maize and Their Differential Responses to Low and High Temperature.

Plant fatty acid desaturases (FADs) catalyze the desaturation of fatty acids in various forms and play important roles in regulating fatty acid composition and maintaining membrane fluidity under temperature stress. A total of 30 FADs were identified from a maize genome, including 13 soluble and 17 membrane-bound FADs, which were further classified into two and five sub-groups, respectively, via phylogenetic analysis. Although there is no evolutionary relationship between the soluble and the membrane-bound FADs, they all harbor a highly conserved FA_desaturase domain, and the types and the distributions of conserved motifs are similar within each sub-group. The transcriptome analysis revealed that genes encoding FADs exhibited different expression profiles under cold and heat stresses. The expression of ZmFAD2.1&2.2, ZmFAD7, and ZmSLD1&3 were significantly up-regulated under cold stress; moreover, the expression of ZmFAD2.1&2.3 and ZmSLD1&3 were obviously down-regulated under heat stress. The co-expression analysis demonstrated close correlation among the transcription factors and the significant responsive FAD genes under cold or heat stress. This study helps to understand the roles of plant FADs in temperature stress responses.

Genome-wide characterization and expression profiling of diacylglycerol acyltransferase genes from maize.

Diacylglycerol acyltransferase (DGAT) catalyzes the only rate-limiting step in the pathway of plant oil (TAG) biosynthesis and is involved in plant development. In this study, five DGAT family members were identified from maize genome database. Phylogenetic analysis classified the ZmDGATs into type-I, II, and III clusters. Conserved functional domain analysis revealed that the proteins encoded by ZmDGAT1 contained conserved MBOAT domains, while two ZmDGAT2-encoding proteins harbored LPLAT domains. qRT-PCR analysis showed that ZmDGAT genes exhibited very high relative expression in developing seeds, especially at the early stage of seed development. Under various abiotic stress conditions, differential responses of ZmDGAT genes were observed. An overall significant induction of ZmDGAT genes under cold stress in leaves and a quick and strong response to osmotic stresses in roots were highlighted. This study provides useful information for understanding the roles of DGATs in oil accumulation and stress responses in higher plants.

Glycerol-3-phosphate dehydrogenase (GPDH) gene family in Zea mays L.: Identification, subcellular localization, and transcriptional responses to abiotic stresses.

Glycerol-3-phosphate dehydrogenase (GPDH) catalyzes the formation of glycerol-3-phosphate, and plays an essential role in glycerolipid metabolism and in response to various stresses in different species. In this study, six ZmGPDH genes were obtained by a thorough search against maize genome, and designated as ZmGPDH1-6, respectively. The structural and evolutionary analyses showed that the ZmGPDHs family had typical conserved domains and similar protein structures as the known GPDHs from other plant species. ZmGPDHs were divided into NAD+-dependent type A form (ZmGPDH1-5) and FAD-dependent type B form (ZmGPDH6) based on their N-terminal sequences. Four full length ZmGPDHs were fused with GFP fusion proteins, and their subcellular localization was determined. ZmGPDH1 and ZmGPDH3 were located to the cytosol and mainly recruited to the surface of endoplasmic reticulum (ER), whereas ZmGPDH4 and ZmGPDH5 were located in the chloroplast. The transcriptional analysis of the ZmGPDHs in different maize tissues revealed relatively high level of transcripts accumulation of ZmGPDHs in roots and early stage developing seeds. Furthermore, we examined the transcriptional responses of the six GPDH genes in maize under various abiotic stresses, including salt, drought, alkali and cold, and significant induction of ZmGPDHs under osmotic stresses was observed. Together, this work will provide useful information for deciphering the roles of GPDHs in plant development and abiotic stress responses.