Increase of isoflavones in the aglycone form in soybeans by targeted crossings of cultivated breeding material (2024)

Isoflavones are naturally occurring group of phytoestrogens. They are mostly found among the members of the Fabaceae family, primarily soybeans and red clover¹. Soybean (Glycine max (L.) Merr.) is an industrial crop with a stable production growth. Soybean is the fourth major crop in the world, following cereal crops, such as maize, wheat and rice, which are the main source of human food. Since soybeans and soy products are the richest source of isoflavones in human diet², the content of isoflavones in new soybean cultivars has become a trait with increased significance. Regarding differences in bioavailability of individual isoflavones, it is important to establish the form in which they are present in the seed³. Soybeans contain twelve different phytoestrogens. Daidzein, genistein and glycitein are the aglycones which can form three glucoside forms, a β-glucoside, a 6″-O-malonyl-glucoside and a 6″-O-acetyl-glucoside⁴. The isoflavone aglycones are absorbed faster and in greater amounts than their glucosides in humans⁵.

Considering the increasing presence of soybean in human diet, the nutritive value of soybean and soy products have become the focus of numerous studies in agronomy, production technology, and medicine. Moreover, the key role of diet in the prevention of some serious pathologies, such as cardiovascular diseases, atherosclerosis and cancer, has been globally recognized. Soybeans have been extensively studied around the world due to their presence in human diet and many health benefits which have been attributed to the consumption of soy and soy products. Special attention has been placed on the research on soybean effects on certain types of cancer. Studies regarding the relationship between soy intake and breast cancer^{6,7,8,9,10,11} showed that increased amounts of soybean in daily diet reduce the risk of breast cancer in Asian women, which is partly attributed to the presence of isoflavones. In addition, it has been confirmed that isoflavones, primarily their aglycones (genistein and daidzein), is involved in a number of biological activities including breast and prostate cancer chemo-preventive activity, and the ability to modify carcinogenesis, e.g. initiation, promotion, and cancer progression^12,13,14,15.

Increasing the content of isoflavones in new varieties of soybean, particularly in their aglycone form, is one of the most important objectives in modern soybean breeding programs. The influence of various factors on isoflavones content and composition in soybeans has been studied by different authors. Content and composition of isoflavones in soybeans largely depend on the soybean genotype. While the origin of soybeans does not have a major impact on these characteristics^16,17, the maturity group i.e. the length of vegetation period does^18,19. According to previous studies, the amount and the composition of isoflavones in soybeans are affected by the environment, primarily by the temperature during the vegetation period and storage^20,21,22,23.

Previous investigations of the inheritance of isoflavone content in soybeans indicate that genetic improvement for these traits should explore the additive genetic variance in superior lines, or the cytoplasmic effect and the epistatic interactions between cytoplasmic and nuclear genes to obtain the maximum gain in selection²⁴. Our previous results showed that F₁ soybean progenies could increase isoflavone content²⁵, and that the content and composition of isoflavones could be passed from the parental genotypes to the hybrids, and therefore utilized for breeding soybean cultivars with desirable traits²⁶.

Since there is not much data on genetics of isoflavones in soybean, elucidation of the mode of inheritance is necessary in order to design an efficient breeding strategy for development of high-isoflavone soybean genotypes. The aim of this study was to investigate the inheritance of total isoflavones and their aglycone form in F₂ generations, as well as prospective increase of the total isoflavone content in genotypes with the low isoflavone content.

Crosses

Based on the isoflavone content in 23 different soybean genotypes from four different maturity groups (00, 0, I, and II), three crosses were made in order to determine the inheritance pattern and increase the content of total isoflavones and their aglycone form. The entire list of isoflavone content in these 23 soybean genotypes can be found as Supplementary TableS1. Genotype with the lowest total isoflavone content (NS-L-146) was crossed with the low- (NS Zenit), mid- (NS Maximus), and high- (NS Virtus) isoflavone content genotypes. Offspring of these crosses was grown the following year as the F₁ generation at the same locality, according to the standard methodology²⁷. The same stands for the F₂ generation.

Isoflavone analysis

Seeds used for in vitro experiments were collected at full maturity. The collected seeds were then extracted by the method of Andlauer, Martena, and Furst²⁸. Powdered soybean seeds (500 mg) were defatted by hexane extraction (2 × 10 ml, 30 min, and subsequent centrifugation, 30 min, 1780 rcf) and then extracted for 2 h with 8 ml methanol/water (4:1, v/v) and centrifuged (30 min, 1780 rcf). Prior to HPLC injection, each extract was filtered using Agilent technologies Teflon filters (0.45 μm, Delaware, Wilmington).

Individual isoflavones were identified and quantified according to the method of Lee et al.⁴, with minor modifications. Separation was achieved on 5 μm Zorbax SB C18 HPLC column (150 × 4.6 mm), with Zorbax SB C18 guard column. Mobile phase consisted of two solvents. Solvent A was 1% (v/v) acetic acid in water and solvent B was 100% acetonitrile. Analysis was conducted under the following conditions: 0–5 min 85% A; 5–44 min from 85 to 65% A; 44–45 min from 65 to 85% A, and 45–50 min 85% A. The column temperature was 25 °C, a solvent flow rate 0.6 ml/min, and injection volume 10 μl. The spectra were collected between 240 and 400 nm by DAD, and components in the eluate were detected at 270 nm.

Isoflavones were identified by comparing the retention times in HPLC chromatograms and UV spectral patterns with those of standard compounds and literature data^4,28 (Figs1 and 2).

Isoflavone concentrations were quantified by external standard (five-point regression curves, r ≥ 0.9997) of daidzein, glycitein and genistein standard compounds. Standard solutions were made by dissolving standard compounds in mixture of methanol/water (4:1, v/v) and linearity was studied for each compound in the range of 0.5 to 50 mg/l (0.5 0.1 2.0 3.0, 5.0, 10.0, 25.0 and 50.0). As only standard phytoestrogen aglycones were used, the content of the corresponding glycoside forms was obtained by calculation. For this purpose, calibration curves of the corresponding aglycone compounds were used and corrections for differences in molecular weight between aglycones and glucosides were applied following the pattern given by Romani et al.²⁹:

Statistical analysis

All sample measurements were done in triplicate. Data were processed by analysis of variance (ANOVA) using statistical program package STATISTICA v. 12.0³⁰. Statistical comparisons between samples were performed by Student’s t-test for independent analysis. Differences were considered significant at p < 0.05. The method proposed by Mahmud and Kramer³¹ was used to estimate broad-sense heritability on a single-plant basis. The method involves the measurement of a character on F2 plants of a single-cross population and on the inbred parents used to form the population. The formula for estimating heritability presented by these investigators is

Where:

h² = heritability, \({\sigma }_{F2}^{2}\) = variance of F₂, \({\sigma }_{P1}^{2}\) = variance of P₁, \({\sigma }_{P2}^{2}\) = variance of P₂

Principle Component Analysis was done in statistical program SPSS 10.0.

Broad-sense heritability for analyzed traits was medium to low. It ranged from 0.20 for acetyl genistin of F₂ NS-L-146 × Virtus, to 0.61 for acetyl glycitin of the same cross (Table1).

Three F₂ populations, obtained from the crosses between the low isoflavone genotype and three other parents (low, medium, high), had the different mode of inheritance of isoflavone content (Fig.3a).

Dominant isoflavone in soybean seed is diadzein, followed by genistein and glycitein. Significant differences in the content of individual isoflavones were not found between populations, except for the content of total glycitein. Malonyl daidzin, followed by daidzein, were the dominant forms of diadzein (Fig.3b). Content of acetyl daidzin in the F₂ population made from NS-L146 × Zenit was higher compared to the both parents. The average content of malonyl daidzin in populations F₂ (NS-L146 × Maximus) and F₂ (NS-L146 × Virtus) low and similar to NS-L-146. Glycitin and malonyl glycitin were the dominant forms of glycitein (Fig.3c). All three populations had a similar, unusually high content of acetyl glycitin, regardless of the parents used for crosses. Furthermore, there were several lines in populations F₂ (NS-L146 × Maximus) and F₂ (NS-L146 × Zenit) with glycitein different from zero. Malonyl glycitin content in populations F₂ (NS-L146 × Maximus) and (NS-L146 × Virtus) exceeded the mid-parent value. In parental genotypes, malonyl genistin was the dominant form of genistein while in F₂ populations genistein was found in the form of malonyl genistin and genistin (Fig.3d). There were no significant differences between F₂ populations regarding their genistein content. Overall, F₂ populations had a similar content of diadzein and genistein (glycon and aglycon form) and population average value was closer to the lower parent. On the other hand, dominance or overdominance of the better parent was observed regarding glycitein content in F₂ populations.

Ratio of aglycone and glycone forms of isoflavones is an important indicator of isoflavone content, due to the faster absorption of the aglycone form in the gastrointestinal tract. Population average indicates that the aglycone ratio showed the mid-parent value mode of inheritance. Even from crosses between parents with zero value aglycone ratio, some progenies had a certain amount of aglycone form of particular isoflavones (Fig.4).

The first two principal components of F₂ lines isoflavone content, explained 68% of the total trait variation from the initial dataset (Fig.5), 43% for principal component I and 25% for principal component II. The first principal component was mostly related to acetyl, 7-glycoside form of diadzein, malonyl glycitin, genistin, and total isoflavones, while the second principal component explained diadzein, glycitein, acetyl genistin content and the aglycone ratio. Although average isoflavone content showed little difference between populations, the principal component analysis divided F₂ lines into two groups. The first group consisted of lines with high isoflavone content (Fig.5, red circle), while the second group consisted of lines with higher aglycone ratio. On the other hand, the lines obtained from the NS-L146 × Zenit cross were grouped on the PCA plot, while the lines from the other two crosses were more disperse, which indicates low breeding potential of the NS-L146 × Zenit cross.

The prospect of increasing total isoflavones in genotypes with low isoflavone content was tested by choosing four parents, based on the tests of isoflavone content in 23 different soybean genotypes.

Mean values of total isoflavone content were the same for all three F₂ populations, regardless of the parents used in the cross. Crossing the parents which had low total isoflavone content (NS-L-146 by Zenit) resulted in the inheritance of total isoflavone content from the better parent (Zenit), which leads to the conclusion that the inheritance is associated with the better parent. However, the same level of total isoflavones was obtained in crosses between the lines with low total isoflavone content and cultivars with medium and high total isoflavone content, such as Maximus and Virtus. Therefore, it was not possible to make the final decision on the mode of inheritance, as this trait is affected by many genes with few individual effects and it is under strong influence of the environmental factors^32,33.

The content of diadzein and diadzein conjugates, as well as genistein and genistein conjugates, followed the same pattern of inheritance as total isoflavones. At the same time, they were the dominant forms of isoflavones in all the analyzed samples. This is in agreement with the previous observations, stating that soybeans and soy foods usually contain similar amounts of genistein and daidzein and a much lower amount of glycitein³⁴. On the other hand, glycitein and glycitein conjugates had the highest variation, which is in agreement with the result of Gutierrez-Gonzales et al.³⁵. The content of glycitein conjugates in F₂ populations is especially interesting. There was no difference in glycitin content between the three populations and the parent with the highest glycitin content (Virtus). Acetyl glycitin content of the three F₂ populations was equally higher, compared to the parent with the highest acetyl glycitin content (Virtus), which is interesting given the fact that F₂ (NS-L-146 × Zenit) derived from the parents with no detected acetyl glycitin content. Low isoflavone levels of the offspring indicate that maternal effect has a significant role in isoflavone inheritance, as showed by Chiari et al.²⁴.

Aglycone ratio of the F₂ populations derived from the cross between NS-L-146, Maximus and Virtus, approximately had an intermediate type of inheritance. However, the cross between NS-L-146 and Zenit, which did not contain free isoflavone forms, resulted in some F₂ lines with a certain amount of free isoflavone forms. Inheritance of aglycone ratio in soybeans is scarcely documented, but our results point out the aglycone ratio as one of the prospective breeding goals in future soybean breeding programs.

The use of Principal Component Analysis to estimate isoflavone content in individual lines provides a clearer view of the variation of this trait in F₂ populations. Most lines from the F₂ (NS-L-146 × Zenit) were grouped on the PC graph, indicating their low variation and thus low significance in breeding for this trait. On the other hand, besides their total isoflavone content, the position of several lines from F₂ (NS-L-146 × Maximus) on the graph confirms high variability for this trait, so they will be used for further breeding programs. Lines from F₂ (NS-L-146 × Virtus) and F₂ (NS-L-146 × Maximus), have a higher aglycone ratio than the better parent. Therefore, they will serve as important starting material in the development of genotypes with increased content of isoflavones in the aglycone form, which is an important objective in modern soybean breeding programs.

According to the obtained results, future soybean breeding for changed isoflavone content and composition is possible. Breeding success is important considering the constant need for improvement i.e. increase of isoflavone content in cultivars with the desired agronomic traits, and the rising importance of isoflavones in human diet. In their choice of parental components, breeding programs should aim at obtaining the highest possible variability, which would be possible by using both the average and high isoflavone genotypes. Further research should focus on detecting the mode of inheritance in order to achieve major breakthrough in breeding for this significant trait.

Authors and Affiliations

1. Soybean Department, Institute of Field and Vegetable Crops, 21000, Novi Sad, Serbia

   Jegor Miladinović,Vuk Đorđević,Svetlana Balešević-Tubić,Kristina Petrović&Marina Ćeran

2. Department of Pharmacy, Faculty of Medicine, University of Novi Sad, 21000, Novi Sad, Serbia

   Jelena Cvejić&Mira Bursać

3. Industrial Crops Department, Institute of Field and Vegetable Crops, 21000, Novi Sad, Serbia

   Dragana Miladinović

Contributions

J.M., V.Đ. and M.Ć. conceived and designed the experiments. V.Đ., K.P., S.B.T. and J.M. carried out experimental determinations. J.C. and M.B. carried out chemical analyses. J.M., D.M. and M.Ć. participated in the data analysis and results interpretation. J.M., V.Đ. and M.Ć. wrote the manuscript and D.M. revised the definitive version. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jegor Miladinović.

Competing Interests

The authors declare no competing interests.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions