Differential Gene Expression Associated with Altered Isoflavone and Fatty Acid Contents in Soybean Mutant Diversity Pool

Identification of candidate genes and genomic prediction of soybean fatty acid components in two soybean populations

Soybean, a source of plant-derived lipids, contains an array of fatty acids essential for health. A comprehensive understanding of the fatty acid profiles in soybean is crucial for enhancing soybean cultivars and augmenting their qualitative attributes. Here, 180 F10 generation recombinant inbred lines (RILs), derived from the cross-breeding of the cultivated soybean variety 'Jidou 12' and the wild soybean 'Y9,' were used as primary experimental subjects. Using inclusive composite interval mapping (ICIM), this study undertook a quantitative trait locus (QTL) analysis on five distinct fatty acid components in the RIL population from 2019 to 2021. Concurrently, a genome-wide association study (GWAS) was conducted on 290 samples from a genetically diverse natural population to scrutinize the five fatty acid components during the same timeframe, thereby aiming to identify loci closely associated with fatty acid profiles. In addition, haplotype analysis and the Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to predict candidate genes. The QTL analysis elucidated 23 stable QTLs intricately associated with the five fatty acid components, exhibiting phenotypic contribution rates ranging from 2.78% to 25.37%. In addition, GWAS of the natural population unveiled 102 significant loci associated with these fatty acid components. The haplotype analysis of the colocalized loci revealed that Glyma.06G221400 on chromosome 6 exhibited a significant correlation with stearic acid content, with Hap1 showing a markedly elevated stearic acid level compared with Hap2 and Hap3. Similarly, Glyma.12G075100 on chromosome 12 was significantly associated with the contents of oleic, linoleic, and linolenic acids, suggesting its involvement in fatty acid biosynthesis. In the natural population, candidate genes associated with the contents of palmitic and linolenic acids were predominantly from the fatty acid metabolic pathway, indicating their potential role as pivotal genes in the critical steps of fatty acid metabolism. Furthermore, genomic selection (GS) for fatty acid components was conducted using ridge regression best linear unbiased prediction based on both random single nucleotide polymorphisms (SNPs) and SNPs significantly associated with fatty acid components identified by GWAS. GS accuracy was contingent upon the SNP set used. Notably, GS efficiency was enhanced when using SNPs derived from QTL mapping analysis and GWAS compared with random SNPs, and reached a plateau when the number of SNP markers exceeded 3,000. This study thus indicates that Glyma.06G221400 and Glyma.12G075100 are genes integral to the synthesis and regulatory mechanisms of fatty acids. It provides insights into the complex biosynthesis and regulation of fatty acids, with significant implications for the directed improvement of soybean oil quality and the selection of superior soybean varieties. The SNP markers delineated in this study can be instrumental in establishing an efficacious pipeline for marker-assisted selection and GS aimed at improving soybean fatty acid components.

QTL mapping reveals key factors related to the isoflavone contents and agronomic traits of soybean (Glycine max)

BACKGROUND

Soybean is a valuable source of edible protein and oil, as well as secondary metabolites that can be used in food products, cosmetics, and medicines. However, because soybean isoflavone content is a quantitative trait influenced by polygenes and environmental interactions, its genetic basis remains unclear.

RESULTS

This study was conducted to identify causal quantitative trait loci (QTLs) associated with soybean isoflavone contents. A mutant-based F2 population (190 individuals) was created by crossing the Korean cultivar Hwanggeum with low isoflavone contents (1,558 µg g-1) and the soybean mutant DB-088 with high isoflavone contents (6,393 µg g-1). A linkage map (3,049 cM) with an average chromosome length of 152 cM was constructed using the 180K AXIOM® SoyaSNP array. Thirteen QTLs related to agronomic traits were mapped to chromosomes 2, 3, 11, 13, 19, and 20, whereas 29 QTLs associated with isoflavone contents were mapped to chromosomes 1, 3, 8, 11, 14, 15, and 17. Notably, the qMGLI11, qMGNI11, qADZI11, and qTI11, which located Gm11_9877690 to Gm11_9955924 interval on chromosome 11, contributed to the high isoflavone contents and explained 11.9% to 20.1% of the phenotypic variation. This QTL region included four candidate genes, encoding β-glucosidases 13, 14, 17-1, and 17-2. We observed significant differences in the expression levels of these genes at various seed developmental stages. Candidate genes within the causal QTLs were functionally characterized based on enriched GO terms and KEGG pathways, as well as the results of a co-expression network analysis. A correlation analysis indicated that certain agronomic traits (e.g., days to flowering, days to maturity, and plant height) are positively correlated with isoflavone content.

CONCLUSIONS

Herein, we reported that the major QTL associated with isoflavone contents was located in the interval from Gm11_9877690 to Gm11_9955924 (78 kb) on chromosome 11. Four β-glucosidase genes were identified that may be involved in high isoflavone contents of soybean DB-088. Thus, the mutant alleles from soybean DB-088 may be useful for marker-assisted selection in developing soybean lines with high isoflavone contents and superior agronomic traits.

Identification of Loci Governing Agronomic Traits and Mutation Hotspots via a GBS-Based Genome-Wide Association Study in a Soybean Mutant Diversity Pool

In this study, we performed a genotyping-by-sequencing analysis and a genome-wide association study of a soybean mutant diversity pool previously constructed by gamma irradiation. A GWAS was conducted to detect significant associations between 37,249 SNPs, 11 agronomic traits, and 6 phytochemical traits. In the merged data set, 66 SNPs on 13 chromosomes were highly associated (FDR p < 0.05) with the following 4 agronomic traits: days of flowering (33 SNPs), flower color (16 SNPs), node number (6 SNPs), and seed coat color (11 SNPs). These results are consistent with the findings of earlier studies on other genetic features (e.g., natural accessions and recombinant inbred lines). Therefore, our observations suggest that the genomic changes in the mutants generated by gamma irradiation occurred at the same loci as the mutations in the natural soybean population. These findings are indicative of the existence of mutation hotspots, or the acceleration of genome evolution in response to high doses of radiation. Moreover, this study demonstrated that the integration of GBS and GWAS to investigate a mutant population derived from gamma irradiation is suitable for dissecting the molecular basis of complex traits in soybeans.

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.

Use of Gamma Radiation for the Genetic Improvement of Underutilized Plant Varieties

Agricultural biodiversity includes many species that have biological variants (natives, ecotypes, races, morphotypes). Their use is restricted to local areas because they do not fulfill the commercial requirements; however, it is well documented that these species are a source of metabolites, proteins, enzymes, and genes. Rescuing and harnessing them through traditional genetic breeding is time-consuming and expensive. Inducing mutagenesis may be a short-time option for its genetic improvement. A review of outstanding research was carried out, in order to become familiar with gene breeding using gamma radiation and its relevance to obtain outstanding agronomic characteristics for underutilized species. An approach was made to the global panorama of the application of gamma radiation in different conventional crop species and in vitro cultivated species, in order to obtain secondary metabolites, as well as molecular tools used for mutation screening. The varied effects of gamma radiation are essentially the result of the individual responses and phenotypic plasticity of each organism. However, even implicit chance can be reduced with specific genetic breeding, environmental adaptation, or conservation objectives.