Quantitative Trait Loci for Seed Protein and Oil Concentration, and Seed Size in Soybean

Genome-Wide Association and RNA-Seq Analyses Reveal a Potential Candidate Gene Related to Oil Content in Soybean Seeds

Soybean is a crucial crop globally, serving as a significant source of unsaturated fatty acids and protein in the human diet. However, further enhancements are required for the related genes that regulate soybean oil synthesis. In this study, 155 soybean germplasms were cultivated under three different environmental conditions, followed by phenotypic identification and genome-wide association analysis using simplified sequencing data. Genome-wide association analysis was performed using SLAF-seq data. A total of 36 QTLs were significantly associated with oil content (-log10(p) > 3). Out of the 36 QTLs associated with oil content, 27 exhibited genetic overlap with previously reported QTLs related to oil traits. Further transcriptome sequencing was performed on extreme high-low oil soybean varieties. Combined with transcriptome expression data, 22 candidate genes were identified (|log2FC| ≥ 3). Further haplotype analysis of the potential candidate genes showed that three potential candidate genes had excellent haplotypes, including Glyma.03G186200, Glyma.09G099500, and Glyma.18G248900. The identified loci harboring beneficial alleles and candidate genes likely contribute significantly to the molecular network's underlying marker-assisted selection (MAS) and oil content.

Across-environment seed protein stability and genetic architecture of seed components in soybean

The recent surge in the plant-based protein market has resulted in high demands for soybean genotypes with improved grain yield, seed protein and oil content, and essential amino acids (EAAs). Given the quantitative nature of these traits, complex interactions among seed components, as well as between seed components and environmental factors and management practices, add complexity to the development of desired genotypes. In this study, the across-environment seed protein stability of 449 genetically diverse plant introductions was assessed, revealing that genotypes may display varying sensitivities to such environmental stimuli. The EAAs valine, phenylalanine, and threonine showed the highest variable importance toward the variation in stability, while both seed protein and oil contents were among the explanatory variables with the lowest importance. In addition, 56 single nucleotide polymorphism (SNP) markers were significantly associated with various seed components. Despite the strong phenotypic Pearson's correlation observed among most seed components, many independent genomic regions associated with one or few seed components were identified. These findings provide insights for improving the seed concentration of specific EAAs and reducing the negative correlation between seed protein and oil contents.

Genetic analysis and QTL mapping of seed hardness trait in a soybean (Glycine max) recombinant inbred line (RIL) population

Seed hardness is a critical quality trait impacting both the suitability of soybeans for consumption and their processing. The primary objective of this study was to explore the genetic foundations underlying seed hardness in soybeans. A 234 recombinant inbred line (RIL) population was evaluated for seed hardness across three years (2015 in Gansu, 2016, and 2017 in Hainan). Notably, the parent varieties, Zhonghuang35 and Jindou21, displayed significant differences in seed hardness. Also, the RIL population exhibited a wide range of genetic variation in seed hardness, with coefficients of variation between 70.53 % and 94.94 %. The frequency distribution of this trait conformed to a relatively normal distribution, making it suitable for QTL analysis. Six QTLs associated with seed hardness were identified with three located on chromosome 2 and three on chromosome 16. The major QTL, qHS-2-1, consistently exhibited the highest percentage of PVE and LOD in Gansu 2015, Hainan 2016, and Hainan 2017, suggesting its central role in determining seed hardness. Further investigation revealed four genes within the qHS-2-1 interval potentially related to seed hardness. GO enrichment analysis provided insights into their functions, including factors such as Glyma.02G307000, a pectin lyase-like superfamily protein, which could influence seed hardness through its role in pectin lyase enzyme activity. Expression analysis of these candidate genes demonstrated significant differences between the two parent varieties, further highlighting their potential role in seed coat hardness. This study offers valuable insights into the genetic mechanisms governing soybean seed coat hardness, providing a foundation for future research and crop improvement efforts.

Identification of the domestication gene GmCYP82C4 underlying the major quantitative trait locus for the seed weight in soybean

KEY MESSAGE

A major quantitative trait locus (QTL) for the hundred-seed weight (HSW) was identified and confirmed in the two distinct soybean populations, and the target gene GmCYP82C4 underlying this locus was identified that significantly associated with soybean seed weight, and it was selected during the soybean domestication and improvement process. Soybean is a major oil crop for human beings and the seed weight is a crucial goal of soybean breeding. However, only a limited number of target genes underlying the quantitative trait loci (QTLs) controlling seed weight in soybean are known so far. In the present study, six loci associated with hundred-seed weight (HSW) were detected in the first population of 573 soybean breeding lines by genome-wide association study (GWAS), and 64 gene models were predicted in these candidate QTL regions. The QTL qHSW_1 exhibits continuous association signals on chromosome four and was also validated by region association study (RAS) in the second soybean population (409 accessions) with wild, landrace, and cultivar soybean accessions. There were seven genes in qHSW_1 candidate region by linkage disequilibrium (LD) block analysis, and only Glyma.04G035500 (GmCYP82C4) showed specifically higher expression in flowers, pods, and seeds, indicating its crucial role in the soybean seed development. Significant differences in HSW trait were detected when the association panels are genotyped by single-nucleotide polymorphisms (SNPs) in putative GmCYP82C4 promoter region. Eight haplotypes were generated by six SNPs in GmCYP82C4 in the second soybean population, and two superior haplotypes (Hap2 and Hap4) of GmCYP82C4 were detected with average HSW of 18.27 g and 18.38 g, respectively. The genetic diversity of GmCYP82C4 was analyzed in the second soybean population, and GmCYP82C4 was most likely selected during the soybean domestication and improvement process, leading to the highest proportion of Hap2 of GmCYP82C4 both in landrace and cultivar subpopulations. The QTLs and GmCYP82C4 identified in this study provide novel genetic resources for soybean seed weight trait, and the GmCYP82C4 could be used for soybean molecular breeding to develop desirable seed weight in the future.

Soybean genetics, genomics, and breeding for improving nutritional value and reducing antinutritional traits in food and feed

Soybean [Glycine max (L.) Merr.] is a globally important crop due to its valuable seed composition, versatile feed, food, and industrial end-uses, and consistent genetic gain. Successful genetic gain in soybean has led to widespread adaptation and increased value for producers, processors, and consumers. Specific focus on the nutritional quality of soybean seed composition for food and feed has further elucidated genetic knowledge and bolstered breeding progress. Seed components are historical and current targets for soybean breeders seeking to improve nutritional quality of soybean. This article reviews genetic and genomic foundations for improvement of nutritionally important traits, such as protein and amino acids, oil and fatty acids, carbohydrates, and specific food-grade considerations; discusses the application of advanced breeding technology such as CRISPR/Cas9 in creating seed composition variations; and provides future directions and breeding recommendations regarding soybean seed composition traits.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

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.

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 major QTLs for soybean seed size and seed weight traits using a RIL population in different environments

INTRODUCTION

The seed weight of soybean [Glycine max (L.) Merr.] is one of the major traits that determine soybean yield and is closely related to seed size. However, the genetic basis of the synergistic regulation of traits related to soybean yield is unclear.

METHODS

To understand the molecular genetic basis for the formation of soybean yield traits, the present study focused on QTLs mapping for seed size and weight traits in different environments and target genes mining.

RESULTS

A total of 85 QTLs associated with seed size and weight traits were identified using a recombinant inbred line (RIL) population developed from Guizao1×B13 (GB13). We also detected 18 environmentally stable QTLs. Of these, qSL-3-1 was a novel QTL with a stable main effect associated with seed length. It was detected in all environments, three of which explained more than 10% of phenotypic variance (PV), with a maximum of 15.91%. In addition, qSW-20-3 was a novel QTL with a stable main effect associated with seed width, which was identified in four environments. And the amount of phenotypic variance explained (PVE) varied from 9.22 to 21.93%. Five QTL clusters associated with both seed size and seed weight were summarized by QTL cluster identification. Fifteen candidate genes that may be involved in regulating soybean seed size and weight were also screened based on gene function annotation and GO enrichment analysis.

DISCUSSION

The results provide a biologically basic reference for understanding the formation of soybean seed size and weight traits.

Identification of major quantitative trait loci and candidate genes for seed weight in soybean

Four major quantitative trait loci for 100-seed weight were identified in a soybean RIL population under five environments, and the most likely candidate genes underlying these loci were identified. Seed weight is an important target of soybean breeding. However, the genes underlying the major quantitative trait loci (QTL) controlling seed weight remain largely unknown. In this study, a soybean population of 300 recombinant inbred lines (RILs) derived from a cross between PI595843 (PI) and WH was used to map the QTL and identify candidate genes for seed weight. The RIL population was genotyped through whole genome resequencing, and phenotyped for 100-seed weight under five environments. A total of 38 QTL were detected, and four major QTL, each explained at least 10% of the variation in 100-seed weight, were identified. Six candidate genes within these four major QTL regions were identified by analyses of their tissue expression patterns, gene annotations, and differential gene expression levels in soybean seeds during four developmental stages between two parental lines. Further sequence variation analyses revealed a C to T substitution in the first exon of the Glyma.19G143300, resulting in an amino acid change between PI and WH, and thus leading to a different predicted kinase domain, which might affect its protein function. Glyma.19G143300 is highly expressed in soybean seeds and encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). Its predicted protein has typical domains of LRR-RLK family, and phylogenetic analyses reveled its similarity with the known LRR-RLK protein XIAO (LOC_Os04g48760), which is involved in controlling seed size. The major QTL and candidate genes identified in this study provide useful information for molecular breeding of new soybean cultivars with desirable seed weight.

Identification of Drought-Tolerance Genes in the Germination Stage of Soybean

Drought stress influences the vigor of plant seeds and inhibits seed germination, making it one of the primary environmental factors adversely affecting food security. The seed germination stage is critical to ensuring the growth and productivity of soybeans in soils prone to drought conditions. We here examined the genetic diversity and drought-tolerance phenotypes of 410 accessions of a germplasm diversity panel for soybean and conducted quantitative genetics analyses to identify loci associated with drought tolerance of seed germination. We uncovered significant differences among the diverse genotypes for four growth indices and five drought-tolerance indices, which revealed abundant variation among genotypes, upon drought stress, and for genotype × treatment effects. We also used 158,327 SNP markers and performed GWAS for the drought-related traits. Our data met the conditions (PCA + K) for using a mixed linear model in TASSEL, and we thus identified 26 SNPs associated with drought tolerance indices for germination stage distributed across 10 chromosomes. Nine SNP sites, including, for example, Gm20_34956219 and Gm20_36902659, were associated with two or more phenotypic indices, and there were nine SNP markers located in or adjacent to (within 500 kb) previously reported drought tolerance QTLs. These SNPs led to our identification of 41 candidate genes related to drought tolerance in the germination stage. The results of our study contribute to a deeper understanding of the genetic mechanisms underlying drought tolerance in soybeans at the germination stage, thereby providing a molecular basis for identifying useful soybean germplasm for breeding new drought-tolerant varieties.

A multiple phenotype imputation method for genetic diversity and core collection in Taiwanese vegetable soybean

Establishment of vegetable soybean (edamame) [Glycine max (L.) Merr.] germplasms has been highly valued in Asia and the United States owing to the increasing market demand for edamame. The idea of core collection (CC) is to shorten the breeding program so as to improve the availability of germplasm resources. However, multidimensional phenotypes typically are highly correlated and have different levels of missing rate, often failing to capture the underlying pattern of germplasms and select CC precisely. These are commonly observed on correlated samples. To overcome such scenario, we introduced the "multiple imputation" (MI) method to iteratively impute missing phenotypes for 46 morphological traits and jointly analyzed high-dimensional imputed missing phenotypes (EC impu ) to explore population structure and relatedness among 200 Taiwanese vegetable soybean accessions. An advanced maximization strategy with a heuristic algorithm and PowerCore was used to evaluate the morphological diversity among the EC impu . In total, 36 accessions (denoted as CC impu ) were efficiently selected representing high diversity and the entire coverage of the EC impu . Only 4 (8.7%) traits showed slightly significant differences between the CC impu and EC impu . Compared to the EC impu , 96% traits retained all characteristics or had a slight diversity loss in the CC impu . The CC impu exhibited a small percentage of significant mean difference (4.51%), and large coincidence rate (98.1%), variable rate (138.76%), and coverage (close to 100%), indicating the representativeness of the EC impu . We noted that the CC impu outperformed the CC raw in evaluation properties, suggesting that the multiple phenotype imputation method has the potential to deal with missing phenotypes in correlated samples efficiently and reliably without re-phenotyping accessions. Our results illustrated a significant role of imputed missing phenotypes in support of the MI-based framework for plant-breeding programs.

Linkage analysis and residual heterozygotes derived near isogenic lines reveals a novel protein quantitative trait loci from a Glycine soja accession

Modern soybean [Glycine max (L.) Merr] cultivars have low overall genetic variation due to repeated bottleneck events that arose during domestication and from selection strategies typical of many soybean breeding programs. In both public and private soybean breeding programs, the introgression of wild soybean (Glycine soja Siebold and Zucc.) alleles is a viable option to increase genetic diversity and identify new sources for traits of value. The objectives of our study were to examine the genetic architecture responsible for seed protein and oil using a recombinant inbred line (RIL) population derived from hybridizing a G. max line ('Osage') with a G. soja accession (PI 593983). Linkage mapping identified a total of seven significant quantitative trait loci on chromosomes 14 and 20 for seed protein and on chromosome 8 for seed oil with LOD scores ranging from 5.3 to 31.7 for seed protein content and from 9.8 to 25.9 for seed oil content. We analyzed 3,015 single F4:9 soybean plants to develop two residual heterozygotes derived near isogenic lines (RHD-NIL) populations by targeting nine SNP markers from genotype-by-sequencing, which corresponded to two novel quantitative trait loci (QTL) derived from G. soja: one for a novel seed oil QTL on chromosome 8 and another for a novel protein QTL on chromosome 14. Single marker analysis and linkage analysis using 50 RHD-NILs validated the chromosome 14 protein QTL, and whole genome sequencing of RHD-NILs allowed us to reduce the QTL interval from ∼16.5 to ∼4.6 Mbp. We identified two genomic regions based on recombination events which had significant increases of 0.65 and 0.72% in seed protein content without a significant decrease in seed oil content. A new Kompetitive allele-specific polymerase chain reaction (KASP) assay, which will be useful for introgression of this trait into modern elite G. max cultivars, was developed in one region. Within the significantly associated genomic regions, a total of eight genes are considered as candidate genes, based on the presence of gene annotations associated with the protein or amino acid metabolism/movement. Our results provide better insights into utilizing wild soybean as a source of genetic diversity for soybean cultivar improvement utilizing native traits.

Ensuring Global Food Security by Improving Protein Content in Major Grain Legumes Using Breeding and 'Omics' Tools

Grain legumes are a rich source of dietary protein for millions of people globally and thus a key driver for securing global food security. Legume plant-based 'dietary protein' biofortification is an economic strategy for alleviating the menace of rising malnutrition-related problems and hidden hunger. Malnutrition from protein deficiency is predominant in human populations with an insufficient daily intake of animal protein/dietary protein due to economic limitations, especially in developing countries. Therefore, enhancing grain legume protein content will help eradicate protein-related malnutrition problems in low-income and underprivileged countries. Here, we review the exploitable genetic variability for grain protein content in various major grain legumes for improving the protein content of high-yielding, low-protein genotypes. We highlight classical genetics-based inheritance of protein content in various legumes and discuss advances in molecular marker technology that have enabled us to underpin various quantitative trait loci controlling seed protein content (SPC) in biparental-based mapping populations and genome-wide association studies. We also review the progress of functional genomics in deciphering the underlying candidate gene(s) controlling SPC in various grain legumes and the role of proteomics and metabolomics in shedding light on the accumulation of various novel proteins and metabolites in high-protein legume genotypes. Lastly, we detail the scope of genomic selection, high-throughput phenotyping, emerging genome editing tools, and speed breeding protocols for enhancing SPC in grain legumes to achieve legume-based dietary protein security and thus reduce the global hunger risk.

Identification of major genomic regions for soybean seed weight by genome-wide association study

The hundred-seed weight (HSW) is an important yield component and one of the principal breeding traits in soybean. More than 250 quantitative trait loci (QTL) for soybean HSW have been identified. However, most of them have a large genomic region or are environmentally sensitive, which provide limited information for improving the phenotype in marker-assisted selection (MAS) and identifying the candidate genes. Here, we utilized 281 soybean accessions with 58,112 single nucleotide polymorphisms (SNPs) to dissect the genetic basis of HSW in across years in the northern Shaanxi province of China through one single-locus (SL) and three multi-locus (ML) genome-wide association study (GWAS) models. As a result, one hundred and fifty-four SNPs were detected to be significantly associated with HSW in at least one environment via SL-GWAS model, and 27 of these 154 SNPs were detected in all (three) environments and located within 7 linkage disequilibrium (LD) block regions with the distance of each block ranged from 40 to 610 Kb. A total of 15 quantitative trait nucleotides (QTNs) were identified by three ML-GWAS models. Combined with the results of different GWAS models, the 7 LD block regions associated with HSW detected by SL-GWAS model could be verified directly or indirectly by the results of ML-GWAS models. Eleven candidate genes underlying the stable loci that may regulate seed weight in soybean were predicted. The significantly associated SNPs and the stable loci as well as predicted candidate genes may be of great importance for marker-assisted breeding, polymerization breeding, and gene discovery for HSW in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01310-y.

Nutritional Diversity in Native Germplasm of Maize Collected From Three Different Fragile Ecosystems of India

Native germplasm resources are adapted to specific ecological niches. They have sustained over generations owing to the preference of local communities for their unique taste, the utility to particular dishes, and the low cost of cultivation. They may help eradicate malnutrition and act as a source for trait-linked genes. The present dataset comprises thirty-three native germplasm of maize collected from Rajasthan, Himachal Pradesh, and Andhra Pradesh states of India with an altitudinal variation of 386–2,028 m. They were evaluated for proximate composition, minerals, nutritional attributes, and antioxidant activity and compared with the standard values reported in the Indian Food Composition Table 2017 (IFCT2017). The nutritional profile showed moisture content in the range of 7.16–10.9%, ash 0.73–1.93%, crude protein 8.68–12.0%, crude fat 3.72–8.03%, dietary fiber 5.21–11.2%, and available carbohydrates 60.6–69.8%. Three accessions, namely, Malan 11 (7.06%), Malan 24 (7.20%), and Yellow Chamba Local 02 (8.03%) exhibited almost double the crude fat content as compared with the values notified in IFCT2017 (3.77). Total sugar content obtained was in the range of 5.00–11.3%, whereas the starch content was found between 50.9 and 64.9%. All the germplasm except Yellow Chamba Local reflected a higher protein content than reported values in IFCT2017 (8.80). Sathi, Safed Chamba Local, and Ragal Makka had nearly 12% protein content. Mineral malnutrition, mainly due to iron (Fe) deficiency, is a worldwide issue to science, humanity, and society. The mineral profile revealed that most germplasm had a higher iron content. Accessions with the iron content of nearly three times of IFCT2017 reported value were identified in germplasm belonging to three states. A negative relationship was observed between the altitude of the sample collection site and available carbohydrate content. In contrast, available carbohydrate showed inverse correlations with dietary fiber, protein, and fat content. The information generated in this study can be utilized to promote these germplasm as nutrifood, nutritional surveillance, labeling, and crop improvement programs.

Genome-Wide Association Study and Genomic Selection for Proteinogenic Methionine in Soybean Seeds

Soybean [Glycine max (L.) Merr.] seeds have an amino acid profile that provides excellent viability as a food and feed protein source. However, low concentrations of an essential amino acid, methionine, limit the nutritional utility of soybean protein. The objectives of this study were to identify genomic associations and evaluate the potential for genomic selection (GS) for methionine content in soybean seeds. We performed a genome-wide association study (GWAS) that utilized 311 soybean accessions from maturity groups IV and V grown in three locations in 2018 and 2019. A total of 35,570 single nucleotide polymorphisms (SNPs) were used to identify genomic associations with proteinogenic methionine content that was quantified by high-performance liquid chromatography (HPLC). Across four environments, 23 novel SNPs were identified as being associated with methionine content. The strongest associations were found on chromosomes 3 (ss715586112, ss715586120, ss715586126, ss715586203, and ss715586204), 8 (ss715599541 and ss715599547) and 16 (ss715625009). Several gene models were recognized within proximity to these SNPs, such as a leucine-rich repeat protein kinase and a serine/threonine protein kinase. Identification of these linked SNPs should help soybean breeders to improve protein quality in soybean seeds. GS was evaluated using k-fold cross validation within each environment with two SNP sets, the complete 35,570 set and a subset of 248 SNPs determined to be associated with methionine through GWAS. Average prediction accuracy (r ²) was highest using the SNP subset ranging from 0.45 to 0.62, which was a significant improvement from the complete set accuracy that ranged from 0.03 to 0.27. This indicated that GS utilizing a significant subset of SNPs may be a viable tool for soybean breeders seeking to improve methionine content.

Identification and Validation of Major QTLs, Epistatic Interactions, and Candidate Genes for Soybean Seed Shape and Weight Using Two Related RIL Populations

Understanding the genetic mechanism underlying seed size, shape, and weight is essential for enhancing soybean cultivars. High-density genetic maps of two recombinant inbred line (RIL) populations, LM6 and ZM6, were evaluated across multiple environments to identify and validate M-QTLs as well as identify candidate genes behind major and stable quantitative trait loci (QTLs). A total of 239 and 43 M-QTLs were mapped by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM) approaches, from which 180 and 18, respectively, are novel QTLs. Twenty-two QTLs including four novel major QTLs were validated in the two RIL populations across multiple environments. Moreover, 18 QTLs showed significant AE effects, and 40 pairwise of the identified QTLs exhibited digenic epistatic effects. Thirty-four QTLs associated with seed flatness index (FI) were identified and reported here for the first time. Seven QTL clusters comprising several QTLs for seed size, shape, and weight on genomic regions of chromosomes 3, 4, 5, 7, 9, 17, and 19 were identified. Gene annotations, gene ontology (GO) enrichment, and RNA-seq analyses of the genomic regions of those seven QTL clusters identified 47 candidate genes for seed-related traits. These genes are highly expressed in seed-related tissues and nodules, which might be deemed as potential candidate genes regulating the seed size, weight, and shape traits in soybean. This study provides detailed information on the genetic basis of the studied traits and candidate genes that could be efficiently implemented by soybean breeders for fine mapping and gene cloning, and for marker-assisted selection (MAS) targeted at improving these traits individually or concurrently.

Fine mapping QTL and mining genes for protein content in soybean by the combination of linkage and association analysis

Soybean is one main source of dietary protein; therefore, improving protein content is an important objective in breeding programs. There is a significant negative correlation between protein and oil content, which influenced mapping quantitative trait locus (QTL) and quantitative trait nucleotides for these two traits. In this study, a linkage map was created with 2232 single-nucleotide polymorphism markers for the four-way recombinant inbred line (FW-RIL) population derived from the cross (Kenfeng 14 × Kenfeng 15) × (Heinong 48 × Kenfeng 19), and then conditional and unconditional QTL analyses were carried out by inclusive complete interval mapping based on the phenotypic data of protein and oil content collected in 10 different environments. As shown in the results of linkage analysis, a total of 85 QTL have been detected. We have performed association analysis using 109,676 markers after quality filtering for FW-RIL, and the results have shown that a total of 60 QTNs were detected. We have performed association analysis using 63,306 markers after quality filtering for resource population, and the results have shown that a total of 123 QTNs were detected. We have combined linkage and association analysis, and there are six QTNs verified by FW-RIL and resource population. We have performed pathway analysis on the genes in these six QTN attenuation regions, and the result shows that a total of four candidate genes are related to the synthesis or metabolism of soybean protein. These findings will facilitate marker-assisted selection and molecular breeding of soybean.

Genome-wide association study identified candidate genes for seed size and seed composition improvement in M. truncatula

Grain legumes are highly valuable plant species, as they produce seeds with high protein content. Increasing seed protein production and improving seed nutritional quality represent an agronomical challenge in order to promote plant protein consumption of a growing population. In this study, we used the genetic diversity, naturally present in Medicago truncatula, a model plant for legumes, to identify genes/loci regulating seed traits. Indeed, using sequencing data of 162 accessions from the Medicago HAPMAP collection, we performed genome-wide association study for 32 seed traits related to seed size and seed composition such as seed protein content/concentration, sulfur content/concentration. Using different GWAS and postGWAS methods, we identified 79 quantitative trait nucleotides (QTNs) as regulating seed size, 41 QTNs for seed composition related to nitrogen (i.e. storage protein) and sulfur (i.e. sulfur-containing amino acid) concentrations/contents. Furthermore, a strong positive correlation between seed size and protein content was revealed within the selected Medicago HAPMAP collection. In addition, several QTNs showed highly significant associations in different seed phenotypes for further functional validation studies, including one near an RNA-Binding Domain protein, which represents a valuable candidate as central regulator determining both seed size and composition. Finally, our findings in M. truncatula represent valuable resources to be exploitable in many legume crop species such as pea, common bean, and soybean due to its high synteny, which enable rapid transfer of these results into breeding programs and eventually help the improvement of legume grain production.

Effect of Planting Date and Cultivar Maturity in Edamame Quality and Harvest Window

Edamame is a food-grade soybean [Glycine max (L.) Merr.] that is harvested immature between the R6 and R7 reproductive stages. To be labeled as a premium product, the edamame market demands large pod size and intense green color. A staggered harvest season is critical for the commercial industry to post-harvest process the crop in a timely manner. Currently, there is little information to assist in predicting the optimum time to harvest edamame when the pods are at their collective largest size and greenest color. The objectives of this study were to assess the impact of cultivar, planting date, and harvest date on edamame color, pod weight, and a newly minted Edamame Harvest Quality Index combining both aforementioned factors. And to predict edamame harvest quality based on phenological stages, thermal units, and planting dates. We observed that pod color and weight depended on the cultivar, planting date, and harvest date combination. Our results also indicated that edamame quality is increased with delayed planting dates and that quality was dependent on harvest date with a quadratic negative response to delaying harvest. Maximum quality depended on cultivar and planting and harvest dates, but it remained stable for an interval of 18-27 days around the peak. Finally, we observed that the number of days between R1 and harvest was consistently identified as a key factor driving edamame quality by both stepwise regression and neural network analysis. These research results will help define a planting and harvest strategy for edamame production in Arkansas and the United States Mid-South.

GmKIX8-1 regulates organ size in soybean and is the causative gene for the major seed weight QTL qSw17-1

Seed weight is one of the most important agronomic traits in soybean for yield improvement and food production. Several quantitative trait loci (QTLs) associated with the trait have been identified in soybean. However, the genes underlying the QTLs and their functions remain largely unknown. Using forward genetic methods and CRISPR/Cas9 gene editing, we identified and characterized the role of GmKIX8-1 in the control of organ size in soybean. GmKIX8-1 belongs to a family of KIX domain-containing proteins that negatively regulate cell proliferation in plants. Consistent with this predicted function, we found that loss-of-function GmKIX8-1 mutants showed a significant increase in the size of aerial plant organs, such as seeds and leaves. Likewise, the increase in organ size is due to increased cell proliferation, rather than cell expansion, and increased expression of CYCLIN D3;1-10. Lastly, molecular analysis of soybean germplasms harboring the qSw17-1 QTL for the big-seeded phenotype indicated that reduced expression of GmKIX8-1 is the genetic basis of the qSw17-1 phenotype.

Genomic regions associated with important seed quality traits in food-grade soybeans

BACKGROUND

The production of soy-based food products requires specific physical and chemical characteristics of the soybean seed. Identification of quantitative trait loci (QTL) associated with value-added traits, such as seed weight, seed protein and sucrose concentration, could accelerate the development of competitive high-protein soybean cultivars for the food-grade market through marker-assisted selection (MAS). The objectives of this study were to identify and validate QTL associated with these value-added traits in two high-protein recombinant inbred line (RIL) populations.

RESULTS

The RIL populations were derived from the high-protein cultivar 'AC X790P' (49% protein, dry weight basis), and two high-yielding commercial cultivars, 'S18-R6' (41% protein) and 'S23-T5' (42% protein). Fourteen large-effect QTL (R2 > 10%) were identified associated with seed protein concentration. Of these QTL, seven QTL were detected in both populations, and eight of them were co-localized with QTL associated with either seed sucrose concentration or seed weight. None of the protein-related QTL was found to be associated with seed yield in either population. Sixteen candidate genes with putative roles in protein metabolism were identified within seven of these protein-related regions: qPro_Gm02-3, qPro_Gm04-4, qPro_Gm06-1, qPro_Gm06-3, qPro_Gm06-6, qPro_Gm13-4 and qPro-Gm15-3.

CONCLUSION

The use of RIL populations derived from high-protein parents created an opportunity to identify four novel QTL that may have been masked by large-effect QTL segregating in populations developed from diverse parental cultivars. In total, we have identified nine protein QTL that were detected either in both populations in the current study or reported in other studies. These QTL may be useful in the curated selection of new soybean cultivars for optimized soy-based food products.

QTL mapping for soybean (Glycine max L.) leaf chlorophyll-content traits in a genotyped RIL population by using RAD-seq based high-density linkage map

BACKGROUND

Different soybean (Glycine max L.) leaf chlorophyll-content traits are considered to be significantly linked to soybean yield. To map the quantitative trait loci (QTLs) of soybean leaf chlorophyll-content traits, an advanced recombinant inbred line (RIL, ZH, Zhonghuang 24 × Huaxia 3) population was adopted to phenotypic data acquisitions for the target traits across six distinct environments (seasons and soybean growth stages). Moreover, the restriction site-associated DNA sequencing (RAD-seq) based high-density genetic linkage map of the RIL population was utilized for QTL mapping by carrying out the composite interval mapping (CIM) approach.

RESULTS

Correlation analyses showed that most traits were correlated with each other under specific chlorophyll assessing method and were regulated both by hereditary and environmental factors. In this study, 78 QTLs for soybean leaf chlorophyll-content traits were identified. Furthermore, 13 major QTLs and five important QTL hotspots were classified and highlighted from the detected QTLs. Finally, Glyma01g15506, Glyma02g08910, Glyma02g11110, Glyma07g15960, Glyma15g19670 and Glyma15g19810 were predicted from the genetic intervals of the major QTLs and important QTL hotspots.

CONCLUSIONS

The detected QTLs and candidate genes may facilitate to gain a better understanding of the hereditary basis of soybean leaf chlorophyll-content traits and may be valuable to pave the way for the marker-assisted selection (MAS) breeding of the target traits.

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 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.

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.

A Combined Linkage and GWAS Analysis Identifies QTLs Linked to Soybean Seed Protein and Oil Content

Soybean is an excellent source of vegetable protein and edible oil. Understanding the genetic basis of protein and oil content will improve the breeding programs for soybean. Linkage analysis and genome-wide association study (GWAS) tools were combined to detect quantitative trait loci (QTL) that are associated with protein and oil content in soybean. Three hundred and eight recombinant inbred lines (RILs) containing 3454 single nucleotide polymorphism (SNP) markers and 200 soybean accessions, including 94,462 SNPs and indels, were applied to identify QTL intervals and significant SNP loci. Intervals on chromosomes 1, 15, and 20 were correlated with both traits, and QTL qPro15-1, qPro20-1, and qOil5-1 reproducibly correlated with large phenotypic variations. SNP loci on chromosome 20 that overlapped with qPro20-1 were reproducibly connected to both traits by GWAS (p < 10-4). Twenty-five candidate genes with putative roles in protein and/or oil metabolisms within two regions (qPro15-1, qPro20-1) were identified, and eight of these genes showed differential expressions in parent lines during late reproductive growth stages, consistent with a role in controlling protein and oil content. The new well-defined QTL should significantly improve molecular breeding programs, and the identified candidate genes may help elucidate the mechanisms of protein and oil biosynthesis.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.

Genome wide association mapping and candidate gene analysis for hundred seed weight in soybean [Glycine max (L.) Merrill]

BACKGROUND

The hundred seed weight (HSW) is one of the yield components of soybean [Glycine max (L.) Merrill] and is especially critical for various soybean food types. In this study, a representative sample consisting of 185 accessions was selected from Northeast China and analysed in three tested environments to determine the quantitative trait nucleotide (QTN) of HSW through a genome-wide association study (GWAS).

RESULT

A total of 24,180 single nucleotide polymorphisms (SNPs) with minor allele frequencies greater than 0.2 and missing data less than 3% were utilized to estimate linkage disequilibrium (LD) levels in the tested association panel. Thirty-four association signals were identified as associated with HSW via GWAS. Among them, nineteen QTNs were novel, and another fifteen QTNs were overlapped or located near the genomic regions of known HSW QTL. A total of 237 genes, derived from 31 QTNs and located near peak SNPs from the three tested environments in 2015 and 2016, were considered candidate genes, were related to plant growth regulation, hormone metabolism, cell, RNA, protein metabolism, development, starch accumulation, secondary metabolism, signalling, and the TCA cycle, some of which have been found to participate in the regulation of HSW. A total of 106 SNPs from 16 candidate genes were significantly associated with HSW in soybean.

CONCLUSIONS

The identified loci with beneficial alleles and candidate genes might be valuable for the molecular network and MAS of HSW.

Genome-wide association study of inflorescence length of cultivated soybean based on the high-throughout single-nucleotide markers

As an important and complex trait, inflorescence length (IL) of soybean [Glycine max (L.) Merr.] significantly affected seed yields. Therefore, elucidating molecular basis of inflorescence architecture, especially for IL, was important for improving soybean yield potentials. Longer IL meaned to have more pod and seed in soybean. Hence, increasing IL and improving yield are targets for soybean breeding. In this study, a association panel, comprising 283 diverse samples, was used to dissect the genetic basis of IL based on genome-wide association analysis (GWAS) and haplotype analysis. GWAS and haplotype analysis were conducted through high-throughout single-nucleotide polymorphisms (SNP) developed by SLAF-seq methodology. A total of 39, 057 SNPs (minor allele frequency ≥ 0.2 and missing data ≤ 10%) were utilized to evaluate linkage disequilibrium (LD) level in the tested association panel. A total of 30 association signals were identified to be associated with IL via GWAS. Among them, 13 SNPs were novel, and another 17 SNPs were overlapped or located near the linked regions of known quantitative trait nucleotide (QTN) with soybean seed yield or yield component. The functional genes, located in the 200-kb genomic region of each peak SNP, were considered as candidate genes, such as the cell division/ elongation, specific enzymes, and signaling or transport of specific proteins. These genes have been reported to participant in the regulation of IL. Ten typical long-IL lines and ten typical short-IL lines were re-sequencing, and then, six SNPs from five genes were obtained based on candidate gene-based association. In addition, 42 haplotypes were defined based on haplotype analysis. Of them, 11 haplotypes were found to regulate long IL (> 14 mm) in soybean. The identified 30 QTN with beneficial alleles and their candidate genes might be valuable for dissecting the molecular mechanisms of IL and further improving the yield potential of soybean.

Genome-Wide Detection of Major and Epistatic Effect QTLs for Seed Protein and Oil Content in Soybean Under Multiple Environments Using High-Density Bin Map

Seed protein and oil content are the two important traits determining the quality and value of soybean. Development of improved cultivars requires detailed understanding of the genetic basis underlying the trait of interest. However, it is prerequisite to have a high-density linkage map for precisely mapping genomic regions, and therefore the present study used high-density genetic map containing 2267 recombination bin markers distributed on 20 chromosomes and spanned 2453.79 cM with an average distance of 1.08 cM between markers using restriction-site-associated DNA sequencing (RAD-seq) approach. A recombinant inbred line (RIL) population of 104 lines derived from a cross between Linhefenqingdou and Meng 8206 cultivars was evaluated in six different environments to identify main- and epistatic-effect quantitative trait loci (QTLs)as well as their interaction with environments. A total of 44 main-effect QTLs for protein and oil content were found to be distributed on 17 chromosomes, and 15 novel QTL were identified for the first time. Out of these QTLs, four were major and stable QTLs, viz., qPro-7-1, qOil-8-3, qOil-10-2 and qOil-10-4, detected in at least two environments plus combined environment with R² values >10%. Within the physical intervals of these four QTLs, 111 candidate genes were screened for their direct or indirect involvement in seed protein and oil biosynthesis/metabolism processes based on gene ontology and annotation information. Based on RNA sequencing (RNA-seq) data analysis, 15 of the 111 genes were highly expressed during seed development stage and root nodules that might be considered as the potential candidate genes. Seven QTLs associated with protein and oil content exhibited significant additive and additive × environment interaction effects, and environment-independent QTLs revealed higher additive effects. Moreover, three digenic epistatic QTLs pairs were identified, and no main-effect QTLs showed epistasis. In conclusion, the use of a high-density map identified closely linked flanking markers, provided better understanding of genetic architecture and candidate gene information, and revealed the scope available for improvement of soybean quality through marker assisted selection (MAS).

Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping

The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high-resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K-SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K-SNP markers. Importantly, the bin map derived from next-generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication-related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome-wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL-based haplotype using whole-genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

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.

Genotyping of Soybean Cultivars With Medium-Density Array Reveals the Population Structure and QTNs Underlying Maturity and Seed Traits

Soybean was domesticated about 5,000 to 6,000 years ago in China. Although genotyping technologies such as genotyping by sequencing (GBS) and high-density array are available, it is convenient and economical to genotype cultivars or populations using medium-density SNP array in genetic study as well as in molecular breeding. In this study, 235 cultivars, collected from China, Japan, USA, Canada and some other countries, were genotyped using SoySNP8k iSelect BeadChip with 7,189 single nucleotide polymorphisms (SNPs). In total, 4,471 polymorphic SNP markers were used to analyze population structure and perform genome-wide association study (GWAS). The most likely K value was 7, indicating this population can be divided into 7 subpopulations, which is well in accordance with the geographic origins of cultivars or accession studied. The LD decay rate was estimated at 184 kb, where r2 dropped to half of its maximum value (0.205). GWAS using FarmCPU detected a stable quantitative trait nucleotide (QTN) for hilum color and seed color, which is consistent with the known loci or genes. Although no universal QTNs for flowering time and maturity were identified across all environments, a total of 30 consistent QTNs were detected for flowering time (R1) or maturity (R7 and R8) on 16 chromosomes, most of them were corresponding to known E1 to E4 genes or QTL region reported in SoyBase (soybase.org). Of 16 consistent QTNs for protein and oil contents, 11 QTNs were detected having antagonistic effects on protein and oil content, while 4 QTNs soly for oil content, and one QTN soly for protein content. The information gained in this study demonstrated that the usefulness of the medium-density SNP array in genotyping for genetic study and molecular breeding.

Quantitative trait loci associated with soybean seed weight and composition under different phosphorus levels

Seed size and composition are important traits in food crops and can be affected by nutrient availability in the soil. Phosphorus (P) is a non-renewable, essential macronutrient, and P deficiency limits soybean (Glycine max) yield and quality. To investigate the associations of seed traits in low- and high-P environments, soybean recombinant inbred lines (RILs) from a cross of cultivars Fiskeby III and Mandarin (Ottawa) were grown under contrasting P availability environments. Traits including individual seed weight, seed number, and intact mature pod weight were significantly affected by soil P levels and showed transgressive segregation among the RILs. Surprisingly, P treatments did not affect seed composition or weight, suggesting that soybean maintains sufficient P in seeds even in low-P soil. Quantitative trait loci (QTLs) were detected for seed weight, intact pods, seed volume, and seed protein, with five significant QTLs identified in low-P environments and one significant QTL found in the optimal-P environment. Broad-sense heritability estimates were 0.78 (individual seed weight), 0.90 (seed protein), 0.34 (seed oil), and 0.98 (seed number). The QTLs identified under low P point to genetic regions that may be useful to improve soybean performance under limiting P conditions.

Molecular mapping and genomics of soybean seed protein: a review and perspective for the future

KEY MESSAGE

Genetic improvement of soybean protein meal is a complex process because of negative correlation with oil, yield, and temperature. This review describes the progress in mapping and genomics, identifies knowledge gaps, and highlights the need of integrated approaches. Meal protein derived from soybean [Glycine max (L) Merr.] seed is the primary source of protein in poultry and livestock feed. Protein is a key factor that determines the nutritional and economical value of soybean. Genetic improvement of soybean seed protein content is highly desirable, and major quantitative trait loci (QTL) for soybean protein have been detected and repeatedly mapped on chromosomes (Chr.) 20 (LG-I), and 15 (LG-E). However, practical breeding progress is challenging because of seed protein content's negative genetic correlation with seed yield, other seed components such as oil and sucrose, and interaction with environmental effects such as temperature during seed development. In this review, we discuss rate-limiting factors related to soybean protein content and nutritional quality, and potential control factors regulating seed storage protein. In addition, we describe advances in next-generation sequencing technologies for precise detection of natural variants and their integration with conventional and high-throughput genotyping technologies. A syntenic analysis of QTL on Chr. 15 and 20 was performed. Finally, we discuss comprehensive approaches for integrating protein and amino acid QTL, genome-wide association studies, whole-genome resequencing, and transcriptome data to accelerate identification of genomic hot spots for allele introgression and soybean meal protein improvement.

Identification of quantitative trait loci underlying seed protein content of soybean including main, epistatic, and QTL × environment effects in different regions of Northeast China

The objective here was to identify QTL underlying soybean protein content (PC), and to evaluate the additive and epistatic effects of the QTLs. A mapping population, consisting of 129 recombinant inbred lines (RILs), was created by crossing 'Dongnong 46' and 'L-100'. Phenotypic data of the parents and RILs were collected for 4 years in three locations of Heilongjiang Province of China. A total of 213 SSR markers were used to construct a genetic linkage map. Eight QTLs, located on seven chromosomes (Chr), were identified to be associated with PC among the 10 tested environments. Of the seven QTLs, five QTLs, qPR-2 (Satt710, on Chr9), qPR-3 (Sat_122, on Chr12), qPR-5 (Satt543, on Chr17), qPR-7 (Satt163, on Chr18), and qPR-8 (Satt614, on Chr20), were detected in six, seven, seven, six, and seven environments, respectively, implying relatively stable QTLs. qPR-3 could explain 3.33%-11.26% of the phenotypic variation across eight tested environments. qPR-5 and qPR-8 explained 3.64%-10.1% and 11.86%-18.40% of the phenotypic variation, respectively, across seven tested environments. Eight QTLs associated with PC exhibited additive and (or) additive × environment interaction effects. The results showed that environment-independent QTLs often had higher additive effects. Moreover, five epistatic pairwise QTLs were identified in the 10 environments.

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

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

Genetic architecture of kernel composition in global sorghum germplasm

Background

Sorghum [Sorghum bicolor (L.) Moench] is an important cereal crop for dryland areas in the United States and for small-holder farmers in Africa. Natural variation of sorghum grain composition (protein, fat, and starch) between accessions can be used for crop improvement, but the genetic controls are still unresolved. The goals of this study were to quantify natural variation of sorghum grain composition and to identify single-nucleotide polymorphisms (SNPs) associated with variation in grain composition concentrations.

Results

In this study, we quantified protein, fat, and starch in a global sorghum diversity panel using near-infrared spectroscopy (NIRS). Protein content ranged from 8.1 to 18.8%, fat content ranged from 1.0 to 4.3%, and starch content ranged from 61.7 to 71.1%. Durra and bicolor-durra sorghum from Ethiopia and India had the highest protein and fat and the lowest starch content, while kafir sorghum from USA, India, and South Africa had the lowest protein and the highest starch content. Genome-wide association studies (GWAS) identified quantitative trait loci (QTL) for sorghum protein, fat, and starch. Previously published RNAseq data was used to identify candidate genes within a GWAS QTL region. A putative alpha-amylase 3 gene, which has previously been shown to be associated with grain composition traits, was identified as a strong candidate for protein and fat variation.

Conclusions

We identified promising sources of genetic material for manipulation of grain composition traits, and several loci and candidate genes that may control sorghum grain composition. This survey of grain composition in sorghum germplasm and identification of protein, fat, and starch QTL contributes to our understanding of the genetic basis of natural variation in sorghum grain nutritional traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3403-x) contains supplementary material, which is available to authorized users.

Proteome rebalancing in transgenic Camelina occurs within the enlarged proteome induced by β-carotene accumulation and storage protein suppression

Oilseed crops are global commodities for their oil and protein seed content. We have engineered the oilseed Camelina sativa to exhibit increased protein content with a slight decrease in oil content. The introduction of a phytoene synthase gene with an RNAi cassette directed to suppress the storage protein 2S albumin resulted in seeds with an 11-24 % elevation in overall protein. The phytoene synthase cassette alone produced enhanced β-carotene content of an average 275 ± 6.10 μg/g dry seed and an overall altered seed composition of 11 % less protein and comparable nontransgenic amounts of both oil and carbohydrates. Stacking an RNAi to suppress the major 2S storage protein resulted in seeds that contain elevated protein and slight decrease in oil and carbohydrate amounts showing that Camelina rebalances its proteome within an enlarged protein content genotype. In both β-carotene enhanced seeds with/without RNAi2S suppression, the seed size was noticeably enlarged compared to nontransgenic counterpart seeds. Metabolic analysis of maturing seeds indicate that the enhanced β-carotene trait had the larger effect than the RNAi2S suppression on the seed metabolome. The use of a GRAS (generally regarded as safe) β-carotene as a visual marker in a floral dip transformation system, such as Camelina, might eliminate the need for costly regulatory and controversial antibiotic resistance markers. β-carotene enhanced RNAi2S suppressed Camelina seeds could be further developed as a rapid heterologous protein production platform in a nonfood crop leveraging its enlarged protein content and visual marker.

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.

The transcriptomic signature of developing soybean seeds reveals the genetic basis of seed trait adaptation during domestication

Cultivated soybean has undergone many transformations during domestication. In this paper we report a comprehensive assessment of the evolution of gene co-expression networks based on the analysis of 40 transcriptomes from developing soybean seeds in cultivated and wild soybean accessions. We identified 2680 genes that are differentially expressed during seed maturation and established two cultivar-specific gene co-expression networks. Through analysis of the two networks and integration with quantitative trait locus data we identified two potential key drivers for seed trait formation, GA20OX and NFYA. GA20OX encodes an enzyme in a rate-limiting step of gibberellin biosynthesis, and NFYA encodes a transcription factor. Overexpression of GA20OX and NFYA enhanced seed size/weight and oil content, respectively, in seeds of transgenic plants. The two genes showed significantly higher expression in cultivated than in wild soybean, and the increases in expression were associated with genetic variations in the promoter region of each gene. Moreover, the expression of GA20OX and NFYA in seeds of soybean accessions correlated with seed weight and oil content, respectively. Our study reveals transcriptional adaptation during soybean domestication and may identify a mechanism of selection by expression for seed trait formation, providing strategies for future breeding practice.