Natural variation of GmFATA1B regulates seed oil content and composition in soybean.

Soybean (Glycine max) produces seeds that are rich in unsaturated fatty acids and is an important oilseed crop worldwide. Seed oil content and composition largely determine the economic value of soybean. Due to natural genetic variation, seed oil content varies substantially across soybean cultivars. Although much progress has been made in elucidating the genetic trajectory underlying fatty acid metabolism and oil biosynthesis in plants, the causal genes for many quantitative trait loci (QTLs) regulating seed oil content in soybean remain to be revealed. In this study, we identified GmFATA1B as the gene underlying a QTL that regulates seed oil content and composition, as well as seed size in soybean. Nine extra amino acids in the conserved region of GmFATA1B impair its function as a fatty acyl-acyl carrier protein thioesterase, thereby affecting seed oil content and composition. Heterogeneously overexpressing the functional GmFATA1B allele in Arabidopsis thaliana increased both the total oil content and the oleic acid and linoleic acid contents of seeds. Our findings uncover a previously unknown locus underlying variation in seed oil content in soybean and lay the foundation for improving seed oil content and composition in soybean.

Transgenic soybean of GsMYB10 shapes rhizosphere microbes to promote resistance to aluminum (Al) toxicity.

Plant resistance genes could affect rhizosphere microbiota, which in turn enhanced plant resistance to stresses. Our previous study found that overexpression of the GsMYB10 gene led to enhanced tolerance of soybean plants to aluminum (Al) toxicity. However, whether GsMYB10 gene could regulate rhizosphere microbiota to mitigate Al toxicity remains unclear. Here, we analyzed the rhizosphere microbiomes of HC6 soybean (WT) and transgenic soybean (trans-GsMYB10) at three Al concentrations, and constructed three different synthetic microbial communities (SynComs), including bacterial, fungal and cross-kingdom (bacteria and fungi) SynComs to verify their role in improving Al tolerance of soybean. Trans-GsMYB10 shaped the rhizosphere microbial communities and harbored some beneficial microbes, such as Bacillus, Aspergillus and Talaromyces under Al toxicity. Fungal and cross-kingdom SynComs showed a more effective role than the bacterial one in resistance to Al stress, and these SynComs helped soybean resist Al toxicity via affecting some functional genes that involved cell wall biosynthesis and organic acid transport etc. Overall, this study reveals the mechanism of soybean functional genes regulating the synergistic resistance of rhizosphere microbiota and plants to Al toxicity, and also highlights the possibility of focusing on the rhizobial microbial community as a potential molecular breeding target to produce crops.

GmABR1 encoding an ERF transcription factor enhances the tolerance to aluminum stress in Arabidopsis thaliana.

The ethylene response factor (ERF) transcription factors, which is one of the largest transcription factor families in plants, are involved in biological and abiotic stress response and play an important role in plant growth and development. In this study, the GmABR1 gene from the soybean inbred line Zhonghuang24 (ZH24)×Huaxia 3 (HX3) was investigated its aluminum (Al) tolerance. GmABR1 protein has a conserved domain AP2, which is located in the nucleus and has transcriptional activation ability. The results of real-time quantitative PCR (qRT-PCR) showed that the GmABR1 gene presented a constitutive expression pattern rich in the root tip, stem and leaf tissues of HX3. After Al stress, the GmABR1 transcript was significantly increased in the roots. The transcripts of GmABR1 in the roots of HX3 treated with 50 µM AlCl3 was 51 times than that of the control. The GmABR1 was spatiotemporally specific with the highest expression levels when Al concentration was 50 µM, which was about 36 times than that of the control. The results of hematoxylin staining showed that the root tips of GmABR1-overexpression lines were stained the lightest, followed by the control, and the root tips of GmABR1 RNAi lines were stained the darkest. The concentrations of Al3+ in root tips were 207.40 µg/g, 147.74 µg/g and 330.65 µg/g in wild type (WT), overexpressed lines and RNAi lines, respectively. When AlCl3 (pH4.5) concentration was 100 µM, all the roots of Arabidopsis were significantly inhibited. The taproot elongation of WT, GmABR1 transgenic lines was 69.6%, 85.6%, respectively. When treated with Al, the content of malondialdehyde (MDA) in leaves of WT increased to 3.03 µg/g, while that of transgenic Arabidopsis increased from 1.66-2.21 µg/g, which was lower than that of WT. Under the Al stress, the Al stress responsive genes such as AtALMT1 and AtMATE, and the genes related to ABA pathway such as AtABI1, AtRD22 and AtRD29A were up-regulated. The results indicated that GmABR1 may jointly regulate plant resistance to Al stress through genes related to Al stress response and ABA response pathways.

GmGSTU23 Encoding a Tau Class Glutathione S-Transferase Protein Enhances the Salt Tolerance of Soybean (Glycine max L.).

Salt stress has a detrimental impact on crop yield, quality, and profitability. The tau-like glutathione transferases (GSTs) represent a significant group of enzymes that play a crucial role in plant stress responses, including salt stress. In this study, we identified a tau-like glutathione transferase family gene from soybean named GmGSTU23. Expression pattern analysis revealed that GmGSTU23 was predominantly expressed in the roots and flowers and exhibited a concentration-time-specific pattern in response to salt stress. Transgenic lines were generated and subjected to phenotypic characterization under salt stress. The transgenic lines exhibited increased salt tolerance, root length, and fresh weight compared to the wild type. Antioxidant enzyme activity and malondialdehyde content were subsequently measured, and the data revealed no significant differences between the transgenic and wild-type plants in the absence of salt stress. However, under salt stress, the wild-type plants exhibited significantly lower activities of SOD, POD, and CAT than the three transgenic lines, whereas the activity of APX and the content of MDA showed the opposite trend. We identified changes in glutathione pools and associated enzyme activity to gain insights into the underlying mechanisms of the observed phenotypic differences. Notably, under salt stress, the transgenic Arabidopsis's GST activity, GR activity, and GSH content were significantly higher than those of the wild type. In summary, our findings suggest that GmGSTU23 mediates the scavenging of reactive oxygen species and glutathione by enhancing the activity of glutathione transferase, thereby conferring enhanced tolerance to salt stress in plants.

MOTHER-OF-FT-AND-TFL1 regulates the seed oil and protein content in soybean.

Soybean is a major crop that produces valuable seed oil and protein for global consumption. Seed oil and protein are regulated by complex quantitative trait loci (QTLs) and have undergone intensive selections during the domestication of soybean. It is essential to identify the major genetic components and understand their mechanism behind seed oil and protein in soybean. We report that MOTHER-OF-FT-AND-TFL1 (GmMFT) is the gene of a classical QTL that has been reported to regulate seed oil and protein content in many studies. Mutation of MFT decreased seeds oil content and weight in both Arabidopsis and soybean, whereas increased expression of GmMFT enhanced seeds oil content and weight. Haplotype analysis showed that GmMFT has undergone selection, which resulted in the extended haplotype homozygosity in the cultivated soybean and the enriching of the oil-favorable allele in modern soybean cultivars. This work unraveled the GmMFT-mediated mechanism regulating seed oil and protein content and seed weight, and revealed a previously unknown function of MFT that provides new insights into targeted soybean improvement and breeding.

Identification and Mapping of QTLs for Sulfur-Containing Amino Acids in Soybean (Glycine max L.).

Soybean is a major source of high-quality protein for humans and animals. The content of sulfur-containing amino acids (SAA) in soybean is insufficient, which has become the main factor limiting soybean nutrition. In this study, we used the high-density genetic maps derived from Guizao 1 and Brazil 13 to evaluate the quantitative trait loci of cysteine (Cys), methionine (Met), SAA, glycinin (7S), ß-conglycinin (11S), ratio of glycinin to ß-conglycinin (RGC), and protein content (PC). In genetic map linkage analysis, the major and stable 44 QTLs were detected, which shared nine bin intervals. Among them, the bin interval (bin157-bin160) on chromosome 5 was detected in multiple environments as a stable QTL, which was linked to 11S, 7S, RGC, and SSA. Based on the analysis of bioinformatics and RNA-sequencing data, 16 differentially expressed genes (DEGs) within these QTLs were selected as candidate genes. These results will help to elucidate the genetic mechanism of soybean SAA-related traits and provide the basis for the gene mining of sulfur-containing amino acids.

GmWRKY21, a Soybean WRKY Transcription Factor Gene, Enhances the Tolerance to Aluminum Stress in Arabidopsis thaliana.

The WRKY transcription factors (TFs) are one of the largest families of TFs in plants and play multiple roles in plant growth and development and stress response. In this study, GmWRKY21 encoding a WRKY transcription factor was functionally characterized in Arabidopsis and soybean. The GmWRKY21 protein containing a highly conserved WRKY domain and a C2H2 zinc-finger structure is located in the nucleus and has the characteristics of transcriptional activation ability. The GmWRKY21 gene presented a constitutive expression pattern rich in the roots, leaves, and flowers of soybean with over 6-fold of relative expression levels and could be substantially induced by aluminum stress. As compared to the control, overexpression of GmWRKY21 in Arabidopsis increased the root growth of seedlings in transgenic lines under the AlCl3 concentrations of 25, 50, and 100 µM with higher proline and lower MDA accumulation. The results of quantitative real-time polymerase chain reaction (qRT-PCR) showed that the marker genes relative to aluminum stress including ALMT, ALS3, MATE, and STOP1 were induced in GmWRKY21 transgenic plants under AlCl3 treatment. The stress-related genes, such as KIN1, COR15A, COR15B, COR47, GLOS3, and RD29A, were also upregulated in GmWRKY21 transgenic Arabidopsis under aluminum stress. Similarly, stress-related genes, such as GmCOR47, GmDREB2A, GmMYB84, GmKIN1, GmGST1, and GmLEA, were upregulated in hair roots of GmWRKY21 transgenic plants. In summary, these results suggested that the GmWRKY21 transcription factor may promote the tolerance to aluminum stress mediated by the pathways regulating the expression of the acidic aluminum stress-responsive genes and abiotic stress-responsive genes.

GsMYB7 encoding a R2R3-type MYB transcription factor enhances the tolerance to aluminum stress in soybean (Glycine max L.).

BACKGROUND: MYB transcription factor (TF) is one of the largest families of TFs in plants and play essential roles in plant growth and development, and is involved in responses to biological and abiotic stress. However, there are few reports on GsMYB7 gene in soybean under aluminum acid stress, and its regulatory mechanism remains unclear. RESULTS: The GsMYB7 protein is localized in the nucleus and has transcriptional activation ability. Quantitative real-time PCR (qRT-PCR) results showed that GsMYB7 held a constitutive expression pattern rich in roots. When AlCl3 concentration was 25 µM, the total root surface area (SA) of GsMYB7 transgenic lines were 34.97% higher than that of wild-type Huachun 6 (HC6). While the accumulation of Al3+ in root tip of transgenic plants after aluminum treatment was 17.39% lower than that of wild-type. RNA-sequencing analysis indicated that over 1181 genes were regulated by GsMYB7 and aluminum stress. Among all the regulated genes, the expression levels of glutathione peroxidase, protein kinase, cytochrome and other genes in the transgenic lines were significantly higher than those in wild type by acidic aluminum stress. The bioinformatics and qRT-PCR results showed that 9 candidate genes were induced under the treatments of acidic aluminum stress which were indirectly and/or directly regulated by GsMYB7. After AlCl3 treatments, the transcripts of these genes in GsMYB7 transgenic seedlings were significantly higher than those of wide-type HC6. CONCLUSIONS: The results suggested that GsMYB7 may enhance soybean tolerance to acidic aluminum stress by regulating the downstream genes.

Silicon-enhanced tolerance to cadmium toxicity in soybean by enhancing antioxidant defense capacity and changing cadmium distribution and transport.

Cadmium (Cd) is a widely distributed heavy metal that is toxic to plants and humans. Although silicon (Si) has been reported to reduce Cd accumulation and toxicity in plants, evidence on the functions of Si and its mechanisms in the possible alleviation of soybean are limited. Therefore, a controlled experiment was conducted to investigate the impacts and mechanisms of Si on Cd retention in soybean. Here, we determined the growth index, Cd distribution, and antioxidant activity systems of Si, as well as expression levels of differentially expressed genes (DEGs) in Si under Cd stress, and conducted RNA-seq analysis. We not only found that Si can significantly promote soybean plant growth, increase plant antioxidant activities, and reduce the Cd translocation factor, but also revealed that a total of 636 DEGs were shared between CK and Cd, CK and Cd + Si, and Cd and Cd + Si. Moreover, several genes were significantly enriched in antioxidant systems and Cd distribution and transport systems. Therefore, the expression status of Si-mediated Cd stress response genes is likely involved in improving oxidative stress and changing Cd uptake and transport, as well as improving plant growth that contributes to Si alleviating Cd toxicity in plants. Moreover, numerous potential target genes were identified for the engineering of Cd-tolerant cultivars in soybean breeding programs.

GmWRKY81 Encoding a WRKY Transcription Factor Enhances Aluminum Tolerance in Soybean.

Aluminum (Al) toxicity is an essential factor that adversely limits soybean (Glycine max (L.) Merr.) growth in acid soils. WRKY transcription factors play important roles in soybean responses to abiotic stresses. Here, GmWRKY81 was screened from genes that were differentially expressed under Al treatment in Al-tolerant soybean Baxi10 and Al-sensitive soybean Bendi2. We found that GmWRKY81 was significantly induced by 20 µM AlCl3 and upregulated by AlCl3 treatment for 2 h. In different tissues, the expression of GmWRKY81 was differentially induced. In 0-1 cm root tips, the expression of GmWRKY81 was induced to the highest level. The overexpression of GmWRKY81 in soybean resulted in higher relative root elongation, root weight, depth, root length, volume, number of root tips and peroxidase activity but lower root average diameter, malonaldehyde and H2O2 contents, indicating enhanced Al tolerance. Moreover, RNA-seq identified 205 upregulated and 108 downregulated genes in GmWRKY81 transgenic lines. Fifteen of these genes that were differentially expressed in both AlCl3-treated and GmWRKY81-overexpressing soybean had the W-box element, which can bind to the upstream-conserved WRKY domain. Overall, the combined functional analysis indicates that GmWRKY81 may improve soybean Al tolerance by regulating downstream genes participating in Al3+ transport, organic acid secretion and antioxidant reactions.

GsERF1 enhances Arabidopsis thaliana aluminum tolerance through an ethylene-mediated pathway.

Ethylene response factor (ERF) transcription factors constitute a subfamily of the AP2/ERF superfamily in plants and play multiple roles in plant growth and development as well as in stress responses. In this study, the GsERF1 gene from the wild soybean BW69 line (an Al-resistant Glycine soja line) was rapidly induced in response to aluminum stress. Quantitative real-time PCR (qRT-PCR) analysis showed that the GsERF1 gene maintained a constitutive expression pattern and was induced in soybean in response to aluminum stress, with increased amounts of transcripts detected in the roots. The putative GsERF1 protein, which contains an AP2 domain, was located in the nucleus and maintained transactivation activity. In addition, under AlCl3 treatment, GsERF1 overexpression increased the relative growth rate of the roots of Arabidopsis and weakened the hematoxylin staining of hairy roots. Ethylene synthesis-related genes such as ACS4, ACS5 and ACS6 were upregulated in GsERF1 transgenic lines compared with the wild type under AlCl3 treatment. Furthermore, the expression levels of stress/ABA-responsive marker genes, including ABI1, ABI2, ABI4, ABI5 and RD29B, in the GsERF1 transgenic lines were affected by AlCl3 treatment, unlike those in the wild type. Taken together, the results indicated that overexpression of GsERF1 may enhance aluminum tolerance of Arabidopsis through an ethylene-mediated pathway and/or ABA signaling pathway, the findings of which lay a foundation for breeding soybean plants tolerant to aluminum stress.

Identification of quantitative trait loci (QTLs) and candidate genes of seed Iron and zinc content in soybean [Glycine max (L.) Merr.].

BACKGROUND: Deciphering the hereditary mechanism of seed iron (Fe) and zinc (Zn) content in soybean is important and sustainable to address the "hidden hunger" that presently affects approximately 2 billion people worldwide. Therefore, in order to detect genomic regions related to soybean seed Fe and Zn content, a recombinant inbred line (RIL) population with 248 lines was assessed in four environments to detect Quantitative Trait Loci (QTLs) related to soybean seed Fe and Zn content. RESULT: Wide variation was found in seed Fe and Zn content in four environments, and genotype, environment, and genotype × environment interactions had significant influences on both the seed Fe and Zn content. A positive correlation was observed between seed Fe content and seed Zn content, and broad-sense heritability (H2) of seed Fe and Zn content were 0.73 and 0.75, respectively. In this study, five QTLs for seed Fe content were detected with 4.57 - 32.71% of phenotypic variation explained (PVE) and logarithm of odds (LOD) scores ranging from 3.60 to 33.79. Five QTLs controlling the seed Zn content were detected, and they individually explained 3.35 to 26.48% of the phenotypic variation, with LOD scores ranging from 3.64 to 20.4. Meanwhile, 409,541 high-quality single-nucleotide variants (SNVs) and 85,102 InDels (except intergenic regions) between two bi-parental lines were identified by whole genome resequencing. A total of 12 candidate genes were reported in one major QTL for seed Fe content and two major QTLs for seed Zn content, with the help of RNA-Seq analysis, gene ontology (GO) enrichment, gene annotation, and bi-parental whole genome sequencing (WGS) data. CONCLUSIONS: Limited studies were performed about microelement of soybean, so these results may play an important role in the biofortification of Fe and Zn and accelerate the development of marker-assisted selection (MAS) for breeding soybeans fortified with iron and zinc.

Identification of an ATP-Binding Cassette Transporter Implicated in Aluminum Tolerance in Wild Soybean (Glycine soja).

The toxicity of aluminum (Al) in acidic soil limits global crop yield. The ATP-binding cassette (ABC) transporter-like gene superfamily has functions and structures related to transportation, so it responds to aluminum stress in plants. In this study, one half-size ABC transporter gene was isolated from wild soybeans (Glycine soja) and designated GsABCI1. By real-time qPCR, GsABCI1 was identified as not specifically expressed in tissues. Phenotype identification of the overexpressed transgenic lines showed increased tolerance to aluminum. Furthermore, GsABCI1 transgenic plants exhibited some resistance to aluminum treatment by ion translocation or changing root components. This work on the GsABCI1 identified the molecular function, which provided useful information for understanding the gene function of the ABC family and the development of new aluminum-tolerant soybean germplasm.

Modulation of evening complex activity enables north-to-south adaptation of soybean.

Soybean, a typical short-day crop, is sensitive to photoperiod, which is a major limiting factor defining its north-to-south cultivation range. The long-juvenile (LJ) trait is controlled primarily by the J locus which has been used for decades by soybean breeders to delay flowering and improve grain yield in tropical regions. The J gene encodes an ortholog of the Arabidopsis Evening Complex (EC) component EARLY FLOWERING 3 (ELF3). To identify modifiers of J, we conducted a forward genetic screen and isolated a mutant (eoj57) that in combination with j has longer flowering delay compared with j single mutant plants. Map-based cloning and genome re-sequencing identified eoj57 (designated as GmLUX2) as an ortholog of the Arabidopsis EC component LUX ARRHYTHMO (LUX). To validate that GmLUX2 is a modifier of J, we used trans-complementation and identified a natural variant allele with a similar phenotype. We also show that GmLUX2 physically interacts with GmELF3a/b and binds DNA, whereas the mutant and natural variant are attenuated in both activities. Transcriptome analysis shows that the GmLUX2-GmELF3a complex co-regulates the expression of several circadian clock-associated genes and directly represses E1 expression. These results provide mechanistic insight into how GmLUX2-GmELF3 controls flowering time via synergistic regulation of gene expression. These novel insights expand our understanding of the regulation of the EC complex, and facilitate the development of soybean varieties adapted for growth at lower latitudes.

Wild Soybean Oxalyl-CoA Synthetase Degrades Oxalate and Affects the Tolerance to Cadmium and Aluminum Stresses.

Acyl activating enzyme 3 (AAE3) was identified as being involved in the acetylation pathway of oxalate degradation, which regulates the responses to biotic and abiotic stresses in various higher plants. Here, we investigated the role of Glycine sojaAAE3 (GsAAE3) in Cadmium (Cd) and Aluminum (Al) tolerances. The recombinant GsAAE3 protein showed high activity toward oxalate, with a Km of 105.10 ± 12.30 µM and Vmax of 12.64 ± 0.34 µmol min-1 mg-1 protein, suggesting that it functions as an oxalyl-CoA synthetase. The expression of a GsAAE3-green fluorescent protein (GFP) fusion protein in tobacco leaves did not reveal a specific subcellular localization pattern of GsAAE3. An analysis of the GsAAE3 expression pattern revealed an increase in GsAAE3 expression in response to Cd and Al stresses, and it is mainly expressed in root tips. Furthermore, oxalate accumulation induced by Cd and Al contributes to the inhibition of root growth in wild soybean. Importantly, GsAAE3 overexpression increases Cd and Al tolerances in A. thaliana and soybean hairy roots, which is associated with a decrease in oxalate accumulation. Taken together, our data provide evidence that the GsAAE3-encoded protein plays an important role in coping with Cd and Al stresses.

High Aluminum Drives Different Rhizobacterial Communities Between Aluminum-Tolerant and Aluminum-Sensitive Wild Soybean.

Aluminum (Al)-resistant plant cultivars can recruit beneficial microbes to alleviate the stresses. However, the mechanism of how rhizobacterial communities strengthen Al tolerance of wild soybean has not been addressed. The aim of this study was to investigate the bacterial community structure in the rhizosphere of Al-tolerant (BW69) and Al-sensitive (W270) wild soybean germplasm subjected to three Al concentrations. We analyzed the rhizobacterial communities of the two genotypes by high-throughput sequencing of 16S rRNA genes. The results showed that high Al stress recruited different rhizobacterial communities between two genotypes. In total, 49 OTUs, such as OTU15 (Gammaproteobacteria_KF-JG30-C25_norank), OTU23 (Mizugakiibacter), and OTU93 (Alkanibacter), were enriched in the rhizosphere of BW69 at the low and high Al concentrations. Moreover, bacterial community in the rhizosphere of BW69 had a more complex co-occurrence network than did W270 at the high Al concentration. Overall, our findings highlighted that high Al concentration magnified the difference in rhizobacterial community structure between two genotypes. However, the lower modularity of the co-occurrence network in rhizosphere of BW69 than W270 under Al stress may cause the rhizobacterial community to be less resistant and more influenced by disturbance. This study emphasizes the possibility of using rhizobacteria as an improved crop breeding or gene to produce crops that are more resistant to the toxicity of heavy metal.

Transcription Factor GmWRKY142 Confers Cadmium Resistance by Up-Regulating the Cadmium Tolerance 1-Like Genes.

Cadmium (Cd) is a widespread pollutant that is toxic to living organisms. Previous studies have identified certain WRKY transcription factors, which confer Cd tolerance in different plant species. In the present study, we have identified 29 Cd-responsive WRKY genes in Soybean [Glycine max (L.) Merr.], and confirmed that 26 of those GmWRKY genes were up-regulated, while 3 were down-regulated. We have also cloned the novel, positively regulated GmWRKY142 gene from soybean and investigated its regulatory mechanism in Cd tolerance. GmWRKY142 was highly expressed in the root, drastically up-regulated by Cd, localized in the nucleus, and displayed transcriptional activity. The overexpression of GmWRKY142 in Arabidopsis thaliana and soybean hairy roots significantly enhanced Cd tolerance and lead to extensive transcriptional reprogramming of stress-responsive genes. ATCDT1, GmCDT1-1, and GmCDT1-2 encoding cadmium tolerance 1 were induced in overexpression lines. Further analysis showed that GmWRKY142 activated the transcription of ATCDT1, GmCDT1-1, and GmCDT1-2 by directly binding to the W-box element in their promoters. In addition, the functions of GmCDT1-1 and GmCDT1-2, responsible for decreasing Cd uptake, were validated by heterologous expression in A. thaliana. Our combined results have determined GmWRKYs to be newly discovered participants in response to Cd stress, and have confirmed that GmWRKY142 directly targets ATCDT1, GmCDT1-1, and GmCDT1-2 to decrease Cd uptake and positively regulate Cd tolerance. The GmWRKY142-GmCDT1-1/2 cascade module provides a potential strategy to lower Cd accumulation in soybean.

Identification and Mapping of Stable QTLs for Seed Oil and Protein Content in Soybean [Glycine max (L.) Merr.].

This research aimed to identify stable quantitative trait loci (QTL) associated with oil and protein content in soybean. A population of 196 recombinant inbred lines (RILs) derived from Huachun 2 × Wayao was used to evaluate these target traits. A high-density genetic linkage map was constructed by using high-throughput genome-wide sequencing technology, which contained 3413 recombination bin markers and spanned 5400.4 cM with an average distance of 1.58 cM between markers. Eighteen stable QTLs controlling oil and protein content were detected. Among them, qOil-11-1 was identified for the first time as a novel QTL, while qOil-5-1, qPro-10-1, and qPro-14-1 were strong and stable QTLs with high log-likelihood (LOD) values. Sixteen differentially expressed genes (DEGs) within these four QTLs were shown to be highly expressed during seed development based on RNA sequencing (RNA-seq) data analysis. Our results may contribute toward gene mining and marker-assisted selection (MAS) for breeding a high-quality soybean in the future.