Identification of Major Quantitative Trait Loci for Seed Oil Content in Soybeans by Combining Linkage and Genome-Wide Association Mapping

Soybean oil is the most widely produced vegetable oil in the world and its content in soybean seed is an important quality trait in breeding programs. More than 100 quantitative trait loci (QTLs) for soybean oil content have been identified. However, most of them are genotype specific and/or environment sensitive. Here, we used both a linkage and association mapping methodology to dissect the genetic basis of seed oil content of Chinese soybean cultivars in various environments in the Jiang-Huai River Valley. One recombinant inbred line (RIL) population (NJMN-RIL), with 104 lines developed from a cross between M8108 and NN1138-2, was planted in five environments to investigate phenotypic data, and a new genetic map with 2,062 specific-locus amplified fragment markers was constructed to map oil content QTLs. A derived F₂ population between MN-5 (a line of NJMN-RIL) and NN1138-2 was also developed to confirm one major QTL. A soybean breeding germplasm population (279 lines) was established to perform a genome-wide association study (GWAS) using 59,845 high-quality single nucleotide polymorphism markers. In the NJMN-RIL population, 8 QTLs were found that explained a range of phenotypic variance from 6.3 to 26.3% in certain planting environments. Among them, qOil-5-1, qOil-10-1, and qOil-14-1 were detected in different environments, and qOil-5-1 was further confirmed using the secondary F₂ population. Three loci located on chromosomes 5 and 20 were detected in a 2-year long GWAS, and one locus that overlapped with qOil-5-1 was found repeatedly and treated as the same locus. qOil-5-1 was further localized to a linkage disequilibrium block region of approximately 440 kb. These results will not only increase our understanding of the genetic control of seed oil content in soybean, but will also be helpful in marker-assisted selection for breeding high seed oil content soybean and gene cloning to elucidate the mechanisms of seed oil content.

Introduction

Soybean, Glycine max (L.) Merr, is an important grain and oil seed crop that is widely grown throughout the world, and provides more than 50% of the world’s oilseed production (Wilson, 2008; Soy Stats, 2011). As one of the determinants of oil production, soybean seed oil is an important trait targeted by many breeding programmers. Seed oil accounts for approximately 20% of the seed weight and it is a heritable quantitative trait in soybean. It is controlled by multiple quantitative trait loci (QTLs)/genes, which have small effects, and is also affected by the environment and its interaction with QTLs (Diers et al., 1992). Although seed oil content can be increased in traditional breeding based on phenotypic data selection (Burton and Brim, 1981), marker-assisted selection (MAS) is an efficient method that is independent of time and location (Allen, 1994; Chapman et al., 2003). Therefore, numerous studies have considered the genetics and QTL composition of seed oil content in soybean.

In the past 30 years, numerous QTLs for seed oil content have been identified in soybean (Diers et al., 1992; Lee et al., 1996; Orf et al., 1999; Csanadi et al., 2001; Chung et al., 2003; Hyten et al., 2004; Panthee et al., 2005; Pathan et al., 2013). However, most of these QTLs have low selection accuracy and have not been used effectively in MAS to breed for high seed oil content in soybean varieties (Collard and Mackill, 2008). It is also difficult to narrow down the QTL interval to the gene level based on a single-linkage mapping experiment. Especially in soybean, the size of the linkage mapping population is typically small. The main reason for this limitation is that the linkage mapping method has a relatively low genome resolution, which only uses the recombination events within the mapping populations (Darvasi et al., 1993; Beavis, 1994).

Genome-wide association (GWAS) is another method that takes advantage of a number of historic recombination events that have occurred within natural populations, and can overcome the limitations of linkage analysis (Rafalski, 2010). GWAS has been successfully applied to crop plants such as rice and maize (Zhao et al., 2011; Yang et al., 2014). In soybean, Mamidi et al. (2014) performed a GWAS to identify significant markers associated with iron deficiency chlorosis trait variations and identified seven major QTLs on seven chromosomes. Zhang et al. (2015) conducted a GWAS to dissect the genetic architecture of some agronomically important traits in early maturity soybean accessions, and detected 27, 6, 18, and 27 loci for days to flowering, days to maturity, duration of flowering-to-maturity, and plant height, respectively. Despite the fact that association studies have advantages in QTL mapping, the population structure (Q) and individual relationships may identify false positive correlations between the markers and the phenotype. Therefore, it is necessary to conduct linkage mapping in order to verify the results of a GWAS for dissecting complex traits.

Compared with traditional markers, single nucleotide polymorphisms (SNPs) have a more abundant DNA variation. Following the completion of the whole genome sequencing of soybean cv. Williams 82 (Schmutz et al., 2010) and the rapid development in sequencing technology, SNP markers are beginning to be used in soybean (Hyten et al., 2008). The next generation sequencing technologies enable researchers to quickly obtain numbers of SNPs throughout the genome, and have been successfully used in soybean for high-density genetic map construction and GWASs (Hwang et al., 2014; Song et al., 2016; Wang et al., 2016; Contreras-Soto et al., 2017). A specific length amplified fragment sequencing (SLAF-seq) technology with greater genotyping accuracy and relatively lower sequencing cost has been developed (Sun et al., 2013) and has been used to create high-density genetic maps in several plants, such as rice (Mao et al., 2015), soybean (Zhang D. et al., 2016; Cao et al., 2017), sesame (Zhang et al., 2013), and kiwifruit (Huang et al., 2013).

Seed oil content has a wide range of genetic variation among soybean accessions (Wilson, 2004) and there is also a strong negative correlation with seed protein content (Wilcox, 1998). Therefore, several QTLs can be related to both seed protein and oil content. Two major QTLs for seed oil/protein content have been mapped in many soybean populations: one was mapped at an interval from 20 to 40 centimorgan (cM) on chromosome (Chr.) 20 of the Soybean GmConsensus 4.0 Map and the other one was mapped at an interval from 10 to 30 cM on Chr. 15¹. Many studies have focused on these two QTLs. The QTL on Chr. 20 was narrowed down to a small region (Nichols et al., 2006; Bolon et al., 2010; Hwang et al., 2014; Vaughn et al., 2014); and the QTL on Chr. 15 was fine mapped at a 535 kb interval (Kim et al., 2016). In contrast, fewer studies have focused on refining the position of the QTL for oil content on Chr. 05, which has been reported many times within a large genomic region (Mansur et al., 1996; Brummer et al., 1997; Orf et al., 1999; Pathan et al., 2013).

In the Jiang-Huai River Valley in China, soybean is grown from approximately June to October. However, little is known about the genetics of the seed oil content in these genotypes. The objectives of our study were to utilize the linkage and germplasm populations in various environments to map QTLs for seed oil content in soybean cultivars based on SNP markers, then dissect the genetic basis of seed oil content in these cultivars, and finally to refine a major QTL for seed oil content on Chr. 05 by combined linkage and GWA mapping.

Materials and Methods

Plant Materials and Plant Growth Conditions

For linkage mapping, two populations were used in this study. A recombinant inbred line (RIL) population was used for the initial QTL detection. The RIL population (NJMN-RIL) consisted of 104 F_2:8 lines developed from a cross between M8108 (a landrace material, seed oil content exceeding 20%) and NN1138-2 (an elite cultivar characterized by high yield, seed oil content ∼ 19.0%). NJMN-RIL and its two parents were grown in five environments: Jiangpu Experimental Station, Nanjing, Jiangsu Province, in 2012 and 2014 (12JP and 14JP); Fengyang Experimental Station, Chuzhou, Anhui Province, in 2012 (12FY); Yancheng Experimental Station, Yangcheng, Jiangsu Province, in 2014 (14YC); and Huaian Experimental Station, Huaian, Jiangsu Province, in 2014 (14HA).

A secondary F₂ population was used to confirm the most stable QTL from the initial mapping results. The F₂ population consisted of 203 individuals developed from a backcross using MN-5 with NN1138-2. MN-5 was a line from the NJMN-RIL population with the M8108 type–homozygous qOil-5-1 region, while the NN1138-2 type–homozygous was the other major QTL region identified in this study. The F₂ population was grown at Jiangpu Experimental Station in 2015.

A total of 279 soybean accessions from the Yangtze-Huai soybean breeding germplasm population were used to perform the association study. The genotype of these soybean accessions can be obtained from our previous study (Li et al., 2016). The soybean accessions were grown in two environments: 13JP and 14JP.

The NJMN-RIL and germplasm populations were grown in a randomized complete block design, with three replications from approximately June to October. A 1 m single row plot of 10 plants per RIL/accession was grown, with a distance of 10 cm between plants and a row spacing of 50 cm. The secondary F₂ population was grown in rows with a length of 2 m, with 50 cm between rows and 20 cm between plants. Field management was performed under normal conditions.

Measurement and Analysis of Phenotypic Data

Seed oil content was measured by an Infratec^TM1241 near infrared analysis (NIR) Grain Analyzer (Foss, Hillerød, Denmark) on a 10% moisture basis using approximately 15–20 g samples. Seed oil content for each RIL and soybean accession was an average of three replications within a planting environment.

The frequency distribution of phenotypic data was calculated using the SPSS Statistics 20.0 software (SPSS Inc., Chicago, IL, United States). The broad-sense heritability (h²) and analysis of variance (ANOVA) tests were conducted using the SAS PROC generalized linear model (GLM) program. The h² of seed oil content in RIL and germplasm population was estimated using the following equation:

$$h^2 = {\sigma }_g^2/({\sigma }_g^2+{\sigma }_{ge}^2/n+{\sigma }_e^2/nr)$$

where σ²_g is the genotypic variance, σ²_ge is the variance of the genotype-by-environment interaction, σ²_e is the error variance, n is the number of environments, and r is the number of replications within an environment (Nyquist and Baker, 1991).

Genotyping

For the NJMN-RIL population, approximately 1 g of fresh leaves obtained from two parents and individuals of each RIL were used to extract the genomic DNA using the cetyltrimethylammonium bromide method (Doyle, 1990). All RIL individuals and the two parents were genotyped using SLAF-seq technology to generate genome-wide SLAF markers. SLAF library construction, high-throughput sequencing, and high-quality SLAF marker acquisition were performed as described by Sun et al. (2013). The soybean reference genome (G. max Wm82.a1) was used to determine the physical position of the high-quality SLAF markers.

For the secondary F₂ population, a total of 31 pairs of simple sequence repeat (SSR) markers based on the soybean SSR database (Song et al., 2010) and four pairs of InDel markers based on the results of Lam et al. (2010) and distributed across the qOil-5-1 region (38–40.5 Mb on the Chr. 05), were used to screen for polymorphisms between MN-5 and NN1138-2. Finally, 10 polymorphic pairs (one InDel marker and nine SSR markers) were used to genotype the F₂ population and then construct a genetic map covering the qOil-5-1 region.

A set of 279 soybean germplasm lines was assembled for GWAS. All lines were genotyped for SNP markers using high-throughput genotyping platforms, as described by Li et al. (2016). Finally, 59,845 SNPs, with minor allele frequencies (MAF) > 0.05, were available for the GWAS.

Construction of Linkage Maps for the NJMN-RIL Population and QTL Linkage Mapping

After genotyping the 104 RILs, all polymorphic SLAF markers were filtered four times and further quality assessed as described by Sun et al. (2013). A SLAF with less than three SNPs and an average depth for each sample above three was considered a high-quality SLAF marker. Then, all high-quality SLAF markers were used for construction of the linkage map. The linkage map was constructed using High Map software (Liu et al., 2014). The markers were divided into linkage groups (LGs) using the single-linkage clustering algorithm at a logarithm of odds (LOD) threshold ≥3.0. The Kosambi mapping function (Kosambi, 1943) was used to calculate the map distances in cM from the recombination frequencies.

The QTLs were detected using WinQTLCart 2.5 software. The window size, the working speed, and the control marker number were set at 10 cM, 1 cM and 5, respectively. Model 6 (standard model) for composite interval mapping was used to identify QTLs in each environment, with a log likelihood of 2.5 set as the threshold of an existing QTL (Wang et al., 2007).

The Q and Linkage Disequilibrium (LD) of the Soybean Breeding Germplasm Population

The Q and the kinship matrix (K) for these 279 soybean accessions were calculated as previously described (Li et al., 2016). The 279 soybean materials could be divided into three major subpopulations. The LD was approximately 480 kb, whereas r² dropped to half its maximum value (Li et al., 2016). The LD blocks were within the interval of a QTL based on an LD analysis using the r² correlation between each marker using the HAPLOVIEW 4.2 software (Barrett et al., 2005). The parameters in the program included MAF (≥0.05) and the integrity of each SNP (≥50%).

Genome-Wide Association Mapping

In this study, the GWAS was performed using the TASSEL 5.0 software (Bradbury et al., 2007). Two statistical models were considered: (1) the GLM model without the Q and K; and (2) the mixed linear model (MLM) with the Q and K, which regarded Q and K as fixed and random effects. Compared with the GLM method, the MLM model can reduce the number of false positive results and was considered more appropriate for this study. To obtain more reliable results, the Bonferroni threshold (P < 1/59845, -log₁₀P > 4.78) was used as the threshold to identify significant association SNPs. The association region was also estimated by extending 480 kb (LD distance) upstream and downstream of the most significantly associated SNP position.

Results

Phenotypic Evaluation of Seed Oil Content

The phenotypic performance is presented in Table 1 and Supplementary Figure S1. For the NJMN-RIL population, in the female parent, M8108 had an average seed oil content of ∼20.5%, and NN1138-2 of ∼19%, with these differences being significant in all environments (P < 0.001). A continuous distribution and transgressive segregation were observed for the trait in the RIL population in all environments. Due to the kurtosis and skewness (absolute value) being less than 1, the phenotypic frequency distributions had a relatively normal distribution. This indicates that the phenotypic data was suitable for QTL mapping. For the germplasm population, the seed oil content ranged from 15.63 to 22.90% and from 14.67 to 22.63% in 2013 and 2014, respectively. The phenotypic frequency distributions of the germplasm population also had a normal distribution.

TABLE 1

TABLE 1. Results of seed oil content in the linkage mapping and soybean germplasm populations.

In addition, the ANOVA results indicated significant differences in genotype, environment, and genotype-by-environment (Supplementary Table S1). The h² of seed oil content in the NJMN-RIL population and germplasm population was high (0.93 and 0.88, respectively), which was consistent with previous studies (Chung et al., 2003; Panthee et al., 2005).

Linkage Maps and QTL Mapping of Seed Oil Content in the NJMN-RIL Population

The RIL population was genotyped using SLAF-seq technology. After being filtered and quality assessed, a total of 2,086 SLAF markers were available to construct genetic maps for the NJMN-RIL population. Finally, 2,062 SLAF markers were grouped into 20 LGs. The total genetic distance of this map was 2054.50 cM. The average distance between adjacent markers was approximately 1 cM. The genetic length of 20 LGs ranged from 40.21 cM (Chr. 17) to 199.77 cM (Chr. 12). The largest LG was Chr. 18, with 338 SLAF markers. The LG with the minimum number of markers was Chr. 05, with 28 SLAF markers. The mean LG length was 102.73 cM and each chromosome contained an average of 103 markers. Detailed map information is provided in Figure 1 and Supplementary Table S2.

FIGURE 1

FIGURE 1. Distribution of markers in 20 linkage groups (LGs) in the NJMN-RIL population. The black bars in each LG represent mapped SLAF-seq markers. Detailed information is presented in Supplementary Table S2.

Based on the genetic map, a total of eight QTLs for seed oil content, distributed on seven chromosomes, were identified in the NJMN-RIL population across the different environments by WinQTLCart 2.5 software, with LOD scores ranging from 2.7 to 9.1 and 6.3 to 26.2% of the phenotypic variation explained by individual QTLs (Table 2). Among them, qOil-5-1 was the most stable QTL detected in all environments and could explain 7.5–26.2% of the phenotypic variation. The qOil-5-1 QTL was mapped at an approximate 6.8 cM (43.4–50.2 cM) interval in all environments (Figure 2A and Table 2). This QTL interval has a corresponding physical distance of 2.5 Mb, from 38.04 to 40.51 Mb on Chr. 05 (Figure 2B). The qOil-10-1 and qOil-14-1 QTLs could also be considered major QTLs because they were detected in four environments. The other five QTLs were detected in a single environment. In addition, the positive alleles of qOil-3-1 and qOil-10-1 came from NN1138-2, while the others came from M8108.

TABLE 2

TABLE 2. Analysis of quantitative trait loci (QTLs) for seed oil content in the NJMN-RIL population.

FIGURE 2

FIGURE 2. Detection and confirmation of qOil-5-1 using NJMN-RIL and secondary F₂ populations by linkage mapping. (A) qOil-5-1 was detected in the NJMN-RIL population in five environments. The physical location of the confidence interval of qOil-5-1 was approximately 2.5 Mb, from 38.04 to 40.51 Mb on chromosome 05. (B) Physical location of the 1- logarithm of odds (LOD) confidence interval of qOil-5-1 on chromosome 05. (C) qOil-5-1 was confirmed in the secondary F2 population. The peak marker was InDel-5-7. The 1- logarithm of odds (LOD) confidence interval of the QTL spanned approximately 1 Mb.

Confirmation of the qOil-5-1 QTL Using a Secondary F₂ Population

The qOil-5-1 QTL was detected in all environments and so could be considered the most stable QTL in the NJMN-RIL population. To confirm and refine the position of this QTL, the secondary F₂ population was used. The seed oil content ranged from 17.50 to 21.90% in this F₂ population (Table 1 and Supplementary Figure S1). We also constructed a genetic map with 10 polymorphic pairs for the F₂ population (Figure 2C). This genetic map covered the qOil-5-1 region (from ∼37.8 to 40.8 Mb on Chr. 05). Based on the phenotype and genetic map, a QTL was identified. The highest peak of the QTL was found at the marker InDel-5-7 [primer sequence: F′ (5′–3′): CAAAAATACATTGAATGGAAGGCAA; R′ (5′–3′): TGACAGAGAGATATTATCCCCAGC], with an LOD score of 6.87, which strongly exceeded the threshold level. The InDel-5-7 marker corresponded to position 38494541-38494952 bp on the reference genome (G. max Wm82.a1). The 1-LOD confidence interval of the QTL spanned approximately 1 Mb, from position 37.85 Mb (BARCSOYSSR_05_1217) to 38.90 Mb (BARCSOYSSR_05_1261) of the physical map, which was consistent with the RIL population mapping results (Figures 2B,C).

Genome-Wide Association Study (GWAS) for Seed Oil Content

In this study, 279 soybean gremplasm lines with 59,845 (MAF > 0.05) SNPs were used to perform the association analysis. The quantile–quantile (QQ) plot showed that the MLM method was appropriate for this study and could reduce the number of false positive results (Supplementary Figure S2). Therefore, the MLM model was used for the association analysis.

There were two association peaks identified on Chr. 05 and on Chr. 20 in 2013; and one association peak on Chr.05 and two association peaks on Chr. 20 were identified in 2014. The phenotypic variation explained by each peak ranged from 8.5 to 12.9% (Table 3 and Figure 3A).

TABLE 3

TABLE 3. Results of a genome-wide association study (GWAS) for seed oil content.

FIGURE 3

FIGURE 3. Genome-wide association study (GWAS) results for seed oil content in the germplasm population. (A) Manhattan plot of the seed oil content in two different environments. The dashed horizontal line depicts a significant threshold level [–log₁₀ (1/59845)]. (B) GWAS results of the 36–42 Mb genomic region of chromosome 05. (C) The distribution of linkage disequilibrium (LD) blocks of the major locus on chromosome 05. All the single nucleotide polymorphisms (SNPs) that had a significant association with seed oil content in the 2 years investigated were located in one LD block from 38.13 to 38.57 Mb.

The loci for seed oil content on Chr. 05 could be detected in both the linkage and association population. Therefore, qOil-5-1 can be considered an important locus for seed oil content in the summer planting soybeans. To reduce the qOil-5-1 region, we focused on the results from the association population (Figure 3B). The peaks on Chr. 05 were located at 38506373 and 38473956 bp in 2013 and 2014, respectively. Based on the LD distance, we extended the region of 480 kb upstream and downstream of the significantly associated SNP position. The LD was estimated from the r² correlation between each marker in this region, with r² values less than 0.2 considered unlinked. This region possessed two large LD blocks, one from 38.13 to 38.57 Mb and the other from 38.59 to 38.95 Mb (Figure 3C). However, all the SNPs showing a significant association with seed oil content in the 2 years were located within the first LD block region. This suggested that the major QTL for seed oil content was located in a small region from 38.13 to 38.57 Mb on Chr. 05.

Discussion

Optimizing the Population and Marker System for Mapping QTLs of Complex Traits

Quantitative trait loci mapping is an effective approach to analyze quantitative traits in plants. Linkage mapping and GWAS are commonly used to map the QTLs of complex quantitative traits. However, the linkage mapping method has a relatively low genome resolution and the GWAS results are affected by the Q and individual relationships. A combination of GWAS and linkage mapping could provide mutual authentication and enable more reliable results to be obtained. This strategy has been successfully applied in many studies (Konishi et al., 2006; Zhang D. et al., 2014; Zhang Y. et al., 2016; Zhao et al., 2015; Sun et al., 2016). Because the soybean cultivars in the Jiang-Huai River Valley usually have a relatively low seed oil content (Zhang J. W. et al., 2014), it is difficult to breed high oil content varieties by traditional breeding in this region. Molecular MAS is a possible alternative. However, little is known about the genetics of oil content in these soybean. In this study, we used 104 RILs and 279 accessions to dissect the genetic basis of seed oil content in soybean, which would provide useful information for MAS to breed high oil content varieties in this region.

The marker density should be considered during QTL mapping, regardless of which method is used. Because of the technology limitations, the genetic maps used in most of the studies mentioned above were constructed for mapping QTLs associated with seed oil content in soybean with only hundreds of restriction fragment length polymorphism, amplified fragment length polymorphism, and/or SSR markers. In addition, the interval distance between adjacent markers was large, which limits the efficiency and accuracy of QTL mapping. In our study, SLAF-seq was used to genotype a total of 104 individuals and the two parents. We constructed a genetic map with 2,062 SLAF markers for linkage mapping, and the average distance between adjacent markers was only 1 cM. For GWAS, 59,845 polymorphism SNPs with MAF >0.05 were used. Previous studies have shown that the LD of cultivated soybean is approximately 150 kb (Lam et al., 2010; Hao et al., 2012; Wen et al., 2014). Hwang et al. (2014) reported that the mean LD (r²) rapidly declined to 0.2 within 360 kb in euchromatic regions and declined to 0.2 at 9600 kb in heterochromatic regions. However, there is a high LD in the soybean genome, implying a limitation on the effectiveness of increasing marker density to improve the resolution. This genetic map and high-number of SNPs can meet the demands of linkage and a GWAS of the complex traits in soybean.

Major QTLs for Seed Oil Content in Soybean

Soybean is a primary source of plant oil for humans (Wilcox, 2004). The breeding of high oil content varieties is an important goal of oil crop breeding programs. The QTLs controlling seed oil content have been reported many times in bi-parental populations (Mansur et al., 1996; Brummer et al., 1997; Kabelka et al., 2004; Qi et al., 2011; Wang et al., 2015). Many loci associated with soybean seed oil content have been detected by GWAS methods (Hwang et al., 2014; Bandillo et al., 2015; Wen et al., 2015). However, most QTLs are environmentally sensitive. Brummer et al. (1997) considered that in molecular-assisted breeding programs, breeders should use QTLs that are stable in a range of different environments. Thus, there is a need to identify stable QTLs that can be conveniently used for MAS and map-based cloning to determine the mechanisms of seed oil content.

In this study, the distribution of seed oil content followed a normal distribution in both the natural soybean population and in RILs, indicating that the genetic control of seed oil content is complex (Supplementary Figure S1). Although seed oil content had a high h², the environmental instability of this trait made it difficult to detect stable QTLs in different environments. However, there were some major loci detected in different environments using one method in our study. In linkage mapping, qOil-10-1 and qOil-14-1 were detected in four environments, had adjacent or physically overlapping QTLs, as reported by previous studies, and might actually be the same QTL (Tajuddin et al., 2003; Wang et al., 2014). In the GWAS, the significant SNPs on Chr. 20 were found to be associated over the 2 years investigated. The peak marker, Gm20_31164168, was similar to the one identified by Vaughn et al. (2014) at 31150279 bp on Chr. 20. Therefore, our results were in good agreement with previous studies, and the QTLs identified can be treated as important targets for the identification of candidate genes involved in the modification of oil in future studies.

Combination of Linkage and Association Mapping to Refine the Stable QTL on Chr. 05

Although the number of QTLs controlling seed oil content have been reported by linkage mapping and/or GWASs in soybean, the confirmation of such QTLs is rare. Here, we confirmed and narrowed down the interval of a stable QTL on Chr. 05 by a combination of linkage and association mapping. In the NJMN-RIL population, qOil-5-1 was mapped at an approximate 6.8 cM interval in all environments. Based on the physical position of the SLAF markers, this genetic map interval has a corresponding physical distance of 2.5 Mb, from 38.04 to 40.51 Mb on Chr. 05, which has adjacent or physically overlapping QTLs, as reported by Mansur et al. (1996), Orf et al. (1999), and Pathan et al. (2013). In the secondary F₂ population, the position of this QTL could be refined at a 3 cM interval between the BARCSOYSSR_05_1226 and BARCSOYSSR_05_1250 markers, and the peak marker corresponded to position 38494541–38494952 bp. These results indicated that the region of qOil-5-1 is from ∼37.85 to 38.90 Mb. In the germplasm population, based on the GWAS results and the partial LD of Chr. 05, we found that all the SNPs that had a significant association with seed oil content over the 2 years investigated were located in one LD block. Taken together, this evidence strongly suggests that this stable QTL for seed oil content was located in a small region from 38.13 to 38.57 Mb on Chr. 05.

Plant oil is typically biosynthesized and accumulated in the cotyledon or endosperm tissues of seeds. The pathways for lipid biosynthesis and oil accumulation have been studied and are known to be controlled by multiple genes. Within the LD block region, there were more than 60 functional genes². Thirteen of those genes were highly expressed in seed (information from Phytozome database³) and may be related to oil biosynthesis or oil accumulation (Table 4). However, the functions of these genes need to be verified in future studies. Finally, such fine-mapping of QTLs and close markers (including SNP, SSR and InDel markers) must be used for MAS and map-based cloning to determine the mechanisms of seed oil content.

TABLE 4

TABLE 4. Genes within the linkage disequilibrium (LD) block that are highly expressed in soybean seed.

Conclusion

We identified some major QTLs for seed oil content and dissected the genetic basis of this trait in soybean using linkage and germplasm populations. A stable QTL was confined to a small interval of approximately 440 kb on Chr. 05 by combining a linkage and genome-wide association mapping approach. A genetic map was constructed using more than 2,062 SLAF markers for the NJMN-RIL population in this study, which can provide a good foundation for analyzing quantitative traits. The results of this study are important for future map-based gene cloning and also provide support for the implementation of MAS for breeding soybean with a high seed oil content in the Jiang-Huai River Valley.

Author Contributions

TZ conceived the research; TZ and JG designed the research; YC, SL, ZW, FC, and JK performed the experiments and analyzed the data. YC drafted the manuscript. TZ revised the paper.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by grants from the National Key R & D Program for Crop Breeding (2016YFD0100304), the Natural Science Foundation of China (31271750, 31571691), the MOE Program for Changjiang Scholars and Innovative Research Team in University (PCSIRT13073), and the Jiangsu JCIC-MCP program.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01222/full#supplementary-material

Footnotes

1. ^ http://www.soybase.org/
2. ^ http://www.soybase.org
3. ^ http://phytozome.net

References

Allen, F. L. (1994). “Usefulness of plant genome mapping to plant-breeding,” in Plant Genome Analysis, ed. P. M. Gresshoff (Boca Raton, FL: CRC Press), 11–18.

Google Scholar

Bandillo, N., Jarquin, D., Song, Q. J., Nelson, R., Cregan, P., Specht, J., et al. (2015). A population structure and genome-wide association analysis on the USDA soybean germplasm collection. Plant Genome 8, 1–13. doi: 10.3835/plantgenome2015.04.0024

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2005). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi: 10.1093/bioinformatics/bth457

PubMed Abstract | CrossRef Full Text | Google Scholar

Beavis, W. D. (1994). “The power and deceit of QTL experiments: lessons from comparative QTL studies,” in Proceedings of the Forty-ninth Annual Corn and Sorghum Research Conference, (Washington, DC: American Seed Trade Association), 250–266.

Google Scholar

Bolon, Y. T., Joseph, B., Cannon, S. B., Graham, M. A., Diers, B. W., Farmer, A. D., et al. (2010). Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean. BMC Plant Biol. 10:41. doi: 10.1186/1471-2229-10-41

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | CrossRef Full Text | Google Scholar

Brummer, E. C., Graef, G. L., Orf, J., Wilcox, J. R., and Shoemaker, R. C. (1997). Mapping QTL for seed protein and oil content in eight soybean populations. Crop Sci. 37, 370–378. doi: 10.1007/BF00224058

PubMed Abstract | CrossRef Full Text | Google Scholar

Burton, J. W., and Brim, C. A. (1981). Recurrent selection in soybeans. III. Selection for increased percent oil in seeds. Crop Sci. 21, 31–34. doi: 10.2135/cropsci1981.0011183X002100010009x

CrossRef Full Text | Google Scholar

Cao, Y. C., Li, S. G., He, X. H., Chang, F. G., Kong, J. J., Gai, J. Y., et al. (2017). Mapping QTLs for plant height and flowering time in a Chinese summer planting soybean RIL population. Euphytica 213:39. doi: 10.1007/s10681-016-1834-8

CrossRef Full Text | Google Scholar

Chapman, A., Pantalone, V. R., Ustun, A., Allen, F. L., Landau-Ellis, D., Trigiano, R. N., et al. (2003). Quantitative trait loci for agronomic and seed quality traits in an F2 and F4:6 soybean population. Euphytica 129, 387–393. doi: 10.1023/A:1022282726117

CrossRef Full Text | Google Scholar

Chung, J., Babka, H. L., Graef, G. L., Staswick, P. E., Lee, D. J., Cregan, P. B., et al. (2003). The seed protein, oil, and yield QTL on soybean linkage group I. Crop Sci. 43, 1053–1067. doi: 10.2135/cropsci2003.1053

PubMed Abstract | CrossRef Full Text | Google Scholar

Collard, B. C. Y., and Mackill, D. J. (2008). Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Philos. Trans. R. Soc. B Biol. Sci. 363, 557–572. doi: 10.1098/rstb.2007.2170

PubMed Abstract | CrossRef Full Text | Google Scholar

Contreras-Soto, R. I., Mora, F., de Oliveira, M. A. R., Higashi, W., Scapim, C. A., and Schuster, I. (2017). A genome-wide association study for agronomic traits in soybean using SNP markers and SNP-based haplotype analysis. PLoS ONE 12:e0171105. doi: 10.1371/journal.pone.0171105

PubMed Abstract | CrossRef Full Text | Google Scholar

Csanadi, G., Vollmann, J., Stift, G., and Lelley, T. (2001). Seed quality QTLs identified in a molecular map of early maturing soybean. Theor. Appl. Genet. 103, 912–919. doi: 10.1007/s001220100621

CrossRef Full Text | Google Scholar

Darvasi, A., Weinreb, A., Minke, V., Weller, J. I., and Soller, M. (1993). Detecting marker-Qtl linkage and estimating Qtl gene effect and map location using a saturated genetic-map. Genetics 134, 943–951.

PubMed Abstract | Google Scholar

Diers, B. W., Keim, P., Fehr, W. R., and Shoemaker, R. C. (1992). Rflp analysis of soybean seed protein and oil content. Theor. Appl. Genet. 83, 608–612. doi: 10.1007/BF00226905

PubMed Abstract | CrossRef Full Text | Google Scholar

Doyle, J. J. (1990). Isolation of plant DNA from fresh tissue. Focus 12, 13–15.

Google Scholar

Hao, D. R., Chao, M. N., Yin, Z. T., and Yu, D. Y. (2012). Genome-wide association analysis detecting significant single nucleotide polymorphisms for chlorophyll and chlorophyll fluorescence parameters in soybean (Glycine max) landraces. Euphytica 186, 919–931. doi: 10.1007/s10681-012-0697-x

CrossRef Full Text | Google Scholar

Huang, S. X., Ding, J., Deng, D. J., Tang, W., Sun, H. H., Liu, D. Y., et al. (2013). Draft genome of the kiwifruit Actinidia chinensis. Nat. Commun. 4:2640. doi: 10.1038/Ncomms3640

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, E. Y., Song, Q. J., Jia, G. F., Specht, J. E., Hyten, D. L., Costa, J., et al. (2014). A genome-wide association study of seed protein and oil content in soybean. BMC Genomics 15:1. doi: 10.1186/1471-2164-15-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyten, D. L., Pantalone, V. R., Sams, C. E., Saxton, A. M., Landau-Ellis, D., Stefaniak, T. R., et al. (2004). Seed quality QTL in a prominent soybean population. Theor. Appl. Genet. 109, 552–561. doi: 10.1007/s00122-004-1661-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyten, D. L., Song, Q., Choi, I. Y., Yoon, M. S., Specht, J. E., Matukumalli, L. K., et al. (2008). High-throughput genotyping with the GoldenGate assay in the complex genome of soybean. Theor. Appl. Genet. 116, 945–952. doi: 10.1007/s00122-008-0726-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kabelka, E. A., Diers, B. W., Fehr, W. R., LeRoy, A. R., Baianu, I. C., You, T., et al. (2004). Putative alleles for increased yield from soybean plant introductions. Crop Sci. 44, 784–791. doi: 10.2135/cropsci2004.7840

CrossRef Full Text | Google Scholar

Kim, M., Schultz, S., Nelson, R. L., and Diers, B. W. (2016). Identification and fine mapping of a soybean seed protein QTL from PI 407788A on chromosome 15. Crop Sci. 56, 219–225. doi: 10.2135/cropsci2015.06.0340

CrossRef Full Text | Google Scholar

Konishi, S., Izawa, T., Lin, S. Y., Ebana, K., Fukuta, Y., Sasaki, T., et al. (2006). An SNP caused loss of seed shattering during rice domestication. Science 312, 1392–1396. doi: 10.1126/science.1126410

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosambi, D. D. (1943). The estimation of map distances from recombination values. Ann. Eugen. 12, 172–175. doi: 10.1111/j.1469-1809.1943.tb02321.x

CrossRef Full Text | Google Scholar

Lam, H. M., Xu, X., Liu, X., Chen, W. B., Yang, G. H., Wong, F. L., et al. (2010). Resequencing of 31 wild and cultivated soybean genomes identifies patterns of genetic diversity and selection. Nat. Genet. 42, 1053–1059. doi: 10.1038/ng.715

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. H., Bailey, M. A., Mian, M. A. R., Carter, T. E., Shipe, E. R., Ashley, D. A., et al. (1996). RFLP loci associated with soybean seed protein and oil content across populations and locations. Theor. Appl. Genet. 93, 649–657. doi: 10.1007/Bf00224058

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L. H., Guo, N., Niu, J. P., Wang, Z. L., Cui, X. X., Sun, J. T., et al. (2016). Loci and candidate gene identification for resistance to Phytophthora sojae via association analysis in soybean [Glycine max (L.) Merr.]. Mol. Genet. Genomics 291, 1095–1103. doi: 10.1007/s00438-015-1164-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D. Y., Ma, C. X., Hong, W. G., Huang, L., Liu, M., Liu, H., et al. (2014). Construction and analysis of high-density linkage map using high-throughput sequencing data. PLoS ONE 9:e98855. doi: 10.1371/journal.pone.0098855

PubMed Abstract | CrossRef Full Text | Google Scholar

Mamidi, S., Lee, R. K., Goos, J. R., and McClean, P. E. (2014). Genome-wide association studies identifies seven major regions responsible for iron deficiency chlorosis in soybean (Glycine max). PLoS ONE 9:e107469. doi: 10.1371/journal.pone.0107469

PubMed Abstract | CrossRef Full Text | Google Scholar

Mansur, L. M., Orf, J. H., Chase, K., Jarvik, T., Cregan, P. B., and Lark, K. G. (1996). Genetic mapping of agronomic traits using recombinant inbred lines of soybean. Crop Sci. 36, 1327–1336. doi: 10.2135/cropsci1996.0011183X003600050042x

CrossRef Full Text | Google Scholar

Mao, D. H., Yu, L., Chen, D. Z., Li, L. Y., Zhu, Y. X., Xiao, Y. Q., et al. (2015). Multiple cold resistance loci confer the high cold tolerance adaptation of Dongxiang wild rice (Oryza rufipogon) to its high-latitude habitat. Theor. Appl. Genet. 128, 1359–1371. doi: 10.1007/s00122-015-2511-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Nichols, D. M., Glover, K. D., Carlson, S. R., Specht, J. E., and Diers, B. W. (2006). Fine mapping of a seed protein QTL on soybean linkage group I and its correlated effects on agronomic traits. Crop Sci. 46, 834–839. doi: 10.2135/cropsci205.05-0168

CrossRef Full Text | Google Scholar

Nyquist, W. E., and Baker, R. J. (1991). Estimation of heritability and prediction of selection response in plant populations. Crit. Rev. Plant Sci. 10, 235–322. doi: 10.1080/07352689109382313

CrossRef Full Text | Google Scholar

Orf, J. H., Chase, K., Jarvik, T., Mansur, L. M., Cregan, P. B., Adler, F. R., et al. (1999). Genetics of soybean agronomic traits: I. Comparison of three related recombinant inbred populations. Crop Sci. 39, 1642–1651. doi: 10.2135/cropsci1999.3961642x

CrossRef Full Text | Google Scholar

Panthee, D. R., Pantalone, V. R., West, D. R., Saxton, A. M., and Sams, C. E. (2005). Quantitative trait loci for seed protein and oil concentration, and seed size in soybean. Crop Sci. 45, 2015–2022. doi: 10.2135/cropsci2004.0720

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathan, S. M., Vuong, T., Clark, K., Lee, J. D., Shannon, J. G., Roberts, C. A., et al. (2013). Genetic mapping and confirmation of quantitative trait loci for seed protein and oil contents and seed weight in soybean. Crop Sci. 53, 765–774. doi: 10.2135/cropsci2012.03.0153

CrossRef Full Text | Google Scholar

Qi, Z. M., Wu, Q., Han, X., Sun, Y. N., Du, X. Y., Liu, C. Y., et al. (2011). Soybean oil content QTL mapping and integrating with meta-analysis method for mining genes. Euphytica 179, 499–514. doi: 10.1007/s10681-011-0386-1

CrossRef Full Text | Google Scholar

Rafalski, J. A. (2010). Association genetics in crop improvement. Curr. Opin. Plant Biol. 13, 174–180. doi: 10.1016/j.pbi.2009.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmutz, J., Cannon, S. B., Schlueter, J., Ma, J. X., Mitros, T., Nelson, W., et al. (2010). Genome sequence of the palaeopolyploid soybean. Nature 463, 178–183. doi: 10.1038/nature08670

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Q. J., Jenkins, J., Jia, G. F., Hyten, D. L., Pantalone, V., Jackson, S. A., et al. (2016). Construction of high resolution genetic linkage maps to improve the soybean genome sequence assembly Glyma1.01. BMC Genomics 17:33. doi: 10.1186/s12864-015-2344-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Q. J., Jia, G. F., Zhu, Y. L., Grant, D., Nelson, R. T., Hwang, E. Y., et al. (2010). Abundance of SSR motifs and development of candidate polymorphic SSR markers (BARCSOYSSR_1.0) in soybean. Crop Sci. 50, 1950–1960. doi: 10.2135/cropsci2009.10.0607

CrossRef Full Text | Google Scholar

Soy Stats. (2011). A Reference Guide to Important Soybean Facts and Figures. St. Louis, MO: American Soybean Association.

Google Scholar

Sun, F. M., Liu, J., Hua, W., Sun, X. C., Wang, X. F., and Wang, H. Z. (2016). Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus. Plant Sci. 252, 388–399. doi: 10.1016/j.plantsci.2016.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X. W., Liu, D. Y., Zhang, X. F., Li, W. B., Liu, H., Hong, W. G., et al. (2013). SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE 8:e58700. doi: 10.1371/journal.pone.0058700

PubMed Abstract | CrossRef Full Text | Google Scholar

Tajuddin, T., Watanabe, S., Yamanaka, N., and Harada, K. (2003). Analysis of quantitative trait loci for protein and lipid contents in soybean seeds using recombinant inbred lines. Breed. Sci. 53, 133–140. doi: 10.1270/jsbbs.53.133

CrossRef Full Text | Google Scholar

Vaughn, J. N., Nelson, R. L., Song, Q. J., Cregan, P. B., and Li, Z. L. (2014). The genetic architecture of seed composition in soybean is refined by genome-wide association scans across multiple populations. G3 4, 2283–2294. doi: 10.1534/g3.114.013433

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Chen, P. Y., Wang, D. C., Shannon, G., Shi, A. N., Zeng, A. L., et al. (2015). Identification of quantitative trait loci for oil content in soybean seed. Crop Sci. 55, 23–34. doi: 10.2135/cropsci2014.04.0280

CrossRef Full Text | Google Scholar

Wang, S., Basten, C. J., and Zeng, Z. B. (2007). Windows QTL Cartographer 2.5. Raleigh, NC: North Carolina state university.

Google Scholar

Wang, W. B., Liu, M. F., Wang, Y. F., Li, X. L., Cheng, S. X., Shu, L. P., et al. (2016). Characterizing two inter-specific bin maps for the exploration of the QTLs/Genes that confer three soybean evolutionary traits. Front. Plant Sci. 7:1248. doi: 10.3389/Fpls.2016.01248

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. Z., Jiang, G. L., Green, M., Scott, R. A., Song, Q. J., Hyten, D. L., et al. (2014). Identification and validation of quantitative trait loci for seed yield, oil and protein contents in two recombinant inbred line populations of soybean. Mol. Genet. Genomics 289, 935–949. doi: 10.1007/s00438-014-0865-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, Z. X., Boyse, J. F., Song, Q. J., Cregan, P. B., and Wang, D. C. (2015). Genomic consequences of selection and genome-wide association mapping in soybean. BMC Genomics 16:671. doi: 10.1186/S12864-015-1872-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, Z. X., Tan, R. J., Yuan, J. Z., Bales, C., Du, W. Y., Zhang, S. C., et al. (2014). Genome-wide association mapping of quantitative resistance to sudden death syndrome in soybean. BMC Genomics 15:809. doi: 10.1186/1471-2164-15-809

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilcox, J. R. (1998). Increasing seed protein in soybean with eight cycles of recurrent selection. Crop Sci. 38, 1536–1540. doi: 10.2135/cropsci1998.0011183X003800060021x

CrossRef Full Text | Google Scholar

Wilcox, J. R. (2004). “World distribution and trade of soybean,” in Soybeans: Improvement, production, and Uses, 3rd Edn, eds H. R. Boerma and J. E. Specht (Madison, WI: American Society of Agronomy), 1–14.

Google Scholar

Wilson, R. F. (2004). “Seed composition,” in Soybeans: Improvement, Production, and Uses, 3rd Edn, eds H. R. Boerma and J. E. Specht (Madison, WI: American Society of Agronomy), 621–677.

Google Scholar

Wilson, R. F. (2008). “Soybean: market driven research needs,” in Genetics and Genomics of Soybean, ed. G. Stacey (New York, NY: Springer), 3–15. doi: 10.1007/978-0-387-72299-3_1

CrossRef Full Text | Google Scholar

Yang, N., Lu, Y. L., Yang, X. H., Huang, J., Zhou, Y., Ali, F., et al. (2014). Genome wide association studies using a new nonparametric model reveal the genetic architecture of 17 agronomic traits in an enlarged maize association panel. PLoS Genet. 10:e1004573. doi: 10.1371/journal.pgen.1004573

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Li, H. Y., Wang, J. S., Zhang, H. Y., Hu, Z. B., Chu, S. S., et al. (2016). High-density genetic mapping identifies new major loci for tolerance to low-phosphorus stress in soybean. Front. Plant Sci. 7:372. doi: 10.3389/Fpls.2016.00372

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Song, H. N., Cheng, H., Hao, D. R., Wang, H., Kan, G. Z., et al. (2014). The acid phosphatase-encoding gene GmACP1 contributes to soybean tolerance to low-phosphorus stress. PLoS Genet. 10:e1004061. doi: 10.1371/journal.pgen.1004061

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. P., Song, Q. J., Cregan, P. B., Nelson, R. L., Wang, X. Z., Wu, J. X., et al. (2015). Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm. BMC Genomics 16:217. doi: 10.1186/S12864-015-1441-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. W., Han, F. X., Sun, J. M., Han, G. Z., Yu, S. X., Yu, F. K., et al. (2014). Genetic variation of protein and fat content in soybean mini core collections. J. Plant Genet. Resour. 15, 405–410.

Zhang, Y., Cui, M., Zhang, J. M., Zhang, L., Li, C. L., Kan, X., et al. (2016). Confirmation and fine mapping of a major QTL for aflatoxin resistance in maize using a combination of linkage and association mapping. Toxins 8:258. doi: 10.3390/Toxins8090258

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. X., Wang, L. H., Xin, H. G., Li, D. H., Ma, C. X., Ding, X., et al. (2013). Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC Plant Biol. 13:141. doi: 10.1186/1471-2229-13-141

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, K., Tung, C. W., Eizenga, G. C., Wright, M. H., Ali, M. L., Price, A. H., et al. (2011). Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2:467. doi: 10.1038/Ncomms1467

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X., Han, Y. P., Li, Y. H., Liu, D. Y., Sun, M. M., Zhao, Y., et al. (2015). Loci and candidate gene identification for resistance to Sclerotinia sclerotiorum in soybean (Glycine max L. Merr.) via association and linkage maps. Plant J. 82, 245–255. doi: 10.1111/tpj.12810

PubMed Abstract | CrossRef Full Text | Google Scholar