Introduction
MATERIALS AND METHODS
Plant materials
Determination of contents of soybean isoflavone components
Statistical analysis
RESULTS AND DISCUSSION
Distribution of total isoflavones and their compositions
Analysis of variance of isoflavone contents
Variation of isoflavone contents by soybean cultivars for their end-use
Effect of growth characteristics on isoflavone content
Possibility to improve isoflavone contents by plant breeding
Introduction
Isoflavones in soybean are called natural phytoestrogens because they have a structure and efficacy similar to estrogen, a female hormone. Several studies have shown that soybean isoflavones have effects on bone health, anticancer and antioxidant activities (Peterson, 1995; Uesugi et al., 2002; Messina, 2014). In particular, consumption of various soybean food products has been recommended, considering the bioactivity of isoflavone in relieving menopausal syndrome in women, and in physiological activations that prevent osteoporosis and reduce blood cholesterol. Recently, as the interest in soybean isoflavone has increased, its products were developed into various types such as dietary supplements, cosmetics, and beverages for health purposes, and are now commercially available in the market.
Soy-isoflavone is a flavonoid compound that exists in the form of 12 chemical structures, including the aglycone forms (daidzein, glycitein and genistein), the glucoside derivatives (daidzin, glycitin, genistin), the glucosides with malonyl-esters (6’’-O-malonyldaidzin, 6’’-O-malonylgenistin, 6’’-O-malonylglycitin) and the glucosides with acetyl-esters (6’’-O-acetyldaidzin, 6’’-O-acetylgenistin, 6’’-O-acetylglycitin) (Yu & McGonigle, 2005). While the aglycone forms of isoflavones exist in small quantities, isoflavones in the glucoside forms are dominant quantitatively. Specifically, the malonylglucosides are the highest proportion of the soy isoflavones (Kim et al., 2018; Ahmad et al., 2017). Isoflavones in the form of glucosides are inactive, and are known to be absorbed more rapidly in the body when converted to the aglycone form through hydrolysis metabolism by intestinal microbiota in the human digestive tract (Izumi et al., 2000). Soy-isoflavone content and bioactivity were mainly investigated in earlier studies focusing on the aglycone forms of isoflavones such as daidzein and genistein (Shin et al., 2006; Omoni & Aluko, 2005). In recent years, researches have been underway to complement the absorption efficiency by converting isoflavones in the form of glucoside into aglycone through microbial enzyme process (Kim et al., 2014; Kim et al., 2010).
A wide range of isoflavone contents in soybean cultivars has been reported in food sources such as soy-milk and soy-paste (Kim et al., 2018; Lee et al., 2016; Park et al., 2012; Choung et al., 2006). It was known that 90% of the total isoflavones in soybean seeds were distributed in cotyledon, while the rest in hypocotyl areas (Tsukamoto et al., 1995). The content of total isoflavones present in soybean seeds varies by cultivars and genetic resources, from less than 300 ㎍/g to more than 3000 ㎍/g. According to Kim et al. (2018), the total isoflavone contents for 172 soybean varieties developed in Korea from 1913 to 2013 were 206.3 to 2666.2 ㎍/g. Park et al. (2012) also found that the total isoflavones ranged from 527.9 to 3436.5 ㎍/g for 106 soybean cultivars released in Korea from 1913 to 2006. Choung et al. (2006) analyzed 43 major Korean soybean cultivars for aglycone isoflavone content indicating 277 to 2002 ㎍/g. The isoflavone content of soybeans grown overseas in the U.S. and other countries also varies, indicating the range of 880 to 4220 ㎍/g, 683 to 2441 ㎍/g, and 570 to 1880 ㎍/g for the U.S, Chinese, and Brazilian cultivars, respectively (Xu & Chang, 2008; Wang & Murphy, 1994; Chung et al., 2008; Genoverse et al., 2005).
The development of soybean cultivars in Korea has been specialized for traditional food processing and cropping system applications such as soy-sauce and paste, tofu, soy- sprouts, cooking with rice and short-term cropping properties. With recent market trends towards health functional property, new soybean cultivars are required to create new demand. The objective of this study was to determine the variation of contents of 12 soybean seed isoflavone components in the aglycone, glucoside, malonylglucoside and acetyglucoside groups and to select for the optimum source of soy-isoflavone in the 44 Korean soybean cultivars. These cultivars were developed by the National Institute of Crop Science (NICS) of the Rural Development Administration (RDA), Republic of Korea. Considering the genetic and environmental factors affecting the contents of the total isoflavones and 12 components among the soybean cultivars, the possibility to improve soy-isoflavones in the soybean breeding programs was reviewed in this study.
MATERIALS AND METHODS
Plant materials
In this study, 44 soybean cultivars were used, which included recently developed new cultivars and the cultivars known to have isoflavone content at the time of their release. These cultivars were classified by their end-uses : 17 cultivars for soy-paste and tofu, 13 cultivars for soy- sprouts, 8 cultivars for cooking with rice, 4 cultivars for cropping system, and 2 for vegetable soybeans (Table 1). They were grown in each of 2016 and 2017 with planting on June 14, 2016 and on June 19, 2017 at the Daegu experimental field site in the Department of Southern Area Crop Science, NICS, RDA in Habin-myeon, Dalseong-gun, Gyeongsangbuk-do, Republic of Korea.
Table 1. List of 44 Korean soybean cultivars classified by their end-uses.
§Seedcoat color : Y yellow, B black, G green, PG pale green
Seeds of individual cultivars were planted in rows 4 m long, 0.7 m apart and with a space of 15 cm between plants. The plants were thinned to two plants/hill to maintain a uniform density at 14 days after planting in the plots. Fertilizers were applied at a ratio of N : P2O5 : K2O = 3 : 3 : 3.4 kg/10a, and other field management was conducted in accordance with the RDA’s standard cultivation protocols. Soybean cultivars were completely randomized in the plots, and two and three replicates for individual cultivars in 2016 and in 2017, respectively were applied to analysis. Number of days to flowering and/or maturity from planting were examined for individual cultivars in each year. The harvested seeds were dried for 3 days, threshed, and then selected for analysis. The trait of 100- seed weight of individual cultivars was measured by three repetitions.
Determination of contents of soybean isoflavone components
Soybean seed samples were crushed to 100-mesh size without breaking the seed coat with vibrating sample mill (T1-100 model, CMT Co Ltd., Tokyo Japan) and stored in a 4 °C refrigerator for isoflavone analysis. Standards of daidzein, glycitein, genistein, daidzin, glycitin and genistin were purchased from Sigma (St. Louis, MO, USA). Other standards of 6"-O-malonyldaidzin, 6"-O-malonylglycitin, 6"- O-malonylgenistin, 6"-O-acetyldaidzin, 6"-O-acetylglycitin and 6"-O-acetylgenistin were used and purchased from Nacalai Tesque INC. (Kyoto, Japan).
For the analysis of soybean isoflavone, 20 mL of 50% methanol was added to 1.0 g of seed pulverized sample and stirred for 24 hours at 250 rpm to elute the isoflavone. After the extraction was completed, the samples were allowed to be left unattended for one hour at room temperature and the supernatant was taken and filtered again with a filter paper (No. 1 Whatman) and the 1 mL of the filtered solution was filtered through a 0.2 μm syringe filter and used as an analysis sample for HPLC.
HPLC analysis was conducted using Dionex Ultimate 3000 (Thermo, USA). HPLC column was LiChrospher RP-C18 (5 μm, 125×4.0 mm). The 10 ㎕ filtrate of sample was injected to HPLC and analyzed with mobile phase (H2O (0.1% acetic acid) and acetonitrile (0.1% acetic acid)) at a flow rate of 1.0 mL/min. The detection wavelength was UV 260 nm and the column temperature was limited to 30°C. Isoflavone compositions in soybean seeds were separated according to the retention time (Fig. 1). To calculate the standard calibration curve for quantitative analysis of isoflavones, each of 12 standard substances was dissolved in a methanol containing 1% dimethylsulfoxide (DMSO) and diluted to prepare for 50, 40, 30, 25, 20, 10 μg/ml, respectively. The standard of isoflavone components for each prepared concentration was calculated by the calibration equation derived from the corresponding peak area through the HPLC chromatogram under the same condition. For all standard reference materials, the standard calibration curves showed high linearity (R2 > 0.997).
Statistical analysis
Analysis of variance (ANOVA) of the isoflavone data across multiple year environments was performed using the general linear model (PROC GLM) in SAS 9.4 (SAS Institute Inc., USA). This linear model had the effect of environment of year, replication, genotype, and environment x genotype interaction. Variance components were determined in ANOVA to estimate heritability in broad-sense. Comparisons between means of cultivars for each isoflavone component were made using LSDs at a 0.05 probability level. Pearson correlation was calculated with PROC CORR in SAS 9.4 (SAS Institute Inc., USA) to determine the relationship among the same pairs of variables of agronomic traits and isoflavone components. Broad-sense heritability (h2) was estimated for the contents of the total isoflavone and its 12 components by using the following formula (Toker, 2004):
phenotypic variance ,
where y, g, and r are number of year, genotype and replication, respectively; , and are the components of variance for genotypes, genotype x year and error, respectively. Broad-sense heritability (h2) was calculated as genotypic variance/phenotypic variance :
RESULTS AND DISCUSSION
Distribution of total isoflavones and their compositions
Table 2 shows the results of the total isoflavone content as well as the mean, minimum and maximum values of the isoflavone contents of the 12 composition groups for all 44 soybean cultivars in each growing year or over both 2016 and 2017. The total isoflavone content of the 44 cultivars averaged at 2935.4 ㎍/g and was in the range of 950.6 to 5226.3 ㎍/g over the two years. The isoflavone content of seeds harvested in 2016 averaged at 2745.0 ㎍/g and ranged from 723.6 to 5202.4 ㎍/g. For 2017, the average was 3061.4 ㎍/g and 963.8 to 5770.0 ㎍/g. The content of isoflavones increased by about 11.5% on average in 2017 as compared to 2016.
Table 2. Total isoflavone content as well as the mean, minimum and maximum values of the isoflavone contents of the 12 composition components for all 44 soybean cultivars in each growing year or over 2016−2017.
The average contents in 2016–2017 were analyzed for the four composition groups of isoflavones such as aglycone, glucoside, malonylglucoside, and acetylglucoside : the average amount of malonylglucoside group was 2437.2 ㎍/g with the highest proportion of isoflavone composition (83.0%) and ranges from 803.9–4333.5 ㎍/g. The other groups of isoflavone were small in contents in order of glucoside groups (average of 277.3 ㎍/g, range of 91.2 to 525.1 ㎍/g) > aglycone groups (average of 159.8 ㎍/g, range of 32.9 to 330.5 ㎍/g) > acetylglucoside group (average of 61.1 ㎍/g, range of 17.1 to 117.8 ㎍/g). The tendency of variation in the contents of isoflavone composition groups was similar in each year. The malonylglucoside group showed the highest level with 78.9–85.5% and the acetylglucoside group had the lowest level with 0.8–4.1%.
When comparing the content by isoflavone type, the malonylglucoside group was the major constituent of soybean isoflavone, and specifically, the highest ratio of malonylgenistin content was observed in the malonylglucoside group, which was almost similar to the previous results (Lee et al., 2018; Kim et al., 2018). The average contents of genistein and daidzein, known to have high bio-activity among the aglycone components, were almost the same level with 59.1 ㎍/g and 62.4 ㎍/g, respectively, accounting for about 2% of the total isoflavone content. In contrast, glycitein had an average of 38.3 ㎍/g for over two years, which was relatively small in the aglycone composition. Choung et al. (2006) also found that the proportions of daidzein and genistein in the aglycone form were almost similar with each other in Korean soybean genetic resources except for the early-matured cultivars. A similar trend was notable with the small difference between daidzin and genistin components for the glucoside and malonylglucoside groups as well as the aglycone except for the acetylglucoside (Table 2). Even in the group of malonylglucoside with the highest rate of isoflavone composition, malonylglycitin was shown to be less than malonylgenistin and malonyldaidzin in the isoflavone composition ratio. Acetylglucoside group was present in very small amounts in each composition and could not be analyzed in the study. In particular, acetyldaidzin was rarely detected in the seed samples produced in 2016.
Analysis of variance of isoflavone contents
Analysis of variance (ANOVA) for the content of total isoflavones and 12 composition components for the 44 soybean cultivars developed by NICS, RDA showed highly significant differences among the cultivars (P < 0.0001) (Table 2). Environmental effects between years were significant for the contents of isoflavone group components except daidzein and genistin. The interactions of genotypes with cropping years were also significant in total isoflavone content and its components. This indicates that variation of soybean isoflavone levels is affected by all genotypic, environmental and g x e interaction effects. According to Hoeck et al. (2000), even a single cultivar might have a 1.7-fold difference in single cultivation conditions.
Heritability is an estimate of the genetic difference contributed to the phenotypic variation. As seen in Table 2, the estimated heritability (h2) in a broad-sense ranged 0.8–0.93 for the total isoflavone and its components except for the lowest acetylglucoside group. Considering variation due to environmental factors over years, the estimates of heritability were high within the set of 44 Korean soybean cultivars. According to ANOVA in this study, genotypic variance contributing to the observed variability of the levels of isoflavone components was greater than the variances of nonheritable components of environment and g x e interactions (data not shown). Overall, this result suggests that isoflavone content in soybean seed is influenced by genetic and environmental factors, and that the magnitude of each contribution of variance components to phenotypical expression is different. In the present study, difference between cultivars due to genetic factors in isoflavone content are large as compared with environmental variation and interactions with cultivars, and hence the heritability can be expected to increase. According to Seker et al. (2014), heritability estimate increased as the difference between the values of genotypic and phenotypic variance in percentage of mean decreased for a target trait.
The heritability estimate of malonylglucoside components was relatively high with the range of 0.87–0.96. Unlike other component groups, heritability of acetylglucoside group was low due to a greatly significant g x e interaction and yearly variation as well as its low detection level of composition. The results were almost similar to the high range of heritability values for the soybean isoflavone content in the former studies (Gutierrez-Gonzalez et al., 2009; Chiari et al., 2004). Gutierrez-Gonzalez et al. (2009) found that heritability values of genistein, daidzein and total isoflavone were 0.83, 0.86 and 0.86, respectively, except glycitein (0.43) in the genetic analysis of isoflavone content using F7-derived RIL population of Essex × PI 437654. In addition to the existing results, high heritability values in this study indicates that selection for high isoflavone on the individual plant base might be allowed in the early breeding generations (Chiari et al., 2004).
The correlations of 44 Korean soybean cultivars for the contents of total isoflavone and three composition groups (aglycone, glucoside and malonylglucoside) between two cropping years of 2016–2017 were highly significant except for acetylglucoside group (r = 0.15ns) (Fig. 2). The correlations for total isoflavone and malonylglucoside group contents were similar each other (r = 0.88****) between 2016 and 2017. On the other hand, correlations for the contents of aglycone and glucoside group were lower than malonylglucoside, but all statistically significant with r = 0.78**** and 0.67****, respectively. This result of high correlation among the soybean cultivars over the two years suggests that some cultivars have relatively consistent level of high or low isoflavone contents controlled predominantly by the genetic factors rather than environmental factors as shown in the heritability results for isoflavone content.
Variation of isoflavone contents by soybean cultivars for their end-use
Table 3 shows ANOVA results for the content of total isoflavone and 12 individual components indicating significant variance in the genotype and g x e interactions for each of four groups of 44 soybean cultivars classified by their end-uses such as soy-paste and tofu, soy-sprout, cooking with rice, and early maturity. Significant differences were also observed between years for the isoflavone composition groups except for glucoside content in the cultivar group of soy-paste and tofu, and aglycone content in the cultivar group of soy-sprout. The total isoflavone content in cultivar group for soy-sprout was higher (3850.4 ㎍/g) than that for the other cultivar groups of soy-paste and tofu (3082.8 ㎍/g), black or green soybean cooked with rice (2345.8 ㎍/g) and early maturity group (1298.6 ㎍/g). There were significant differences among the cultivar groups in the contents for total isoflavone and four composition groups (Fig. 3). The variation of total isoflavone content in soybean cultivars grouped by end-use was also consistent with the previous report (Park et al., 2006). The trend of variations of average content of each isoflavone composition group was the same as that of total isoflavone content in rankings of soy-sprout > soy-paste and tofu > cooked with rice > early maturity group (Fig. 3). The content of acetyldaidzin belonging to acetylglucoside, which was relatively low in content of isoflavone composition, was low or almost not detected in all cultivars regardless of their end-uses (Table 3).
Table 3. Mean, minimum and maximum values and ANOVA results for the content of total isoflavone and 12 individual components for four groups of 44 soybean cultivars classified by their end-uses.
Distribution of soybean cultivars by their average content of total isoflavone over two years is shown in Table 4 and Fig. 4, respectively. The total isoflavone content of ‘Sowonkong’, a soybean cultivar for soy-sprout, was the highest (5226.3 ㎍/g) and was about 5.5 times higher than that of ‘Taecheong’, a black soybean cultivar with green cotyledon and the lowest isoflavone content (950.6 ㎍/g). In the cultivar group for soy-paste and tofu, the average contents of total isoflavone for ‘Daepung’, ‘Daepung2ho’, ‘Saegeum’, ‘Uram’, and ‘Jinpung’ were higher than 4000 ㎍/g. However, ‘Taekwangkong’, ‘Saedanbaek’, ‘Danbaegkong’ in the same cultivar group were relatively low with less than 1500 ㎍/g. In particular, these cultivars at the upper and lower levels of isoflavone content were stable without significant annual variations. In the small seeded cultivar group for soy-sprout, in addition to ‘Sowonkong’, such cultivars as ‘Sinhwa’, ‘Sorog’ and ‘Haewon’ were high in isoflavone content with ≥ 4480 ㎍/g despite their annual variation. The colored soybean cultivars with black or green seed coats were in the range of 950.6 to 3411.2 ㎍/g in total isoflavone. Early-matured soybean cultivars for cropping system or green bean vegetables only had 1000 to 1500 ㎍/g except for ‘Chamol (2050.8 ㎍/g)’, and generally had lower isoflavones than the other cultivars of soy-paste, tofu, and soy-sprouts. The early-matured cultivars were expected to have a short growth period, which would have prevented isoflavone synthesis due to high temperature from flowering to maturity (Park et al., 2012; Tsukamoto et al., 1995).
Table 4. Distribution of 44 soybean cultivars by average content of total isoflavone over 2016−2017.
Effect of growth characteristics on isoflavone content
Since the soybean cultivar group for soy-sprout, which had a 100-seed weight of about 12 g or less, had a higher content of isoflavones on average than other cultivars, so the relationship between seed size and isoflavone content was analyzed (Table 5). Seed weight and total isoflavone content of 44 soybean cultivars were significantly negatively correlated in 2016 and 2017, respectively (r = -0.47** and -0.49**). This means that the smaller the seed size, the higher the isoflavone content. However, isoflavone contents of some cultivars with small seed size were exceptionally lower than the average cultivars within the group : the isoflavone contents of ‘Sohwang’ for soy-sprout and ‘Socheongja’, a small black soybean with green cotyledon were 1895 ㎍/g and 2320 ㎍/g, respectively.
Table 5. Correlations of growth characteristics with the average contents of total isoflavone and 12 individual components for 44 soybean cultivars during 2016 and 2017.
†fd-flowering date, pd-planting date, md-maturing date
The correlations of seed weight with the contents of individual isoflavone component groups in 2016 and 2017 were varied. They were significant with r = -0.50** and -0.52** for glucoside group and r = -0.45** and -0.49*** for malonylglucoside group. In the case of the aglycone group, its correlation with seed weight varied over years, and was significant for all three components of genistein, daidzein and glycitein in 2016 but only daidzein in 2017. The number of days to maturity after sowing 44 cultivars was positively correlated with the total isoflavone content, although the correlation was not high (r = 0.4**–0.44**) (Table 5). However, the number of days from sowing to maturity was not related to isoflavone content when considering 38 cultivars except for the early-matured six cultivars (r = 0.09ns). It was difficult to find the correlation between the two traits because 80% of the soybean cultivars except the early-matured had a small variation in the range of 120–130 days of growth, while the variation of the isoflavone content was relatively large.
The number of days of growth from flowering to maturing did not affect the isoflavone content variation (Table 5). This was contrary to the results of Yun et al. (2006) who reported that soybean isoflavone contents increased as the periods from flowering to maturity were longer in the analysis of Korean soybean cultivars of the early-, mid- and late-matured and Chinese cultivars. These results indicated no specific trend for the relationship between the length of growth period and the isoflavone content of soybean, and it may vary depending on the characteristics of the cultivars and resources to be analyzed.
Possibility to improve isoflavone contents by plant breeding
It is important to improve the utilization of soybean genetic resources with high isoflavone contents as the consumer’s interest in bioactive components increases. From the study and various reports, the soybean cultivars developed in Korea have a large genetic variation and high heritability values in regards to isoflavone content. Therefore, it is surmised that selection of putative germplasm and improvement of elite cultivars for soy-isoflavone contents will be effective for Korean soybean. In particular, recent molecular genetic approaches based on QTL mapping have been actively conducted to promote understanding of genome- based isoflavone content and composition components (Han et al., 2016; Meng et al., 2016; Smallwood et al., 2014; Zeng et al., 2009). It will be allowed in the near future to suggest the possibility to develop soybean cultivars for specialized purposes that enhance isoflavone contents by combining genome-based technology with conventional breeding schemes. Some positive effects will be expected such as the increase of consumption of domestic soybeans and substitution of the foreign imports by diversifying usability of functional materials in Korean soybean as well as traditional end-use of food processing.