CALCIUM-DEPENDENT PROTEIN KINASE38 regulates flowering time and common cutworm resistance in soybean.

Photoperiod-sensitive plants such as soybean (Glycine max) often face threats from herbivorous insects throughout their whole growth period and especially during flowering; however, little is known about the relationship between plant flowering and insect resistance. Here, we used gene editing, multiple omics, genetic diversity and evolutionary analyses to confirm that the calcium-dependent protein kinase GmCDPK38 plays a dual role in coordinating flowering time regulation and insect resistance of soybean. Haplotype 2 (Hap2)-containing soybeans flowered later and were more resistant to the common cutworm (Spodoptera litura Fabricius) than those of Hap3. gmcdpk38 mutants with Hap3 knocked out exhibited similar flowering and resistance phenotypes as Hap2. Knocking out GmCDPK38 altered numerous flowering- and resistance-related phosphorylated proteins, genes, and metabolites. For example, the S-adenosylmethionine synthase GmSAMS1 was post-translationally upregulated in the gmcdpk38 mutants. GmCDPK38 has abundant genetic diversity in wild soybeans and was likely selected during soybean domestication. We found that Hap2 was mostly distributed at low latitudes and had a higher frequency in cultivars than in wild soybeans, while Hap3 was widely selected at high latitudes. Overall, our results elucidated that the two distinct traits (flowering time and insect resistance) are mediated by GmCDPK38.

Genetic dissection reveals the complex architecture of amino acid composition in soybean seeds.

KEY MESSAGE: Five loci related to soybean protein and amino acid contents were colocated by performing linkage mapping and GWAS. The haplotype analysis showed that Glyma.08G109100 may be useful to improve the soybean seed composition. Soybean (Glycine max (L.) Merr.) seeds are good protein sources. Although genetic variation is abundant, natural variation in seed amino acids and their derived traits is lacking across soybean accessions. Here, we determined the contents of protein and 17 amino acids, obtained 36 derived traits based on the protein and total amino acid contents, and derived 34 traits based on seven amino acid family groups. Furthermore, we performed a linkage analysis of the contents of 17 amino acids and 73 amino acid-derived traits based on the recombinant inbred line (RIL)-derived Kefeng No. 1 × Nannong 1138-2. Six hundred thirty-nine quantitative trait loci (QTLs) were identified, explaining 6.07-39.00% of the phenotypic variation. Among these loci, five were detected in diverse soybean accessions using a genome-wide association study. A network analysis revealed that some loci that were significantly associated with multiple amino acids were tightly linked on chromosome 8 based on linkage disequilibrium values, which also further confirmed the results of the correlation analysis among amino acid traits. Through a combination of a genome-wide association study, linkage analysis, qRT-PCR, and genomic polymorphism comparison, Glyma.08G109100 on chromosome 8, which may affect amino acid contents, was selected. The haplotype analysis showed that Hap-T of Glyma.08G109100 may be useful to improve the contents of protein and 16 amino acids in soybean. This study provides new insights into the genetic basis of the amino acid composition in soybean seeds and may facilitate marker-based breeding of soybean with improved nutritional value.

GmFtsH25 overexpression increases soybean seed yield by enhancing photosynthesis and photosynthates.

Increasing plant photosynthetic capacity is a promising approach to boost yields, but it is particularly challenging in C3 crops, such as soybean (Glycine max (L.) Merr.). Here, we identified GmFtsH25, encoding a member of the filamentation temperature-sensitive protein H protease family, as a major gene involved in soybean photosynthesis, using linkage mapping and a genome-wide association study. Overexpressing GmFtsH25 resulted in more grana thylakoid stacks in chloroplasts and increased photosynthetic efficiency and starch content, while knocking out GmFtsH25 produced the opposite phenotypes. GmFtsH25 interacted with photosystem I light harvesting complex 2 (GmLHCa2), and this interaction may contribute to the observed enhanced photosynthesis. GmFtsH25 overexpression lines had superior yield traits, such as yield per plant, compared to the wild type and knockout lines. Additionally, we identified an elite haplotype of GmFtsH25, generated by natural mutations, which appears to have been selected during soybean domestication. Our study sheds light on the molecular mechanism by which GmFtsH25 modulates photosynthesis and provides a promising strategy for improving the yields of soybean and other crops.

Ectopic expression of GmRNF1a encoding a soybean E3 ubiquitin ligase affects Arabidopsis silique development and dehiscence.

MAIN CONCLUSION: A soybean E3 ubiquitin ligase, GmRNF1a, may affect pod dehiscence and seed development through MADS family genes. These results would be useful for the study of soybean pod and seed development. Pod dehiscence is one of the critical causes of yield loss in cultivated soybeans, and it is of great significance to understand the molecular mechanisms underlying pod dehiscence in soybeans. In this study, we identified a new RING family member of the E3 ubiquitin ligase, GmRNF1a, which was observed to interact with the MADS-box protein GmAGL1 to regulate siliques dehiscence. Tissue-specific gene expression analysis revealed that GmRNF1a was mainly expressed in flowers and pods in soybean. The subcellular localization assay showed the nuclear and cytoplasmic localization of GmRNF1a. In addition, it was found that GmRNF1a exhibits higher promoter activity in soybean hairy roots as well as in Arabidopsis leaves, flowers, and siliques. Heterologous expression of GmRNF1a in Arabidopsis showed that the transgenic Arabidopsis siliques had a faster maturation rate and cracked earlier than the wild-type plants. The functional and nucleotide diversity analysis suggests that GmRNF1a might play an important role in pod maturation and dehiscence and has been strongly selected for during soybean domestication.

Downregulation of a gibberellin 3ß-hydroxylase enhances photosynthesis and increases seed yield in soybean.

Seed yield, determined mainly by seed numbers and seed weight, is the primary target of soybean breeding. Identifying the genes underlying yield-related traits is of great significance. Through joint linkage mapping and a genome-wide association study for 100-seed weight, we cloned GmGA3ox1, a gene encoding gibberellin 3ß-hydroxylase, which is the key enzyme in the gibberellin synthesis pathway. Genome resequencing identified a beneficial GmGA3ox1 haplotype contributing to high seed weight, which was further confirmed by soybean transformants. CRISPR/Cas9-generated gmga3ox1 mutants showed lower seed weight, but promoted seed yield by increasing seed numbers. The gmga3ox1 mutants reduced gibberellin biosynthesis while enhancing photosynthesis. Knockout of GmGA3ox1 resulted in the upregulation of numerous photosynthesis-related genes, particularly the GmRCA family encoding ribulose-1,5-bispho-sphate carboxylase-oxygenase (Rubisco) activases. The basic leucine zipper transcription factors GmbZIP97 and GmbZIP159, which were both upregulated in the gmga3ox1 mutants and induced by the gibberellin synthesis inhibitor uniconazole, could bind to the promoter of GmRCAß and activate its expression. Analysis of genomic sequences with over 2700 soybean accessions suggested that GmGA3ox1 is being gradually utilized in modern breeding. Our results elucidated the important role of GmGA3ox1 in soybean yield. These findings reveal important clues for future high-yield breeding in soybean and other crops.

Artificial selection on GmOLEO1 contributes to the increase in seed oil during soybean domestication.

Increasing seed oil content is one of the most important breeding goals for soybean due to a high global demand for edible vegetable oil. However, genetic improvement of seed oil content has been difficult in soybean because of the complexity of oil metabolism. Determining the major variants and molecular mechanisms conferring oil accumulation is critical for substantial oil enhancement in soybean and other oilseed crops. In this study, we evaluated the seed oil contents of 219 diverse soybean accessions across six different environments and dissected the underlying mechanism using a high-resolution genome-wide association study (GWAS). An environmentally stable quantitative trait locus (QTL), GqOil20, significantly associated with oil content was identified, accounting for 23.70% of the total phenotypic variance of seed oil across multiple environments. Haplotype and expression analyses indicate that an oleosin protein-encoding gene (GmOLEO1), colocated with a leading single nucleotide polymorphism (SNP) from the GWAS, was significantly correlated with seed oil content. GmOLEO1 is predominantly expressed during seed maturation, and GmOLEO1 is localized to accumulated oil bodies (OBs) in maturing seeds. Overexpression of GmOLEO1 significantly enriched smaller OBs and increased seed oil content by 10.6% compared with those of control seeds. A time-course transcriptomics analysis between transgenic and control soybeans indicated that GmOLEO1 positively enhanced oil accumulation by affecting triacylglycerol metabolism. Our results also showed that strong artificial selection had occurred in the promoter region of GmOLEO1, which resulted in its high expression in cultivated soybean relative to wild soybean, leading to increased seed oil accumulation. The GmOLEO1 locus may serve as a direct target for both genetic engineering and selection for soybean oil improvement.

A cation diffusion facilitator, GmCDF1, negatively regulates salt tolerance in soybean.

Salt stress is one of the major abiotic factors that affect the metabolism, growth and development of plants, and soybean [Glycine max (L.) Merr.] germination is sensitive to salt stress. Thus, to ensure the successful establishment and productivity of soybeans in saline soil, the genetic mechanisms of salt tolerance at the soybean germination stage need to be explored. In this study, a population of 184 recombinant inbred lines (RILs) was utilized to map quantitative trait loci (QTLs) related to salt tolerance. A major QTL related to salt tolerance at the soybean germination stage named qST-8 was closely linked with the marker Sat_162 and detected on chromosome 8. Interestingly, a genome-wide association study (GWAS) identified several single nucleotide polymorphisms (SNPs) significantly associated with salt tolerance in the same genetic region on chromosome 8. Resequencing, bioinformatics and gene expression analyses were implemented to identify the candidate gene Glyma.08g102000, which belongs to the cation diffusion facilitator (CDF) family and was named GmCDF1. Overexpression and RNA interference of GmCDF1 in soybean hairy roots resulted in increased sensitivity and tolerance to salt stress, respectively. This report provides the first demonstration that GmCDF1 negatively regulates salt tolerance by maintaining K+-Na+ homeostasis in soybean. In addition, GmCDF1 affected the expression of two ion homeostasis-associated genes, salt overly sensitive 1 (GmSOS1) and Na+/H+ exchanger 1 (GmNHX1), in transgenic hairy roots. Moreover, a haplotype analysis detected ten haplotypes of GmCDF1 in 31 soybean genotypes. A candidate-gene association analysis showed that two SNPs in GmCDF1 were significantly associated with salt tolerance and that Hap1 was more sensitive to salt stress than Hap2. The results demonstrated that the expression level of GmCDF1 was negatively correlated with salt tolerance in the 31 soybean accessions (r = -0.56, P < 0.01). Taken together, these results not only indicate that GmCDF1 plays a negative role in soybean salt tolerance but also help elucidate the molecular mechanisms of salt tolerance and accelerate the breeding of salt-tolerant soybean.

Genome-wide association mapping for protein, oil and water-soluble protein contents in soybean.

As a globally important legume crop, soybean provides excellent sources of protein and oil for human and livestock nutrition. Improving seed protein and oil contents has always been an important objective in soybean breeding. Water-soluble protein plays a significant role in the processing and efficacy of soybean protein. Here, a genome-wide association study (GWAS) of seed compositions (protein, oil, and water-soluble protein contents) was conducted using 211 diverse soybean accessions genotyped with a 355 K SoySNP array. Three, four, and five QTLs were identified related to the protein, oil, and water-soluble protein contents, respectively. Furthermore, five QTLs (qPC-15-1, qOC-8-1, qOC-12-1, qOC-20-1 and qWSPC-8-1) were detected in multiple environments. Analysis of the favorable alleles for oil and water-soluble protein contents showed that qOC-8-1 (qWSPC-8-1) exerted inverse effects on oil and water-soluble protein synthesis. Relative expression analysis suggested that Glyma.15G049200 in qPC-15-1 affects protein synthesis and Glyma.08G107800 in qOC-8-1 and qWSPC-8-1 might be involved in oil and water-soluble protein synthesis, producing opposite effects. The candidate genes and significant SNPs detected in the present study will allow a deeper understanding of the genetic basis for the regulation of protein, oil and water-soluble protein contents and provide important information that could be utilized in marker-assisted selection for soybean quality improvement.

Linkage and association study discovered loci and candidate genes for glycinin and ß-conglycinin in soybean (Glycine max L. Merr.).

KEY MESSAGE: Linkage mapping and GWAS identified 67 QTLs related to soybean glycinin, ß-conglycinin and relevant traits. Polymorphisms of the candidate gene Gy1 promoter were associated with the glycinin content in soybean. The major components of storage proteins in soybean seeds are glycinin and ß-conglycinin, which play important roles in determining protein nutrition and soy food processing properties. Increasing the protein content while improving the ratio of glycinin to ß-conglycinin is substantially important for soybean protein improvement. To investigate the genetic mechanism of storage protein subunits, 184 recombinant inbred lines (RILs) derived from a cross of Kefeng No. 1 and Nannong 1138-2 and 211 diverse soybean cultivars were used to detect loci related to glycinin (11S), ß-conglycinin (7S), the sum of glycinin and ß-conglycinin (SGC), and the ratio of glycinin to ß-conglycinin (RGC). Sixty-seven QTLs and 11 hot genomic regions were identified as affecting the four traits. One genetic region (q10-1) on chromosome 10 was associated with multiple traits by both linkage and association analysis. Eight genes in 11 hot genomic regions might be related to soybean protein subunit. The candidate gene analysis showed that polymorphisms in Gy1 promoters were significantly correlated with the 11S content. The QTLs and candidate genes identified in the present study allow for further understanding the genetic basis of 11S and 7S regulation and provide useful information for marker-assisted selection (MAS) in soybean quality improvement.

Association analysis of GmMAPKs and functional characterization of GmMMK1 to salt stress response in soybean.

Salinity is one of the major abiotic constraints affecting the growth and yield of plants including soybean. In this context, the previous studies have documented the role of the mitogen-activated protein kinase (MAPK) cascade in the regulation of salt signaling in model plants. However, there is not a systematic analysis of salt-related MAPKs in soybean. Hence, in this study, we identified a total of 32 GmMAPKs via., genome-wide reanalysis of the MAPK family using the soybean genome v4.0. Based on the transcriptome datasets in the public database, we observed that GmMAPKs are induced by different abiotic stresses, especially salt stress. Furthermore, based on the candidate gene association mapping and haplotype analysis of the GmMAPKs, we identified a salt-related MAPK member, GmMMK1. GmMMK1 possesses significant sequence variations, which affect salt tolerance in soybean at the germination stage. Besides, the overexpression of the GmMMK1 in soybean hairy roots has a significant negative effect on the root growth, leading to increased sensitivity of the GmMMK1-OE plants to salt stress. Moreover, the heterologous expression of the GmMMK1 in Arabidopsis has been also observed to have a negative effect on the germination and root growth under salt stress. The transcriptome analysis and yeast two-hybrid screening showed that hormone signaling and the homeostasis of reactive oxygen species are involved in the GmMMK1 regulation network. In conclusion, the results of this work demonstrated that GmMMK1 is an important negative regulator of the salt stress response, and provides better insights for understanding the role of the MAPKs in soybean salt signaling.

Genetic dissection of yield-related traits via genome-wide association analysis across multiple environments in wild soybean (Glycine soja Sieb. and Zucc.).

MAIN CONCLUSION: A total of 41 SNPs were identified as significantly associated with five yield-related traits in wild soybean populations across multiple environments, and the candidate gene GsCID1 was found to be associated with seed weight. These results may facilitate improvements in cultivated soybean. Crop-related wild species contain new sources of genetic diversity for crop improvement. Wild soybean (Glycine soja Sieb. and Zucc.) is the progenitor of cultivated soybean [Glycine max (L.) Merr.] and can be used as an essential genetic resource for yield improvements. In this research, using genome-wide association study (GWAS) in 96 out of 113 wild soybean accessions with 114,090 single nucleotide polymorphisms (SNPs) (with minor allele frequencies ≤ 0.05), SNPs associated with five yield-related traits were identified across multiple environments. In total, 41 SNPs were significantly associated with the traits in two or more environments (significance threshold P ≤ 8.76 × 10-6), with 29, 7, 3, and 2 SNPs detected for 100-seed weight (SW), maturity time (MT), seed yield per plant (SY) and flowering time (FT), respectively. BLAST search against the Glycine soja W05 reference genome was performed, 20 candidate genes were identified based on these 41 significant SNPs. One candidate gene, GsCID1 (Glysoja.04g010563), harbored two significant SNPs-AX-93713187, with a non-synonymous mutation, and AX-93713188, with a synonymous mutation. GsCID1 was highly expressed during seed development based on public information resources. The polymorphisms in this gene were associated with SW. We developed a derived cleaved amplified polymorphic sequence (dCAPS) marker for GsCID1 that was highly associated with SW and was validated as a functional marker. In summary, the revealed SNPs/genes are useful for understanding the genetic architecture of yield-related traits in wild soybean, which could be used as a potential exotic resource to improve cultivated soybean yields.

The acid phosphatase-encoding gene GmACP1 contributes to soybean tolerance to low-phosphorus stress.

Phosphorus (P) is essential for all living cells and organisms, and low-P stress is a major factor constraining plant growth and yield worldwide. In plants, P efficiency is a complex quantitative trait involving multiple genes, and the mechanisms underlying P efficiency are largely unknown. Combining linkage analysis, genome-wide and candidate-gene association analyses, and plant transformation, we identified a soybean gene related to P efficiency, determined its favorable haplotypes and developed valuable functional markers. First, six major genomic regions associated with P efficiency were detected by performing genome-wide associations (GWAs) in various environments. A highly significant region located on chromosome 8, qPE8, was identified by both GWAs and linkage mapping and explained 41% of the phenotypic variation. Then, a regional mapping study was performed with 40 surrounding markers in 192 diverse soybean accessions. A strongly associated haplotype (P = 10(-7)) consisting of the markers Sat_233 and BARC-039899-07603 was identified, and qPE8 was located in a region of approximately 250 kb, which contained a candidate gene GmACP1 that encoded an acid phosphatase. GmACP1 overexpression in soybean hairy roots increased P efficiency by 11-20% relative to the control. A candidate-gene association analysis indicated that six natural GmACP1 polymorphisms explained 33% of the phenotypic variation. The favorable alleles and haplotypes of GmACP1 associated with increased transcript expression correlated with higher enzyme activity. The discovery of the optimal haplotype of GmACP1 will now enable the accurate selection of soybeans with higher P efficiencies and improve our understanding of the molecular mechanisms underlying P efficiency in plants.

Association analysis for detecting significant single nucleotide polymorphisms for phosphorus-deficiency tolerance at the seedling stage in soybean [Glycine max (L) Merr].

Tolerance to low-phosphorus soil is a desirable trait in soybean cultivars. Previous quantitative trait locus (QTL) studies for phosphorus-deficiency tolerance were mainly derived from bi-parental segregating populations and few reports from natural population. The objective of this study was to detect QTLs that regulate phosphorus-deficiency tolerance in soybean using association mapping approach. Phosphorus-deficiency tolerance was evaluated according to five traits (plant shoot height, shoot dry weight, phosphorus concentration, phosphorus acquisition efficiency and use efficiency) comprising a conditional phenotype at the seedling stage. Association mapping of the conditional phenotype detected 19 SNPs including 13 SNPs that were significantly associated with the five traits across two years. A novel cluster of SNPs, including three SNPs that consistently showed significant effects over two years, that associated with more than one trait was detected on chromosome 3. All favorable alleles, which were determined based on the mean of conditional phenotypic values of each trait over the two years, could be pyramided into one cultivar through parental cross combination. The best three cross combinations were predicted with the aim of simultaneously improving phosphorus acquisition efficiency and use efficiency. These results will provide a thorough understanding of the genetic basis of phosphorus deficiency tolerance in soybean.

Identification of novel loci for salt stress at the seed germination stage in soybean.

Salt tolerance in soybean [Glycine max (L.) Merr.] at the seed germination stage is a critical determinant of stable stand establishment in saline soil. This study examined one population of 184 recombinant inbred lines (RILs, F7:11) derived from a cross between Kefeng1 and Nannong1138-2 and one natural population consisting of 196 soybean landraces. A total of 11 quantitative trait loci (QTLs) and 22 simple sequence repeat (SSR) loci associated with three salt tolerance indices were detected by linkage and association mapping. The SSR marker Sat_162 was found to be closely linked to the co-localized QTLs at a site 792,811 bp from the gene Glyma08g12400.1, which was verified in response to salt stress at the germination stage. Five SSR markers, Satt201, BE475343, CSSR306, Satt664 and Satt567, were co-associated with two of the salt tolerance indices, and two SSR markers, Satt156 and Satt636, were co-associated with all three salt tolerance indices. Furthermore, elite alleles and their carrier materials were identified by analyzing alleles at the loci associated with these salt tolerance indices. These results may be beneficial for the future breeding of soybean salt tolerance at the germination stage using marker-assisted selection and molecular pyramiding breeding.

Association mapping of soybean seed germination under salt stress.

Soil salinity is a serious threat to agriculture sustainability worldwide. Seed germination is a critical phase that ensures the successful establishment and productivity of soybeans in saline soils. However, little information is available regarding soybean salt tolerance at the germination stage. The objective of this study was to identify the genetic mechanisms of soybean seed germination under salt stress. One natural population consisting of 191 soybean landraces was used in this study. Soybean seeds produced in four environments were used to evaluate the salt tolerance at their germination stage. Using 1142 single-nucleotide polymorphisms (SNPs), the molecular markers associated with salt tolerance were detected by genome-wide association analysis. Eight SNP-trait associations and 13 suggestive SNP-trait associations were identified using a mixed linear model and the TASSEL 4.0 software. Eight SNPs or suggestive SNPs were co-associated with two salt tolerance indices, namely (1) the ratio of the germination index under salt conditions to the germination index under no-salt conditions (ST-GI) and (2) the ratio of the germination rate under salt conditions to the germination rate under no-salt conditions (ST-GR). One SNP (BARC-021347-04042) was significantly associated with these two traits (ST-GI and ST-GR). In addition, nine possible candidate genes were located in or near the genetic region where the above markers were mapped. Of these, five genes, Glyma08g12400.1, Glyma08g09730.1, Glyma18g47140.1, Glyma09g00460.1, and Glyma09g00490.3, were verified in response to salt stress at the germination stage. The SNPs detected could facilitate a better understanding of the genetic basis of soybean salt tolerance at the germination stage, and the marker BARC-021347-04042 could contribute to future breeding for soybean salt tolerance by marker-assisted selection.

Use of single nucleotide polymorphisms and haplotypes to identify genomic regions associated with protein content and water-soluble protein content in soybean.

KEY MESSAGE: Four major SPC-specific loci were identified, and these accounted for 8.5-15.1 % of the phenotypic variation, thus explaining why certain soybean varieties have a high PC but a low SPC. Water-soluble protein content (SPC) is a critical factor in both food quality and the production of isolated soybean proteins. However, few data are available regarding the genetic control and the mechanisms contributing to elevated SPC. In this study, a soybean collection of 192 accessions from a wide geographic range was used to identify genomic regions associated with soybean protein content (PC) and SPC using an association mapping approach employing 1,536 SNP makers and 232 haplotypes. The diverse panel revealed a large genetic variation in PC and SPC. Association mapping was performed using three methods to minimize false-positive associations. This resulted in 4/8 SNPs and 3/6 haplotypes that were significantly associated with soybean PC/SPC in two or more environments based on the mixed model. An SNP that was highly significantly associated with PC, BARC-021267-04016, was localized 0.28 cM away from a published glycinin gene, G7, and was detected across all four environments. Four major SPC-specific loci, BARC-029149-06088, BARC-018023-02499, BARC-041663-08059 and haplotype 15 (hp15), were stably identified on chromosomes five and eight and explained 8.5-15.1 % of the phenotypic variation. Moreover, a glutelin type-B 2-like gene was identified on chromosome eight and may be related to soybean protein solubility. These markers, which are located in previously reported QTL, reconfirmed previous findings and may be important targets for the identification of protein-related genes. These novel SNPs and haplotypes are important for further understanding the genetic basis of PC and SPC. In addition, by comparing the correlation and genetic loci between PC and SPC, we provide new insights into why certain soybean varieties have a high protein content but a low SPC.

Association mapping of yield-related traits and SSR markers in wild soybean (Glycine soja Sieb. and Zucc.).

Wild soybean, the progenitor of cultivated soybean, is an important gene pool for ongoing soybean breeding efforts. To identify yield-enhancing quantitative trait locus (QTL) or gene from wild soybean, 113 wild soybeans accessions were phenotyped for five yield-related traits and genotyped with 85 simple sequence repeat (SSR) markers to conduct association mapping. A total of 892 alleles were detected for the 85 SSR markers, with an average 10.49 alleles; the corresponding PIC values ranged from 0.07 to 0.92, with an average 0.73. The genetic diversity of each SSR marker ranged from 0.07 to 0.93, with an average 0.75. A total of 18 SSR markers were identified for the five traits. Two SSR markers, sct_010 and satt316, which are associated with the yield per plant were stably expressed over two years at two experimental locations. Our results suggested that association mapping can be an effective approach for identifying QTL from wild soybean.

Determination of the genetic architecture of seed size and shape via linkage and association analysis in soybean (Glycine max L. Merr.).

Seed-size traits, which are controlled by multiple genes in soybean, play an important role in determining seed yield, quality and appearance. However, the molecular mechanisms controlling the size of soybean seeds remain unclear, and little research has been done to investigate these mechanisms. In this study, we performed a genetic analysis to determine the genetic architecture of soybean seed size and shape via linkage and association analyses. We used 184 recombinant inbred lines (RILs) and 219 cultivated soybean accessions to evaluate seed length, seed width and seed height as seed-size traits, and their ratios of these values as seed-shape traits. Our results showed that all six traits had high heritability ranging from 92.46 to 98.47 %. Linkage analysis in the RILs identified 12 quantitative traits loci (QTLs), with five of these QTLs being associated with seed size, five with seed shape and two with the two first principal components of our principal component analysis (PCA). Association analysis in the 219 accessions detected 41 single nucleotide polymorphism (SNP)-trait associations, with 20 of these SNPs being associated with seed-size traits, seven with seed-shape traits and 14 with the two first principal components of our PCA. This analysis reveals that seed-size and seed-shape may be controlled by different genetic factors. Our results provide a greater understanding of phenotypic structure and genetic architecture of soybean seed, and the QTLs detected in this study form a basis for future fine mapping, quantitative trait gene cloning and molecular breeding in soybean.

QTL Mapping for Protein and Sulfur-Containing Amino Acid Contents Using a High-Density Bin-Map in Soybean (Glycine max L. Merr.).

Soybean provides essential protein and amino acids for humans and animals, while sulfur-containing amino acids (SAA), including methionine (Met) and cysteine (Cys), are very limited. In this study, we constructed a high-density bin-map with 3420 bin markers using 676 857 SNPs of a recombinant-inbred line (RIL) population derived from a cross between Kefeng no. 1 and Nannong 1138-2. Quantitative trait loci (QTL) mapping was performed for Cys, Met, SAA, and the protein content using this high-density bin-map. Twenty-five QTLs linked to these four traits were identified, and four genomic regions located on chromosomes (Chr) 07, 08, 15, and 20 were overlapped by multiple QTLs. Among them, bin 115-124 located on Chr 15 was associated with all four traits and was a novel locus with a high LOD value. These findings will provide a basis for nutritional quality improvement using marker-assisted selection breeding and clarify the genetic mechanisms of SAA and protein in soybean.

Identification of Loci and Candidate Genes Responsible for Pod Dehiscence in Soybean via Genome-Wide Association Analysis Across Multiple Environments.

Pod dehiscence (shattering) is the main cause of serious yield loss during the soybean mechanical harvesting process. A better understanding of the genetic architecture and molecular mechanisms of pod dehiscence is of great significance for soybean breeding. In this study, genome-wide association analysis (GWAS) with NJAU 355K SoySNP array was performed to detect single nucleotide polymorphisms (SNPs) associated with pod dehiscence in an association panel containing 211 accessions across five environments. A total of 163 SNPs were identified as significantly associated with pod dehiscence. Among these markers, 136 SNPs identified on chromosome 16 were located in the known QTL qPDH1. One, one, three, eleven, three, one, three, three and one SNPs were distributed on chromosome 1, 4, 6, 8, 9, 11, 17, 18, and 20, respectively. Favorable SNPs and six haplotypes were identified based on ten functional SNPs; among those Hap2 and Hap3 were considered as optimal haplotypes. In addition, based on GWAS results, the candidate gene Glyma09g06290 was identified. Quantitative real-time PCR (qRT-PCR) results and polymorphism analysis suggested that Glyma09g06290 might be involved in pod dehiscence. Furthermore, a derived cleaved amplified polymorphic sequences (dCAPS) marker for Glyma09g06290 was developed. Overall, the loci and genes identified in this study will be helpful in breeding soybean accessions resistant to pod dehiscence.