QTL Analysis of Major Agronomic Traits in Soybean

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

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

Identification of superior and rare haplotypes to optimize branch number in soybean

KEY MESSAGE

Using the integrated approach in the present study, we identified eleven significant SNPs, seven stable QTLs and 20 candidate genes associated with branch number in soybean. Branch number is a key yield-related quantitative trait that directly affects the number of pods and seeds per soybean plant. In this study, an integrated approach with a genome-wide association study (GWAS) and haplotype and candidate gene analyses was used to determine the detailed genetic basis of branch number across a diverse set of soybean accessions. The GWAS revealed a total of eleven SNPs significantly associated with branch number across three environments using the five GWAS models. Based on the consistency of the SNP detection in multiple GWAS models and environments, seven genomic regions within the physical distance of ± 202.4 kb were delineated as stable QTLs. Of these QTLs, six QTLs were novel, viz., qBN7, qBN13, qBN16, qBN18, qBN19 and qBN20, whereas the remaining one, viz., qBN12, has been previously reported. Moreover, 11 haplotype blocks, viz., Hap4, Hap7, Hap12, Hap13A, Hap13B, Hap16, Hap17, Hap18, Hap19A, Hap19B and Hap20, were identified on nine different chromosomes. Haplotype allele number across the identified haplotype blocks varies from two to five, and different branch number phenotype is regulated by these alleles ranging from the lowest to highest through intermediate branching. Furthermore, 20 genes were identified underlying the genomic region of ± 202.4 kb of the identified SNPs as putative candidates; and six of them showed significant differential expression patterns among the soybean cultivars possessing contrasting branch number, which might be the potential candidates regulating branch number in soybean. The findings of this study can assist the soybean breeding programs for developing cultivars with desirable branch numbers.

Identification of a Branch Number Locus in Soybean Using BSA-Seq and GWAS Approaches

The determination of the soybean branch number plays a pivotal role in plant morphogenesis and yield components. This polygenic trait is subject to environmental influences, and despite its significance, the genetic mechanisms governing the soybean branching number remain incompletely understood. To unravel these mechanisms, we conducted a comprehensive investigation employing a genome-wide association study (GWAS) and bulked sample analysis (BSA). The GWAS revealed 18 SNPs associated with the soybean branch number, among which qGBN3 on chromosome 2 emerged as a consistently detected locus across two years, utilizing different models. In parallel, a BSA was executed using an F2 population derived from contrasting cultivars, Wandou35 (low branching number) and Ruidou1 (high branching number). The BSA results pinpointed a significant quantitative trait locus (QTL), designated as qBBN1, located on chromosome 2 by four distinct methods. Importantly, both the GWAS and BSA methods concurred in co-locating qGBN3 and qBBN1. In the co-located region, 15 candidate genes were identified. Through gene annotation and RT-qPCR analysis, we predicted that Glyma.02G125200 and Glyma.02G125600 are candidate genes regulating the soybean branch number. These findings significantly enhance our comprehension of the genetic intricacies regulating the branch number in soybeans, offering promising candidate genes and materials for subsequent investigations aimed at augmenting the soybean yield. This research represents a crucial step toward unlocking the full potential of soybean cultivation through targeted genetic interventions.

A loss-of-function mutant allele of a glycosyl hydrolase gene has been co-opted for seed weight control during soybean domestication

The resultant DNA from loss-of-function mutation can be recruited in biological evolution and development. Here, we present such a rare and potential case of "to gain by loss" as a neomorphic mutation during soybean domestication for increasing seed weight. Using a population derived from a chromosome segment substitution line of Glycine max (SN14) and Glycine soja (ZYD06), a quantitative trait locus (QTL) of 100-seed weight (qHSW) was mapped on chromosome 11, corresponding to a truncated β-1, 3-glucosidase (βGlu) gene. The novel gene hsw results from a 14-bp deletion, causing a frameshift mutation and a premature stop codon in the βGlu. In contrast to HSW, the hsw completely lost βGlu activity and function but acquired a novel function to promote cell expansion, thus increasing seed weight. Overexpressing hsw instead of HSW produced large soybean seeds, and surprisingly, truncating hsw via gene editing further increased the seed size. We further found that the core 21-aa peptide of hsw and its variants acted as a promoter of seed size. Transcriptomic variation in these transgenic soybean lines substantiated the integration hsw into cell and seed size control. Moreover, the hsw allele underwent selection and expansion during soybean domestication and improvement. Our work cloned a likely domesticated QTL controlling soybean seed weight, revealed a novel genetic variation and mechanism in soybean domestication, and provided new insight into crop domestication and breeding, and plant evolution.

Genetic Diversity Analysis of Soybean Collection Using Simple Sequence Repeat Markers

Soybean [Glycine max (L.) Merr.] is a nutrient-rich crop that offers a sustainable source of dietary protein and edible oil. Determining the level of genetic diversity and relationships between various genetic resources involved in breeding programs is very important in crop improvement strategies. This study evaluated 100 soybean accessions with diverse origins for 10 important agronomic traits, including plant height (PH), an important plant adaptation-related trait impacting yield, in conditions in southeastern Kazakhstan for 2 years. The comparison of different groups of PH (tall, middle, and short) using a t-test suggested that the group of plants with the tallest PH provided a higher yield (p < 0.001) in relatively dry field conditions. The genetic diversity of the accessions was estimated using 25 simple sequence repeat (SSR) markers previously known to be associated with plant height. The results showed a significant variation among different groups of origin for all measured agronomic traits, as well as high genetic diversity, with the PIC (polymorphism information content) varying from 0.140 to 0.732, with an average of 0.524. Nei's diversity index ranged between 0.152 and 0.747, with an average of 0.526. The principal coordinate analysis (PCoA) of the studied soybean collection showed that Kazakhstan accessions were genetically distant from European, East Asian, and North American cultivars. Twelve out of twenty-five SSR markers demonstrated significant associations with ten studied agronomic traits, including PH (p < 0.05). Six SSRs with pleiotropic effects for studied traits were selected, and their haplotypes with phenotypic effects were generated for each soybean accession. The obtained results can be used in soybean improvement programs, including molecular-assisted breeding projects.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

The genetic basis of shoot architecture in soybean

Shoot architecture refers to the three-dimensional body plan of the above ground organs of the plant. The patterning of this body plan results from the tight genetic control of the size and maintenance of meristems, the initiation of axillary growth, and the timing of developmental phase transition. Variation in shoot architecture can result in dramatic differences in plant productivity and/or grain yield due to their effects on light interception, photosynthetic efficiency, response to agronomic inputs, and environmental adaptation. The fine-tuning of shoot architecture has consequently been of great interest to plant breeders, driving the need for deeper understanding of the genes and molecular mechanisms governing these traits. In soybean, the world's most important oil and protein crop, major components of shoot architecture include stem growth habit, plant height, branch angle, branch number, leaf petiole angle, and the size and shape of leaves. Key genes underlying some of these traits have been identified to integrate hormonal, developmental, and environmental signals modulating the growth and orientation of shoot organs. Here we summarize the current knowledge and recent advances in the understanding of the genetic control of these important architectural traits in soybean.

Identification of quantitative trait loci controlling soybean seed protein and oil content

Soybean is a major source of seed protein and oil globally with an average composition of 40% protein and 20% oil in the seed. The goal of this study was to identify quantitative trait loci (QTL) conferring seed protein and oil content utilizing a population constructed by crossing an above average protein content line, PI 399084 to another line that had a low protein content value, PI 507429, both from the USDA soybean germplasm collection. The recombinant inbred line (RIL) population, PI 507429 x PI 399084, was evaluated in two replications over four years (2018-2021); the seeds were analyzed for seed protein and oil content using near-infrared reflectance spectroscopy. The recombinant inbred lines and the two parents were re-sequenced using genotyping by sequencing. A total of 12,761 molecular markers, which came from genotyping by sequencing, the SoySNP6k BeadChip and selected simple sequence repeat (SSR) markers from known protein QTL chromosomal regions were used for mapping. One QTL was identified on chromosome 2 explaining up to 56.8% of the variation for seed protein content and up to 43% for seed oil content. Another QTL identified on chromosome 15 explained up to 27.2% of the variation for seed protein and up to 41% of the variation for seed oil content. The protein and oil QTLs of this study and their associated molecular markers will be useful in breeding to improve nutritional quality in soybean.

Genetic mapping and functional genomics of soybean seed protein

Soybean is an utterly important crop for high-quality meal protein and vegetative oil. Soybean seed protein content has become a key factor in nutrients for livestock feed as well as human dietary consumption. Genetic improvement of soybean seed protein is highly desired to meet the demands of rapidly growing world population. Molecular mapping and genomic analysis in soybean have identified many quantitative trait loci (QTL) underlying seed protein content control. Exploring the mechanisms of seed storage protein regulation will be helpful to achieve the improvement of protein content. However, the practice of breeding higher protein soybean is challenging because soybean seed protein is negatively correlated with seed oil content and yield. To overcome the limitation of such inverse relationship, deeper insights into the property and genetic control of seed protein are required. Recent advances of soybean genomics have strongly enhanced the understandings for molecular mechanisms of soybean with better seed quality. Here, we review the research progress in the genetic characteristics of soybean storage protein, and up-to-date advances of molecular mappings and genomics of soybean protein. The key factors underlying the mechanisms of the negative correlation between protein and oil in soybean seeds are elaborated. We also briefly discuss the future prospects of breaking the bottleneck of the negative correlation to develop high protein soybean without penalty of oil and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01373-5.

Identification of Quantitative Trait Locus and Candidate Genes for Drought Tolerance in a Soybean Recombinant Inbred Line Population

With global warming and regional decreases in precipitation, drought has become a problem worldwide. As the number of arid regions in the world is increasing, drought has become a major factor leading to significant crop yield reductions and food crises. Soybean is a crop that is relatively sensitive to drought. It is also a crop that requires more water during growth and development. The aim of this study was to identify the quantitative trait locus (QTL) that affects drought tolerance in soybean by using a recombinant inbred line (RIL) population from a cross between the drought-tolerant cultivar 'Jindou21' and the drought-sensitive cultivar 'Zhongdou33'. Nine agronomic and physiological traits were identified under drought and well-watered conditions. Genetic maps were constructed with 923,420 polymorphic single nucleotide polymorphism (SNP) markers distributed on 20 chromosomes at an average genetic distance of 0.57 centimorgan (cM) between markers. A total of five QTLs with a logarithm of odds (LOD) value of 4.035-8.681 were identified on five chromosomes. Under well-watered conditions and drought-stress conditions, one QTL related to the main stem node number was located on chromosome 16, accounting for 17.177% of the phenotypic variation. Nine candidate genes for drought resistance were screened from this QTL, namely Glyma.16G036700, Glyma.16G036400, Glyma.16G036600, Glyma.16G036800, Glyma.13G312700, Glyma.13G312800, Glyma.16G042900, Glyma.16G043200, and Glyma.15G100700. These genes were annotated as NAC transport factor, GATA transport factor, and BTB/POZ-MATH proteins. This result can be used for molecular marker-assisted selection and provide a reference for breeding for drought tolerance in soybean.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

Physiological breeding for yield improvement in soybean: solar radiation interception-conversion, and harvest index

KEY MESSAGE

Efficiency of light interception, Radiation use efficiency and harvest index can be used as targets to improve grain yield potential in soybean. Grain yield (GY) production can be expressed as the result of three main efficiencies: light interception (Ei), radiation use (RUE), and harvest index (HI). Although dissecting GY through these three efficiencies is not entirely new, there is a lack of knowledge about the phenotypic variation, the genetic architecture, and the relative contribution of these three efficiencies on GY in soybean. This knowledge gap coupled with laborious phenotyping prevents the active consideration of these efficiencies into breeding programs. This study aims to reveal the phenotypic variation, heritability, genetic relationships, genetic architecture, and genomic prediction for Ei, RUE, and HI in soybean. We evaluated a maturity control panel of 383 Recombinant Inbred Lines (RILs) selected from the soybean nested association mapping (SoyNAM) population. Dry matter ground measured along with canopy coverage (CC) from UAS imagery were collected in three environments. Light interception was modeled through a logistic curve using CC as a proxy. The total above-ground biomass collected during the growing season and its respective cumulative light intercepted were used to derive RUE through linear models fitting. Additive-genetic correlations, genome-wide association (GWA) and whole-genome regressions (WGR) were performed to evaluate the relationship between traits, their association with genomic regions, and the feasibility of predicting these efficiencies with genomic information. Correlation analyses considered three groups: the entire data set, and the high- and low-yielding RILs to determine association as a function of the GY. Our results revealed moderate to high phenotypic variation for Ei, RUE, and HI with ranges of 8.5%, 1.1 g MJ-1, and 0.2, respectively. Additive-genetic correlation revealed a strong relationship of GY with HI and moderate with RUE and Ei when whole data set was considered, but negligible contribution of HI on GY when just the top 100 was analyzed. The GWA analyses showed that Ei is associated with three SNPs; two of them located on chromosome 7 and one on chromosome 11 with no previous quantitative trait loci (QTLs) reported for these regions. RUE is associated with four SNPs on chromosomes 1, 7, 11, and 18. Some of these QTLs are novel, while others are previously documented for plant architecture and chlorophyll content. Two SNPs positioned on chromosome 13 and 15 with previous QTLs reported for plant height and seed set, weight and abortion were associated with HI. WGR showed high predictive ability for Ei, RUE, and HI with maximum correlation ranging between 0.75 and 0.80. Future improvements in GY can be expected through strategies prioritizing Ei for short-term results when using high yielding germplasm and RUE for medium- and long-term outcomes. This work is a pioneer attempt to integrate traditional physiological traits into the breeding process in the context of physiological breeding.

Identification of quantitative trait loci underlying five major agronomic traits of soybean in three biparental populations by specific length amplified fragment sequencing (SLAF-seq)

Flowering time, plant height, branch number, node numbers of main stem and pods per plant are important agronomic traits related to photoperiodic sensitivity, plant type and yield of soybean, which are controlled by multiple genes or quantitative trait loci (QTL). The main purpose of this study is to identify new QTL for five major agronomic traits, especially for flowering time. Three biparental populations were developed by crossing cultivars from northern and central China. Specific loci amplified fragment sequencing (SLAF-seq) was used to construct linkage map and QTL mapping was carried out. A total of 10 QTL for flowering time were identified in three populations, some of which were related to E1 and E2 genes or the other reported QTL listed in Soybase. In the Y159 population (Xudou No.9 × Kenfeng No.16), QTL for flowering time on chromosome 4, qFT4_1 and qFT4_2 were new. Compared with the QTL reported in Soybase, 1 QTL for plant height (PH), 3 QTL for branch number (BR), 5 QTL for node numbers of main stem, and 3 QTL for pods per plant were new QTL. Major E genes were frequently detected in different populations indicating that major the E loci had a great effect on flowering time and adaptation of soybean. Therefore, in order to further clone minor genes or QTL, it may be of great significance to carefully select the genotypes of known loci. These results may lay a foundation for fine mapping and clone of QTL/genes related to plant-type, provided a basis for high yield breeding of soybean.

Novel QTL and Meta-QTL Mapping for Major Quality Traits in Soybean

Isoflavone, protein, and oil are the most important quality traits in soybean. Since these phenotypes are typically quantitative traits, quantitative trait locus (QTL) mapping has been an efficient way to clarify their complex and unclear genetic background. However, the low-density genetic map and the absence of QTL integration limited the accurate and efficient QTL mapping in previous researches. This paper adopted a recombinant inbred lines (RIL) population derived from 'Zhongdou27'and 'Hefeng25' and a high-density linkage map based on whole-genome resequencing to map novel QTL and used meta-analysis methods to integrate the stable and consentaneous QTL. The candidate genes were obtained from gene functional annotation and expression analysis based on the public database. A total of 41 QTL with a high logarithm of odd (LOD) scores were identified through composite interval mapping (CIM), including 38 novel QTL and 2 Stable QTL. A total of 660 candidate genes were predicted according to the results of the gene annotation and public transcriptome data. A total of 212 meta-QTL containing 122 stable and consentaneous QTL were mapped based on 1,034 QTL collected from previous studies. For the first time, 70 meta-QTL associated with isoflavones were mapped in this study. Meanwhile, 69 and 73 meta-QTL, respectively, related to oil and protein were obtained as well. The results promote the understanding of the biosynthesis and regulation of isoflavones, protein, and oil at molecular levels, and facilitate the construction of molecular modular for great quality traits in soybean.

Fine mapping of a QTL locus (QNFSP07-1) and analysis of candidate genes for four-seeded pods in soybean

Soybean [Glycine max (L.) Merr.] is an important grain and oil crop in the world, and it is the main source of high-quality protein. The number of four-seeded pods is a quantitative trait in soybean and is closely related to yield in terms of breeding. Therefore, it is of great significance to study the inheritance of four-seed pods and to excavate related genes for improving soybean yield. In this study, individuals with high ratio of four-seed pods which from chromosome segment substitution lines (CSSLs) that can be stably inherited were selected as the parent, and Suinong 14 (SN14) was used as recurrent parent to construct secondary mapping population via marker-assisted selection. From 2006 to 2017, QTL analysis was performed using secondary mapping populations, and the initial QTL mapping interval was 0.67 Mb and was located on Gm07. Based on the initial QTL mapping results, individuals that were heterozygous at the interval (36,116,118-37,399,738 bp) were screened in 2018, and the heterozygous individuals were subjected to inbreeding to obtain 13 F3 populations, with a target interval of 321 kb. Gene annotation was performed on the fine mapping interval, and 27 genes were obtained. Among 27 genes, Glyma.07G200900 and Glyma.07G201200 were identified as candidate genes. qRT-PCR was used to measure the expression of the 2 candidate genes at different developmental stages of soybean, and the expression levels of the 2 candidate genes in terms of cell division (axillary buds, COTs, EMs) were higher than those in terms of cell expansion (MM, LM), and these genes play a positive regulatory role in the formation of four-seeded pods. Haplotype analysis of 2 candidate genes which shows that Glyma.07G201200 has two excellent haplotypes, and the significance level between the two excellent haplotypes at p < 0.05. Those results provide the information for gene map-based cloning and molecular marker-assisted breeding of the number of four-seeded pod in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01265-6.

QTL for Main Stem Node Number and Its Response to Plant Densities in 144 Soybean FW-RILs

Although the main stem node number of soybean [Glycine max (L.) Merr. ] is an important yield-related trait, there have been limited studies on the effect of plant density on the identification of quantitative trait loci (QTL) for main stem node number (MSNN). To address this issue, here, 144 four-way recombinant inbred lines (FW-RILs) derived from Kenfeng 14, Kenfeng 15, Heinong 48, and Kenfeng 19 were used to identify QTL for MSNN with densities of 2.2 × 105 (D1) and 3 × 105 (D2) plants/ha in five environments by linkage and association studies. As a result, the linkage and association studies identified 40 and 28 QTL in D1 and D2, respectively, indicating the difference in QTL in various densities. Among these QTL, five were common in the two densities; 36 were singly identified for response to density; 12 were repeatedly identified by both response to density and phenotype of two densities. Thirty-one were repeatedly detected across various methods, densities, and environments in the linkage and association studies. Among the 24 common QTL in the linkage and association studies, 15 explained a phenotypic variation of more than 10%. Finally, Glyma.06G094400, Glyma.06G147600, Glyma.19G160800.1, and Glyma.19G161100 were predicted to be associated with MSNN. These findings will help to elucidate the genetic basis of MSNN and improve molecular assistant selection in high-yield soybean breeding.

Quantitative Trait Locus Analysis of Protein and Oil Content in Response to Planting Density in Soybean (Glycine max [L.] Merri.) Seeds Based on SNP Linkage Mapping

Soybean varieties suitable for high planting density allow greater yields. However, the seed protein and oil contents, which determine the value of this crop, can be influenced by planting density. Thus, it is important to understand the genetic basis of the responses of different soybean genotypes to planting density. In this study, we quantified the protein and oil contents in a four-way recombinant inbred line (FW-RIL) soybean population under two planting densities and the response to density. We performed quantitative trait locus (QTL) mapping using a single nucleotide polymorphism (SNP) linkage map generated by inclusive composite interval mapping. We identified 14 QTLs for protein content and 17 for oil content at a planting density of 2.15 × 10⁵ plant/ha (D1) and 14 QTLs for protein content and 20 for oil content at a planting density 3.0 × 10⁵ plant/ha (D2). Among the QTLs detected, two oil-content QTLs was detected at both plant densities. In addition, we identified 38 QTLs for the responses of protein and oil contents to planting density. Of the QTLs detected, 70 were identified in previous studies, while 33 were newly identified. Fourty-five QTLs accounted for over 10% of the phenotypic variation of the corresponding trait, based on 23 QTLs at a marker interval distance of ~600 kb detected under different densities and with the responses to density difference. Pathway analysis revealed four candidate genes involved in protein and oil biosynthesis/metabolism. These results improve our understanding of the genetic underpinnings of protein and oil biosynthesis in soybean, laying the foundation for enhancing protein and oil contents and increasing yields in soybean.

Identification of QTNs and Their Candidate Genes for 100-Seed Weight in Soybean (Glycine max L.) Using Multi-Locus Genome-Wide Association Studies

100-seed weight (100-SW) in soybeans is a yield component trait and controlled by multiple genes with different effects, but limited information is available for its quantitative trait nucleotides (QTNs) and candidate genes. To better understand the genetic architecture underlying the trait and improve the precision of marker-assisted selection, a total of 43,834 single nucleotide polymorphisms (SNPs) in 250 soybean accessions were used to identify significant QTNs for 100-SW in four environments and their BLUP values using six multi-locus and one single-locus genome-wide association study methods. As a result, a total of 218 significant QTNs were detected using multi-locus methods, whereas eight QTNs were identified by a single-locus method. Among 43 QTNs or QTN clusters identified repeatedly across various environments and/or approaches, all of them exhibited significant trait differences between their corresponding alleles, 33 were found in the genomic region of previously reported QTLs, 10 were identified as new QTNs, and three (qHSW-4-1, qcHSW-7-3, and qcHSW-10-4) were detected in all the four environments. The number of seed weight (SW) increasing alleles for each accession ranged from 8 (18.6%) to 36 (83.72%), and three accessions (Yixingwuhuangdou, Nannong 95C-5, and Yafanzaodou) had more than 35 SW increasing alleles. Among 36 homologous seed-weight genes in Arabidopsis underlying the above 43 stable QTNs, more importantly, Glyma05g34120, GmCRY1, and GmCPK11 had known seed-size/weight-related genes in soybean, and Glyma07g07850, Glyma10g03440, and Glyma10g36070 were candidate genes identified in this study. These results provide useful information for genetic foundation, marker-assisted selection, genomic prediction, and functional genomics of 100-SW.

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

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

Exploring the genetic base of the soybean germplasm from Africa, America and Asia as well as mining of beneficial allele for flowering and seed weight

Genetic diversity is the foundation for any breeding program. The present study analyzed the genetic base of 163 soybean genotypes from three continents viz. Africa, America and Asia using 68 trait-linked simple sequence repeats (SSR) markers. The average number of alleles among the germplasm from the three continents followed the trend as Asia (9) > America (8) > Africa (7). Similar trends were observed for gene diversity (0.76 > 0.74 > 0.71) and polymorphism information content (PIC) (0.73 > 0.71 > 0.68). These findings revealed that soybean germplasm from Asia has wider genetic base followed by America, and least in Africa. The 163 genotypes were grouped into 4 clusters by phylogenetic analysis, whereas model-based population structure analysis also divided them into 4 subpopulations comprising 80.61% pure lines and 19.39% admixtures. The genotypes from Africa were easily distinguished from those of other two continents using phylogenetic analysis, indicating important role of geographyical differentiation for this genetic variability. Our results indicated that soybean germplasm has moved from Asia to America, and from America to Africa. Analysis of molecular variance (AMOVA) showed 8.41% variation among the four subpopulations, whereas 63.12% and 28.47% variation existed among and within individuals in the four subpopulations, respectively. Based on the association mapping, a total of 21 SSR markers showed significant association with days to flowering (DoF) and 100-seed weight (HSW). Two markers Satt365 and Satt581 on chromosome 6 and 10, respectively, showed pleiotropic effect or linkage on both traits. Genotype A50 (Gakuran Daizu/PI 506679) from Japan has 8 out of the 13 beneficial alleles for increased HSW. The diverse genotypes, polymorphic SSR markers and desirable alleles identified for DoF and HSW will be used in future breeding programs to improve reproductive, yield and quality traits.

Linkage Analysis and Multi-Locus Genome-Wide Association Studies Identify QTNs Controlling Soybean Plant Height

Plant height is an important target for soybean breeding. It is a typical quantitative trait controlled by multiple genes and is susceptible to environmental influences. Here, we carried out phenotypic analysis of 156 recombinant inbred lines derived from "Dongnong L13" and "Henong 60" in nine environments at four locations over 6 years using interval mapping and inclusive composite interval mapping methods. We performed quantitative trait locus (QTL) analysis by applying pre-built simple-sequence repeat maps. We detected 48 QTLs, including nine significant QTLs detected by multiple methods and in multiple environments. Meanwhile, genotyping of all lines using the SoySNP660k BeadChip produced 54,836 non-redundant single-nucleotide polymorphism (SNP) genotypes. We used five multi-locus genome-wide association analysis methods to locate 10 quantitative trait nucleotides (QTNs), four of which overlap with previously located QTLs. Five candidate genes related to plant height are predicted to lie within 200 kb of these four QTNs. We identified 19 homologous genes in Arabidopsis, two of which may be associated with plant height. These findings further our understanding of the multi-gene regulatory network and genetic determinants of soybean plant height, which will be important for breeding high-yielding soybean.

Combined Linkage Mapping and BSA to Identify QTL and Candidate Genes for Plant Height and the Number of Nodes on the Main Stem in Soybean

Soybean is one of the most important food and oil crops in the world. Plant height (PH) and the number of nodes on the main stem (NNMS) are quantitative traits closely related to soybean yield. In this study, we used 208 chromosome segment substitution lines (CSSL) populations constructed using "SN14" and "ZYD00006" for quantitative trait locus (QTL) mapping of PH and NNMS. Combined with bulked segregant analysis (BSA) by extreme materials, 8 consistent QTLs were identified. According to the gene annotation of the QTL interval, a total of 335 genes were obtained. Five of which were associated with PH and NNMS, potentially representing candidate genes. RT-qPCR of these 5 candidate genes revealed two genes with differential relative expression levels in the stems of different materials. Haplotype analysis showed that different single nucleotide polymorphisms (SNPs) between the excellent haplotypes in Glyma.04G251900 and Glyma.16G156700 may be the cause of changes in these traits. These results provide the basis for research on candidate genes and marker-assisted selection (MAS) in soybean breeding.

Mapping developmental QTL for plant height in soybean [Glycine max (L.) Merr.] using a four-way recombinant inbred line population

Plant height (PH) is an important trait in soybean, as taller plants may have higher yields but may also be at risk for lodging. Many genes act jointly to influence PH throughout development. To map the quantitative trait loci (QTL) controlling PH, we used the unconditional variable method (UVM) and conditional variable method (CVM) to analyze PH data for a four-way recombinant inbred line (FW-RIL) population derived from the cross of (Kenfeng14 × Kenfeng15) × (Heinong48 × Kenfeng19). We identified 7, 8, 16, 19, 15, 27, 17, 27, 22, and 24 QTL associated with PH at 10 developmental stages, respectively. These QTL mapped to 95 genomic regions. Among these QTL, 9 were detected using UVM and CVM, and 89 and 66 were only detected by UVM or CVM, respectively. In total, 36 QTL controlling PH were detected at multiple developmental stages and these made unequal contributions to genetic variation throughout development. Among 19 novel regions discovered in our study, 7 could explain over 10% of the phenotypic variation and contained only one single QTL. The unconditional and conditional QTL detected here could be used in molecular design breeding across the whole developmental procedure.

Genome-wide associations and epistatic interactions for internode number, plant height, seed weight and seed yield in soybean

BACKGROUND

Breeding programs benefit from information about marker-trait associations for many traits, whether the goal is to place those traits under active selection or to maintain them through background selection. Association studies are also important for identifying accessions bearing potentially useful alleles by characterizing marker-trait associations and allelic states across germplasm collections. This study reports the results of a genome-wide association study and evaluation of epistatic interactions for four agronomic and seed-related traits in soybean.

RESULTS

Using 419 diverse soybean accessions, together with genotyping data from the SoySNP50K Illumina Infinium BeadChip, we identified marker-trait associations for internode number (IN), plant height (PH), seed weight (SW), and seed yield per plant (SYP). We conducted a genome-wide epistatic study (GWES), identifying candidate genes that show evidence of SNP-SNP interactions. Although these candidate genes will require further experimental validation, several appear to be involved in developmental processes related to the respective traits. For IN and PH, these include the Dt1 determinacy locus (a soybean meristematic transcription factor), as well as a pectinesterase gene and a squamosa promoter binding gene that in other plants are involved in cell elongation and the vegetative-to-reproductive transition, respectively. For SW, candidate genes include an ortholog of the AP2 gene, which in other species is involved in maintaining seed size, embryo size, seed weight and seed yield. Another SW candidate gene is a histidine phosphotransfer protein - orthologs of which are involved in cytokinin-mediated seed weight regulating pathways. The SYP association loci overlap with regions reported in previous QTL studies to be involved in seed yield.

CONCLUSIONS

This study further confirms the utility of GWAS and GWES approaches for identifying marker-trait associations and interactions within a diverse germplasm collection.

Integration of lodging resistance QTL in soybean

Poor lodging resistance could limit increases in soybean yield. Previously, a considerable number of observations of quantitative trait loci (QTL) for lodging resistance have been reported by independent studies. The integration of these QTL into a consensus map will provide further evidence of their usefulness in soybean improvement. To improve informative QTL in soybean, a mapping population from a cross between the Harosoy and Clark cultivars, which inherit major U.S. soybean genetic backgrounds, was used along with previous mapping populations to identify QTL for lodging resistance. Together with 78 QTL for lodging collected from eighteen independent studies, a total of 88 QTL were projected onto the soybean consensus map. A total of 16 significant QTL clusters were observed; fourteen of them were confirmed in either two or more mapping populations or a single population subjected to different environmental conditions. Four QTL (one on chromosome 7 and three on 10) were newly identified in the present study. Further, meta-analysis was used to integrate QTL across different studies, resulting in two significant meta-QTL each on chromosomes 6 and 19. Our results provide deeper knowledge of valuable lodging resistance QTL in soybean, and these QTL could be used to increase lodging resistance.

QTL fine-mapping of soybean (Glycine max L.) leaf type associated traits in two RILs populations

BACKGROUND

The different leaf type associated traits of soybean (Glycine max L.) including leaf area, leaf length, leaf width, leaf shape and petiole length are considered to be associated with seed yield. In order to identify quantitative trait loci (QTLs) affecting leaf type traits, two advanced recombinant inbred line (RIL, ZH, Zhonghuang 24 × Huaxia 3; GB, Guizao 1 × Brazil 13) populations were introduced to score phenotypic values in plants across nine different environments (years, seasons, locations and soybean growth stages). Two restriction site-associated DNA sequencing (RAD-seq) based high-density genetic linkage maps with an average distance of 1.00 centimorgan (cM) between adjacent bin markers were utilized for QTL fine mapping.

RESULTS

Correlation analysis showed that most of the traits were correlated with each other and regulated both by hereditary and environmental factors. A total of 190 QTLs were identified for leaf type associated traits in the two populations, of which 14 loci were found to be environmentally stable. Moreover, these detected QTLs were categorized into 34 QTL hotspots, and four important QTL hotspots with phenotypic variance ranging from 3.89-23.13% were highlighted. Furthermore, Glyma04g05840, Glyma19g37820, Glyma14g07140 and Glyma19g39340 were predicted in the intervals of the stable loci and important QTL hotspots for leaf type traits by adopting Gene Ontology (GO) enrichment analysis.

CONCLUSIONS

Our findings of the QTLs and the putative genes will be beneficial to gain new insights into the genetic basis for soybean leaf type traits and may further accelerate the breeding process for reasonable leaf type soybean.

Genetic map construction and QTL analysis of leaf-related traits in soybean under monoculture and relay intercropping

Soybean (Glycine max L.) is an important food and oil crop widely planted by intercropping in southwest China. The shade caused by intercropping changes plant growth traits, such as soybean leaf and dry mass, thereby reducing yields. To improve the yield and elucidate the genetic mechanism of the leaf-related traits in intercropped soybeans, we measured the F6:7-8 recombinant inbred lines (RILs) derived from the cross of 'Nandou 12' and 'Jiuyuehuang' for six leaf-related traits under monoculture and relay intercropping in 2015 and 2016. We found 6366 single-nucleotide polymorphisms (SNPs) markers that covered the whole genome of soybean distributed in 20 linkage groups, which spanned 2818.67 cM with an average interval of 0.44 cM between adjacent markers. Nineteen quantitative trait loci (QTLs) were detected in two environments in 2 years. Three candidate genes associated to leaf-related traits were found according to gene expression and GO enrichment analyses. These results revealed the susceptibility of leaf phenotype to shading and helped elucidate the mechanisms that control leaf-related traits.

GmBRC1 is a Candidate Gene for Branching in Soybean (Glycine max (L.) Merrill)

Branch number is one of the main factors affecting the yield of soybean (Glycine max (L.)). In this study, we conducted a genome-wide association study combined with linkage analysis for the identification of a candidate gene controlling soybean branching. Five quantitative trait nucleotides (QTNs) were associated with branch numbers in a soybean core collection. Among these QTNs, a linkage disequilibrium (LD) block qtnBR6-1 spanning 20 genes was found to overlap a previously identified major quantitative trait locus qBR6-1. To validate and narrow down qtnBR6-1, we developed a set of near-isogenic lines (NILs) harboring high-branching (HB) and low-branching (LB) alleles of qBR6-1, with 99.96% isogenicity and different branch numbers. A cluster of single nucleotide polymorphisms (SNPs) segregating between NIL-HB and NIL-LB was located within the qtnBR6-1 LD block. Among the five genes showing differential expression between NIL-HB and NIL-LB, BRANCHED1 (BRC1; Glyma.06G210600) was down-regulated in the shoot apex of NIL-HB, and one missense mutation and two SNPs upstream of BRC1 were associated with branch numbers in 59 additional soybean accessions. BRC1 encodes TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORS 1 and 2 transcription factor and functions as a regulatory repressor of branching. On the basis of these results, we propose BRC1 as a candidate gene for branching in soybean.

GmBRC1 is a Candidate Gene for Branching in Soybean (Glycine max (L.) Merrill)

Branch number is one of the main factors affecting the yield of soybean (Glycine max (L.)). In this study, we conducted a genome-wide association study combined with linkage analysis for the identification of a candidate gene controlling soybean branching. Five quantitative trait nucleotides (QTNs) were associated with branch numbers in a soybean core collection. Among these QTNs, a linkage disequilibrium (LD) block qtnBR6-1 spanning 20 genes was found to overlap a previously identified major quantitative trait locus qBR6-1. To validate and narrow down qtnBR6-1, we developed a set of near-isogenic lines (NILs) harboring high-branching (HB) and low-branching (LB) alleles of qBR6-1, with 99.96% isogenicity and different branch numbers. A cluster of single nucleotide polymorphisms (SNPs) segregating between NIL-HB and NIL-LB was located within the qtnBR6-1 LD block. Among the five genes showing differential expression between NIL-HB and NIL-LB, BRANCHED1 (BRC1; Glyma.06G210600) was down-regulated in the shoot apex of NIL-HB, and one missense mutation and two SNPs upstream of BRC1 were associated with branch numbers in 59 additional soybean accessions. BRC1 encodes TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORS 1 and 2 transcription factor and functions as a regulatory repressor of branching. On the basis of these results, we propose BRC1 as a candidate gene for branching in soybean.

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.

Meta-analysis and transcriptome profiling reveal hub genes for soybean seed storage composition during seed development

Soybean is an important crop providing edible oil and protein source. Soybean oil and protein contents are quantitatively inherited and significantly affected by environmental factors. In this study, meta-analysis was conducted based on soybean physical maps to integrate quantitative trait loci (QTLs) from multiple experiments in different environments. Meta-QTLs for seed oil, fatty acid composition, and protein were identified. Of them, 11 meta-QTLs were located on hot regions for both seed oil and protein. Next, we selected 4 chromosome segment substitution lines with different seed oil and protein contents to characterize their 3 years of phenotype selection in the field. Using strand-specific RNA-sequencing analysis, we profile the time-course transcriptome patterns of soybean seeds at early maturity, middle maturity, and dry seed stages. Pairwise comparison and K-means clustering analysis revealed 7,482 differentially expressed genes and 45 expression patterns clusters. Weighted gene coexpression network analysis uncovered 46 modules of gene expression patterns. The 2 most significant coexpression networks were visualized, and 7 hub genes were identified that were involved in soybean oil and seed storage protein accumulation processes. Our results provided a transcriptome dataset for soybean seed development, and the candidate hub genes represent a foundation for further research.

Quantitative Trait Locus Mapping of Flowering Time and Maturity in Soybean Using Next-Generation Sequencing-Based Analysis

Soybean (Glycine max L.) is a major legume crop that is mainly distributed in temperate regions. The adaptability of soybean to grow at relatively high latitudes is attributed to natural variations in major genes and quantitative trait loci (QTLs) that control flowering time and maturity. Identification of new QTLs and map-based cloning of candidate genes are the fundamental approaches in elucidating the mechanism underlying soybean flowering and adaptation. To identify novel QTLs/genes, we developed two F8:10 recombinant inbred lines (RILs) and evaluated the traits of time to flowering (R1), maturity (R8), and reproductive period (RP) in the field. To rapidly and efficiently identify QTLs that control these traits, next-generation sequencing (NGS)-based QTL analysis was performed. This study demonstrates that only one major QTL on chromosome 4 simultaneously controls R1, R8, and RP traits in the Dongnong 50 × Williams 82 (DW) RIL population. Furthermore, three QTLs were mapped to chromosomes 6, 11, and 16 in the Suinong 14 × Enrei (SE) RIL population. Two major pleiotropic QTLs on chromosomes 4 and 6 were shown to affect flowering time, maturity, and RP. A QTL influencing RP was identified on chromosome 11, and QTL on chromosome 16 was associated with time to flowering responses. All these QTLs contributed to soybean maturation. The QTLs identified in this study may be utilized in fine mapping and map-based cloning of candidate genes to elucidate the mechanisms underlying flowering and soybean adaptation to different latitudes and to breed novel soybean cultivars with optimal yield-related traits.

Identification of QTLs related to the vertical distribution and seed-set of pod number in soybean [Glycine max (L.) Merri]

Pod number is an important factor that influences yield in soybean. Here, we used two associated recombinant inbred line (RIL) soybean populations, RIL3613 (containing 134 lines derived from Dongnong L13 × Heihe 36) and RIL6013 (composed of 156 individuals from Dongnong L13 × Henong 60), to identify quantitative trait loci (QTLs) regulating the vertical distribution and quantity of seeds and seed pods. The numbers of pods were quantified in the upper, middle, and lower sections of the plant, as well as in the plants as a whole, and QTLs regulating these spatial traits were mapped using an inclusive complete interval mapping method. A total of 21 and 26 QTLs controlling pod-number-related traits were detected in RIL3613 and RIL6013, respectively, which explained 1.25-11.6698% and 0.0001-7.91% of the phenotypic variation. A total of 34 QTLs were verified by comparison with previous research, were identified in both populations, or were found to regulate multiple traits, indicating their authenticity. These results enhance our understanding of the vertical distribution of pod-number-related traits and support molecular breeding for seed yield.

Construction of high-density genetic map and QTL mapping of yield-related and two quality traits in soybean RILs population by RAD-sequencing

BACKGROUND

One of the overarching goals of soybean breeding is to develop lines that combine increased yield with improved quality characteristics. High-density-marker QTL mapping can serve as an effective strategy to identify novel genomic information to facilitate crop improvement. In this study, we genotyped a recombinant inbred line (RIL) population (Zhonghuang 24 × Huaxia 3) using a restriction-site associated DNA sequencing (RAD-seq) approach. A high-density soybean genetic map was constructed and used to identify several QTLs that were shown to influence six yield-related and two quality traits.

RESULTS

A total of 47,472 single-nucleotide polymorphisms (SNPs) were detected for the RILs that were integrated into 2639 recombination bin units, with an average distance of 1.00 cM between adjacent markers. Forty seven QTLs for yield-related traits and 13 QTLs for grain quality traits were found to be distributed on 16 chromosomes in the 2 year studies. Among them, 18 QTLs were stable, and were identified in both analyses. Twenty six QTLs were identified for the first time, with a single QTL (qNN19a) in a 56 kb region explaining 32.56% of phenotypic variation, and an additional 10 of these were novel, stable QTLs. Moreover, 8 QTL hotpots on four different chromosomes were identified for the correlated traits.

CONCLUSIONS

With RAD-sequencing, some novel QTLs and important QTL clusters for both yield-related and quality traits were identified based on a new, high-density bin linkage map. Three predicted genes were selected as candidates that likely have a direct or indirect influence on both yield and quality in soybean. Our findings will be helpful for understanding common genetic control mechanisms of co-localized traits and to select cultivars for further analysis to predictably modulate soybean yield and quality simultaneously.

Meta-Analyses of QTLs Associated with Protein and Oil Contents and Compositions in Soybean [Glycine max (L.) Merr.] Seed

Soybean [Glycine max (L.) Merr.] is a valuable and nutritious crop in part due to the high protein meal and vegetable oil produced from its seed. Soybean producers desire cultivars with both elevated seed protein and oil concentrations as well as specific amino acid and fatty acid profiles. Numerous studies have identified quantitative trait loci (QTLs) associated with seed composition traits, but validation of these QTLs has rarely been carried out. In this study, we have collected information, including genetic location and additive effects, on each QTL for seed contents of protein and oil, as well as amino acid and fatty acid compositions from over 80 studies. Using BioMercator V. 4.2, a meta-QTL analysis was performed with genetic information comprised of 175 QTLs for protein, 205 QTLs for oil, 156 QTLs for amino acids, and 113 QTLs for fatty acids. A total of 55 meta-QTL for seed composition were detected on 6 out of 20 chromosomes. Meta-QTL possessed narrower confidence intervals than the original QTL and candidate genes were identified within each meta-QTL. These candidate genes elucidate potential natural genetic variation in genes contributing to protein and oil biosynthesis and accumulation, providing meaningful information to further soybean breeding programs.

Genotyping-by-Sequencing-Based Investigation of the Genetic Architecture Responsible for a ∼Sevenfold Increase in Soybean Seed Stearic Acid

Soybean oil is highly unsaturated but oxidatively unstable, rendering it nonideal for food applications. Until recently, the majority of soybean oil underwent partial chemical hydrogenation, which produces trans fats as an unavoidable consequence. Dietary intake of trans fats and most saturated fats are conclusively linked to negative impacts on cholesterol levels and cardiovascular health. Two major soybean oil breeding targets are: (1) to reduce or eliminate the need for chemical hydrogenation, and (2) to replace the functional properties of partially hydrogenated soybean oil. One potential solution is the elevation of seed stearic acid, a saturated fat which has no negative impacts on cardiovascular health, from 3 to 4% in typical cultivars to > 20% of the seed oil. We performed QTL analysis of a population developed by crossing two mutant lines, one with a missense mutation affecting a stearoyl-acyl-carrier protein desaturase gene resulting in ∼11% seed stearic acid crossed to another mutant, A6, which has 24-28% seed stearic acid. Genotyping-by-sequencing (GBS)-based QTL mapping identified 21 minor and major effect QTL for six seed oil related traits and plant height. The inheritance of a large genomic deletion affecting chromosome 14 is the basis for largest effect QTL, resulting in ∼18% seed stearic acid. This deletion contains SACPD-C and another gene(s); loss of both genes boosts seed stearic acid levels to ≥ 18%. Unfortunately, this genomic deletion has been shown in previous studies to be inextricably correlated with reduced seed yield. Our results will help inform and guide ongoing breeding efforts to improve soybean oil oxidative stability.

Comparison of Genetic Diversity between Chinese and American Soybean (Glycine max (L.)) Accessions Revealed by High-Density SNPs

Soybean is one of the most important economic crops for both China and the United States (US). The exchange of germplasm between these two countries has long been active. In order to investigate genetic relationships between Chinese and US soybean germplasm, 277 Chinese soybean accessions and 300 US soybean accessions from geographically diverse regions were analyzed using 5,361 SNP markers. The genetic diversity and the polymorphism information content (PIC) of the Chinese accessions was higher than that of the US accessions. Population structure analysis, principal component analysis, and cluster analysis all showed that the genetic basis of Chinese soybeans is distinct from that of the USA. The groupings observed in clustering analysis reflected the geographical origins of the accessions; this conclusion was validated with both genetic distance analysis and relative kinship analysis. F_ST-based and EigenGWAS statistical analysis revealed high genetic variation between the two subpopulations. Analysis of the 10 loci with the strongest selection signals showed that many loci were located in chromosome regions that have previously been identified as quantitative trait loci (QTL) associated with environmental-adaptation-related and yield-related traits. The pattern of diversity among the American and Chinese accessions should help breeders to select appropriate parental accessions to enhance the performance of future soybean cultivars.

SNP-SNP Interaction Analysis on Soybean Oil Content under Multi-Environments

Soybean oil content is one of main quality traits. In this study, we used the multifactor dimensionality reduction (MDR) method and a soybean high-density genetic map including 5,308 markers to identify stable single nucleotide polymorphism (SNP)—SNP interactions controlling oil content in soybean across 23 environments. In total, 36,442,756 SNP-SNP interaction pairs were detected, 1865 of all interaction pairs associated with soybean oil content were identified under multiple environments by the Bonferroni correction with p <3.55×10⁻¹¹. Two and 1863 SNP-SNP interaction pairs detected stable across 12 and 11 environments, respectively, which account around 50% of total environments. Epistasis values and contribution rates of stable interaction (the SNP interaction pairs were detected in more than 2 environments) pairs were detected by the two way ANOVA test, the available interaction pairs were ranged 0.01 to 0.89 and from 0.01 to 0.85, respectively. Some of one side of the interaction pairs were identified with previously research as a major QTL without epistasis effects. The results of this study provide insights into the genetic architecture of soybean oil content and can serve as a basis for marker-assisted selection breeding.

Marker-assisted breeding for transgressive seed protein content in soybean [Glycine max (L.) Merr]

KEY MESSAGE

After two cycles of marker-assisted breeding on three loci, lines with transgressive segregation of 8.22-9.32 % protein content were developed based on four original soybean parents with 35.35-44.83 % protein content. Marker-assisted breeding has been an innovative approach in conventional breeding, which is to be further demonstrated, especially for quantitative traits. A study on continuous transgressive breeding for seed protein content (SPC) in soybean using marker-assisted procedures is reported here. The SPC of the recombinant inbred line (RIL) population XG varied in 38.04-47.54 % under five environments with P 1 of 35.35 %, P 2 of 44.34 % and total heritability of 89.11 %. A transgressive segregant XG30 with SPC 45.53 % was selected for further improvement. The linkage mapping of XG showed its genetic constitution composed of five additive QTL (32.16 % of phenotypic variation or PV) and two pairs of epistatic QTL (2.96 % PV) using 400 SSR markers with the remnant heritability 53.99 % attributed to the undetected collective of minor QTL. Another transgressive segregant WT133 with SPC 48.39 % was selected from the RIL population WT (44.83 % SPC for both parents). XG30 and WT133 were genotyped on the three major additive QTL (Prot-08-1, Prot-14-1 and Prot-19-2) as A 2 A 2 B 2 B 2 L 1 L 1 and A 1 A 1 B 1 B 1 L 2 L 2 , respectively. From WT133×XG30, surprising transgressive progenies were obtained, among which the recombinants with all three positive alleles A 2 _B 2 _L 2 _ performed the highest SPC, especially that of Prot-08-1. The five F 2-derived superior families showed their means higher than the high parent value in F 2:3 and F 2:4 and more transgressive effect in F 2:5:6, with the highest as high as 54.15 %, or 4.82 and 9.32 % more than WT133 and its original high parent, respectively. This study demonstrated the efficiency of marker-assisted procedure in breeding for transgressive segregation of quantitative trait.

Genomic resources in mungbean for future breeding programs

Among the legume family, mungbean (Vigna radiata) has become one of the important crops in Asia, showing a steady increase in global production. It provides a good source of protein and contains most notably folate and iron. Beyond the nutritional value of mungbean, certain features make it a well-suited model organism among legume plants because of its small genome size, short life-cycle, self-pollinating, and close genetic relationship to other legumes. In the past, there have been several efforts to develop molecular markers and linkage maps associated with agronomic traits for the genetic improvement of mungbean and, ultimately, breeding for cultivar development to increase the average yields of mungbean. The recent release of a reference genome of the cultivated mungbean (V. radiata var. radiata VC1973A) and an additional de novo sequencing of a wild relative mungbean (V. radiata var. sublobata) has provided a framework for mungbean genetic and genome research, that can further be used for genome-wide association and functional studies to identify genes related to specific agronomic traits. Moreover, the diverse gene pool of wild mungbean comprises valuable genetic resources of beneficial genes that may be helpful in widening the genetic diversity of cultivated mungbean. This review paper covers the research progress on molecular and genomics approaches and the current status of breeding programs that have developed to move toward the ultimate goal of mungbean improvement.

Multi-environment mapping and meta-analysis of 100-seed weight in soybean

100-Seed weight (100-SW) of soybean is an important but complicated quantitative trait to yield. This study was focus on the quantitative trait loci (QTLs) of soybean 100-SW from 2006 to 2010, using recombination inbred lines population that was derived from a cross between Charleston and Dongnong 594. A total of 23 QTLs for 100-SW were detected in the linkage group C2, D1a, F, G and O. Nine QTLs were identified by composite interval mapping including one QTL with the minimum confidence interval (CI) of 1.3 cM, while 14 QTLs by multiple interval mapping. Furthermore, 94 reported QTLs of 100-SW were integrated with our QTL mapping results using BioMercator. As a result, 15 consensus QTLs and their corresponding markers were identified. The minimum CI was reduced to 1.52 cM by the combination of meta-analysis. These findings may merit fine-mapping of these QTL in soybean.