Next Article in Journal
Drought and Heat Stress in Cool-Season Food Legumes in Sub-Tropical Regions: Consequences, Adaptation, and Mitigation Strategies
Next Article in Special Issue
Development of an Aus-Derived Nested Association Mapping (Aus-NAM) Population in Rice
Previous Article in Journal
Wetland Restoration with Hydrophytes: A Review
Previous Article in Special Issue
Characterization of the Genetic Diversity Present in a Diverse Sesame Landrace Collection Based on Phenotypic Traits and EST-SSR Markers Coupled With an HRM Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 6
10.3390/plants10061037
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Dong-Gun Kim
1,†,
Jae-Il Lyu
1,†,
You-Jin Lim
2,
Jung-Min Kim
1,
Nguyen-Ngoc Hung
1,
Seok-Hyun Eom
2,
Sang-Hoon Kim
1,
Jin-Baek Kim
1,
Chang-Hyu Bae
3,* and
Soon-Jae Kwon
1,*
1
Advanced Radiation Technology Institute, Korea Atomic Energy Research Institute, Jeongup 56212, Korea
2
Department of Horticultural Biotechnology, Institute of Life Sciences & Resources, Kyung Hee University, Yongin 17104, Korea
3
Department of Life Resources, Graduate School, Sunchon National University, Suncheon 57922, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2021, 10(6), 1037; https://doi.org/10.3390/plants10061037
Submission received: 8 April 2021 / Revised: 18 May 2021 / Accepted: 18 May 2021 / Published: 21 May 2021
(This article belongs to the Special Issue Genetic Resources and Crop Improvement)
Browse Figures
Versions Notes

Abstract

:
Soybean seeds are consumed worldwide owing to their nutritional value and health benefits. In this study we investigated the metabolic properties of 208 soybean mutant diversity pool (MDP) lines by measuring the isoflavone and fatty acid contents of the seed. The total isoflavone content (TIC) ranged from 0.88 mg/g to 7.12 mg/g and averaged 3.08 mg/g. The proportion of oleic acid among total fatty acids (TFA) ranged from 0.38% to 24.66% and averaged 11.02%. Based on the TIC and TFA among the 208 MDP lines, we selected six lines with altered isoflavone content and six lines with altered oleic acid content compared with those of the corresponding wild-types for measuring gene expression. Each of twelve genes from the isoflavone and fatty acid biosynthesis pathways were analyzed at three different seed developmental stages. Isoflavone biosynthetic genes, including CHI1A, IFS1, and IFS2, showed differences in stages and expression patterns among individuals and wild-types, whereas MaT7 showed consistently higher expression levels in three mutants with increased isoflavone content at stage 1. Expression patterns of the 12 fatty acid biosynthetic genes were classifiable into two groups that reflected the developmental stages of the seeds. The results will be useful for functional analysis of the regulatory genes involved in the isoflavone and fatty acid biosynthetic pathways in soybean.

1. Introduction

Soybean is an important oilseed crop and is among the most widely used and healthy edible oils. Soybean oil is low in saturated fat, rich in essential fatty acids, an excellent source of vitamin E, and contains no cholesterol. The oil quality is primarily a function of the fatty acid composition. The ratio and amount of saturated and unsaturated fatty acids determine the physical, chemical, and nutritional values of the oil [1]. The average fatty acid composition of commercial soybean oil is approximately 12% palmitic (16:0), 4% stearic (18:0), 23% oleic (18:1), 53% linoleic (18:2), and 8% linolenic (18:3) acids [2]. These are the predominant fatty acids in soybean oil, but the ideal fatty acid composition depends on the specific use for which it is intended. Thus, modification of the fatty acid composition of soybean oil is beneficial for food and industrial applications.
Soybean plant growth and development are traditionally classified into six vegetative and eight reproductive (R) stages. R1 and R2 are stages of flowering, R3 and R4 encompass pod formation, and seed development occurs in stages R5–R8 [3,4]. The final seed composition of secondary metabolites, including isoflavones, is strongly influenced by environmental stresses during stages R5–R7 [5,6,7]. Seed development in stage R5 is characterized by a rapid increase in weight and nutrient accumulation, which continues until R6 [3]. Seeds in the R6 stage fill the pod cavity but are still immature [4]. Typical of the R7 stage, the seed coat color begins to change from green to either tan or yellow, depending on the cultivar [8]. At this stage, accumulation of dry weight ceases and the seed has attained physiological maturity [9]. Isoflavones accumulate in the seeds during the advanced stages of seed maturation [10,11], and isoflavone contents are strongly influenced by water availability during this period.
Isoflavones are polyphenolic secondary plant metabolites found in seedlings, flowers, and roots, and are especially abundant in seeds and leaves of soybean. Within seeds, diverse tissues have the ability to synthesize isoflavones [12]. Multigenic responses to abiotic stimuli influence soybean development and are highly variable throughout the plant, and among organs and tissues, with respect to the environment [10,13]. Genistein, daidzein, and glycitein, the known soybean isoflavones, are synthesized by a branch of the phenylpropanoid pathway. This extended metabolic route is also involved in the synthesis of other important compounds in plants, such as tannins, lignins, lignans, anthocyanins, flavones, flavonols, and the soybean phytoalexins, glyceollins, which are pterocarpans that possess antimicrobial activities [14]. The precursor in the pathway is the amino acid L-phenylalanine, which in the initial step is stripped of its amine group to produce cinnamic acid catalyzed by phenylalanine ammonia lyase (PAL). In the second and third reactions, cinnamate 4-hydroxylase (C4H) and 4-coumarate CoA ligase (4CL) convert cinnamic acid into p-coumaryol CoA. The first critical enzyme required for flavonoid synthesis is chalcone synthase (CHS), which is a multigene family in soybean, although not all copies are expressed in seeds at detectable levels. Other important enzymes in the pathway for isoflavone synthesis are chalcone isomerase (CHI), which converts chalcones to flavanones, and chalcone reductase (CHR), which is required for daidzein and glycitein formation. However, the enzyme that specifically differentiates isoflavone-producing plant species from those with no isoflavone content is isoflavone synthase (IFS), an endoplasmic reticulum (ER)-associated cytochrome P450 monooxygenase, that catalyzes 2,3-aryl ring migration of flavanones to their corresponding isoflavones [15,16,17]. In the soybean genome, IFS is present in two copies, IFS1 and IFS2 that differ by several amino acids. Both enzymes convert naringenin and liquiritigenin to genistein and daidzein, respectively. Despite their homology, IFS1 and IFS2 are differentially regulated at the transcriptional level. For instance, although both proteins contribute to the isoflavone content in the seed [18], expression of IFS2 increases at advanced stages of seed development, whereas IFS1 transcription remains relatively constant [10,12]. Moreover, only IFS2 is induced in soybean hypocotyls and transgenic roots in response to pathogen attack [19].
Three types of fatty acid metabolic enzymes, namely stearoyl-acyl carrier protein-desaturases (encoded by GmSACPD genes), omega-6 desaturases (GmFAD2 genes), and omega-3 desaturases (GmFAD3 genes), largely determine the relative degree of unsaturated fatty acids and the contents of the C18 fatty acids stearate (18:0), oleate (18:1), linoleate (18:2), and linolenate (18:3) in vegetative and seed lipids. The Δ9 stearoyl-acyl carrier protein-desaturases are soluble enzymes localized to the stroma fraction of plastids that insert the first double bond to stearoyl-ACP (18:0-ACP) to produce oleoyl (18:1Δ9)-ACP. Thus, Δ9-stearoyl-ACP-desaturases perform a crucial step in C18 fatty acid biosynthesis because perturbation of SACPD expression and/or enzyme activity may modulate the relative cellular content of both stearate and oleate. Three alleles of SACPD have been identified and characterized from soybean. The soybean ω-6 oleate fatty acid desaturases (FAD2 genes) are microsomal enzymes that initiate the primary route of polyunsaturated lipid biosynthesis by catalyzing the first extraplastidal desaturation to convert 18:1 esterified fatty acids to phosphatidylcholine to α-18:2 fatty acids [20]. Four ω-6 desaturase genes comprise the soybean FAD2 gene family [21], consisting of GmFAD2-1 and GmFAD2-2 and their alleles [20,22,23], GmFAD2-3 [24], and GmFAD6 [20,22].
Radiation breeding is an effective method of inducing mutations in seeds and other propagative materials, such as pollen, whole plants, or tissue-cultured calli [25,26], and directly produces mutant varieties [27]. With regard to agronomic traits, gamma radiation frequently alters flowering, maturation traits, seed coat color, chloroplast number, and biomass yield in soybean [28,29]. Gamma radiation directly produces mutant varieties without the need to go through the otherwise lengthy and laborious process of conventional breeding. For that reason, mutation breeding has been used in various crops and ornamental plants, and has proven to be a promising means of producing new genetic variants [30,31].
Previously, we successfully constructed a soybean mutant diversity pool (MDP), which exhibited dramatic changes in agronomic and morphological traits [30]. In the present study, gamma-irradiated mutants were selected by screening for changes in isoflavone and fatty acid content using MDP lines. In addition, the relative expression level of genes associated with isoflavone and fatty acid biosynthesis was evaluated.

2. Results

2.1. Seed Isoflavone Content of 208 Soybean MDP Lines

First, we assessed the total isoflavone content (TIC) in seeds of the 208 MDP lines. The TIC ranged from 0.88 mg/g (KAS360-22) to 7.12 mg/g (DB-088) with an average of 3.08 mg/g for the 208 MDP lines. With regard to TIC, the seven wild-type cultivars were ranked, in descending order, as DP (4.95 mg/g), PD (2.31 mg/g), BS (2.26 mg/g), 94Seori (1.97 mg/g), HK (1.85 mg/g), DB (1.03 mg/g), and KAS360-22 (0.88 mg/g) (Supplementary Table S1, Figure 1B). For the 60 DP-mutants and 64 DB-mutants, the TIC ranged from 1.47 to 7.08 mg/g and from 1.20 to 7.12 mg/g, respectively. The TIC for the three mutant populations derived from PD, BS, and HK were relatively lower and ranged from 1.14 to 4.07 mg/g, 1.05 to 1.51 mg/g, and 1.08 to 3.21 mg/g, respectively. Detailed information on the isoflavone contents for each mutant-population is summarized in Table 1. We selected nine MDP lines that exhibited either increased isoflavone content (IIC: DB-088, 7.12 mg/g; DP-084, 7.08 mg/g; and HK-17, 3.21 mg/g) or decreased isoflavone content (DIC: DB-064, 1.34 mg/g; DP-093, 2.51 mg/g; and HK25-165, 1.18 mg/g), and the corresponding wild-type cultivars (DB, DP, and HK) for gene expression analysis. Also, variations of the isoflavone content and fatty acid content levels of the selected mutants are shown in Figure 1A.

2.2. Seed Fatty Acid Content of 208 Soybean Lines

For investigation of the fatty acid composition of the seeds, we measured the concentrations of five saturated/unsaturated fatty acids, comprising palmitic (16:0), stearic (18:0), oleic (18:1), linoleic (18:2), and linolenic (18:3) acids, by GC-MS analysis. The proportion of palmitic acid ranged from 12.42% to 21% of the total fatty acid (TFA) composition in the 208 MDP lines (Supplementary Table S2). Interestingly, KAS360-22 showed a comparatively high percentage stearic acid content (22.99%), which was approximately 11-fold higher than the average percentage for the 208 MDP lines (2.08%). However, the KAS360-22-W mutant derived from KAS360-22 exhibited a similar percentage (3.66%) to that of other mutant lines. Linoleic acid, a major fatty acid compound in soybean seeds, constituted more than 50% of the TFA composition in the 208 MDP lines. The proportion of oleic acid in the 208 MDP lines ranged from 0.38% to 15.43% in the DB-mutant population, 1.54% to 19.83% in the DP-mutant population, and 0.41% to 24.66% in the HK-mutant population (Table 2). With regard to the percentage oleic acid, the wild-type cultivars were ranked, in descending order, as HK (18.52%), PD (16.95%), 94Seori (14.98%), DP (6.05%), DB (5.83%), BS (4.57%), and KAS360-22 (2.25%) (Figure 1C). The proportion of oleic acid differed significantly among the 208 MDP lines and thus represents variation useful for genetic engineering. Therefore, for further genetic analyses, we focused on the change in oleic acid content among the MDP lines. After screening the oleic acid content among the 208 MDP lines, we selected six mutants that exhibited either increased oleic acid proportions (IOC: DB-075, 15.43%; DP-056, 19.68%; and HK-30, 24.66%) or decreased oleic acid proportions (DOC: DB-041, 0.38%; DP-184, 3.81%; and HK-37, 0.41%), and the corresponding wild-type cultivars (DB, DP, and HK).

2.3. Differential Expression of Isoflavone Biosynthesis Genes during Seed Development

In the present study, we analyzed the expression patterns of 12 genes involved in isoflavone biosynthesis in soybean using qRT-PCR (Figure 2). The IIC (DB-088, DP-084, and HK-17) and DIC mutants (DB-064, DP-093, and HK25-165), and the corresponding wild-type cultivars (DB, DP, and HK) were selected for analysis of the differential expression of isoflavone biosynthetic genes at seed developmental stages 1, 2, and 3. The expression level of upstream genes in the flavanone biosynthesis pathway, comprising three genes acting upstream of CHI (4CL, CHS1, and CHS7), CHI1A, and two genes leading to the synthesis of isoflavones or flavonols from flavanone (IFS1 and IFS2, F3H, MaT1, MaT3, MaT7, UGT1, and UGT9), was evaluated in the nine selected lines. The 4CL gene showed significant differences in expression level between IIC and DIC lines, which was consistently higher at stage 1 in the DB-, DP-, and HK-mutants. CHS1 and CHS7 exhibited similar expression patterns between the DB- and DP-mutants, whereas the HK-mutants showed a different expression pattern compared with those of the DB- and DP-mutants. CHI1A, IFS1, and IFS2 showed differences in stages and expression pattern depending on the individual or wild-type cultivar. For example, the IIC mutant lines DB-088, DP-084, and HK-17 showed differential expression of IFS2 at stage 2, stage 1, and stage 3, respectively, based on the maximum relative expression values. These results indicated that IFS2 was involved in isoflavone accumulation, but the regulation timing differed owing to the genetic background. In contrast, MaT7, an isoflavone-specific malonyltransferase, consistently showed higher relative expression levels at stage 1 in the three IIC lines, and thus may be the main contributor to isoflavone biosynthesis in soybean seeds. Interestingly, UGT1 expression was not detected in HK-mutants and HK.

2.4. Differential Expression of Fatty Acid Biosynthesis Gene during Seed Development

Based on the oleic acid content in seeds of the MDP lines, we analyzed the differential expression of fatty acid biosynthetic genes at three seed developmental stages. The qRT-PCR analyses were conducted with primers that targeted genes up- or downstream of oleic acid (18:1) biosynthesis (Figure 3). To confirm the interaction between oleic acid content and related gene expression, the expression of 12 major structural genes, including FAD2-1A, FAD2-1B, FAD2-2A, FAD2-2B, FAD2-2C, FAD2-2D, SACPD-A, SACPD-B, SACPD-C, ACT1A, ACT1B, and FAD6, was analyzed at three seed developmental stages for nine MDP lines. The IOC mutant DB-075 showed consistently higher expression levels of expression stage (ES) I genes than those of DB and DB-041. In contrast, ES III genes showed a similar or decreased expression level in DB-075 compared with those of DB and DB-041 at stage 3. At stage 2, DB-075 showed consistently high expression levels for 11 major structural genes, except FAD6, compared with those of DB and DB-041. Expression patterns differed somewhat between the IOC mutant DP-056 and the DOC mutant DP-184; for example, FAD2-2C and SACPD-C expression differed significantly at stage 2, but no significant differences between DP-056 and DP-184 were observed for other genes analyzed. The HK-30 mutant showed higher relative expression levels of ES III genes compared with those of HK and HK-37; these expression patterns differed from those observed for the DB- and DP-mutants.

3. Discussion

In this study, we investigated the metabolic properties, represented by the contents of isoflavones and five fatty acids, of the seeds of 208 MDP lines. The TIC ranged from 0.88 mg/g (KAS360-22) to 7.12 mg/g (DB-088) and averaged 3.08 mg/g. The TIC of the seven wild-type cultivars were ranked, in descending order, as DP (4.95 mg/g), PD (2.31 mg/g), BS (2.26 mg/g), 94Seori (1.97 mg/g), HK (1.85 mg/g), DB (1.03 mg/g), and KAS360-22 (0.88 mg/g) (Table 1, Supplementary Table S1). A previous study reported that, among five Korean soybean cultivars analyzed, ‘Saegeum’ showed the highest TIC (1.22 mg/g) in the seeds [32]. This content is similar to that of DB observed in the present study and was 5.9-times lower than that of the DB-088 line. These results indicate that the DB-088 line shows the highest seed TIC reported for the Korean region. Park et al. [33] reported that seeds of the cultivar ‘Bosek’ showed the highest TIC (3.43 mg/g) among 106 Korean soybean cultivars surveyed, but this content was lower than that of the DB-088 line. In addition, Choi et al. [34] reported a TIC that ranged from 1.827 to 5.777 mg/g in 49 Korean soybean accessions, and Azam et al. [35] reported a TIC that ranged from 0.745 to 5.253 mg/g in 1168 China soybean accessions. These results confirmed that the TIC of DB-088 was approximately 1.2- to 9.5-times higher than that of other accessions.
Soybean oil with high oleic acid content has the potential to improve the overall oil composition profile because oleic acid is an important component of soybean oil. Therefore, many researchers have attempted to develop soybeans with the desired fatty acid phenotypes through breeding [36]. In the present study, we measured the contents in the seed of five saturated/unsaturated fatty acids, consisting of palmitic (16:0), stearic (18:0), oleic (18:1), linoleic (18:2), and linolenic (18:3) acids, by GC-MS analysis. The proportion of oleic acid in 208 MDP lines ranged from 0.38% to 15.43% in the DB-mutant population, 1.54% to 19.83% in the DP-mutant population, and 0.41% to 24.66% in the HK-mutant population. The oleic acid content ranged from 0.38% to 24.66% and averaged 11.02% (Table 2, Supplementary Table S2). Our MDP lines identified a variety of oleic acids. Therefore, it was confirmed that the oleic acid of the mutant lines increased in a range from approximately 1.3 to 5.3-times compared to the wild-type. In previous studies, Rahman et al. [37] conducted X-ray irradiation of soybean seeds and selected the M23 line, which showed an approximately two-fold increase in oleic acid content in the seed compared with that of the wild-type. This mutant line showed 46.1% oleic acid content in the seed, but also exhibited a parallel decrease in linoleic acid content. In another study, the increase in oleic acid content of the transgenic varieties with FAD-2 downregulation was approximately 4-times higher than that of soybean wild-type [38]. As an extreme case, Lee et al. [39] developed the cultivar ‘Hosim’ from the cross 17D × S08-14788, which shows a high proportion of oleic acid (79%) in the seed. These results confirmed that the MDP lines were approximately 3.2 to 207-times lower than that of ‘Hosim’ cultivar. In addition, Choung [40] analyzed the fatty acid composition of 563 soybean accessions and reported the proportions 13.3–14.3% oleic acid, 38.2–61.5% linoleic acid, and 5.0–11.5% linolenic acid. Similarly, Song et al. [41] reported the oleic acid content of 379 soybean accessions (128 Korean, 214 Chinese, and 37 USA accessions), which ranged from 15.3% to 56.0% (mean 28.2%). Commodity soybean oil typically contains 20–25% monounsaturated oleic acid but recently lines producing higher proportions of oleic acid have entered commercial production [33]. A variety of genetic and biotechnological approaches to increase the oleic acid content in soybean seeds are currently being explored [42].
We selected six mutant lines that showed altered TIC for analysis of the expression patterns of isoflavone biosynthetic genes. The genes CHI1A, IFS1, and IFS2 showed differences in stages and expression patterns depending on the individual or the wild-type cultivar. The genes IFS1 and IFS2 encode proteins that differ by 14 amino acids. The IFS1 and IFS2 proteins convert the flavanones naringenin and liquiritigenin to the corresponding isoflavones genistein and daidzein, respectively [16]. The expression patterns of the IFS gene family were consistent with isoflavonoid accumulation in the seed, which indicated there is a close relationship between expression of these genes and metabolite accumulation in the seed [9]. In this study, the IFS2 relative expression level was significantly higher in mutant lines that exhibited increased isoflavone accumulation compared with lines that showed decreased isoflavone content in the seed. Also, the MaT7 gene showed consistently higher expression levels at stage 1 in the three IIC lines and may be the predominant contributor to isoflavone biosynthesis in the seeds, while the MaT3 gene only increased in DP-084 during stage 2. Previous studies have shown the diverse subcellular localization of MaTs, such as the cytosol [43], the ER [44], and the nucleus and cytoplasm [14]. In addition, the soybean genome contains a number of MaT homologues because of the genome’s paleopolyploid nature [45]. The primary subcellular localization of the protein encoded by MaT1 in the ER and the protein encoded by MaT3 in the cytosol might reflect their functions in malonylation following the synthesis of isoflavone glucosides. Similar to other flavonoids, isoflavone glucosides are considered to be synthesized on the cytosolic side of the ER and are subsequently modified by MaTs [46]. We found that expression patterns of the major isoflavone structural genes in selected mutants were dependent on the seed developmental stages and thus were also cultivar-specific. So, in order to more clearly define the gene expression patterns, we displayed these relationships in the phenylalanine pathway (Figure 4).
The expression level of 12 major structural genes in fatty acid biosynthesis was analyzed in nine selected lines, comprising three IOC mutant lines (DB-075, DP-056, and HK-30), three DOC mutant lines (DB-041, DP-184, and HK-37), and the corresponding cultivars (DB, DP, and HK) at seed developmental stages 1 to 3. In the study, the 12 genes were classifiable into two groups: (i) ES I: highly accumulated in stage 1, consisting of ACT1A, ACT1B, FAD2-2A, FAD2-2B, FAD2-2C, and FAD6; and (ii) ES III: highly accumulated in stage 3, comprising SACPD-A, SACPD-B, SACPD-C, FAD2-1A, FAD2-1B, and FAD2-2D (Figure 5). Previously, we developed mutant soybean populations by gamma irradiation of the cultivars ‘Danbaek’ and ‘Daepung’ and evaluated the linolenic acid content of the seed in 78 and 154 M9 mutant progenies. The selected mutant line showed 33.9%–67.7% higher linolenic acid content compared with that of the wild-type cultivars, and increased expression levels of the fatty acid desaturation enzyme (FAD) gene during seed development. In addition, the linolenic acid content was associated with a significant increase in the expression levels of FAD3C and FAD3D in the ER [47]. According to soybean RNA-seq data, FAD gene expression patterns vary depending on the tissue and developmental stage [48]. In the present study, the microsomal omega-6 desaturase genes FAD2-1A and FAD2-1B showed increased expression levels among the selected DB-, DP-, and HK-mutants during seed development compared with those of the wild-type cultivars. Similarly, Robert et al [49] reported the gene expression of eight FAD genes in three seed development stages. Among genes, FAD2-1A and FAD2-1B increased their expression from 22 DAF to 28 DAF stages. However, it decreased at the 35 DAF stages. In addition, the two microsomal FAD2-1 desaturases FAD2-1A and FAD2-1B were mainly expressed in developing seeds [23]. Thus, FAD2-1A and FAD2-1B are considered to play an important role in controlling the oleic acid level in soybean seed development. The seed-specific expression of FAD2-1 and FAD2-2 plays a role in desaturation of 18:1 in soybean seeds [20]. The majority of acyl carrier proteins and fatty acid elongases were expressed at a higher level during the early maturity stage (stage 2) and subsequently the expression level declined at stage 3. Similarly, soybean also reported that FAD expression patterns varied depending on the tissue and developmental stage [49] and in Arabidopsis, the expression level of fatty acid-associated genes is elevated during seed maturation [50].

4. Materials and Methods

4.1. Plant Materials

Soybean seeds were irradiated with gamma-rays generated using a 60Co gamma-irradiator (150 TBq capacity; ACEL, Ottawa, ON, Canada) at the Korea Atomic Energy Research Institute (KAERI). The irradiated seeds and control seeds were sown at the Breeding Research Farm at the KAERI. Seeds from the M1 plants were harvested individually and carried forward to the M2 generation. Genetically fixed mutant lines (M12 generation) were selected for desirable agronomic traits, such as grain yield, growth type, climate adaptability, and morphological phenotypes (Figure 6). For analysis of isoflavones and fatty acids, 208 soybean MDP lines, consisting of one landrace (KAS360-22) and six representative Korean soybean cultivars (‘94Seori’, ‘Bangsa’ [BS], ‘Paldal’ [PD], ‘Danbaek’ [DB], ‘Daepung’ [DP], and ‘Hwangkeum’ [HK]), were cultivated from 2008 [30]. Ultimately, we selected 15 MDP lines for gene expression analysis: six isoflavone-altered lines, comprising two DB-mutants (DB-088 and DB-064), two DP-mutants (DP-084 and DP-093), and two HK-mutants (HK-17 and HK25-165); six fatty acid-altered lines, comprising two DB-mutants (DB-075 and DB-041), two DP-mutants (DP-056 and DP-184), and two HK-mutants (HK-30 and HK-37); and the corresponding wild-type cultivars (DB, DP, and HK).

4.2. Isoflavone Extraction and Quantification

Lyophilized whole soybean seeds were finely ground with a mortar. Each ground sample (7 mg) was immersed in 1 mL of 58% (v/v) aqueous acetonitrile. The mixture was sonicated for 30 min and centrifuged at 13,000 rpm for 5 min, after which the supernatant was retained. The pellet was resuspended in an equal volume of solvent, and both retained supernatants were combined and diluted with distilled water. The extract volume was adjusted to 4 mL for each extraction from a 7-mg freeze-dried seed sample. The diluted extracts were filtered through a 0.45-μm syringe filter (Futecs Co., Ltd., Daejeon, Korea) and used for reversed-phase high performance liquid chromatography (HPLC) analysis.
Extracts were analyzed using reversed-phase HPLC (Waters 2695 Alliance HPLC; Waters Inc., Milford, MA, USA) with an octadecylsilane column (Prontosil 120–5-C18-ace-EPS 5.0 μm (250 × 4.6 mm; Bischoff, Leonberg, Germany). The flow rate of the mobile phase was 1.0 mL/min and the sample injection volume was 5 μL. The mobile phase was a combination of (A) water with 0.1% formic acid and (B) acetonitrile with 0.1% formic acid. Gradient elution was performed by adding 15% of solvent B at the initial running time and increasing the concentration to 34% over 60 min. Peaks were monitored at 254 nm using a Waters 996 photodiode array detector (Waters Inc.). Twelve isoflavone standards were purchased from Sigma-Aldrich (St. Louis, MO, USA) and used for quantification of the isoflavones from soybean seeds in the HPLC analysis.

4.3. Fatty Acid Analysis

For gas chromatography–mass spectrometry (GC-MS) analysis, fatty acids were extracted as described by Ryu et al. [51] with the following modifications. A powdered freeze-dried seed sample (0.5 g) was extracted in 1 mL n-hexane for 4 h, then 0.1 mL of 2 N potassium hydroxide in methanol was added. After centrifugation for 5 min at 3000× g, the collected supernatant was filtered using a 0.45-μm syringe filter. The fatty acid composition was analyzed using a GC-MS (Plus-2010, Shimadzu, Japan) instrument equipped with a HP-88 capillary column (J&W Scientific, Folsom, CA, USA, 60 m × 0.25 mm × 0.25 m) under the following conditions: ionization voltage, 70 eV; mass scan range, 50–450 mass units; injector temperature, 230 °C; detector temperature, 230 °C; injection volume, 1 L; split ratio, 1:30; carrier gas, helium; and flow rate, 1.7 mL/min. The column temperature program specified an isothermal temperature of 40 °C for 5 min increasing to 180 °C at the rate of 5 °C/min, then a subsequent increase to 28 °C at the rate of 1 °C/min. We identified the substances present in the extracts in accordance with their retention time and with reference to a mass spectral database (NIST 62 Library).

4.4. RNA Isolation and cDNA Synthesis

Immature seeds from the 15 selected MDP lines were collected at stage 1 (length 4 to <7 mm; R5e, DAF20), stage 2 (7–10 mm; R5L, DAF30), and stage 3 (11–14 mm; R6, DAF40) in accordance with a previous report [13] with some modifications (Supplementary Figure S1). Total RNA was isolated from seeds with TRIzol Reagent in accordance with the manufacturer’s protocol (Invitrogen, Carlsbad, CA, USA). The RNA concentration and quality were measured using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) prior to DNase digestion. For each sample, 15 μg total RNA was digested in a volume of 20 μL using the Invitrogen DNA-free Kit (Life Technologies, Grand Island, NY, USA) to remove genomic DNA contamination following the manufacturer’s instructions. After DNase I digestion, the RNA concentration was determined using a NanoDrop ND-1000 spectrophotometer. First-strand cDNA synthesis was performed using 1 μg DNase-treated total RNA in a 20-μL reaction using the SuperScript III First-Strand Synthesis SuperMix kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions.

4.5. Quantitative Real-Time PCR

Transcript levels of selected genes involved in the phenylpropanoid and fatty acid pathway were quantified by quantitative real-time PCR (qRT-PCR) analysis. Relative expression level was calculated using the 2−ΔΔCt comparative threshold method [52]. Primer specificity was confirmed by blasting each primer sequence against soybean genome sequences lodged in the Phytozome database (http://www.phytozome.net/, last accessed on 19 May 2021) using the BLASTN algorithm. The qRT-PCR reactions were performed in 96-well plates using the CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The iTaq Universal SYBR Green Supermix (Bio-Rad) was used for real-time cDNA quantification. A 10 pmol primer concentration and 3 μL of prepared cDNA were used in a final volume of 20 μL per reaction. The PCR protocol was as follows: 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s, 50 °C for 15 s, and 72 °C for 30 s. The results were normalized to the constitutive expression level of ELF1B, which was selected as an internal reference gene owing to its expression stability. Gene-specific primers used for qRT-PCR analyses are listed in Supplementary Tables S3 and S4.

5. Conclusions

In the present study, we analyzed the metabolic properties, including the isoflavones and five fatty acid contents, of 208 MDP lines. The genetically fixed mutant lines that showed significantly increased or decreased isoflavone and fatty acid contents were selected from the DB-, DP-, and HK- mutant population. The lines were selected to analyze the differential expression of isoflavones and fatty acid biosynthetic genes at three seed developmental stages. Isoflavone biosynthetic genes, including CHI1A, IFS1, and IFS2, showed differences in stages and expression patterns depending on the individual or wild-type cultivar, whereas MaT7 showed consistently higher expression levels in seeds at stage 1. The fatty acid biosynthetic genes were classifiable into two groups based on the developmental stages of the seeds. Our results can serve as a foundation for future functional analysis of the regulatory genes involved in the isoflavone and fatty acid biosynthetic pathways.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10061037/s1, Figure S1: Soybean seed developmental stages. Stage 1, length 4 to <7 mm; stage 2, length 7–10 mm; and stage 3, length 11–14 mm, Table S1: Total isoflavone content in the seeds of 208 soybean MDP lines, Table S2: Fatty acid content in the seeds of 208 soybean MDP lines, Table S3: Primer pairs used for quantitative real-time PCR analysis (isoflavone biosynthesis genes), Table S4: Primer pairs used for quantitative real-time PCR analysis (fatty acid biosynthesis genes).

Author Contributions

Conceptualization, D.-G.K. and J.-I.L.; methodology, D.-G.K., J.-I.L. and S.-J.K.; software, N.-N.H.; validation, J.-B.K., C.-H.B. and S.-J.K.; investigation, J.-M.K.; resources, S.-H.K.; data curation, Y.-J.L. and S.-H.E.; writing—original draft preparation, D.-G.K.; writing—review and editing, J.-I.L. and S.-J.K.; visualization, S.-H.K. and J.-B.K.; supervision, C.-H.B. and S.-J.K.; funding acquisition, J.-B.K. and S.-J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the research program of KAERI, Republic of Korea and the Radiation Technology R&D Program (NRF-2017M2A2A6A05018538) through the National Research Foundation of Korea funded by the Ministry of Science and ICT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gunstone, F.; Norris, F. The Biosynthesis and Metabolism of Fatty Acids and Lipids. In Lipids in Foods: Chemistry Biochemistry and Technology; Pergamon: New York, NY, USA, 1983. [Google Scholar]
  2. Wilson, R.F. Seed composition. Soybeans Improv. Prod. Uses 2004, 16, 621–677. [Google Scholar]
  3. Fehr, W.; Caviness, C.; Burmood, D.; Pennington, J. Stage of development descriptions for soybeans, Glycine Max (L.) Merrill 1. Crop Sci. 1971, 11, 929–931. [Google Scholar] [CrossRef]
  4. Dhaubhadel, S.; Marsolais, F. Transcriptomics of legume seed: Soybean a model grain legume. In Seed Development: OMICS Technologies toward Improvement of Seed Quality and Crop Yield; Springer: Dordrecht, The Netherlands, 2012; pp. 129–142. [Google Scholar]
  5. Brevedan, R.; Egli, D. Short periods of water stress during seed filling, leaf senescence, and yield of soybean. Crop Sci. 2003, 43, 2083–2088. [Google Scholar] [CrossRef]
  6. Gayler, K.R.; Sykes, G.E. Effects of nutritional stress on the storage proteins of soybeans. Plant Physiol. 1985, 78, 582–585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lozovaya, V.V.; Lygin, A.V.; Ulanov, A.V.; Nelson, R.L.; Daydé, J.; Widholm, J.M. Effect of temperature and soil moisture status during seed development on soybean seed isoflavone concentration and composition. Crop Sci. 2005, 45, 1934–1940. [Google Scholar] [CrossRef] [Green Version]
  8. Pedersen, P. Soybean growth stages. In Soybean Growth and Development; Iowa State University Cooperative Extension Publication PM: Ames, IA, USA, 1945. [Google Scholar]
  9. McWilliams, D.; Berglund, D.R.; Endres, G. Soybean Growth and Management Quick Guide; North Dakota State University: Fargo, ND, USA, 1999. [Google Scholar]
  10. Dhaubhadel, S.; Gijzen, M.; Moy, P.; Farhangkhoee, M. Transcriptome analysis reveals a critical role of CHS7 and CHS8 genes for isoflavonoid synthesis in soybean seeds. Plant Physiol. 2007, 143, 326–338. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, J.A.; Chung, I.M. Change in isoflavone concentration of soybean (Glycine max L.) seeds at different growth stages. J. Sci. Food Agric. 2007, 87, 496–503. [Google Scholar] [CrossRef]
  12. Dhaubhadel, S.; McGarvey, B.D.; Williams, R.; Gijzen, M. Isoflavonoid biosynthesis and accumulation in developing soybean seeds. Plant Mol. Biol. 2003, 53, 733–743. [Google Scholar] [CrossRef]
  13. Devi, M.A.; Kumar, G.; Giridhar, P. Effect of biotic and abiotic elicitors on isoflavone biosynthesis during seed development and in suspension cultures of soybean (Glycine max L.). 3 Biotech 2020, 10, 98. [Google Scholar] [CrossRef] [PubMed]
  14. Yu, O.; McGonigle, B. Metabolic engineering of isoflavone biosynthesis. Adv. Agron. 2005, 86, 147–190. [Google Scholar]
  15. Akashi, T.; Aoki, T.; Ayabe, S.-i. Cloning and functional expression of a cytochrome P450 cDNA encoding 2-hydroxyisoflavanone synthase involved in biosynthesis of the isoflavonoid skeleton in licorice. Plant Physiol. 1999, 121, 821–828. [Google Scholar] [CrossRef] [Green Version]
  16. Jung, W.; Yu, O.; Lau, S.-M.C.; O’Keefe, D.P.; Odell, J.; Fader, G.; McGonigle, B. Identification and expression of isoflavone synthase, the key enzyme for biosynthesis of isoflavones in legumes. Nat. Biotechnol. 2000, 18, 208–212. [Google Scholar] [CrossRef]
  17. Steele, C.L.; Gijzen, M.; Qutob, D.; Dixon, R.A. Molecular characterization of the enzyme catalyzing the aryl migration reaction of isoflavonoid biosynthesis in soybean. Arch. Biochem. Biophys. 1999, 367, 146–150. [Google Scholar] [CrossRef]
  18. Cheng, H.; Yu, O.; Yu, D. Polymorphisms of IFS1 and IFS2 gene are associated with isoflavone concentrations in soybean seeds. Plant Sci. 2008, 175, 505–512. [Google Scholar] [CrossRef]
  19. Subramanian, S.; Hu, X.; Lu, G.; Odelland, J.T.; Yu, O. The promoters of two isoflavone synthase genes respond differentially to nodulation and defense signals in transgenic soybean roots. Plant Mol. Biol. 2004, 54, 623–639. [Google Scholar] [CrossRef]
  20. Heppard, E.P.; Kinney, A.J.; Stecca, K.L.; Miao, G.-H. Developmental and growth temperature regulation of two different microsomal [omega]-6 desaturase genes in soybeans. Plant Physiol. 1996, 110, 311–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Schlueter, J.A.; Vasylenko-Sanders, I.F.; Deshpande, S.; Yi, J.; Siegfried, M.; Roe, B.A.; Schlueter, S.D.; Scheffler, B.E.; Shoemaker, R.C. The FAD2 gene family of soybean: Insights into the structural and functional divergence of a paleopolyploid genome. Crop Sci. 2007, 47, S-14–s-26. [Google Scholar] [CrossRef]
  22. Bachlava, E.; Dewey, R.E.; Burton, J.W.; Cardinal, A.J. Mapping and comparison of quantitative trait loci for oleic acid seed content in two segregating soybean populations. Crop Sci. 2009, 49, 433–442. [Google Scholar] [CrossRef] [Green Version]
  23. Tang, G.Q.; Novitzky, W.P.; Carol Griffin, H.; Huber, S.C.; Dewey, R.E. Oleate desaturase enzymes of soybean: Evidence of regulation through differential stability and phosphorylation. Plant J. 2005, 44, 433–446. [Google Scholar] [CrossRef]
  24. Li, L.; Wang, X.; Gai, J.; Yu, D. Molecular cloning and characterization of a novel microsomal oleate desaturase gene from soybean. J. Plant Physiol. 2007, 164, 1516–1526. [Google Scholar] [CrossRef] [PubMed]
  25. Ali, H.; Ghori, Z.; Sheikh, S.; Gul, A. Effects of gamma radiation on crop production. In Crop Production and Global Environmental Issues; Springer: Cham, Switzerland, 2015; pp. 27–78. [Google Scholar]
  26. Oladosu, Y.; Rafii, M.Y.; Abdullah, N.; Hussin, G.; Ramli, A.; Rahim, H.A.; Miah, G.; Usman, M. Principle and application of plant mutagenesis in crop improvement: A review. Biotechnol. Biotechnol. Equip. 2016, 30, 1–16. [Google Scholar] [CrossRef] [Green Version]
  27. Beyaz, R.; Yildiz, M. The use of gamma irradiation in plant mutation breeding. In Plant Engineering; Juri’c, S., Ed.; IntechOpen: London, UK, 2017; pp. 33–46. [Google Scholar]
  28. Ha, B.-K.; Lee, K.J.; Velusamy, V.; Kim, J.-B.; Kim, S.H.; Ahn, J.-W.; Kang, S.-Y.; Kim, D.S. Improvement of soybean through radiation-induced mutation breeding techniques in Korea. Plant Genet. Resour. 2014, 12, S54–S57. [Google Scholar] [CrossRef]
  29. Lee, K.J.; Kim, J.-B.; Kim, S.H.; Ha, B.-K.; Lee, B.-M.; Kang, S.-Y.; Kim, D.S. Alteration of seed storage protein composition in soybean [Glycine max (L.) Merrill] mutant lines induced by γ-irradiation mutagenesis. J. Agric. Food Chem. 2011, 59, 12405–12410. [Google Scholar] [CrossRef] [PubMed]
  30. Kim, D.-G.; Lyu, J.I.; Lee, M.-K.; Kim, J.M.; Hung, N.N.; Hong, M.J.; Kim, J.-B.; Bae, C.-H.; Kwon, S.-J. Construction of Soybean Mutant Diversity Pool (MDP) Lines and an Analysis of Their Genetic Relationships and Associations Using TRAP Markers. Agronomy 2020, 10, 253. [Google Scholar] [CrossRef] [Green Version]
  31. Song, H.; Kang, S. Application of natural variation and induced mutation in breeding and functional genomics: Papers for International Symposium; Current Status and Future of Plant Mutation Breeding. Korean J. Breed. Sci 2003, 35, 24–34. [Google Scholar]
  32. Yun, G.; Hwang, S.; Hong, S.; Hong, U.; Yun, H.; NICS, R.; Kim, H. Comparison of Growth Characteristics and Isoflavones Contents of Major Soybean Cultivars in Mid Part of Korea. J. Korean Soc. Int. Agric. 2016, 28, 496–504. [Google Scholar] [CrossRef]
  33. Park, K.-H.; Piao, X.-M.; Jang, E.-K.; Yoo, Y.-E.; Hwang, T.-Y.; Kim, S.-L.; Jong, J.-H.; Shin, H.-M.; Kim, H.-S. Variation of Isoflavone Contents in Korean Soybean Cultivars Released from 1913 to 2006. Korean J. Breed. Sci. 2012, 44, 149–159. [Google Scholar]
  34. Choi, Y.-M.; Yoon, H.; Lee, S.; Ko, H.-C.; Shin, M.-J.; Lee, M.-C.; Oh, S.; Desta, K.T. Comparison of Isoflavone Composition and Content in Seeds of Soybean (Glycine max (L.) Merrill) Germplasms with Different Seed Coat Colors and Days to Maturity. Korean J. Plant Resour. 2020, 33, 558–577. [Google Scholar]
  35. Azam, M.; Zhang, S.; Abdelghany, A.M.; Shaibu, A.S.; Feng, Y.; Li, Y.; Tian, Y.; Hong, H.; Li, B.; Sun, J. Seed isoflavone profiling of 1168 soybean accessions from major growing ecoregions in China. Food Res. Int. 2020, 130, 108957. [Google Scholar] [CrossRef]
  36. Thapa, R.; Carrero-Colon, M.; Crowe, M.; Gaskin, E.; Hudson, K. Novel FAD2–1A alleles confer an elevated oleic acid phenotype in soybean seeds. Crop Sci. 2016, 56, 226–231. [Google Scholar] [CrossRef]
  37. Rahman, S.M.; Takagi, Y.; Kubota, K.; Miyamoto, K.; Kawakita, T. High oleic acid mutant in soybean induced by X-ray irradiation. Biosci. Biotechnol. Biochem. 1994, 58, 1070–1072. [Google Scholar] [CrossRef]
  38. Spear, J.D.; Fehr, W.R.; Schnebly, S.R. Stability of Oleate Concentration in Soybean Lines Containing the High-Oleate Transgene DP-305423-1. Crop Sci. 2013, 53, 888–893. [Google Scholar] [CrossRef]
  39. Lee, J.-D.; Kim, M.; Kulkarni, K.P.; Song, J.T. Agronomic Traits and Fatty Acid Composition of High–Oleic Acid Cultivar Hosim. Plant Breed. Biotechnol. 2018, 6, 44–50. [Google Scholar] [CrossRef] [Green Version]
  40. Choung, M.-G. Variation of oil contents and fatty acid compositions in Korean soybean germplasms. Korean J. Crop Sci. 2006, 51, 139–145. [Google Scholar]
  41. Song, H.; Jang, E.; Son, E.; Lee, K.; Hwang, S.; Yun, G.; Lee, J.; Kim, H. Evaluation of oil contents and fatty acid compositions and selection of useful resources in soybean germplasm accessions. Korean J. Breed. Sci. 2015, 47, 209–218. [Google Scholar] [CrossRef]
  42. Clemente, T.E.; Cahoon, E.B. Soybean oil: Genetic approaches for modification of functionality and total content. Plant Physiol. 2009, 151, 1030–1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Achnine, L.; Huhman, D.V.; Farag, M.A.; Sumner, L.W.; Blount, J.W.; Dixon, R.A. Genomics-based selection and functional characterization of triterpene glycosyltransferases from the model legume Medicago truncatula. Plant J. 2005, 41, 875–887. [Google Scholar] [CrossRef] [Green Version]
  44. Zhao, J.; Huhman, D.; Shadle, G.; He, X.-Z.; Sumner, L.W.; Tang, Y.; Dixon, R.A. MATE2 mediates vacuolar sequestration of flavonoid glycosides and glycoside malonates in Medicago truncatula. Plant Cell 2011, 23, 1536–1555. [Google Scholar] [CrossRef] [Green Version]
  45. Shoemaker, R.C.; Schlueter, J.; Doyle, J.J. Paleopolyploidy and gene duplication in soybean and other legumes. Curr. Opin. Plant Biol. 2006, 9, 104–109. [Google Scholar] [CrossRef]
  46. Funaki, A.; Waki, T.; Noguchi, A.; Kawai, Y.; Yamashita, S.; Takahashi, S.; Nakayama, T. Identification of a highly specific isoflavone 7-O-glucosyltransferase in the soybean (Glycine max (L.) Merr.). Plant Cell Physiol. 2015, 56, 1512–1520. [Google Scholar] [CrossRef] [Green Version]
  47. Hong, M.J.; Jang, Y.E.; Kim, D.G.; Kim, J.M.; Lee, M.K.; Kim, J.B.; Eom, S.H.; Ha, B.K.; Lyu, J.I.; Kwon, S.J. Selection of mutants with high linolenic acid contents and characterization of fatty acid desaturase 2 and 3 genes during seed development in soybean (Glycine max). J. Sci. Food Agric. 2019, 99, 5384–5391. [Google Scholar] [CrossRef]
  48. Severin, A.J.; Woody, J.L.; Bolon, Y.-T.; Joseph, B.; Diers, B.W.; Farmer, A.D.; Muehlbauer, G.J.; Nelson, R.T.; Grant, D.; Specht, J.E. RNA-Seq Atlas of Glycine max: A guide to the soybean transcriptome. BMC Plant Biol. 2010, 10, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Upchurch, R.G.; Ramirez, M.E. Gene Expression Profiles of Soybeans with Mid-Oleic Acid Seed Phenotype. J. Am. Oil Chem. Soc. 2010, 87, 857–864. [Google Scholar] [CrossRef]
  50. Ruuska, S.A.; Girke, T.; Benning, C.; Ohlrogge, J.B. Contrapuntal networks of gene expression during Arabidopsis seed filling. Plant Cell 2002, 14, 1191–1206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ryu, J.; Kwon, S.-J.; Ahn, J.-W.; Jo, Y.D.; Kim, S.H.; Jeong, S.W.; Lee, M.K.; Kim, J.-B.; Kang, S.-Y. Phytochemicals and antioxidant activity in the kenaf plant (Hibiscus cannabinus L.). J. Plant Biotechnol. 2017, 44, 191–202. [Google Scholar] [CrossRef] [Green Version]
  52. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Plants 10 01037 g001 550
Figure 1. Changes in phytochemical traits (isoflavone and fatty acid contents) of 15 selected MDP lines compared with six Korean cultivars of soybean. (A) Box plots of the phenotypic distributions for three wild-type cultivars and their selected MDP lines (six isoflavone-altered lines comprising DB-088, DB-064, DP-084, DP-093, HK-17, and HK25-165; and six oleic acid-altered lines comprising DB-075, DB-041, DP-056, DP-184, HK-30, and HK-37; represented by gray symbols). The data are presented as log2-based mean values for individual lines. Distribution of the phytochemical traits among the six Korean cultivars is presented for (B) isoflavone (mg/g dry weight) and (C) oleic acid (%) contents.
Figure 1. Changes in phytochemical traits (isoflavone and fatty acid contents) of 15 selected MDP lines compared with six Korean cultivars of soybean. (A) Box plots of the phenotypic distributions for three wild-type cultivars and their selected MDP lines (six isoflavone-altered lines comprising DB-088, DB-064, DP-084, DP-093, HK-17, and HK25-165; and six oleic acid-altered lines comprising DB-075, DB-041, DP-056, DP-184, HK-30, and HK-37; represented by gray symbols). The data are presented as log2-based mean values for individual lines. Distribution of the phytochemical traits among the six Korean cultivars is presented for (B) isoflavone (mg/g dry weight) and (C) oleic acid (%) contents.
Plants 10 01037 g001
Plants 10 01037 g002 550
Figure 2. Relative expression level of isoflavone biosynthesis genes in seeds of soybean ‘Danbaek-mutants’ (Black line, DB; red line, DB-088; blue line, DB-064), ‘Daepung-mutants’ (Black line, DP; red line, DP-084; blue line, DB-093), and ‘Hwangkeum-mutants’ (Black line, HK; red line, HK-17; blue line, HK25-165) MDP lines at three developmental stages. Expression levels were estimated using quantitative real-time PCR analysis. The relative expression level was normalized to the internal reference gene ELF1B.
Figure 2. Relative expression level of isoflavone biosynthesis genes in seeds of soybean ‘Danbaek-mutants’ (Black line, DB; red line, DB-088; blue line, DB-064), ‘Daepung-mutants’ (Black line, DP; red line, DP-084; blue line, DB-093), and ‘Hwangkeum-mutants’ (Black line, HK; red line, HK-17; blue line, HK25-165) MDP lines at three developmental stages. Expression levels were estimated using quantitative real-time PCR analysis. The relative expression level was normalized to the internal reference gene ELF1B.
Plants 10 01037 g002
Plants 10 01037 g003 550
Figure 3. Relative expression level of fatty acid biosynthesis genes in seeds of soybean ‘Danbaek-mutants’ (Black line, DB; red line, DB-075; blue line, DB-041), ‘Daepung-mutants’ (Black line, DP; red line, DP-056; blue line, DB-184), and ‘Hwangkeum-mutants’ (Black line, HK; red line, HK-30; blue line, HK-37) MDP lines at three developmental stages. Expression levels were estimated using quantitative real-time PCR analysis. The relative expression level was normalized to the internal reference gene ELF1B.
Figure 3. Relative expression level of fatty acid biosynthesis genes in seeds of soybean ‘Danbaek-mutants’ (Black line, DB; red line, DB-075; blue line, DB-041), ‘Daepung-mutants’ (Black line, DP; red line, DP-056; blue line, DB-184), and ‘Hwangkeum-mutants’ (Black line, HK; red line, HK-30; blue line, HK-37) MDP lines at three developmental stages. Expression levels were estimated using quantitative real-time PCR analysis. The relative expression level was normalized to the internal reference gene ELF1B.
Plants 10 01037 g003
Plants 10 01037 g004 550
Figure 4. Transcript levels of the 12 major structural genes in the isoflavone biosynthetic pathway. For each gene, the transcript level for six selected mutant lines and three wild-type cultivars analyzed in this study are displayed as a heatmap based on the normalized Ct values.
Figure 4. Transcript levels of the 12 major structural genes in the isoflavone biosynthetic pathway. For each gene, the transcript level for six selected mutant lines and three wild-type cultivars analyzed in this study are displayed as a heatmap based on the normalized Ct values.
Plants 10 01037 g004
Plants 10 01037 g005 550
Figure 5. Transcript levels of the 12 major structural genes in the fatty acid biosynthetic pathway. For each gene, the transcript level for six selected mutant lines and three wild-type cultivars analyzed in this study are displayed as a heatmap based on the normalized Ct values.
Figure 5. Transcript levels of the 12 major structural genes in the fatty acid biosynthetic pathway. For each gene, the transcript level for six selected mutant lines and three wild-type cultivars analyzed in this study are displayed as a heatmap based on the normalized Ct values.
Plants 10 01037 g005
Plants 10 01037 g006 550
Figure 6. Schematic representation of the development of the 208 MDP lines. In total, 1000 irradiated seeds of each of seven soybean cultivars were sown in a research field of KAERI. Seeds from the 1695 M2 plants were harvested and first selected for 523 mutant lines in the M5 generation based on agricultural traits such as grain yield, growth type, and climate adaptability. In a second selection phase in the M12 generation, 208 genetically fixed mutants were selected to eliminate redundant phenotypes from the 523 mutant lines (additional information is described in [30]).
Figure 6. Schematic representation of the development of the 208 MDP lines. In total, 1000 irradiated seeds of each of seven soybean cultivars were sown in a research field of KAERI. Seeds from the 1695 M2 plants were harvested and first selected for 523 mutant lines in the M5 generation based on agricultural traits such as grain yield, growth type, and climate adaptability. In a second selection phase in the M12 generation, 208 genetically fixed mutants were selected to eliminate redundant phenotypes from the 523 mutant lines (additional information is described in [30]).
Plants 10 01037 g006
Table
Table 1. Variation of isoflavone contents (mg/g dry weight) in seeds of 208 soybean MDP lines.
Table 1. Variation of isoflavone contents (mg/g dry weight) in seeds of 208 soybean MDP lines.
LinesValuesTIC x
KAS360-22
(N = 2)		Minimum0.88
Maximum2.02
Mean1.45
SD y0.81
CV (%) z55.59
94seori
(N = 5)		Minimum1.27
Maximum2.42
Mean1.82
SD0.44
CV (%)24.22
Bangsa
(N = 6)		Minimum1.05
Maximum2.26
Mean1.36
SD0.47
CV (%)34.48
Paldal
(N = 16)		Minimum1.14
Maximum4.07
Mean2.45
SD0.91
CV (%)37.31
Danbaek
(N = 65)		Minimum1.03
Maximum7.12
Mean3.43
SD1.46
CV (%)42.55
Daepung
(N = 61)		Minimum1.59
Maximum5.04
Mean3.46
SD1.06
CV (%)30.75
Hwangkeum
(N = 53)		Minimum1.08
Maximum3.21
Mean2.05
SD0.53
CV (%)26.02
z CV: Coefficient of variance, y SD: Standard deviation, x TIC: Total isoflavone content.
Table
Table 2. Variation of fatty acid contents (%) in seeds of 208 soybean MDP lines.
Table 2. Variation of fatty acid contents (%) in seeds of 208 soybean MDP lines.
LinesValuesPA xSA wOA vLA uALA t
KAS360-22
(N = 2)		Minimum15.373.662.2550.595.28
Maximum18.8922.9912.1359.009.84
Mean17.1313.337.1954.807.56
SD y2.4913.676.995.953.22
CV (%) z14.53102.5897.1710.8542.65
94seori
(N = 5)		Minimum15.262.229.5515.266.91
Maximum17.923.8314.9863.909.90
Mean16.803.0511.4560.278.43
SD1.180.652.293.231.13
CV (%)7.0421.3320.005.3713.43
Bangsa
(N = 6)		Minimum17.572.814.1158.839.14
Maximum19.033.399.1664.9611.58
Mean18.263.116.4861.5310.63
SD0.570.232.022.150.88
CV (%)3.157.4631.133.538.32
Paldal
(N = 16)		Minimum14.512.344.3055.517.71
Maximum17.943.7216.9566.1611.02
Mean15.763.0210.2961.879.06
SD0.870.412.972.500.98
CV (%)5.5213.7028.874.0510.81
Danbaek
(N = 65)		Minimum12.460.990.3857.535.97
Maximum20.644.1215.4368.6914.00
Mean16.132.609.3462.799.14
SD1.770.632.792.011.28
CV (%)10.9524.2329.873.2014.05
Daepung
(N = 61)		Minimum13.050.001.5456.515.62
Maximum20.283.2319.8373.4012.10
Mean16.241.4810.0363.898.35
SD1.711.044.554.031.45
CV (%)10.5270.2545.376.3117.40
Hwangkeum
(N = 53)		Minimum12.420.000.4153.111.00
Maximum21.007.4924.6674.7013.28
Mean16.191.2015.0561.406.15
SD1.701.496.034.452.12
CV (%)10.50123.9640.067.2534.38
z CV: Coefficient of variance, y SD: Standard deviation, x PA: Palmitic acid, w SA: Stearic acid, v OA: Oleic acid, u LA: Linoleic acid, t ALA: a-Linolenic acid.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, D.-G.; Lyu, J.-I.; Lim, Y.-J.; Kim, J.-M.; Hung, N.-N.; Eom, S.-H.; Kim, S.-H.; Kim, J.-B.; Bae, C.-H.; Kwon, S.-J. Differential Gene Expression Associated with Altered Isoflavone and Fatty Acid Contents in Soybean Mutant Diversity Pool. Plants 2021, 10, 1037. https://doi.org/10.3390/plants10061037

AMA Style

Kim D-G, Lyu J-I, Lim Y-J, Kim J-M, Hung N-N, Eom S-H, Kim S-H, Kim J-B, Bae C-H, Kwon S-J. Differential Gene Expression Associated with Altered Isoflavone and Fatty Acid Contents in Soybean Mutant Diversity Pool. Plants. 2021; 10(6):1037. https://doi.org/10.3390/plants10061037

Chicago/Turabian Style

Kim, Dong-Gun, Jae-Il Lyu, You-Jin Lim, Jung-Min Kim, Nguyen-Ngoc Hung, Seok-Hyun Eom, Sang-Hoon Kim, Jin-Baek Kim, Chang-Hyu Bae, and Soon-Jae Kwon. 2021. "Differential Gene Expression Associated with Altered Isoflavone and Fatty Acid Contents in Soybean Mutant Diversity Pool" Plants 10, no. 6: 1037. https://doi.org/10.3390/plants10061037

APA Style

Kim, D. -G., Lyu, J. -I., Lim, Y. -J., Kim, J. -M., Hung, N. -N., Eom, S. -H., Kim, S. -H., Kim, J. -B., Bae, C. -H., & Kwon, S. -J. (2021). Differential Gene Expression Associated with Altered Isoflavone and Fatty Acid Contents in Soybean Mutant Diversity Pool. Plants, 10(6), 1037. https://doi.org/10.3390/plants10061037

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop