Identification of the Genomic Region Underlying Seed Weight per Plant in Soybean (Glycine max L. Merr.) via High-Throughput Single-Nucleotide Polymorphisms and a Genome-Wide Association Study

Seed weight per plant (SWPP) of soybean (Glycine max (L.) Merr.), a complicated quantitative trait controlled by multiple genes, was positively associated with soybean seed yields. In the present study, a natural soybean population containing 185 diverse accessions primarily from China was used to analyze the genetic basis of SWPP via genome-wide association analysis (GWAS) based on high-throughput single-nucleotide polymorphisms (SNPs) generated by the Specific Locus Amplified Fragment Sequencing (SLAF-seq) method. A total of 33,149 SNPs were finally identified with minor allele frequencies (MAF) > 5% which were present in 97% of all the genotypes. Twenty association signals associated with SWPP were detected via GWAS. Among these signals, eight SNPs were novel loci, and the other twelve SNPs were overlapped or located in the linked genomic regions of the reported QTL from SoyBase database. Several genes belonging to the categories of hormone pathways, RNA regulation of transcription in plant development, ubiquitin, transporting systems, and other metabolisms were considered as candidate genes associated with SWPP. Furthermore, nine genes from the flanking region of Gm07:19488264, Gm08:15768591, Gm08:15768603, or Gm18:23052511 were significantly associated with SWPP and were stable among multiple environments. Nine out of 18 haplotypes from nine genes showed the effect of increasing SWPP. The identified loci along with the beneficial alleles and candidate genes could be of great value for studying the molecular mechanisms underlying SWPP and for improving the potential seed yield of soybean in the future.

Introduction

Seed weight per plant, a complicated and agronomically important quantitative trait, was significantly related with yield in soybean (Glycine max L. Merr.) (Chen et al., 2007; Liu et al., 2016). SWPP is controlled by multiple genes or quantitative trait loci (QTL). An additive effect dominates the inheritance pattern of SWPP (Chen et al., 2007). The development of cultivars with suitable SWPP has been an important breeding object of many soybean breeders because SWPP was an important soybean yield component. SWPP is influenced by the environment or genotype by environment interactions, and this trait performs differently in different environments. Hence, breeding soybean cultivars with suitable SWPP via traditional methods requires evaluation in multiple environments over several years, which is expensive, time consuming, and labor-intensive.

The advances in molecular marker technologies have enabled the efficient elucidation of the genetic architecture of soybean SWPP. To date, fewer than 50 QTL, located on chromosome (Chr.) 3 (linkage group, LG N), Chr.4 (LG C1), Chr.5 (LG A1), Chr.6 (LG C2), Chr.7 (LG M), Chr.8 (LG A2), Chr.9 (LG K), Chr.10 (LG O), Chr.18 (LG G), and Chr.19 (LG L), have been reported in the SoyBase databank¹. Of these identified QTL, some QTL, located on Chr. 3 (LG N), Chr.4 (LG C1), and Chr.6 (LG C2), were verified by many studies (Chen et al., 2007; Kuroda et al., 2013; Yao et al., 2015). The genetic maps used in most of these studies had incomplete coverage of the soybean genome, with large gaps in the maps. Thus, some QTL were difficult to directly apply for marker assisted selection (MAS) of SWPP.

Genome-wide association studies (GWAS), based on high-density markers and natural populations, have been regarded as an efficient alternative to linkage analysis with more extensive recombination events and shorter LD segments. Thus, the resolution and accuracy of marker-phenotype associations could be further increased via GWAS compared with the conventional QTL mapping of segregating populations (Li et al., 2015). To date, GWAS has been widely utilized to elucidate the genetic basis of many complex traits in some crops, including rice (Huang et al., 2010), maize (Weng et al., 2011), wheat (Raman et al., 2010), and barley (Cockram et al., 2010). In soybean, the genetic architecture of some important traits, such as protein (Hwang et al., 2014), fatty acid content (Li et al., 2015), and SCN resistance (Han et al., 2015; Zhao et al., 2017), has been well dissected. Moreover, the rapid development of next-generation sequencing technology and SNP genotyping technology has propelled much of the practicability of GWAS. In previous studies, Yan et al. (2017), Contreras-Soto et al. (2017), and Zhang et al. (2016) identified seventeen, seven, and forty-eight SNPs, respectively, which were all significantly associated with 100-seed weight of soybean (HSW). However, to date, only few studies were conducted to identify QTL underlying SWPP based on high-throughput sequencing technology. Han et al. (2016) detected rs18976374 (located on Chr.16) significantly related with SWPP by using a bar coded multiplex sequencing approach with an Illumina Genome Analyzer II. Liu et al. (2016) found one SNP (ss244932137 located on Chr.3) associated with SWPP through GWAS strategy based on an association panel of 138 cultivars genotyped by Illumina SoySNP6KiSelect BeadChip.

The aims of the present study were to understand better the genetic architecture of SWPP via GWAS based on 185 tested accessions and 33,149 SNPs and to analyze the potential candidate genes that might regulate soybean SWPP in associated genomic regions with peak SNPs.

Materials and Methods

Soybean Germplasms and Field Trials

A total of 185 diverse soybean accessions (Supplementary Table S1), including landraces and elite cultivars, were applied to evaluate the variation of SWPP and for the subsequent sequencing analysis. All samples were grown in Harbin in 2015 and 2016, Gongzhuling in 2015 and 2016, and Shenyang in 2015 and 2016. Field experiments were performed with single row plots of 3 m in length with 0.65 m between rows. A randomized complete block design was used with three replications in each tested environment. After reaching full maturity, a total of 10 randomly selected plants per row in each plot were weighed and used to evaluate SWPP.

DNA Isolation and Sequencing Analysis

Total DNA from the fresh leaves of a single test sample was extracted by the CTAB method according to Han et al. (2015). The isolated high-quality DNA was partly sequenced via specific locus amplified fragment sequencing (SLAF-seq) methodology (Sun et al., 2013). Soybean reference genome (Version:Glyma.Wm82.a2) was preliminarily analyzed and digested enzyme Mse I (EC 3.1.21.4) and Hae III (EC: 3.1.21.4) (Thermo Fisher Scientific Inc., Waltham, MA, United States) were used to generate more than 50,000 sequencing tags (approximately 300–500 bp in length) of all tested samples. The obtained tags were evenly distributed among the unique genomic regions of the 20 soybean chromosomes. The sequencing libraries of each tested accession were built based on the sequencing tags. A barcode method combined with the Illumina Genome Analyzer II system (Illumina Inc., San Diego, CA, United States) was used to generate the 45-bp sequence reads at both ends of the sequencing tags from each accession library. Short Oligonucleotide Alignment Program 2 (SOAP2) software was used to align raw paired-end reads to the soybean reference genome (Version:Glyma.Wm82.a2). The SLAF groups were designed based on the raw reads, which mapped to the same unique genomic positions. Approximately 58,000 high-quality SLAF tags were acquired from each tested accession. The SNPs were called as such based on MAF ≥ 0.05. If the minor allele depth or the total depth of the sample was larger than 1/3, then the genotype was considered heterozygous.

For thirty lines, a genome resequencing with 10-fold in depth was conducted on an Illumina HiSeq 2000 sequencer. Paired-end resequencing reads were mapped to the soybean Williams 82 reference genome (Version:Glyma.Wm82.a2) with BWA (Version: 0.6.1-r104) (Zhou et al., 2015) using the default parameters. SAMtools48 (Version: 0.1.18) software (Zhou et al., 2015) was used for converting mapping results into the BAM format and to filter the unmapped and non-unique reads. Duplicated reads were filtered with the Picard package (picard.sourceforge.net, Version:1.87) (Zhou et al., 2015). The BEDtools (Version: 2.17.0) (Zhou et al., 2015) coverageBed program was applied to compute the coverage of sequence alignments. A sequence was defined as absent when the coverage was lower than 90% and present when the coverage was higher than 90%. SNP detection was performed by the Genome Analysis Toolkit (GATK, version 2.4-7-g5e89f01) and SAMtools (Zhou et al., 2015). Only the SNPs detected by both methods could be analyzed further. SNPs with allele frequencies lower than 1% in the population were discarded. SNP annotation was performed based on the soybean genome (Version:Glyma.Wm82.a2) using the package ANNOVAR (Version: 2013-08-23) (Zhou et al., 2015).

Population Structure Evaluation and Linkage Disequilibrium (LD) Analysis

The population structure analysis of the natural group was conducted based on the PCA programs in the GAPIT software package (Lipka et al., 2012). The LD block was evaluated across the soybean genome based on SNPs with MAF ≥ 0.05 and missing data ≤ 10% by using squared allele frequency correlations (r²) in TASSEL version 3.0 (Bradbury et al., 2007). In contrast to the GWAS, missing SNP genotypes were not imputed with the major allele prior to LD analysis. The parameters in the software programs were set according to the MAF (≥0.05) and the integrity of each SNP (≥80%).

Genome-Wide Association Analysis

The SWPP association signals were identified based on 33,149 SNPs (Supplementary Table S2) from 185 tested samples with CMLM in GAPIT (Lipka et al., 2012). The P value was calculated using the Bonferroni method with α ≤ 0.05 (≤2.58 × 10^-4) and was used as the threshold to declare whether a significant association signal existed (Holm, 1979).

Prediction of Candidate Genes Controlling SWPP

According to the studies of Hwang et al. (2014), Han et al. (2015), and Zhao et al. (2017), the average LD decay distance of soybean genome (Version:Glyma.Wm82.a2) was approximately 200 kbp. Thus, candidate genes located in the 200-kbp genomic region of each peak SNP were classified and annotated underlying the soybean reference genome Williams 82².

Haplotype Analysis of Candidate Genes

Based on the genome annotation, SNPs were classified in exonic regions (overlapping with a coding exon), splicing sites (within 2 bp of a splicing junction), 5′UTRs and 3′UTRs, intronic regions (overlapping with an intron), upstream and downstream regions (within a 1 kb region upstream or downstream from the transcription start site), and intergenic regions. SNPs in coding exons were further grouped into synonymous SNPs (did not cause amino acid changes) and nonsynonymous SNPs (caused amino acid changes). The variation happened in these regions (except for intergenic regions) of candidate genes in thirty lines generated from genome re-sequencing data which were analyzed using the General Linear Model (GLM) method in TASSEL version 3.0 (Bradbury et al., 2007) to identify related SNPs and haplotypes. Significant SNPs affecting the SWPP were claimed when the test statistics reached P < 0.01.

Results

Statistical and Variation Analysis of SWPP

The SWPP of the 185 tested soybean accessions, grown in multiple locations over years, was determined. The mean results of the SWPP analysis showed that the effects of genotype, environment, and genotype by environment interactions were significant. The skewness and kurtosis of SWPP across the six environments were both less than 1, indicating continuous variation and a near normal distribution (Figure 1). Therefore, the SWPP in the present study was suitable for the subsequent GWAS.

FIGURE 1

FIGURE 1. Distribution of seed weight per plant (SWPP) among 185 soybean accessions.

Specific-Locus Amplified Fragment Sequencing (SLAF-seq) and Genotyping

The selective population contained 185 diverse accessions, primarily collected from China. The genomic DNA from all tested accessions was extracted and partially sequenced based on the SLAF-seq approach. For each tested sample, a mean of 49,571 high-quality tags (or SLAFs) was scanned from 153 million paired-end reads with a 45-bp read length and 6.51-fold sequencing depth. A total of 33,149 SNPs with MAF ≥ 0.05 were generated from the high-quality tags and subsequently used to perform GWAS for SWPP. The obtained SNPs were evenly distributed among the 20 soybean chromosomes, resulting in a marker density of 28.7 kbp per SNP. Chr. 6 and Chr. 11 included the most and least numbers of SNPs, which was 3,556 and 638, respectively (Figure 2).

FIGURE 2

FIGURE 2. Distribution of SNP markers across 20 soybean chromosomes.

Extent of Linkage Disequilibrium (LD) and Analysis of Population Structure

The average distance of LD decay was analyzed to characterize the mapping resolution for GWAS, and LD decayed differently among all the 20 chromosomes. Accordingly, the mean LD decay of the panel was evaluated at 214 kbp when r² dropped to half of its maximum value (Figure 3A). To scan the population stratification of the association panel, principal component analysis and kinship analysis were conducted based on all 33,149 SNP markers. The first three PCs explained 13.83% of the genetic variation. A drastic decline in the genetic variation appeared at PC2 (Figure 3B). However, analysis of the variation of the first 10 PCs revealed an inflection point at PC3 (Figure 3C). Thus, these results suggested that mainly the first three PCs dominated the population structure. Additionally, the heatmap of kinship matrix revealed low levels of genetic relatedness among the 185 tested samples (Figure 3D).

FIGURE 3

FIGURE 3. Linkage disequilibrium (LD), principal component, and kinship analyses of soybean genetic data. (A) LD decay of the genome-wide association study (GWAS) population. (B) The first three principal components of more than 30,000 SNPs used in the GWAS. (C) Population structure of soybean germplasm collection reflected by principal components. (D) A heatmap of the kinship matrix of the 185 soybean accessions calculated from the same SNPs.

Quantitative Trait Nucleotide (QTN) Associated With SWPP Evaluated by GWAS

A total of 20 association signals, distributed on 12 of the 20 soybean chromosomes, were detected by CMLM in the present study (Figure 4 and Table 1). Among these signals, only one QTN (Gm18:23052511 located on Chr.18) was identified in the linked region of a known SWPP QTL. However, as a specific member of seed weight, another 11 SNPs were overlapped or located in the linked region of a known seed weight QTL from the SoyBase databank, particularly for 100-seed weight (Table 2). The remaining 8 association signals were regarded as novel loci, which were first reported for SWPP in the present study (Table 2). The average SWPP for all tested samples with two different alleles were compared (Table 1), and the results demonstrated that the SWPP among these samples were so different that the appropriate alleles might be effectively applied in the marker assisted selection (MAS) of soybean cultivars with suitable SWPP.

FIGURE 4

FIGURE 4. Manhattan plot of association mapping of the SWPP in soybean.

TABLE 1

TABLE 1. Peak SNP associated with SWPP and the evaluation of beneficial alleles.

TABLE 2

TABLE 2. Significant SNPs and candidate genes associated with SWPP.

Prediction of Candidate Genes Controlling SWPP

The genes located in the 200-kbp genomic region of each peak SNP of the identified loci were considered as candidate genes, consistent with a mean LD decay distance of 214 kbp for the entire association panel. Approximately 126 candidate genes were identified. Among these genes, 26 genes had no functional annotation, and one gene had unknown function domains. The remaining 99 genes were classified into different groups by MAPMAN (Thimm et al., 2004) to determine clearly the potential functions. A total of 16 categories were scanned, including co-factor and vitamin metabolism, misc, RNA regulation of transcription, DNA synthesis/chromatin structure, protein synthesis/modification/degradation, signaling, development, transport, amino acid metabolism, hormone metabolism, abiotic stress, cell wall, major CHO metabolism, nucleotide metabolism, other groups, and unassigned genes (Supplementary Figure S1).

The factors that influenced SWPP were almost the same as those that influence seed weight, and seed size was the main component for determining seed weight. Thus, the factors that could regulate the mechanism of seed size might be important elements in directly or indirectly adjusting SWPP (Wang et al., 2017). In some plants, several genes controlling seed size/weight have been identified, including the genes associated with hormone metabolism, such as auxin and cytokinin (Schruff et al., 2006; Li et al., 2013), various transcription factors associated with RNA regulation of transcription in plant development and maturation, such as TTG2, AP2, MINI3, and C2H2 (Garcia et al., 2005; Ohto et al., 2009; Costa et al., 2010; Yin et al., 2010; Cao et al., 2016), the genes associated with protein modification or degradation (particularly ubiquitin ligase genes) (Li et al., 2008; Xia et al., 2013; Wang et al., 2017), and regulatory genes dominating transport systems and other metabolic processes associated with plant growth and development, such as ABC transporters and amino acid metabolism (Wu et al., 2007; Kim et al., 2010; Less et al., 2010; Hildebrandt et al., 2015). Among these genes detected by GWAS in the present study, those associated with the RNA regulation of transcription, including different transcription factor families, such as AP2 and C2H2, were scanned and identified as the functional genes for SWPP, including Glyma.02G159600, Glyma.02G173900, Glyma.07G157600, Glyma.07G157800, Glyma.07G157900, Glyma.08G195300, Glyma.08G195500, Glyma.08G195600, Glyma.08G195700, Glyma.08G195800, Glyma.17G127400, Glyma.17G127700, Glyma.19G057700, and Glyma.19G067900 located on Chr.2, Chr.7, Chr.8, Chr.17, and Chr.19. Two brassinosteroid response factor genes (Glyma.10G129700 and Glyma.17G127100 with distances of 34.31 and 39.83 kbp from peak SNPs Gm10:34747895 and Gm17:10196755, respectively) were identified as candidate genes. Similarly, Glyma.10G130000, an auxin response factor gene located 44.84 kbp from SNP Gm10:34946492 on Chr.10, Glyma.15G240600, an ethylene factor gene located 31.54 kbp from SNP Gm15:45663581 on Chr.15, and Glyma.17G127500, a gibberellin factor gene located 15.58 kbp from SNP Gm17:10196755 on Chr.17, were all considered as the genes controlling SWPP. Three E3 ubiquitin ligase genes (Glyma.08G195200, Glyma.08G195400, and Glyma.08G195900 with distances of 40.42, 17.31, and 15.42 kbp, respectively, from the peak SNP Gm08:15768591 located on Chr.8), which are associated with protein modification or degradation, might also affect SWPP. An additional 10 genes belonging to transport systems (such as ABC transporters) and other main metabolic processes (such as amino acid metabolism, N-metabolism, and vitamin metabolism), including Glyma.02G159900, Glyma.08G196200, Glyma.12G163700, Glyma.12G163800, Glyma.14G207200, Glyma.14G207300, Glyma.14G207700, Glyma.14G207800, Glyma.17G127800, and Glyma.18G143700, were also selected and might also contribute to SWPP.

To predict the possible roles of candidate genes associated with SWPP, haplotype analysis of the 99 genes was performed. A total of 578 SNPs in 99 candidate genes were found among the thirty soybean lines (MAF > 0.1) through genome re-sequencing. Finally, 44 SNPs from nine genes were significantly associated with SWPP among multiple environments (Figure 5 and Supplementary Table S3). Glyma.07G157900, with only one SNP, was not shown in the figure (Figure 5). Two haplotypes were identified for each of the nine genes and the SWPP between each pair of haplotypes showed significant or highly significant difference (Figure 6). Glyma.07G157900, from Gm07:19488264, was detected under the environments of Shenyang in 2016 and the average at Shenyang in 2015 and 2016, which was located in the genomic region of the known loci, “seed weight 12-4” (Csanadi et al., 2001). A total of seven genes from Gm08:15768591 and Gm08:15768603 were detected, which overlapped the region of the known loci, “seed weight 22-1” (Zhang et al., 2004) and “seed weight 35-1” (Han et al., 2012). Of the seven genes, Glyma.08G195200 and Glyma.08G195900 were detected under the environments of Harbin in 2016 and the average at Harbin in 2015 and 2016. Glyma.08G195300 and Glyma.08G195400 were screened under Shenyang in 2015, Harbin in 2016, and the average at Harbin in 2015 and 2016. Glyma.08G195500 and Glyma.08G195700 were detected under Harbin in 2015, Harbin in 2016, and the average at Harbin in 2015 and 2016. Glyma.08G195600 was scanned under Shenyang in 2015, the average at Shenyang in 2015 and 2016, Harbin in 2016, and the average at Harbin in 2015 and 2016. Glyma.18G143700 was detected from Gm18:23052511 that overlapped the region of known loci, “SWPP 6-6” (Yao et al., 2015) and was detected under the condition of Shenyang in 2015, Shenyang in 2016, and the average at Shenyang in 2015 and 2016. These genes and beneficial haplotypes might be valuable for MAS in regulating SWPP of soybean.

FIGURE 5

FIGURE 5. Candidate gene-based association. Gene-based association analysis of candidate genes with SNPs that were significantly correlated to SWPP. HB, Harbin; SY, Shenyang; Ave, Average seed weight per plant in 2015 and 2016.

FIGURE 6

FIGURE 6. Haplotypes analysis of genes that related to SWPP. The ^∗ and ^∗∗ suggested significance of ANOVA at p < 0.05 and p < 0.01, respectively. HB, Harbin; GZ, Gongzhuling; SY, Shenyang; Ave, Average seed weight per plant in 2015 and 2016.

Discussion

Seed weight per plant, controlled by multiple genes, was an important component of seed yield in soybean. Thus far, some QTL have been identified based on linkage analysis from the SoyBase databank. However, few studies dissecting the genetic basis of SWPP in soybean via GWAS in combination with high-throughput SNPs and diverse accessions across multiple environments have been discussed, and even fewer candidate genes have been reported. In the present study, 185 soybean accessions were widely collected from China and used to conduct GWAS via high-throughput SNPs and diverse environments. In total, twenty SNPs were identified among twelve different soybean chromosomes, and these SNPs might have value for further breeding cultivars with appropriate SWPP.

In addition to the QTL associated with SWPP in the SoyBase databank, the genes associated with seed weight were referred for more accurately identifying the genomic region underlying the SWPP of soybean via the twenty SNPs used in the present study. Among the twenty SNPs, eight loci, including Gm02:27731352 located on Chr.2, Gm03:25789435 and Gm03:26388034 located on Chr.3, Gm12:24101071 located on Chr.12, Gm14:47326193 located on Chr.14, Gm15:45663581 located on Chr.15, and Gm19:19848537 and Gm19:19848932 located on Chr.19, were novel genes first reported in the present study. The genomic region of Chr.19, which included two close loci (Gm19:19848537 and Gm19:19848932), might be a primary region associated with SWPP. Additionally, another twelve SNPs were overlapped or located near known QTL (Table 2). Among these SNPs, one SNP (Gm18:23052511 located on Chr.18) was identified in the same region reported by Mian et al. (1996), and this polymorphism was also located near the locus named SWPP 6-6, which had been reported by Yao et al. (2015). A set of seven SNPs, including Gm02:17792241 on Chr.2, Gm08:15768591 and Gm08:15768603 on Chr.8, Gm10:34747895 and Gm10:34946492 on Chr.10, Gm12:31782089 on Chr.12, and Gm20:29716176 on Chr.20, were located in the same genomic region previously reported by Han et al. (2012) by using three RIL populations, and furthermore, among these seven SNPs, two close loci (Gm08:15768591 and Gm08:15768603 on Chr.8), Gm12:31782089 and Gm20:29716176 were also stably identified near previously reported loci (Sebolt et al., 2000; Zhang et al., 2004; Kato et al., 2014). Additionally, four genomic regions (Gm07:19488264 on Chr.7, Gm14:24306475 on Chr.14, Gm17:10196755 on Chr.17, and Gm19:11082560 on Chr.19) were also identified in previous studies (Csanadi et al., 2001; Hoeck et al., 2003; Teng et al., 2009). Moreover, some SNPs from candidate genes based on gene-based association analysis were found close to the loci that were verified in the present study by GWAS and previous studies. The main loci were Gm07:19488264, Gm08:15768591, Gm08:15768603, and Gm18:23052511. These four loci might be major genomic regions containing candidate genes associated with SWPP.

Currently, GWAS had been an effective method to acquire and confirm genes with a suitable LD block (Zhao et al., 2017). In the present study, a total of 99 candidate genes were identified in 200-kbp genomic regions based on the twenty association signals by GWAS with an LD block of approximately 214 kbp in length. Among these genes, those involved in hormone pathways, RNA regulation of transcription in plant development and maturation (TTG2, AP2, MINI3, and C2H2, for example), ubiquitin, transporting systems, and other metabolic processes associated with plant growth (ABC transporters and amino acid metabolism, for example), were considered as the important factors in determining the regulation of SWPP. Therefore, fourteen genes (Glyma.02G159600, Glyma.02G173900, Glyma.07G157600, Glyma.07G157800, Glyma.07G157900, Glyma.08G195300, Glyma.08G195500, Glyma.08G195600, Glyma.08G195700, Glyma.08G195800, Glyma.17G127400, Glyma.17G127700, Glyma.19G057700, and Glyma.19G067900) mainly belonging to transcription factor families of AP2 and C2H2 were proposed as responsible for the SWPP of soybean. Another five novel genes associated with the pathways of brassinosteroid, auxin, ethylene, and gibberellin were also regarded as the candidate genes, including Glyma.10G129700, Glyma.10G130000, Glyma.15G240600, Glyma.17G127100, and Glyma.17G127500. As E3 ubiquitin ligase genes catalyze the ubiquitination of proteins in regulating the growth of plants associated with SWPP (Xia et al., 2013; Yao et al., 2017), Glyma.08G195200, Glyma.08G195400, and Glyma.08G195900, belonging to the E3 ubiquitin family, might act as important factors in controlling SWPP. Additionally, in Arabidopsis, the regulators of amino acid metabolism and ABC transporters, which participate in auxin transport, played a key role in regulating plant growth and seed maturation (Wu et al., 2007; Less et al., 2010). Thus, Glyma.02G159900, Glyma.08G196200, Glyma.12G163700, Glyma.12G163800, Glyma.14G207200, Glyma.14G207300, Glyma.14G207700, Glyma.14G207800, Glyma.17G127800, and Glyma.18G143700, which belong to transport systems (ABC transporters, for example) and other main metabolic processes (amino acid metabolism, for example), might be novel genes associated with SWPP. To further accurately detect the candidate genes controlling SWPP, haplotype analysis of candidate genes was performed. As a result, nine genes and 15 beneficial haplotypes were detected by gene-based association analysis. Definitive function of all candidate genes would be discussed and verified in further studies.

Author Contributions

YJ and XZ conceived the study and contributed to population development. JW contributed to phenotypic evaluation. WT and LQ contributed to genotyping. YH and WL contributed to the experimental design and drafting the manuscript. All authors contributed to and approved the final manuscript.

Funding

This study was financially supported by the National Key R&D Project for Crop Breeding (2016YFD0100304), the National Supporting Project (2014BAD22B01), the Youth Leading Talent Project of the Ministry of Science and Technology in China (2015RA228), the Chinese National Natural Science Foundation (31471517 and 31671717), the Heilongjiang Provincial Project (GX17B002), the ‘Academic Backbone’ Project of Northeast Agricultural University (15XG04), the “Young Talents” Project of Northeast Agricultural University (14QC27), and the Postdoctoral Fund in Heilongjiang Province (LBH-Z15017 and LBH-Q17015).

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 study was conducted in the Key Laboratory of Soybean Biology of the Chinese Education Ministry, Soybean Research and Development Center (CARS), and the Key Laboratory of Northeastern Soybean Biology and Breeding/Genetics of the Chinese Agriculture Ministry.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01392/full#supplementary-material

FIGURE S1 | Functional classifications of the candidate genes of seed weight per plant in soybean.

TABLE S1 | List of 185 soybean accessions.

TABLE S2 | Information of 33,149 SNPs.

TABLE S3 | Haplotype analysis of candidate genes.

Abbreviations

CMLM, compressed mixed linear model; GWAS, genome-wide association study; LD, linkage disequilibrium; MAF, minor allele frequency; PCA, principle component analysis; QTL, quantitative trait locus; SLAF-seq, specific-locus amplified fragment sequencing; SNPs, single-nucleotide polymorphisms; SWPP, seed weight per plant.

Footnotes

1. ^www.soybase.org
2. ^http://www.soybase.org

References

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

Cao, H., Huang, P., Zhang, L., Shi, Y., Sun, D., Yan, Y., et al. (2016). Characterization of 47 Cys2 -His2 zinc finger proteins required for the development and pathogenicity of the rice blast fungus Magnaporthe oryzae. New Phytol. 211, 1035–1051. doi: 10.1111/nph.13948

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Q., Zhang, Z., Liu, C., Xin, D., Qiu, H., Shan, D., et al. (2007). QTL analysis of major agronomic traits in soybean. Agricult. Sci. China 6, 399–405. doi: 10.1016/S1671-2927(07)60062-5

CrossRef Full Text | Google Scholar

Cockram, J., White, J., Zuluaga, D. L., Smith, D., Comadran, J., Macaulay, M., et al. (2010). Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome. Proc. Natl. Acad. Sci. U.S.A. 107, 21611–21616. doi: 10.1073/pnas.1010179107

PubMed Abstract | CrossRef Full Text | Google Scholar

Contreras-Soto, R. I., Mora, F., de Oliveira, M. A., 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

Costa, J. H., Mota, E. F., Cambursano, M. V., Lauxmann, M. A., de Oliveira, L. M., Silva Lima Mda, G., et al. (2010). Stress-induced co-expression of two alternative oxidase (VuAox1 and 2b) genes in Vigna unguiculata. J. Plant Physiol. 167, 561–570. doi: 10.1016/j.jplph.2009.11.001

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

Garcia, D., Fitz Gerald, J. N., and Berger, F. (2005). Maternal control of integument cell elongation and zygotic control of endosperm growth are coordinated to determine seed size in Arabidopsis. Plant Cell 17, 52–60. doi: 10.1105/tpc.104.027136

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Li, D., Zhu, D., Li, H., Li, X., Teng, W., et al. (2012). QTL analysis of soybean seed weight across multi-genetic backgrounds and environments. Theor. Appl. Genet. 125, 671–683. doi: 10.1007/s00122-012-1859-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Zhao, X., Cao, G., Wang, Y., Li, Y., Liu, D., et al. (2015). Genetic characteristics of soybean resistance to HG type 0 and HG type 1.2.3.5.7 of the cyst nematode analyzed by genome-wide association mapping. BMC Genomics 16:598. doi: 10.1186/s12864-015-1800-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Zhao, X., Liu, D., Li, Y., Lightfoot, D. A., Yang, Z., et al. (2016). Domestication footprints anchor genomic regions of agronomic importance in soybeans. New Phytol. 209, 871–884. doi: 10.1111/nph.13626

PubMed Abstract | CrossRef Full Text | Google Scholar

Hildebrandt, T. M., Nunes Nesi, A., Araujo, W. L., and Braun, H. P. (2015). Amino acid catabolism in plants. Mol. Plant 8, 1563–1579. doi: 10.1016/j.molp.2015.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoeck, J. A., Fehr, W. R., Shoemaker, R. C., Welke, G. A., Johnson, S. L., and Cianzio, S. R. (2003). Molecular Marker Analysis of Seed Size in Soybean. Crop Sci. 43:68. doi: 10.2135/cropsci2003.0068

CrossRef Full Text | Google Scholar

Holm, S. (1979). A simple sequentially rejective multiple test procedure. Scand. J. Stat. 6, 65–70.

Google Scholar

Huang, X., Wei, X., Sang, T., Zhao, Q., Feng, Q., Zhao, Y., et al. (2010). Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 42, 961–967. doi: 10.1038/ng.695

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, E. Y., Song, Q., Jia, G., 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-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, S., Sayama, T., Fujii, K., Yumoto, S., Kono, Y., Hwang, T. Y., et al. (2014). A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds. Theor. Appl. Genet. 127, 1365–1374. doi: 10.1007/s00122-014-2304-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. Y., Henrichs, S., Bailly, A., Vincenzetti, V., Sovero, V., Mancuso, S., et al. (2010). Identification of an ABCB/P-glycoprotein-specific inhibitor of auxin transport by chemical genomics. J. Biol. Chem. 285, 23309–23317. doi: 10.1074/jbc.M110.105981

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuroda, Y., Kaga, A., Tomooka, N., Yano, H., Takada, Y., Kato, S., et al. (2013). QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields. Ecol. Evol. 3, 2150–2168. doi: 10.1002/ece3.606

PubMed Abstract | CrossRef Full Text | Google Scholar

Less, H., Angelovici, R., Tzin, V., and Galili, G. (2010). Principal transcriptional regulation and genome-wide system interactions of the Asp-family and aromatic amino acid networks of amino acid metabolism in plants. Amino Acids 39, 1023–1028. doi: 10.1007/s00726-010-0566-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Nie, X., Tan, J. L., and Berger, F. (2013). Integration of epigenetic and genetic controls of seed size by cytokinin in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 110, 15479–15484. doi: 10.1073/pnas.1305175110

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Zheng, L., Corke, F., Smith, C., and Bevan, M. W. (2008). Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana. Genes Dev. 22, 1331–1336. doi: 10.1101/gad.463608

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y. H., Reif, J. C., Ma, Y. S., Hong, H. L., Liu, Z. X., Chang, R. Z., et al. (2015). Targeted association mapping demonstrating the complex molecular genetics of fatty acid formation in soybean. BMC Genomics 16:841. doi: 10.1186/s12864-015-2049-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipka, A. E., Tian, F., Wang, Q., Peiffer, J., Li, M., Bradbury, P. J., et al. (2012). GAPIT: genome association and prediction integrated tool. Bioinformatics 28, 2397–2399. doi: 10.1093/bioinformatics/bts444

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Li, H., Fan, X., Huang, W., Yang, J., Wen, Z., et al. (2016). Phenotypic characterization and genetic dissection of nine agronomic traits in Tokachi nagaha and its derived cultivars in soybean (Glycine max (L.) Merr.). Plant Sci. 256, 72–86. doi: 10.1016/j.plantsci.2016.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Mian, M. A., Bailey, M. A., Tamulonis, J. P., Shipe, E. R., Cater, T. E., Parrot, W. A., et al. (1996). Molecular markers associated with seed weight in two soybean populations. Theor. Appl. Genet. 93, 1011–1016. doi: 10.1007/BF00230118

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohto, M. A., Floyd, S. K., Fischer, R. L., Goldberg, R. B., and Harada, J. J. (2009). Effects of APETALA2 on embryo, endosperm, and seed coat development determine seed size in Arabidopsis. Sex. Plant Reprod. 22, 277–289. doi: 10.1007/s00497-009-0116-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Raman, H., Stodart, B., Ryan, P. R., Delhaize, E., Emebiri, L., Raman, R., et al. (2010). Genome-wide association analyses of common wheat (Triticum aestivum L.) germplasm identifies multiple loci for aluminium resistance. Genome 53, 957–966. doi: 10.1139/G10-058

PubMed Abstract | CrossRef Full Text | Google Scholar

Schruff, M. C., Spielman, M., Tiwari, S., Adams, S., Fenby, N., and Scott, R. J. (2006). The AUXIN RESPONSE FACTOR 2 gene of Arabidopsis links auxin signalling, cell division, and the size of seeds and other organs. Development 133, 251–261. doi: 10.1242/dev.02194

PubMed Abstract | CrossRef Full Text | Google Scholar

Sebolt, A. M., Shoemaker, R. C., and Diers, B. W. (2000). Analysis of a quantitative trait locus allele from wild soybean that increases seed protein concentration in soybean. Crop Sci. 40, 1438–1444. doi: 10.2135/cropsci2000.4051438x

CrossRef Full Text | Google Scholar

Sun, X., Liu, D., Zhang, X., Li, W., Liu, H., Hong, W., 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

Teng, W., Han, Y., Du, Y., Sun, D., Zhang, Z., Qiu, L., et al. (2009). QTL analyses of seed weight during the development of soybean (Glycine max L. Merr.). Heredity 102, 372–380. doi: 10.1038/hdy.2008.108

PubMed Abstract | CrossRef Full Text | Google Scholar

Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313x.2004.02016.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. L., Tang, M. Q., Chen, S., Zheng, X. F., Mo, H. X., Li, S. J., et al. (2017). Down-regulation of BnDA1, whose gene locus is associated with the seeds weight, improves the seeds weight and organ size in Brassica napus. Plant Biotechnol. J. 15, 1024–1033. doi: 10.1111/pbi.12696

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, J., Xie, C., Hao, Z., Wang, J., Liu, C., Li, M., et al. (2011). Genome-wide association study identifies candidate genes that affect plant height in Chinese elite maize (Zea mays L.) inbred lines. PLoS One 6:e29229. doi: 10.1371/journal.pone.0029229

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G., Lewis, D. R., and Spalding, E. P. (2007). Mutations in Arabidopsis multidrug resistance-like ABC transporters separate the roles of acropetal and basipetal auxin transport in lateral root development. Plant Cell 19, 1826–1837. doi: 10.1105/tpc.106.048777

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, T., Li, N., Dumenil, J., Li, J., Kamenski, A., Bevan, M. W., et al. (2013). The ubiquitin receptor DA1 interacts with the E3 ubiquitin ligase DA2 to regulate seed and organ size in Arabidopsis. Plant Cell 25, 3347–3359. doi: 10.1105/tpc.113.115063

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, L., Hofmann, N., Li, S., Ferreira, M. E., Song, B., Jiang, G., et al. (2017). Identification of QTL with large effect on seed weight in a selective population of soybean with genome-wide association and fixation index analyses. BMC Genomics 18:529. doi: 10.1186/s12864-017-3922-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, D., Liu, Z. Z., Zhang, J., Liu, S. Y., Qu, J., Guan, S. Y., et al. (2015). Analysis of quantitative trait loci for main plant traits in soybean. Genet. Mol. Res. 14, 6101–6109. doi: 10.4238/2015.June.8.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, W., Wang, L., Wang, J., Ma, F., Yang, Y., Wang, C., et al. (2017). VpPUB24, a novel gene from Chinese grapevine, Vitis pseudoreticulata, targets VpICE1 to enhance cold tolerance. J. Exp. Bot. 68, 2933–2949. doi: 10.1093/jxb/erx136

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, Z., Meng, F., Song, H., Wang, X., Xu, X., and Yu, D. (2010). Expression quantitative trait loci analysis of two genes encoding rubisco activase in soybean. Plant Physiol. 152, 1625–1637. doi: 10.1104/pp.109.148312

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Song, Q., Cregan, P. B., and Jiang, G. L. (2016). Genome-wide association study, genomic prediction and marker-assisted selection for seed weight in soybean (Glycine max). Theor. Appl. Genet. 129, 117–130. doi: 10.1007/s00122-015-2614-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W. K., Wang, Y. J., Luo, G. Z., Zhang, J. S., He, C. Y., Wu, X. L., et al. (2004). QTL mapping of ten agronomic traits on the soybean (Glycine max L. Merr.) genetic map and their association with EST markers. Theor. Appl. Genet. 108, 1131–1139. doi: 10.1007/s00122-003-1527-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X., Teng, W., Li, Y., Liu, D., Cao, G., Li, D., et al. (2017). Loci and candidate genes conferring resistance to soybean cyst nematode HG type 2.5.7. BMC Genomics 18:462. doi: 10.1186/s12864-017-3843-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., Jiang, Y., Wang, Z., Gou, Z., Lyu, J., Li, W., et al. (2015). Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat. Biotechnol. 33, 408–414. doi: 10.1038/nbt.3096

PubMed Abstract | CrossRef Full Text | Google Scholar