QTL Mapping of Isoflavone, Oil and Protein Contents in Soybean (Glycine max L. Merr.)

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.

The challenge of breeding for reduced off-flavor in faba bean ingredients

The growing interest in plant protein sources, such as pulses, is driven by the necessity for sustainable food production and climate change mitigation strategies. Faba bean (Vicia faba L.) is a promising protein crop for temperate climates, owing to its remarkable yield potential (up to 8 tonnes ha-1 in favourable growing conditions) and high protein content (~29% dry matter basis). Nevertheless, the adoption of faba bean protein in plant-based products that aim to resemble animal-derived counterparts is hindered by its distinctive taste and aroma, regarded as "off-flavors". In this review, we propose to introduce off-flavor as a trait in breeding programs by identifying molecules involved in sensory perception and defining key breeding targets. We discuss the role of lipid oxidation in producing volatile and non-volatile compounds responsible for the beany aroma and bitter taste, respectively. We further investigate the contribution of saponin, tannin, and other polyphenols to bitterness and astringency. To develop faba bean varieties with diminished off-flavors, we suggest targeting genes to reduce lipid oxidation, such as lipoxygenases (lox) and fatty acid desaturases (fad), and genes involved in phenylpropanoid and saponin biosynthesis, such as zero-tannin (zt), chalcone isomerase (chi), chalcone synthase (chs), β-amyrin (bas1). Additionally, we address potential challenges, including the need for high-throughput phenotyping and possible limitations that could arise during the genetic improvement process. The breeding approach can facilitate the use of faba bean protein in plant-based food such as meat and dairy analogues more extensively, fostering a transition toward more sustainable and climate-resilient diets.

Genetic analysis of protein content and oil content in soybean by genome-wide association study

Soybean seed protein content (PC) and oil content (OC) have important economic value. Detecting the loci/gene related to PC and OC is important for the marker-assisted selection (MAS) breeding of soybean. To detect the stable and new loci for PC and OC, a total of 320 soybean accessions collected from the major soybean-growing countries were used to conduct a genome-wide association study (GWAS) by resequencing. The PC ranged from 37.8% to 46.5% with an average of 41.1% and the OC ranged from 16.7% to 22.6% with an average of 21.0%. In total, 23 and 29 loci were identified, explaining 3.4%-15.4% and 5.1%-16.3% of the phenotypic variations for PC and OC, respectively. Of these, eight and five loci for PC and OC, respectively, overlapped previously reported loci and the other 15 and 24 loci were newly identified. In addition, nine candidate genes were identified, which are known to be involved in protein and oil biosynthesis/metabolism, including lipid transport and metabolism, signal transduction, and plant development pathway. These results uncover the genetic basis of soybean protein and oil biosynthesis and could be used to accelerate the progress in enhancing soybean PC and OC.

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.

MOTHER-OF-FT-AND-TFL1 regulates the seed oil and protein content in soybean

Soybean is a major crop that produces valuable seed oil and protein for global consumption. Seed oil and protein are regulated by complex quantitative trait loci (QTLs) and have undergone intensive selections during the domestication of soybean. It is essential to identify the major genetic components and understand their mechanism behind seed oil and protein in soybean. We report that MOTHER-OF-FT-AND-TFL1 (GmMFT) is the gene of a classical QTL that has been reported to regulate seed oil and protein content in many studies. Mutation of MFT decreased seeds oil content and weight in both Arabidopsis and soybean, whereas increased expression of GmMFT enhanced seeds oil content and weight. Haplotype analysis showed that GmMFT has undergone selection, which resulted in the extended haplotype homozygosity in the cultivated soybean and the enriching of the oil-favorable allele in modern soybean cultivars. This work unraveled the GmMFT-mediated mechanism regulating seed oil and protein content and seed weight, and revealed a previously unknown function of MFT that provides new insights into targeted soybean improvement and breeding.

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.

Fine-Mapping and Functional Analyses of a Candidate Gene Controlling Isoflavone Content in Soybeans Seed

Isoflavones, one of the most important secondary metabolites produced by soybeans (Glycine max (L.) Merr.), are important for a variety of biological processes, and are beneficial for human health. To identify genetic loci underlying soybean isoflavone content, a mapping population containing 119 F_5:18 recombinant inbred lines, derived by crossing soybean cultivar "Zhongdou27" with "Dongong8004," was used. We identified 15 QTLs associated with isoflavone contents. A novel loci, qISO19-1, was mapped onto soybean chromosome 19 and was fine-mapped to a 62.8 kb region using a BC₂F₂ population. We considered GmMT1 as a candidate gene for the qISO19-1 locus due to the significant positive correlation recovered between its expression level and isoflavone content in the seeds of 43 soybean germplasms. Overexpression of GmMT1 in Arabidopsis and soybean cultivars increased isoflavone contents. Transgenic soybeans overexpressing GmMT1 also exhibited improved resistance to pathogenic infection, while transgenic Arabidopsis resisted salt and drought stress.

Protein Biofortification in Lentils (Lens culinaris Medik.) Toward Human Health

Lentil (Lens culinaris Medik.) is a nutritionally dense crop with significant quantities of protein, low-digestible carbohydrates, minerals, and vitamins. The amino acid composition of lentil protein can impact human health by maintaining amino acid balance for physiological functions and preventing protein-energy malnutrition and non-communicable diseases (NCDs). Thus, enhancing lentil protein quality through genetic biofortification, i.e., conventional plant breeding and molecular technologies, is vital for the nutritional improvement of lentil crops across the globe. This review highlights variation in protein concentration and quality across Lens species, genetic mechanisms controlling amino acid synthesis in plants, functions of amino acids, and the effect of antinutrients on the absorption of amino acids into the human body. Successful breeding strategies in lentils and other pulses are reviewed to demonstrate robust breeding approaches for protein biofortification. Future lentil breeding approaches will include rapid germplasm selection, phenotypic evaluation, genome-wide association studies, genetic engineering, and genome editing to select sequences that improve protein concentration and quality.

Identification of Candidate Genes and Genomic Selection for Seed Protein in Soybean Breeding Pipeline

Soybean is a primary meal protein for human consumption, poultry, and livestock feed. In this study, quantitative trait locus (QTL) controlling protein content was explored via genome-wide association studies (GWAS) and linkage mapping approaches based on 284 soybean accessions and 180 recombinant inbred lines (RILs), respectively, which were evaluated for protein content for 4 years. A total of 22 single nucleotide polymorphisms (SNPs) associated with protein content were detected using mixed linear model (MLM) and general linear model (GLM) methods in Tassel and 5 QTLs using Bayesian interval mapping (IM), single-trait multiple interval mapping (SMIM), single-trait composite interval mapping maximum likelihood estimation (SMLE), and single marker regression (SMR) models in Q-Gene and IciMapping. Major QTLs were detected on chromosomes 6 and 20 in both populations. The new QTL genomic region on chromosome 6 (Chr6_18844283-19315351) included 7 candidate genes and the Hap.X ^AA at the Chr6_19172961 position was associated with high protein content. Genomic selection (GS) of protein content was performed using Bayesian Lasso (BL) and ridge regression best linear unbiased prediction (rrBULP) based on all the SNPs and the SNPs significantly associated with protein content resulted from GWAS. The results showed that BL and rrBLUP performed similarly; GS accuracy was dependent on the SNP set and training population size. GS efficiency was higher for the SNPs derived from GWAS than random SNPs and reached a plateau when the number of markers was >2,000. The SNP markers identified in this study and other information were essential in establishing an efficient marker-assisted selection (MAS) and GS pipelines for improving soybean protein content.

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.

Quantitative Trait Locus Mapping of Seed Vigor in Soybean under -20 °C Storage and Accelerated Aging Conditions via RAD Sequencing

Due to its fast deterioration, soybean (Glycine max L.) has an inherently poor seed vigor. Vigor loss occurring during storage is one of the main obstacles to soybean production in the tropics. To analyze the genetic background of seed vigor, soybean seeds of a recombinant inbred line (RIL) population derived from the cross between Zhonghuang24 (ZH24, low vigor cultivar) and Huaxia3hao (HX3, vigorous cultivar) were utilized to identify the quantitative trait loci (QTLs) underlying the seed vigor under -20 °C conservation and accelerated aging conditions. According to the linkage analysis, multiple seed vigor-related QTLs were identified under both -20 °C and accelerated aging storage. Two major QTLs and eight QTL hotspots localized on chromosomes 3, 6, 9, 11, 15, 16, 17, and 19 were detected that were associated with seed vigor across two storage conditions. The indicators of seed vigor did not correlate well between the two aging treatments, and no common QTLs were detected in RIL populations stored in two conditions. These results indicated that deterioration under accelerated aging conditions was not reflective of natural aging at -20 °C. Additionally, we suggest 15 promising candidate genes that could possibly determine the seed vigor in soybeans, which would help explore the mechanisms responsible for maintaining high seed vigor.

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Identification of SNP markers associated with soybean fatty acids contents by genome-wide association analyses

Composition of fatty acids (FAs) in soybean seed is important for the quality and uses of soybean oil. Using gas chromatography, we have measured soybean FAs profiles of 621 soybean accessions (maturity groups I through IV) grown in five different environments; Columbus, OH (2015), Wooster, OH (2014 and 2015), Plymouth, NC (2015), and Urbana, IL (2015). Using publicly available SoySNP50K genotypic data and the FA profiles from this study, a genome-wide association analysis was completed with a compressed mixed linear model to identify 43 genomic regions significantly associated with a fatty acid at a genome wide significance threshold of 5%. Among these regions, one and three novel genomic regions associated with palmitic acid and stearic acid, respectively, were identified across all five environments. Additionally, nine novel environment-specific FA-related genomic regions were discovered providing new insights into the genetics of soybean FAs. Previously reported FA-related loci, such as FATB1a, SACPD-C, and KASIII, were also confirmed in this study. Our results will be useful for future functional studies and marker-assisted breeding for soybean FAs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01216-1.

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 mapping for protein, oil and water-soluble protein contents in soybean

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

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.

The Effects of Domestication on Secondary Metabolite Composition in Legumes

Legumes are rich in secondary metabolites, such as polyphenols, alkaloids, and saponins, which are important defense compounds to protect the plant against herbivores and pathogens, and act as signaling molecules between the plant and its biotic environment. Legume-sourced secondary metabolites are well known for their potential benefits to human health as pharmaceuticals and nutraceuticals. During domestication, the color, smell, and taste of crop plants have been the focus of artificial selection by breeders. Since these agronomic traits are regulated by secondary metabolites, the basis behind the genomic evolution was the selection of the secondary metabolite composition. In this review, we will discuss the classification, occurrence, and health benefits of secondary metabolites in legumes. The differences in their profiles between wild legumes and their cultivated counterparts will be investigated to trace the possible effects of domestication on secondary metabolite compositions, and the advantages and drawbacks of such modifications. The changes in secondary metabolite contents will also be discussed at the genetic level to examine the genes responsible for determining the secondary metabolite composition that might have been lost due to domestication. Understanding these genes would enable breeding programs and metabolic engineering to produce legume varieties with favorable secondary metabolite profiles for facilitating adaptations to a changing climate, promoting beneficial interactions with biotic factors, and enhancing health-beneficial secondary metabolite contents for human consumption.

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.

Overcoming collinearity in path analysis of soybean [Glycine max (L.) Merr.] grain oil content

Path analysis allows understanding the direct and indirect effects among traits. Multicollinearity in correlation matrices may cause a bias in path analysis estimates. This study aimed to: a) understand the correlation among soybean traits and estimate their direct and indirect effects on gain oil content; b) verify the efficiency of ridge path analysis and trait culling to overcome colinearity. Three different matrices with different levels of collinearity were obtained by trait culling. Ridge path analysis was performed on matrices with strong collinearity; otherwise, a traditional path analysis was performed. The same analyses were run on a simulated dataset. Trait culling was applied to matrix R originating the matrices R1 and R2. Path analysis for matrices R1 and R2 presented a high determination coefficient (0.856 and 0.832, respectively) and low effect of the residual variable (0.379 and 0.410 respectively). Ridge path analysis presented low determination coefficient (0.657) and no direct effects greater than the effects of the residual variable (0.585). Trait culling was more effective to overcome collinearity. Mass of grains, number of nodes, and number of pods are promising for indirect selection for oil content.

Genome-wide association study of soybean seed germination under drought stress

Drought stress, which is increasing with climate change, is a serious threat to agricultural sustainability worldwide. Seed germination is an essential growth phase that ensures the successful establishment and productivity of soybean, which can lose substantial productivity in soils with water deficits. However, only limited genetic information is available about how germinating soybean seeds may exert drought tolerance. In this study, we examined the germinating seed drought-tolerance phenotypes and genotypes of a panel of 259 released Chinese soybean cultivars panel. Based on 4616 Single-Nucleotide Polymorphisms (SNPs), we conducted a mixed-linear model GWAS that identified a total of 15 SNPs associated with at least one drought-tolerance index. Notably, three of these SNPs were commonly associated with two drought-tolerance indices. Two of these SNPs are positioned upstream of genes, and 11 of them are located in or near regions where QTLs have been previously mapped by linkage analysis, five of which are drought-related. The SNPs detected in this study can both drive hypothesis-driven research to deepen our understanding of genetic basis of soybean drought tolerance at the germination stage and provide useful genetic resources that can facilitate the selection of drought stress traits via genomic-assisted selection.

Quantitative trait loci analysis of seed oil content and composition of wild and cultivated soybean

BACKGROUND

Soybean oil is a major source of edible oil, and the domestication of wild soybean has resulted in significant changes in oil content and composition. Extensive efforts have been made to identify genetic loci that are related to soybean oil traits. The objective of this study was to identify quantitative trait loci (QTLs) related to soybean seed oil and compare the fatty acid composition between wild and cultivated soybean.

RESULTS

Using the specific-locus amplified fragment sequencing (SLAF-seq) method, a total of 181 recombinant inbred lines (RILs) derived from a cross between wild soybean ZYD00463 (Glycine soja) and cultivated soybean WDD01514 (Glycine max) were genotyped. Finally, a high-density genetic linkage map comprising 11,398 single-nucleotide polymorphism (SNP) markers on 20 linkage groups (LGs) was constructed. Twenty-four stable QTLs for seed oil content and composition were identified by model-based composite interval mapping (CIM) across multiple environments. Among these QTLs, 23 overlapped with or were adjacent to previously reported QTLs. One QTL, qPA10_1 (5.94-9.98 Mb) on Chr. Ten is a novel locus for palmitic acid. In the intervals of stable QTLs, some interesting genes involved in lipid metabolism were detected.

CONCLUSIONS

We developed 181 RILs from a cross between wild soybean ZYD00463 and cultivated soybean WDD01514 and constructed a high-density genetic map using the SLAF-seq method. We identified 24 stable QTLs for seed oil content and compositions, which includes qPA10_1 on Chr. 10, a novel locus for palmitic acid. Some interesting genes in the QTL regions were also detected. Our study will provide useful information for scientists to learn about genetic variations in lipid metabolism between wild and cultivated soybean.

Genome Wide Association Study and Genomic Selection of Amino Acid Concentrations in Soybean Seeds

Soybean is a major source of protein for human consumption and animal feed. Releasing new cultivars with high nutritional value is one of the major goals in soybean breeding. To achieve this goal, genome-wide association studies of seed amino acid contents were conducted based on 249 soybean accessions from China, US, Japan, and South Korea. The accessions were evaluated for 15 amino acids and genotyped by sequencing. Significant genetic variation was observed for amino acids among the accessions. Among the 231 single nucleotide polymorphisms (SNPs) significantly associated with variations in amino acid contents, fifteen SNPs localized near 14 candidate genes involving in amino acid metabolism. The amino acids were classified into two groups with five in one group and seven amino acids in the other. Correlation coefficients among the amino acids within each group were high and positive, but the correlation coefficients of amino acids between the two groups were negative. Twenty-five SNP markers associated with multiple amino acids can be used to simultaneously improve multi-amino acid concentration in soybean. Genomic selection analysis of amino acid concentration showed that selection efficiency of amino acids based on the markers significantly associated with all 15 amino acids was higher than that based on random markers or markers only associated with individual amino acid. The identified markers could facilitate selection of soybean varieties with improved seed quality.

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

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

Identification and characterization of a fast-neutron-induced mutant with elevated seed protein content in soybean

KEY MESSAGE

Protein content of soybean is critical for utility of soybean meal. A fast-neutron-induced deletion on chromosome 12 was found to be associated with increased protein content. Soybean seed composition affects the utility of soybean, and improving seed composition is an essential breeding goal. Fast neutron radiation introduces genomic mutations resulting in novel variation for traits of interest. Two elite soybean lines were irradiated with fast neutrons and screened for altered seed composition. Twenty-three lines with altered protein, oil, or sucrose content were selected based on near-infrared spectroscopy data from five environments and yield tested at five locations. Mutants with significantly increased protein averaged 19.1-36.8 g kg-1 more protein than the parents across 10 environments. Comparative genomic hybridization (CGH) identified putative mutations in a mutant, G15FN-12, that has 36.8 g kg-1 higher protein than the parent genotype, and whole genome sequencing (WGS) of the mutant has confirmed these mutations. An F2:3 population was developed from G15FN-12 to determine association between genomic changes and increased protein content. Bulked segregant analysis of the population using the SoySNP50K BeadChip identified a CGH- and WGS-confirmed deletion on chromosome 12 to be responsible for elevated protein content. The population was genotyped using a KASP marker designed at the mutation region, and significant association (P < 0.0001) between the deletion on chromosome 12 and elevated protein content was observed and confirmed in the F3:4 generation. The F2 segregants homozygous for the deletion averaged 27 g kg-1 higher seed protein and 8 g kg-1 lower oil than homozygous wild-type segregants. Mutants with altered seed composition are a new resource for gene function studies and provide elite materials for genetic improvement of seed composition.

Dissecting the Genetic Architecture of Seed Protein and Oil Content in Soybean from the Yangtze and Huaihe River Valleys Using Multi-Locus Genome-Wide Association Studies

Soybean is a globally important legume crop that provides a primary source of high-quality vegetable protein and oil. Seed protein and oil content are two valuable quality traits controlled by multiple genes in soybean. In this study, the restricted two-stage multi-locus genome-wide association analysis (RTM-GWAS) procedure was performed to dissect the genetic architecture of seed protein and oil content in a diverse panel of 279 soybean accessions from the Yangtze and Huaihe River Valleys in China. We identified 26 quantitative trait loci (QTLs) for seed protein content and 23 for seed oil content, including five associated with both traits. Among these, 39 QTLs corresponded to previously reported QTLs, whereas 10 loci were novel. As reported previously, the QTL on chromosome 20 was associated with both seed protein and oil content. This QTL exhibited opposing effects on these traits and contributed the most to phenotype variation. From the detected QTLs, 55 and 51 candidate genes were identified for seed protein and oil content, respectively. Among these genes, eight may be promising candidate genes for improving soybean nutritional quality. These results will facilitate marker-assisted selective breeding for soybean protein and oil content traits.

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

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.

Fine-mapping of QTLs for individual and total isoflavone content in soybean (Glycine max L.) using a high-density genetic map

KEY MESSAGE

Fifteen stable QTLs were identified using a high-density soybean genetic map across multiple environments. One major QTL, qIF5-1, contributing to total isoflavone content explained phenotypic variance 49.38, 43.27, 46.59, 45.15 and 52.50%, respectively. Soybeans (Glycine max L.) are a major source of dietary isoflavones. To identify novel quantitative trait loci (QTL) underlying isoflavone content, and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population comprised of 196 F7:8-10 recombinant inbred lines (RILs, Huachun 2 × Wayao) was utilized to evaluate individual and total isoflavone content in plants grown in four different environments in Guangdong. A high-density genetic linkage map containing 3469 recombination bin markers based on 0.2 × restriction site-associated DNA tag sequencing (RAD-seq) technology was used to finely map QTLs for both individual and total isoflavone contents. Correlation analyses showed that total isoflavone content, and that of five individual isoflavone, was significantly correlated across the four environments. Based on the high-density genetic linkage map, a total of 15 stable quantitative trait loci (QTLs) associated with isoflavone content across multiple environments were mapped onto chromosomes 02, 05, 07, 09, 10, 11, 13, 16, 17, and 19. Further, one of them, qIF5-1, localized to chromosomes 05 (38,434,171-39,045,620 bp) contributed to almost all isoflavone components across all environments, and explained 6.37-59.95% of the phenotypic variance, especially explained 49.38, 43.27, 46.59, 45.15 and 52.50% for total isoflavone. The results obtained in the present study will pave the way for a better understanding of the genetics of isoflavone accumulation and reveals the scope available for improvement of isoflavone content through marker-assisted selection.

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.

A Genome-Wide Association Study for Agronomic Traits in Soybean Using SNP Markers and SNP-Based Haplotype Analysis

Mapping quantitative trait loci through the use of linkage disequilibrium (LD) in populations of unrelated individuals provides a valuable approach for dissecting the genetic basis of complex traits in soybean (Glycine max). The haplotype-based genome-wide association study (GWAS) has now been proposed as a complementary approach to intensify benefits from LD, which enable to assess the genetic determinants of agronomic traits. In this study a GWAS was undertaken to identify genomic regions that control 100-seed weight (SW), plant height (PH) and seed yield (SY) in a soybean association mapping panel using single nucleotide polymorphism (SNP) markers and haplotype information. The soybean cultivars (N = 169) were field-evaluated across four locations of southern Brazil. The genome-wide haplotype association analysis (941 haplotypes) identified eleven, seventeen and fifty-nine SNP-based haplotypes significantly associated with SY, SW and PH, respectively. Although most marker-trait associations were environment and trait specific, stable haplotype associations were identified for SY and SW across environments (i.e., haplotypes Gm12_Hap12). The haplotype block 42 on Chr19 (Gm19_Hap42) was confirmed to be associated with PH in two environments. These findings enable us to refine the breeding strategy for tropical soybean, which confirm that haplotype-based GWAS can provide new insights on the genetic determinants that are not captured by the single-marker approach.

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.

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.

Detecting the QTL-allele system of seed isoflavone content in Chinese soybean landrace population for optimal cross design and gene system exploration

KEY MESSAGE

Utilizing an innovative GWAS in CSLRP, 44 QTL 199 alleles with 72.2 % contribution to SIFC variation were detected and organized into a QTL-allele matrix for cross design and gene annotation. The seed isoflavone content (SIFC) of soybeans is of great importance to health care. The Chinese soybean landrace population (CSLRP) as a genetic reservoir was studied for its whole-genome quantitative trait loci (QTL) system of the SIFC using an innovative restricted two-stage multi-locus genome-wide association study procedure (RTM-GWAS). A sample of 366 landraces was tested under four environments and sequenced using RAD-seq (restriction-site-associated DNA sequencing) technique to obtain 116,769 single nucleotide polymorphisms (SNPs) then organized into 29,119 SNP linkage disequilibrium blocks (SNPLDBs) for GWAS. The detected 44 QTL 199 alleles on 16 chromosomes (explaining 72.2 % of the total phenotypic variation) with the allele effects (92 positive and 107 negative) of the CSLRP were organized into a QTL-allele matrix showing the SIFC population genetic structure. Additional differentiation among eco-regions due to the SIFC in addition to that of genome-wide markers was found. All accessions comprised both positive and negative alleles, implying a great potential for recombination within the population. The optimal crosses were predicted from the matrices, showing transgressive potentials in the CSLRP. From the detected QTL system, 55 candidate genes related to 11 biological processes were χ (2)-tested as an SIFC candidate gene system. The present study explored the genome-wide SIFC QTL/gene system with the innovative RTM-GWAS and found the potentials of the QTL-allele matrix in optimal cross design and population genetic and genomic studies, which may have provided a solution to match the breeding by design strategy at both QTL and gene levels in breeding programs.

Quantitative Trait Loci (QTL) Mapping for Glycinin and β-Conglycinin Contents in Soybean (Glycine max L. Merr.)

Compared to β-conglycinin, glycinin contains 3-4 times the methionine and cysteine (sulfur-containing amino acids), accounting for approximately 40 and 30%, respectively, of the total storage protein in soybean. Increasing the soybean storage protein content while improving the ratio of glycinin to β-conglycinin is of great significance for soybean breeding and soy food products. The objective of this study is to analyze the genetic mechanism regulating the glycinin and β-conglycinin contents of soybean by using a recombinant inbred line (RIL) population derived from a cross between Kefeng No. 1 and Nannong 1138-2. Two hundred and twenty-one markers were used to map quantitative trait loci (QTLs) for glycinin (11S) and β-conglycinin (7S) contents, the ratio of glycinin to β-conglycinin (RGC), and the sum of glycinin and β-conglycinin (SGC). A total of 35 QTLs, 3 pairs of epistatic QTLs, and 5 major regions encompassing multiple QTLs were detected. Genes encoding the subunits of β-conglycinin were localized to marker intervals sat_418-satt650 and sat_196-sat_303, which are linked to RGC and SGC; marker sat_318, associated with 11S, 7S, and SGC, was located near Glyma10g04280 (Gy4), which encodes a subunit of glycinin. These results, which take epistatic interactions into account, will improve our understanding of the genetic basis of 11S and 7S contents and will lay a foundation for marker-assisted selection (MAS) breeding of soybean and improving the quality of soybean products.

Evaluation of genetic variation among Brazilian soybean cultivars through genome resequencing

BACKGROUND

Soybean [Glycine max (L.) Merrill] is one of the most important legumes cultivated worldwide, and Brazil is one of the main producers of this crop. Since the sequencing of its reference genome, interest in structural and allelic variations of cultivated and wild soybean germplasm has grown. To investigate the genetics of the Brazilian soybean germplasm, we selected soybean cultivars based on the year of commercialization, geographical region and maturity group and resequenced their genomes.

RESULTS

We resequenced the genomes of 28 Brazilian soybean cultivars with an average genome coverage of 14.8X. A total of 5,835,185 single nucleotide polymorphisms (SNPs) and 1,329,844 InDels were identified across the 20 soybean chromosomes, with 541,762 SNPs, 98,922 InDels and 1,093 CNVs that were exclusive to the 28 Brazilian cultivars. In addition, 668 allelic variations of 327 genes were shared among all of the Brazilian cultivars, including genes related to DNA-dependent transcription-elongation, photosynthesis, ATP synthesis-coupled electron transport, cellular respiration, and precursors of metabolite generation and energy. A very homogeneous structure was also observed for the Brazilian soybean germplasm, and we observed 41 regions putatively influenced by positive selection. Finally, we detected 3,880 regions with copy-number variations (CNVs) that could help to explain the divergence among the accessions evaluated.

CONCLUSIONS

The large number of allelic and structural variations identified in this study can be used in marker-assisted selection programs to detect unique SNPs for cultivar fingerprinting. The results presented here suggest that despite the diversification of modern Brazilian cultivars, the soybean germplasm remains very narrow because of the large number of genome regions that exhibit low diversity. These results emphasize the need to introduce new alleles to increase the genetic diversity of the Brazilian germplasm.

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.

Mapping isoflavone QTL with main, epistatic and QTL × environment effects in recombinant inbred lines of soybean

Soybean (Glycine max (L.) Merr.) isoflavone is important for human health and plant defense system. To identify novel quantitative trait loci (QTL) and epistatic QTL underlying isoflavone content in soybean, F5:6, F5:7 and F5:8 populations of 130 recombinant inbred (RI) lines, derived from the cross of soybean cultivar 'Zhong Dou 27' (high isoflavone) and 'Jiu Nong 20' (low isoflavone), were analyzed with 95 new SSR markers. A new linkage map including 194 SSR markers and covering 2,312 cM with mean distance of about 12 cM between markers was constructed. Thirty four QTL for both individual and total seed isoflavone contents of soybean were identified. Six, seven, ten and eleven QTL were associated with daidzein (DZ), glycitein (GC), genistein (GT) and total isoflavone (TI), respectively. Of them 23 QTL were newly identified. The qTIF_1 between Satt423 and Satt569 shared the same marker Satt569 with qDZF_2, qGTF_1 and qTIF_2. The qGTD2_1 between Satt186 and Satt226 was detected in four environments and explained 3.41%-10.98% of the phenotypic variation. The qGTA2_1, overlapped with qGCA2_1 and detected in four environments, was close to the previously identified major QTL for GT, which were responsible for large a effects. QTL (qDZF_2, qGTF_1 and qTIF_2) between Satt144-Satt569 were either clustered or pleiotropic. The qGCM_1, qGTM_1 and qTIM_1 between Satt540-Sat_244 explained 2.02%-9.12% of the phenotypic variation over six environments. Moreover, the qGCE_1 overlapped with qGTE_1 and qTIE_1, the qTIH_2 overlapped with qGTH_1, qGCI_1 overlapped with qDZI_1, qTIL_1 overlapped with qGTL_1, and qTIO_1 overlapped with qGTO_1. In this study, some of unstable QTL were detected in different environments, which were due to weak expression of QTL, QTL by environment interaction in the opposite direction to a effects, and/or epistasis. The markers identified in multi-environments in this study could be applied in the selection of soybean cultivars for higher isoflavone content and in the map-based gene cloning.

Expression quantitative trait loci infer the regulation of isoflavone accumulation in soybean (Glycine max L. Merr.) seed

BACKGROUND

Mapping expression quantitative trait loci (eQTL) of targeted genes represents a powerful and widely adopted approach to identify putative regulatory variants. Linking regulation differences to specific genes might assist in the identification of networks and interactions. The objective of this study is to identify eQTL underlying expression of four gene families encoding isoflavone synthetic enzymes involved in the phenylpropanoid pathway, which are phenylalanine ammonia-lyase (PAL; EC 4.3.1.5), chalcone synthase (CHS; EC 2.3.1.74), 2-hydroxyisoflavanone synthase (IFS; EC1.14.13.136) and flavanone 3-hydroxylase (F3H; EC 1.14.11.9). A population of 130 recombinant inbred lines (F5:11), derived from a cross between soybean cultivar 'Zhongdou 27' (high isoflavone) and 'Jiunong 20' (low isoflavone), and a total of 194 simple sequence repeat (SSR) markers were used in this study. Overlapped loci of eQTLs and phenotypic QTLs (pQTLs) were analyzed to identify the potential candidate genes underlying the accumulation of isoflavone in soybean seed.

RESULTS

Thirty three eQTLs (thirteen cis-eQTLs and twenty trans-eQTLs) underlying the transcript abundance of the four gene families were identified on fifteen chromosomes. The eQTLs between Satt278-Sat_134, Sat_134-Sct_010 and Satt149-Sat_234 underlie the expression of both IFS and CHS genes. Five eQTL intervals were overlapped with pQTLs. A total of eleven candidate genes within the overlapped eQTL and pQTL were identified.

CONCLUSIONS

These results will be useful for the development of marker-assisted selection to breed soybean cultivars with high or low isoflavone contents and for map-based cloning of new isoflavone related genes.

Problems with and a system to eliminate single-primer PCR product contamination in simple sequence repeat molecular marker-assisted selection in soybean

Polymerase chain reaction (PCR) provides a foundation for simple sequence repeat molecular marker-assisted selection (SSR MAS) in soybean. This PCR system and its various conditions have been optimized by many researchers. However, current research on the optimization of the PCR system focuses on double-primer PCR products. We compared single- and double-SSR primer PCR products from 50 soybean samples and found that the use of single-PCR primers in the reaction system can lead to amplified fragments of portions of the SSR primers in the PCR process, resulting in both false-positives and fragment impurity of double-primer PCR amplification, inconvenient for subsequent analysis. We used "single-primer PCR correction" to eliminate interference caused by single-primer nonspecific PCR amplification and improve PCR quality. Using this method, the precision and success rates of SSR MAS in soybean can be increased.