Genetic control of soybean seed isoflavone content: importance of statistical model and epistasis in complex traits

Cytochrome GmGLY1 is Involved in the Biosynthesis of Glycitein in Soybean

Isoflavones, the major secondary metabolites of interest due to their benefits to both human and plant health, are exclusively produced by legumes. In this study, we profiled the isoflavone content in dry seeds from 211 soybean [Glycine max (L.) Merr.] accessions grown across five environments. Broad and discernible phenotypic variations were observed among accessions, regions, and years of growth. Twenty-six single-nucleotide polymorphisms (SNPs) associated with the sum of glycitein (GLE), glycitin (GL), 6″-O-acetylglycitin (AGL), and 6″-O-malonylglycitin (MGL) contents were detected in multiple environments via a genome-wide association study (GWAS). These SNPs were located on chromosome 11 (8,148,438 bp to 8,296,956 bp, renamed qGly11-01). Glyma.11g108300 (GmGLY1), a gene that encodes a P450 family protein, was identified via sequence variation analysis, functional annotation, weighted gene coexpression network analysis (WGCNA), and expression profile analysis of candidate gene, and hairy roots transformation in soybean. Overexpression of GmGLY1 increased the glycitein content (GLC) in soybean hairy roots and transgenic seeds, while CRISPR/Cas9-generated mutants exhibited decreased GLC and increased daidzein content (DAC). Haplotype analysis revealed that GmGLY1 allelic variations significantly affect the GLC accumulation. These findings enhance our understanding of genes influencing GLC in soybean and may guide breeding for lines with high and stable GLC.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

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

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

Identification of Potential Gene Regulatory Pathways Affecting the Ratio of Four-Seed Pod in Soybean

The number of four-seed pods is one of the most important agronomic traits affected by gene and environment that can potentially improve soybean (Glycine max) yield. However, the gene regulatory network that affects the ratio of four-seed pod (the ratio of the number of four-seed pods to the total number of pods in each individual plant) is yet unclear. Here, we performed bulked segregant RNA sequencing (BSR-seq) on a series of recombinant inbred lines (RILs) derived from hybrid progenies between Heinong 48 (HN48), a cultivar with a high ratio of four-seed pod, and Henong 64 (HN64), a cultivar with a low ratio of four-seed pod. Two tissues, flower bud and young pod, at two different growth stages, R1 and R3, were analyzed under the ratios of four-seed pod at less than 10% and greater than 30%, respectively. To identify the potential gene regulation pathways associated with the ratio of soybean four-seed pod, we performed differentially expressed analysis on the four bulked groups. A differentially expressed gene (DEG) encoding a photosystem II 5-kDa protein had the function of participating in the energy conversion of photosynthesis. In addition, 79 common DEGs were identified at different developmental stages and under different ratios of four-seed pod. Among them, four genes encoding calcium-binding proteins and a WRKY transcription factor were enriched in the plant-pathogen interaction pathway, and they showed a high level of expression in roots. Moreover, 10 DEGs were identified in the reported quantitative trait locus (QTL) interval of four-seed pod, and two of them were significantly enriched in the pentose and glucuronate interconversion pathway. These findings provide basic insights into the understanding of the underlying gene regulatory network affected by specific environment and lay the foundation for identifying the targets that affect the ratio of four-seed pod in soybean.

Identification of a candidate gene associated with isoflavone content in soybean seeds using genome-wide association and linkage mapping

Isoflavone, a secondary metabolite produced by Glycine max (L.) Merr. (soybean), is valuable for human and plant health. The genetic architecture of soybean isoflavone content remains unclear, however, despite several mapping studies. We generated genomic data for 200 soybean cultivars and 150 recombinant inbred lines (RILs) to localize putative loci associated with soybean seed isoflavone content. Using a genome-wide association study (GWAS), we identified 87 single-nucleotide polymorphisms (SNPs) that were significantly associated with isoflavone concentration. Using linkage mapping, we identified 37 quantitative trait loci (QTLs) underlying the content of four isoflavones found in the RILs. A major locus on chromosome 8 (qISO8-1) was co-located by both the GWAS and linkage mapping. qISO8-1 was fine mapped to a 99.5-kb region, flanked by SSR_08_1651 and SSR_08_1656, in a BC2 F5 population. GmMPK1, encoding a mitogen-activated protein kinase, was identified as the causal gene in qISO8-1, and two natural GmMPK1 polymorphisms were significantly associated with isoflavone content. Overexpression of GmMPK1 in soybean hairy roots resulted in increased isoflavone concentrations. Overexpressing GmMPK1 in transgenic soybeans had greater resistance to Phytophthora root rot, suggesting that GmMPK1 might increase soybean resistance to biotic stress by influencing isoflavone content. Our results not only increase our understanding of the genetic architecture of soybean seed isoflavone content, but also provide a framework for the future marker-assisted breeding of high isoflavone content in soybean cultivars.

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

Isoflavones are a group of phytoestrogens, naturally-occurring substances important for their role in human health. Legumes, particularly soybeans (Glycine max (L.) Merr.), are the richest source of isoflavones in human diet. Since there is not much current data on genetics of isoflavones in soybean, particularly in the aglycone form, elucidation of the mode of inheritance is necessary in order to design an efficient breeding strategy for the development of high-isoflavone soybean genotypes. Based on the isoflavone content in 23 samples of soybeans 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. Genotype with the lowest total isoflavone content (NS-L-146) was crossed with the low- (NS Zenit), medium (NS Maximus), and high- (NS Virtus) isoflavone genotypes. There were no significant differences in the total isoflavone content (TIF) between F2 populations, and there was no transgression among genotypes within the populations. Each genotype within all three populations had a higher TIF value than the lower parent (NS-L-146), while genotypes with a higher TIF value than the better parent were found only in the NS-L-146 × NS Zenit cross. However, significant differences in the aglycone ratio (ratio of aglycone to glycone form of isoflavones) were found between the populations. The highest aglycone ratio was found in the NS-L-146 × NS Maximus cross. The results indicate that the genetic improvement for the trait is possible.

Is Chickpea a Potential Substitute for Soybean? Phenolic Bioactives and Potential Health Benefits

Legume seeds are rich sources of protein, fiber, and minerals. In addition, their phenolic compounds as secondary metabolites render health benefits beyond basic nutrition. Lowering apolipoprotein B secretion from HepG2 cells and decreasing the level of low-density lipoprotein (LDL)-cholesterol oxidation are mechanisms related to the prevention of cardiovascular diseases (CVD). Likewise, low-level chronic inflammation and related disorders of the immune system are clinical predictors of cardiovascular pathology. Furthermore, DNA-damage signaling and repair are crucial pathways to the etiology of human cancers. Along CVD and cancer, the prevalence of obesity and diabetes is constantly increasing. Screening the ability of polyphenols in inactivating digestive enzymes is a good option in pre-clinical studies. In addition, in vivo studies support the role of polyphenols in the prevention and/or management of diabetes and obesity. Soybean, a well-recognized source of phenolic isoflavones, exerts health benefits by decreasing oxidative stress and inflammation related to the above-mentioned chronic ailments. Similar to soybeans, chickpeas are good sources of nutrients and phenolic compounds, especially isoflavones. This review summarizes the potential of chickpea as a substitute for soybean in terms of health beneficial outcomes. Therefore, this contribution may guide the industry in manufacturing functional foods and/or ingredients by using an undervalued feedstock.

Comprehensive RNA sequencing and co-expression network analysis to complete the biosynthetic pathway of coumestrol, a phytoestrogen

Coumestrol (CMS), a coumestan isoflavone, plays key roles in nodulation through communication with rhizobia, and has been used as phytoestrogens for hormone replacement therapy in humans. Because CMS content is controlled by multiple genetic factors, the genetic basis of CMS biosynthesis has remained unclear. We identified soybean genotypes with consistently high (Daewonkong) or low (SS0903-2B-21-1-2) CMS content over 2 years. We performed RNA sequencing of leaf samples from both genotypes at developmental stage R7, when CMS levels are highest. Within the phenylpropanoid biosynthetic pathway, 41 genes were tightly connected in a functional co-expression gene network; seven of these genes were differentially expressed between two genotypes. We identified 14 candidate genes involved in CMS biosynthesis. Among them, seven were annotated as encoding oxidoreductases that may catalyze the transfer of electrons from daidzein, a precursor of CMS. Two of the other genes, annotated as encoding a MYB domain protein and a MLP-like protein, may increase CMS accumulation in response to stress conditions. Our results will help to complete our understanding of the CMS biosynthetic pathway, and should facilitate development of soybean cultivars with high CMS content that could be used to promote the fitness of plants and human beings.

Dense genotyping-by-sequencing linkage maps of two Synthetic W7984×Opata reference populations provide insights into wheat structural diversity

Wheat (Triticum aestivum) genetic maps are a key enabling tool for genetic studies. We used genotyping-by-sequencing-(GBS) derived markers to map recombinant inbred line (RIL) and doubled haploid (DH) populations from crosses of W7984 by Opata, and used the maps to explore features of recombination control. The RIL and DH populations, SynOpRIL and SynOpDH, were composed of 906 and 92 individuals, respectively. Two high-density genetic linkage framework maps were constructed of 2,842 and 2,961 cM, harboring 3,634 and 6,580 markers, respectively. Using imputation, we added 43,013 and 86,042 markers to the SynOpRIL and SynOpDH maps. We observed preferential recombination in telomeric regions and reduced recombination in pericentromeric regions. Recombination rates varied between subgenomes, with the D genomes of the two populations exhibiting the highest recombination rates of 0.26-0.27 cM/Mb. QTL mapping identified two additive and three epistatic loci associated with crossover number. Additionally, we used published POPSEQ data from SynOpDH to explore the structural variation in W7984 and Opata. We found that chromosome 5AS is missing from W7984. We also found 2,332 variations larger than 100 kb. Structural variants were more abundant in distal regions, and overlapped 9,196 genes. The two maps provide a resource for trait mapping and genomic-assisted breeding.

Acid phosphatase gene GmHAD1 linked to low phosphorus tolerance in soybean, through fine mapping

KEY MESSAGE

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress. Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23-13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80-13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4-16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

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

KEY MESSAGE

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

Detecting the QTL-allele system conferring flowering date in a nested association mapping population of soybean using a novel procedure

KEY MESSAGE

The RTM-GWAS was chosen among five procedures to identify DTF QTL-allele constitution in a soybean NAM population; 139 QTLs with 496 alleles accounting for 81.7% of phenotypic variance were detected. Flowering date (days to flowering, DTF) is an ecological trait in soybean, closely related to its ability to adapt to areas. A nested association mapping (NAM) population consisting of four RIL populations (LM, ZM, MT and MW with M8206 as their common parent) was established and tested for their DTF under five environments. Using restriction-site-associated DNA sequencing the population was genotyped with SNP markers. The restricted two-stage multi-locus (RTM) genome-wide association study (GWAS) (RTM-GWAS) with SNP linkage disequilibrium block (SNPLDB) as multi-allele genomic markers performed the best among the five mapping procedures with software publicly available. It identified the greatest number of quantitative trait loci (QTLs) (139) and alleles (496) on 20 chromosomes covering almost all of the QTLs detected by four other mapping procedures. The RTM-GWAS provided the detected QTLs with highest genetic contribution but without overflowing and missing heritability problems (81.7% genetic contribution vs. heritability of 97.6%), while SNPLDB markers matched the NAM population property of multiple alleles per locus. The 139 QTLs with 496 alleles were organized into a QTL-allele matrix, showing the corresponding DTF genetic architecture of the five parents and the NAM population. All lines and parents comprised both positive and negative alleles, implying a great potential of recombination for early and late DTF improvement. From the detected QTL-allele system, 126 candidate genes were annotated and χ 2 tested as a DTF candidate gene system involving nine biological processes, indicating the trait a complex, involving several biological processes rather than only a handful of major genes.

Dicamba-Tolerant Soybeans (Glycine max L.) MON 87708 and MON 87708 × MON 89788 Are Compositionally Equivalent to Conventional Soybean

Herbicide-tolerant crops can expand both tools for and timing of weed control strategies. MON 87708 soybean has been developed through genetic modification and confers tolerance to the dicamba herbicide. As part of the safety assessment conducted for new genetically modified (GM) crop varieties, a compositional assessment of MON 87708 was performed. Levels of key soybean nutrients and anti-nutrients in harvested MON 87708 were compared to levels of those components in a closely related non-GM variety as well as to levels measured in other conventional soybean varieties. From this analysis, MON 87708 was shown to be compositionally equivalent to its comparator. A similar analysis conducted for a stacked trait product produced by conventional breeding, MON 87708 × MON 89788, which confers tolerance to both dicamba and glyphosate herbicides, reached the same conclusion. These results are consistent with other results that demonstrate no compositional impact of genetic modification, except in those cases where an impact was an intended outcome.

An R2R3-type MYB transcription factor, GmMYB29, regulates isoflavone biosynthesis in soybean

Isoflavones comprise a group of secondary metabolites produced almost exclusively by plants in the legume family, including soybean [Glycine max (L.) Merr.]. They play vital roles in plant defense and have many beneficial effects on human health. Isoflavone content is a complex quantitative trait controlled by multiple genes, and the genetic mechanisms underlying isoflavone biosynthesis remain largely unknown. Via a genome-wide association study (GWAS), we identified 28 single nucleotide polymorphisms (SNPs) that are significantly associated with isoflavone concentrations in soybean. One of these 28 SNPs was located in the 5'-untranslated region (5'-UTR) of an R2R3-type MYB transcription factor, GmMYB29, and this gene was thus selected as a candidate gene for further analyses. A subcellular localization study confirmed that GmMYB29 was located in the nucleus. Transient reporter gene assays demonstrated that GmMYB29 activated the IFS2 (isoflavone synthase 2) and CHS8 (chalcone synthase 8) gene promoters. Overexpression and RNAi-mediated silencing of GmMYB29 in soybean hairy roots resulted in increased and decreased isoflavone content, respectively. Moreover, a candidate-gene association analysis revealed that 11 natural GmMYB29 polymorphisms were significantly associated with isoflavone contents, and regulation of GmMYB29 expression could partially contribute to the observed phenotypic variation. Taken together, these results provide important genetic insights into the molecular mechanisms underlying isoflavone biosynthesis in soybean.

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

KEY MESSAGE

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

Exploiting Phenylpropanoid Derivatives to Enhance the Nutraceutical Values of Cereals and Legumes

Phenylpropanoids are a diverse chemical class with immense health benefits that are biosynthesized from the aromatic amino acid L-phenylalanine. This article reviews the progress for accessing variation in phenylpropanoids in germplasm collections, the genetic and molecular basis of phenylpropanoid biosynthesis, and the development of cultivars dense in seed-phenylpropanoids. Progress is also reviewed on high-throughput assays, factors that influence phenylpropanoids, the site of phenylpropanoids accumulation in seed, Genotype × Environment interactions, and on consumer attitudes for the acceptance of staple foods rich in phenylpropanoids. A paradigm shift was noted in barley, maize, rice, sorghum, soybean, and wheat, wherein cultivars rich in phenylpropanoids are grown in Europe and North and Central America. Studies have highlighted some biological constraints that need to be addressed for development of high-yielding cultivars that are rich in phenylpropanoids. Genomics-assisted breeding is expected to facilitate rapid introgression into improved genetic backgrounds by minimizing linkage drag. More research is needed to systematically characterize germplasm pools for assessing variation to support crop genetic enhancement, and assess consumer attitudes to foods rich in phenylpropanoids.

Glycone-rich Soy Isoflavone Extracts Promote Estrogen Receptor Positive Breast Cancer Cell Growth

Due to the association of hormone replacement therapy (HRT) with breast cancer risk, estrogenically active soy isoflavones are considered as an HRT alternative to alleviate menopausal symptoms. However, several recent reports challenged the health benefits of soy isoflavones and associated them with breast cancer promotion. While glyconic isoflavones are the major constituents of soybean seeds, due to their low cell permeability, they are considered to be biologically inactive. The glyconic isoflavones may exert their effects on membrane-bound estrogen receptors or could be converted to aglycones by extracellular β-glucosidases. Therefore, we hypothesized that despite their low cell permeability, soybean cultivars with high glyconic isoflavones may promote breast cancer cell growth. To test this, composition and estrogenic activity of isoflavones from 54 commercial soybean cultivars were determined. Soybean seeds produced in identical climate and growth conditions were used to minimize the effects of extraneous factors on isoflavone profile and concentrations. The glyconic daidzin concentration negatively correlated with genistin and with other aglycones. Relative to control, isoflavone extracts from 51 cultivars were estrogenic and promoted the growth of estrogen receptor positive (ER+) breast cancer cell line MCF-7 from 1.14 to 4.59 folds and other three cultivars slightly reduced the growth. Among these, extracts from three cultivars were highly estrogenic and promoted MCF-7 cell growth by 2.59-4.64 folds (P<0.005). Among six isoflavones, daidzin was positively associated with MCF-7 cell growth (P<0.005, r = 0.13966), whereas the negative correlation between genistin and MCF-7 cell growth was nearly significant (P≤0.0562, r = -0.026141). Furthermore, in drug interaction studies daidzin-rich isoflavone extracts antagonized tamoxifen, an ER inhibitor. Taken together, our results suggest that the glyconic daidzin-rich soy isoflavone extracts may exert estrogenic effects and promote ER+ breast cancer growth.

QTL Location and Epistatic Effect Analysis of 100-Seed Weight Using Wild Soybean (Glycine soja Sieb. & Zucc.) Chromosome Segment Substitution Lines

Increasing the yield of soybean (Glycine max L. Merrill) is a main aim of soybean breeding. The 100-seed weight is a critical factor for soybean yield. To facilitate genetic analysis of quantitative traits and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population consisting of 194 chromosome segment substitution lines (CSSLs) was developed. In these lines, different chromosomal segments of the Chinese cultivar Suinong 14 were substituted into the genetic background of wild soybean (Glycine soja Sieb. & Zucc.) ZYD00006. Based on these CSSLs, a genetic map covering the full genome was generated using 121 simple sequence repeat (SSR) markers. In the quantitative trait loci (QTL) analysis, twelve main effect QTLs (qSW-B1-1/2/3, qSW-D1b-1/2, qSW-D2-1/2, qSW-G-1/2/3, qSW-M-2 and qSW-N-2) underlying 100-seed weight were identified in 2011 and 2012. The epistatic effects of pairwise interactions between markers were analyzed in 2011 and 2012. The results clearly demonstrated that these CSSLs could be used to identify QTLs, and that an epistatic analysis was able to detect several sites with important epistatic effects on 100-seed weight. Thus, we identified loci that will be valuable for improving soybean 100-seed weight. These results provide a valuable foundation for identifying the precise location of genes of interest, and for designing cloning and marker-assisted selection breeding strategies targeting the 100-seed weight of soybean.

Evaluation of genetic variation among Brazilian soybean cultivars through genome resequencing

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genetical genomics of Populus leaf shape variation

BACKGROUND

Leaf morphology varies extensively among plant species and is under strong genetic control. Mutagenic screens in model systems have identified genes and established molecular mechanisms regulating leaf initiation, development, and shape. However, it is not known whether this diversity across plant species is related to naturally occurring variation at these genes. Quantitative trait locus (QTL) analysis has revealed a polygenic control for leaf shape variation in different species suggesting that loci discovered by mutagenesis may only explain part of the naturally occurring variation in leaf shape. Here we undertook a genetical genomics study in a poplar intersectional pseudo-backcross pedigree to identify genetic factors controlling leaf shape. The approach combined QTL discovery in a genetic linkage map anchored to the Populus trichocarpa reference genome sequence and transcriptome analysis.

RESULTS

A major QTL for leaf lamina width and length:width ratio was identified in multiple experiments that confirmed its stability. A transcriptome analysis of expanding leaf tissue contrasted gene expression between individuals with alternative QTL alleles, and identified an ADP-ribosylation factor (ARF) GTPase (PtARF1) as a candidate gene for regulating leaf morphology in this pedigree. ARF GTPases are critical elements in the vesicular trafficking machinery. Disruption of the vesicular trafficking function of ARF by the pharmacological agent Brefeldin A (BFA) altered leaf lateral growth in the narrow-leaf P. trichocarpa suggesting a molecular mechanism of leaf shape determination. Inhibition of the vesicular trafficking processes by BFA interferes with cycling of PIN proteins and causes their accumulation in intercellular compartments abolishing polar localization and disrupting normal auxin flux with potential effects on leaf expansion.

CONCLUSIONS

In other model systems, ARF proteins have been shown to control the localization of auxin efflux carriers, which function to establish auxin gradients and apical-basal cell polarity in developing plant organs. Our results support a model where PtARF1 transcript abundance changes the dynamics of endocytosis-mediated PIN localization in leaf cells, thus affecting lateral auxin flux and subsequently lamina leaf expansion. This suggests that evolution of differential cellular polarity plays a significant role in leaf morphological variation observed in subgenera of genus Populus.

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

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

Heterosis and combining ability estimates in isoflavone content using different parental soybean accessions: wild soybean, a valuable germplasm for soybean breeding

Isoflavone, a group of secondary metabolites in soybean, is beneficial to human health. Improving isoflavone content in soybean seeds has become one of the most important breeding objectives. However, the narrow genetic base of soybean cultivars hampered crop improvement. Wild soybean is an extraordinarily important gene pool for soybean breeding. In order to select an optimal germplasm for breeding programs to increase isoflavone concentration, 36 F1 soybean progenies from different parental accessions (cultivars, wild, Semi-wild and Interspecific) with various total isoflavone (TIF) concentration (High, Middle, Low) were analyzed for their isoflavone content. Results showed that male parents, except for Cultivars, showed positive GCA effects. In particular, wild soybean had higher positive GCA effects for TIF concentration. Both MP and BP heterosis value declined in the hybrid in which male parents were wild soybean, semi-wild soybean, interspecific offspring and cultivar in turn. In general, combining ability and heterosis in hybrids which had relative higher TIF concentration level parents showed better performance than those which had lower TIF concentration level parents. These results indicated characteristics of isoflavone content were mainly governed by additive type of gene action, and wild relatives could be utilized for breeding of soybean cultivars with this trait. A promising combination was found as the best potential hybrid for isoflavone content improvement.

Profile of the contents of different forms of soybean isoflavones and the effect of germination time on these compounds and the physical parameters in soybean sprouts

The aim of this study was to evaluate the profile of the contents of different forms of soybean isoflavones and the effect of germination time on these compounds and the physical parameters in BRS 284 soybean sprouts. Soybean seeds were germinated for 168 h, and the sprouts were collected every 24 h. The physical parameters and contents of different forms of isoflavones of the seeds and soybean sprouts were evaluated, and the data were subjected to regression analysis. The soybean seeds contained 26.0% β-glucosides, 72.9% malonylglucosides and 1.2% aglycones. The yield of soybean sprouts was 632.4%. The effect of germination time was quadratic on the length, moisture and on the daidzin, genistin and genistein content; linear on the fresh weight and on the malonyldaidzin content. The dry matter and malonylglycitin content was constant, and glycitin and glycitein were not detected in the soybean sprouts.

Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max

BACKGROUND

Quantitative trait locus (QTL) mapping is an efficient approach to discover the genetic architecture underlying complex quantitative traits. However, the low density of molecular markers in genetic maps has limited the efficiency and accuracy of QTL mapping. In this study, specific length amplified fragment sequencing (SLAF-seq), a new high-throughput strategy for large-scale SNP discovery and genotyping based on next generation sequencing (NGS), was employed to construct a high-density soybean genetic map using recombinant inbred lines (RILs, Luheidou2×Nanhuizao, F5:8). With this map, the consistent QTLs for isoflavone content across various environments were identified.

RESULTS

In total, 23 Gb of data containing 87,604,858 pair-end reads were obtained. The average coverage for each SLAF marker was 11.20-fold for the female parent, 12.51-fold for the male parent, and an average of 3.98-fold for individual RILs. Among the 116,216 high-quality SLAFs obtained, 9,948 were polymorphic. The final map consisted of 5,785 SLAFs on 20 linkage groups (LGs) and spanned 2,255.18 cM in genome size with an average distance of 0.43 cM between adjacent markers. Comparative genomic analysis revealed a relatively high collinearity of 20 LGs with the soybean reference genome. Based on this map, 41 QTLs were identified that contributed to the isoflavone content. The high efficiency and accuracy of this map were evidenced by the discovery of genes encoding isoflavone biosynthetic enzymes within these loci. Moreover, 11 of these 41 QTLs (including six novel loci) were associated with isoflavone content across multiple environments. One of them, qIF20-2, contributed to a majority of isoflavone components across various environments and explained a high amount of phenotypic variance (8.7%-35.3%). This represents a novel major QTL underlying isoflavone content across various environments in soybean.

CONCLUSIONS

Herein, we reported a high-density genetic map for soybean. This map exhibited high resolution and accuracy. It will facilitate the identification of genes and QTLs underlying essential agronomic traits in soybean. The novel major QTL for isoflavone content is useful not only for further study on the genetic basis of isoflavone accumulation, but also for marker-assisted selection (MAS) in soybean breeding in the future.

Quantitative trait loci analysis of individual and total isoflavone contents in soybean seeds

Soybean isoflavones play diverse roles in human health, including cancers, osteoporosis, heart disease, menopausal symptoms and pabulums. The objective of this study was to identify the quantitative trait loci (QTL) associated with the isoflavones daidzein (DC), genistein (GeC), glycitein (GlC) and total isoflavone contents (TIC) in soybean seeds. A population of 184 F2:10 recombinant inbred lines derived from a 'Xiaoheidou' x 'GR8836' cross was planted in pot and field conditions to evaluate soybean isoflavones. Twenty-one QTL were detected by composite interval mapping. Several QTL were associated with the traits for DC, GeC, GlC and TIC only. QDGeGlTIC4_1 and QDGlTIC12_1 are reported first in this study and were associated with the DC, GeC, GlC and TIC traits simultaneously. The QTL identified have potential value for marker-assisted selection to develop soybean varieties with desirable isoflavone content.

Analytical method evaluation and discovery of variation within maize varieties in the context of food safety: transcript profiling and metabolomics

Profiling techniques such as microarrays, proteomics, and metabolomics are used widely to assess the overall effects of genetic background, environmental stimuli, growth stage, or transgene expression in plants. To assess the potential regulatory use of these techniques in agricultural biotechnology, we carried out microarray and metabolomic studies of 3 different tissues from 11 conventional maize varieties. We measured technical variations for both microarrays and metabolomics, compared results from individual plants and corresponding pooled samples, and documented variations detected among different varieties with individual plants or pooled samples. Both microarray and metabolomic technologies are reproducible and can be used to detect plant-to-plant and variety-to-variety differences. A pooling strategy lowered sample variations for both microarray and metabolomics while capturing variety-to-variety variation. However, unknown genomic sequences differing between maize varieties might hinder the application of microarrays. High-throughput metabolomics could be useful as a tool for the characterization of transgenic crops. However, researchers will have to take into consideration the impact on the detection and quantitation of a wide range of metabolites on experimental design as well as validation and interpretation of results.

Quantitative trait loci for seed isoflavone contents in 'MD96-5722' by 'Spencer' recombinant inbred lines of soybean

Isoflavones from soybeans [ Glycine max (L.) Merr.] have a significant impact on human health to reduce the risk of several major diseases. Breeding soybean for high isoflavone content in the seed is possible through marker-assisted selection (MAS) which can be based on quantitative trait loci (QTL). The objective of this study was to identify QTL controlling isoflavone content in a set of 'MD96-5722' by 'Spencer' recombinant inbred line (RIL) populations of soybean. Wide variations were found for seed concentrations of daidzein, glycitein, genistein, and total isoflavones among the RIL populations. Three QTL were identified on three different linkage groups (LG) represented by three different chromosomes (Chr). One QTL that controlled daidzein content was identified on LG A1 (Chr 5), and two QTL that underlay glycitein content were identified on LG K (Chr 9) and LG B2 (Chr 14). Identified QTL could be functional in developing soybean with preferable isoflavone concentrations in the seeds through MAS.

Mapping QTL, epistasis and genotype × environment interaction of antioxidant activity, chlorophyll content and head formation in domesticated lettuce (Lactuca sativa)

Fruits and vegetables are rich sources of antioxidants in human diets and their intake is associated with chronic disease prevention. Lettuce (Lactuca sativa L.) is a common vegetable in diets worldwide, but its nutritional content is relatively low. To elucidate the genetic basis of antioxidant content in lettuce, we measured the oxygen radical absorbance capacity (ORAC) and chlorophyll (Chl) content as a proxy of β-carotene in an F(8) recombinant inbred line (RIL) in multiple production cycles at two different production sites. Plants were phenotyped at the open-leaf stage to measure genetic potential (GP) or at market maturity (MM) to measure the influence of head architecture ('head' or 'open'). Main effect quantitative trait loci (QTL) were identified at MM (three Chl and one ORAC QTL) and GP (two ORAC QTL). No main effect QTL for Chl was detected at GP, but epistatic interaction was identified in one pair of marker intervals for each trait at GP. Interactions with environment were also detected for both main and epistatic effects (two for main effect, and one for epistatic effect). Main effect QTL for plant architecture and nutritional traits at MM colocated to a single genomic region. Chlorophyll contents and ORAC values at MM were significantly higher and Chl a to Chl b ratios were lower in 'open' types compared to 'head' types. The nutritional traits assessed for GP showed a significant association with plant architecture suggesting pleiotropic effects or closely linked genes. Taken together, the antioxidant and chlorophyll content of lettuce is controlled by complex mechanisms and participating alleles change depending on growth stage and production environment.

Major locus and other novel additive and epistatic loci involved in modulation of isoflavone concentration in soybean seeds

Seeds of soybean [Glycine max (L.) Merr.] accumulate more isoflavones than any tissue of any plant species. In other plant parts, isoflavones are usually released to counteract the effects of various biotic and abiotic stresses. Because of the benefits to the plant and positive implications that consumption may have on human health, increasing isoflavones is a goal of many soybean breeding programs. However, altering isoflavone levels through marker-assisted selection (MAS) has been impractical due to the small and often environmentally variable contributions that each individual quantitative trait locus (QTL) has on total isoflavones. In this study, we developed a Magellan × PI 437654 F(7)-RIL population to construct a highly saturated non-redundant linkage map that encompassed 451 SNP and SSR molecular markers and used it to locate genomic regions that govern accumulation of isoflavones in the seeds of soybean. Five QTLs were found that contribute to the concentration of isoflavones, having single or multiple additive effects on isoflavone component traits. We also validated a major locus which alone accounted for up to 10% of the phenotypic variance for glycitein, and 35-37% for genistein, daidzein and the sum of all three soybean isoflavones. This QTL was consistently associated with increased concentration of isoflavones across different locations, years and crosses. It was the most important QTL in terms of net increased amounts of all isoflavone forms. Our results suggest that this locus would be an excellent candidate to target for MAS. Also, several minor QTLs were identified that interacted in an additive-by-additive epistatic manner, to increase isoflavone concentration.

Compositional variability in conventional and glyphosate-tolerant soybean (Glycine max L.) varieties grown in different regions in Brazil

The compositions of a diverse range of commercially available conventional and genetically modified (GM; glyphosate-tolerant) soybean varieties from maturity groups 8 and 5, respectively, grown in the northern and southern soybean regions of Brazil during the 2007-2008 and 2008-2009 growing seasons were compared. Compositional analyses included measurement of essential macro- and micronutrients, antinutrients, and selected secondary metabolites in harvested seed as well as measurement of proximates in both forage and harvested seed. Statistical comparisons utilized a mixed analysis of variance model to evaluate the relative contributions of growing season, soybean growing region, production site, phenotype (GM or conventional), and variety. The study highlighted extensive variability in the overall data set particularly for components such as fatty acids, vitamin E, and isoflavones. There were few differences between the GM and non-GM populations, and most of the variability in the data set could be attributed to regional and variety differences. Overall, the results were consistent with the expanding literature on the lack of any meaningful impact of transgene insertion on crop composition.

Stability in the composition equivalence of grain from insect-protected maize and seed from glyphosate-tolerant soybean to conventional counterparts over multiple seasons, locations, and breeding germplasms

Insect-protected maize MON 810 and Roundup Ready soybean 40-3-2 represent major milestones in the adoption of genetically modified (GM) crops to enhance agricultural productivity. This study provides an assessment of the compositional stability of these products over multiple seasons, multiple germplasms, and diverse geographies encompassing North, Central, and South America and Europe. The compositional assessment evaluated levels of proximates in MON 810 and proximates, antinutrients, and isoflavones in 40-3-2. The means and range values for component levels in the GM crops and their conventional comparators were consistently similar to each other within each corresponding year from 2000 to 2009. To our knowledge, this study represents the first meta-analysis of comparative composition assessments of GM products. This approach, combined with graphical approaches, provided an effective summary of the overall data set and confirmed the continued compositional equivalence of these important crops to their conventional counterparts over time.

Metabolic profiling of strawberry (Fragaria x ananassa Duch.) during fruit development and maturation

Strawberry (Fragaria × ananassa Duch), a fruit of economic and nutritional importance, is also a model species for fleshy fruits and genomics in Rosaceae. Strawberry fruit quality at different harvest stages is a function of the fruit's metabolite content, which results from physiological changes during fruit growth and ripening. In order to investigate strawberry fruit development, untargeted (GC-MS) and targeted (HPLC) metabolic profiling analyses were conducted. Principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were employed to explore the non-polar and polar metabolite profiles from fruit samples at seven developmental stages. Different cluster patterns and a broad range of metabolites that exerted influence on cluster formation of metabolite profiles were observed. Significant changes in metabolite levels were found in both fruits turning red and fruits over-ripening in comparison with red-ripening fruits. The levels of free amino acids decreased gradually before the red-ripening stage, but increased significantly in the over-ripening stage. Metabolite correlation and network analysis revealed the interdependencies of individual metabolites and metabolic pathways. Activities of several metabolic pathways, including ester biosynthesis, the tricarboxylic acid cycle, the shikimate pathway, and amino acid metabolism, shifted during fruit growth and ripening. These results not only confirmed published metabolic data but also revealed new insights into strawberry fruit composition and metabolite changes, thus demonstrating the value of metabolomics as a functional genomics tool in characterizing the mechanism of fruit quality formation, a key developmental stage in most economically important fruit crops.

Assessing the natural variability in crop composition

The number of evaluations of the nutrient composition of food and feed crops has increased over the past 15years due to the introduction of new crops using the tools of modern biotechnology. The composition of these crops has been extensively compared with conventional (non-transgenic) controls as an integral part of the comparative safety assessment process. Following guidelines outlined in the Organization of Economic Co-operation and Development (OECD) Consensus Documents, most of these studies have incorporated field trials at multiple geographies and a diverse range of commercially available varieties/hybrids that are analyzed to understand natural variability in composition due to genetic and environmental influences. Using studies conducted in the US, Argentina and Brazil over multiple growing seasons, this report documents the effect of geography, growing season, and genetic background on soybean composition where fatty acids and isoflavones were shown to be particularly variable. A separate investigation of 96 different maize hybrids grown at three locations in the US demonstrated that levels of free amino acids, sugars/polyols, and molecules associated with stress response can vary to a greater degree than that observed for more abundant components. The International Life Sciences Institute (ILSI) crop composition database has proven to be an important resource for collecting and disseminating nutrient composition data to promote a further understanding of the variability that occurs naturally in crops used for food and feed.

Intricate environment-modulated genetic networks control isoflavone accumulation in soybean seeds

BACKGROUND

Soybean (Glycine max [L] Merr.) seed isoflavones have long been considered a desirable trait to target in selection programs for their contribution to human health and plant defense systems. However, attempts to modify seed isoflavone contents have not always produced the expected results because their genetic basis is polygenic and complex. Undoubtedly, the extreme variability that seed isoflavones display over environments has obscured our understanding of the genetics involved.

RESULTS

In this study, a mapping population of RILs with three replicates was analyzed in four different environments (two locations over two years). We found a total of thirty-five main-effect genomic regions and many epistatic interactions controlling genistein, daidzein, glycitein and total isoflavone accumulation in seeds. The use of distinct environments permitted detection of a great number of environment-modulated and minor-effect QTL. Our findings suggest that isoflavone seed concentration is controlled by a complex network of multiple minor-effect loci interconnected by a dense epistatic map of interactions. The magnitude and significance of the effects of many of the nodes and connections in the network varied depending on the environmental conditions. In an attempt to unravel the genetic architecture underlying the traits studied, we searched on a genome-wide scale for genomic regions homologous to the most important identified isoflavone biosynthetic genes. We identified putative candidate genes for several of the main-effect and epistatic QTL and for QTL reported by other groups.

CONCLUSIONS

To better understand the underlying genetics of isoflavone accumulation, we performed a large scale analysis to identify genomic regions associated with isoflavone concentrations. We not only identified a number of such regions, but also found that they can interact with one another and with the environment to form a complex adaptable network controlling seed isoflavone levels. We also found putative candidate genes in several regions and overall we advanced the knowledge of the genetics underlying isoflavone synthesis.

Differential expression of isoflavone biosynthetic genes in soybean during water deficits

Numerous environmental factors influence isoflavone accumulation and have long hampered their genetic dissection. Temperature and water regimes are two of the most significant abiotic factors. However, while the effects of temperature have been widely studied, little is known about how water scarcity might affect isoflavone concentration in seeds. Studies have shown that accumulation of isoflavones is promoted by well-watered conditions, but the molecular basis remains elusive. The length and severity of the water stress required to induce changes are also still unknown. In the present work, several intensities of water stress were evaluated at various critical stages for soybean [Glycine max (L.) Merr.] seed development, in both field and controlled environments. The results suggested that only long-term progressive drought, spanning most of the seed developmental stages, significantly decreased isoflavone content in seeds. The reduction is proportional to the intensity of the stress and appears to occur in a genotype-dependent manner. However, regardless of water regime, isoflavone compounds were mainly accumulated in the later seed developmental stages. Transcripts of the most important genes for isoflavone biosynthesis were also quantified from samples collected at key seed developmental stages under well-watered and long-term water deficit conditions. Expression of CHS7, CHS8 and IFS2 correlated with isoflavone accumulation under well-watered conditions. Interestingly, we found that the two isoflavone synthase genes in soybean (IFS1 and IFS2) showed different patterns of expression. The abundance of IFS1 transcripts was maintained at a constant rate, whereas IFS2 was down-regulated and highly correlated with isoflavone accumulation under both water deficit and well-watered conditions, suggesting IFS2 as a main contributor to isoflavone diminution under drought.