Soybean tocopherol concentrations are affected by crop management

Geographic distributions and the regionalization of soybean seed compositions across China

As one of major food crops, soybean is grown over a broad ecological region in China with considerable variations in environmental conditions, and the seed compositions of soybeans are diverse among different regions. To clarify the spatial patterns of soybean seed compositions, crude oil, protein, and 11 categories of functional components were quantified in 1792 soybean samples collected from a vast range of soybean planting regions across China spanning from 2010 to 2017. The Kriging interpolation maps presented a clear north-to-south (high latitude to low latitude) increasing trend in contents of crude protein and dietary fiber and decreasing trend in contents of crude oil, phospholipids, saponins, and carotenoids. Soybeans with high-level of total oligosaccharide were concentrated in the central region. Based on the geographical distribution of soybean nutritional components, weather conditions, and cultivation systems, the soybean production areas in China were divided into three regions and 10 subregions. This study highlights the geographic distribution of soybean nutritional compositions and provides scientific evidence for guiding the construction of high-quality edible soybean production bases in China.

Genome-wide association studies reveal novel QTLs, QTL-by-environment interactions and their candidate genes for tocopherol content in soybean seed

Genome-wide association studies (GWAS) is an efficient method to detect quantitative trait locus (QTL), and has dissected many complex traits in soybean [Glycine max (L.) Merr.]. Although these results have undoubtedly played a far-reaching role in the study of soybean biology, environmental interactions for complex traits in traditional GWAS models are frequently overlooked. Recently, a new GWAS model, 3VmrMLM, was established to identify QTLs and QTL-by-environment interactions (QEIs) for complex traits. In this study, the GLM, MLM, CMLM, FarmCPU, BLINK, and 3VmrMLM models were used to identify QTLs and QEIs for tocopherol (Toc) content in soybean seed, including δ-Tocotrienol (δ-Toc) content, γ-Tocotrienol (γ-Toc) content, α-Tocopherol (α-Toc) content, and total Tocopherol (T-Toc) content. As a result, 101 QTLs were detected by the above methods in single-environment analysis, and 57 QTLs and 13 QEIs were detected by 3VmrMLM in multi-environment analysis. Among these QTLs, some QTLs (Group I) were repeatedly detected three times or by at least two models, and some QTLs (Group II) were repeatedly detected only by 3VmrMLM. In the two Groups, 3VmrMLM was able to correctly detect all known QTLs in group I, while good results were achieved in Group II, for example, 8 novel QTLs were detected in Group II. In addition, comparative genomic analysis revealed that the proportion of Glyma_max specific genes near QEIs was higher, in other words, these QEIs nearby genes are more susceptible to environmental influences. Finally, around the 8 novel QTLs, 11 important candidate genes were identified using haplotype, and validated by RNA-Seq data and qRT-PCR analysis. In summary, we used phenotypic data of Toc content in soybean, and tested the accuracy and reliability of 3VmrMLM, and then revealed novel QTLs, QEIs and candidate genes for these traits. Hence, the 3VmrMLM model has broad prospects and potential for analyzing the genetic structure of complex quantitative traits in soybean.

Identification of Genomic Loci and Candidate Genes Related to Seed Tocopherol Content in Soybean

Soybean seeds are primary sources of natural tocopherols used by the food and pharmaceutical industries, owing to their beneficial impacts on human health. Selection for higher tocopherol contents in seeds along with other desirable traits is an important goal in soybean breeding. In order to identify the genomic loci and candidate genes controlling tocopherol content in soybean seeds, the bulked-segregant analysis technique was performed using a natural population of soybean consisting of 1525 accessions. We constructed the bulked-segregant analysis based on 98 soybean accessions that showed extreme phenotypic variation for the target trait, consisting of 49 accessions with extremely-high and 49 accessions with extremely-low tocopherol content. A total of 144 variant sites and 109 predicted genes related to tocopherol content were identified, in which a total of 83 genes were annotated by the gene ontology functions. Furthermore, 13 enriched terms (p < 0.05) were detected, with four of them found to be highly enriched: response to lipid, response to abscisic acid, transition metal ion transmembrane transporter activity, and double-stranded DNA binding. Especially, six candidate genes were detected at 41.8−41.9 Mb genomic hotspots on chromosome 5 based on ANNOtate VARiation analysis. Among the genes, only Glyma.05G243400 carried a non-synonymous mutation that encodes a “translation elongation factor EF1A or initiation factor IF2gamma family protein” was identified. The haplotype analysis confirmed that Glyma.05G243400 exhibited highly significant variations in terms of tocopherol content across multiple experimental locations, suggesting that it can be the key candidate gene regulating soybean seed tocopherols. The present findings provide novel gene resources related to seed tocopherols for further validation by genome editing, functional characterization, and genetic improvement targeting enhanced tocopherol composition in soybean molecular breeding.

QTL and Candidate Genes for Seed Tocopherol Content in ‘Forrest’ by ‘Williams 82’ Recombinant Inbred Line (RIL) Population of Soybean

Soybean seeds are rich in secondary metabolites which are beneficial for human health, including tocopherols. Tocopherols play an important role in human and animal nutrition thanks to their antioxidant activity. In this study, the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306) was used to map quantitative trait loci (QTL) for seed α-tocopherol, β-tocopherol, δ -tocopherol, γ-tocopherol, and total tocopherol contents in Carbondale, IL over two years. Also, the identification of the candidate genes involved in soybean tocopherols biosynthetic pathway was performed. A total of 32 QTL controlling various seed tocopherol contents have been identified and mapped on Chrs. 1, 2, 5, 6, 7, 8, 9, 10, 12, 13, 16, 17, and 20. One major and novel QTL was identified on Chr. 6 with an R² of 27.8, 9.9, and 6.9 for δ-tocopherol, α-tocopherol, and total tocopherol content, respectively. Reverse BLAST analysis of the genes that were identified in Arabidopsis allowed the identification of 37 genes involved in soybean tocopherol pathway, among which 11 were located close to the identified QTLs. The tocopherol cyclase gene (TC) Glyma.06G084100 is located close to the QTLs controlling δ-tocopherol (R² = 27.8), α-tocopherol (R² = 9.96), and total-tocopherol (R² = 6.95). The geranylgeranyl diphosphate reductase (GGDR) Glyma.05G026200 gene is located close to a QTL controlling total tocopherol content in soybean (R² = 4.42). The two methylphytylbenzoquinol methyltransferase (MPBQ-MT) candidate genes Glyma.02G002000 and Glyma.02G143700 are located close to a QTL controlling δ-tocopherol content (R² = 3.57). The two γ-tocopherol methyltransferase (γ-TMT) genes, Glyma.12G014200 and Glyma.12G014300, are located close to QTLs controlling (γ+ß) tocopherol content (R² = 8.86) and total tocopherol (R² = 5.94). The identified tocopherol seed QTLs and candidate genes will be beneficial in breeding programs to develop soybean cultivars with high tocopherol contents.

Soybean Callus—A Potential Source of Tocopherols

In vitro cultures have been used as an effective means to achieve a high level of secondary metabolites in various plants, including soy. In this study, the contents of α-, γ-, and δ- tocopherol were quantified in soybean callus, and their amounts were compared to those of soybeans cultivated using the conventional tillage system with three weed controls (respectively without herbicide and with two variants of herbicide). Soybean callus was produced using Murashige and Skoog 1962 (MS) medium supplemented with 0.1 mg/L 6-Benzylaminopurine (BAP) and 0. 1 mg/L Thidiazuron (TDZ). The highest amount of fresh callus was obtained from soybeans from the conventional tillage system with second weed control (S-metolachlor 960 g/L, imazamox 40 g/L, and propaquizafop 100 g/L) respectively 13,652.4 ± 1177.62 mg. The analyzed tocopherols were in much higher content in soy dry callus than the soybean seeds (5.63 µg/g compared with the 0.35 α-toco in soybean, 47.57 µg/g compared with 18.71 µg/g γ-toco or, 5.56 µg/g compared with 1.74 µg/g β-toco). The highest content of the three analyzed tocopherols was γ -tocopherol, both in callus and soybeans. Furthermore, the data showed that herbicides used in soybean culture significantly influenced both the in vitro callus production and the tocopherol callus content (p ˂ 0.05). Altogether, soybean callus can be an important source of tocopherols, and herbicides significantly influence in vitro callus production and the tocopherol callus content.

Identification of Loci and Candidate Genes Analyses for Tocopherol Concentration of Soybean Seed

Tocopherol (Toc) occurs in soybean seeds and is extracted together with the soybean oil. Toc is utilized as an antioxidant in food and an additive in animal feed. A total of 180 representative accessions and 144 recombinant inbred lines (RILs) from the cross of 'Hefeng 25' and 'OAC Bayfield' were selected to evaluate individuals and total Toc concentrations in soybean seeds. The 180 soybean samples were sequenced by the approach of Specific Locus Amplified Fragment Sequencing (SLAF-seq). A total of 22,611 single nucleotide polymorphisms (SNPs) were developed. Nineteen quantitative trait nucleotides (QTNs) were identified associated with individual or total-Toc based on genome-wide association analysis (GWAS). Among them, three QTNs located near known QTLs, and 16 were novel. Eighteen QTLs and nine eQTLs were also detected by linkage mapping. The QTN rs9337368 on Chr.02 was colocalized according to the linkage mapping of the RILs and genome-wide association analysis and regarded as a stable locus for mining the candidate genes in association with Toc. A total of 42 candidate genes near the 200 kbp flanking region of this identified locus were found. Upon a gene-based association, 11 SNPs from five genes out of the 42 candidates were detected. Expression level analysis of five candidate genes revealed that two genes were significantly related to Toc content. The identified loci, along with the candidate genes, might be valuable for increasing the Toc concentration in soybean seeds and improving the nutritional value of soybean oil.

Identification of quantitative trait loci for increased α-tocopherol biosynthesis in wild soybean using a high-density genetic map

BACKGROUND

Soybean is one of the most important crop sources of tocopherols (Toc). However, the content of α-Toc, an isoform with the highest vitamin E activity in humans, is low in most cultivars. With the aim of broadening genetic variability, we performed quantitative trait locus (QTL) analysis for a high seed α-Toc trait detected in a wild soybean and characterized the sequence polymorphisms and expression profiles of γ-tocopherol methyltransferase (γ-TMT) genes as potential candidates.

RESULTS

A recombinant inbred line population was developed from a cross between the low α-Toc breeding line TK780 and the high α-Toc wild accession B04009. The α-Toc content in seeds correlated strongly with the ratio of α-Toc to γ-Toc contents. QTL analysis using a high-density map constructed with 7710 single nucleotide polymorphisms (SNPs) generated by restriction site-associated DNA sequencing detected six QTLs involved in α-Toc biosynthesis. Of these, three in chromosomes (Chr) 9, 11, and 12 produced consistent effects during a 2-year trial. B04009 allele at QTLs in Chr9 and Chr12 and TK780 allele at the QTL in Chr11 each promoted the conversion of γ-Toc to α-Toc, which elevated the seed α-Toc content. SNPs and indels were detected between the parents in three γ-TMT genes (γ-TMT1, γ-TMT2, and γ-TMT3) co-located in the QTLs in Chr9 and Chr12, of which some existed in the cis-regulatory elements associated with seed development and functions. In immature cotyledons, γ-TMT3 was expressed at higher levels in B04009 than TK780, irrespective of two thermal conditions tested, whereas the expression of γ-TMT2 was markedly upregulated under higher temperatures, particularly in B04009.

CONCLUSIONS

We identified QTLs consistently controlling α-Toc biosynthesis in wild soybean seeds in 2-year trials. The QTL on Chr9 had been previously identified in soybean, whereas the QTLs on Chr11 and Chr12 were novel. Further molecular dissections and characterization of the QTLs may facilitate the use of high α-Toc alleles from wild soybean in soybean breeding and an understanding of the molecular mechanisms underlying α-Toc biosynthesis in soybean seeds.

Factors Affecting Tocopherol Concentrations in Soybean Seeds

Soybean seeds contain several health-beneficial compounds, including tocopherols, which are used by the nutraceutical and functional food industries. Soybean tocopherol concentrations are, however, highly variable. Large differences observed in tocopherol concentrations among soybean genotypes together with the relatively simple biosynthetic pathway involving few genes support the feasibility of selecting for high-tocopherol soybean. Tocopherol concentrations are also highly influenced by environmental factors and field management. Temperature during seed filling and soil moisture appear to be the main factors affecting tocopherol concentrations; other factors such as soil fertility and solar radiation also affect concentrations and composition. Field management decisions including seeding date, row spacing, irrigation, and fertilization also affect tocopherols. Knowledge of factors affecting soybean tocopherols is essential to develop management strategies that will lead to the production of seeds with consistent target concentrations that will meet the needs of the nutraceutical and functional food industries.

Natural diversity of seed α-tocopherol ratio in wild soybean (Glycine soja) germplasm collection

Alpha-tocopherol is one of four tocopherol isoforms and has the highest vitamin E activity in humans. Most cultivated soybean seeds contain γ-tocopherol as the predominant form, and the ratio of α-tocopherol content to total tocopherol content (α-tocopherol ratio) is <10%. Three soybean accessions from Eastern Europe have α-tocopherol ratios of >20%. This higher content is likely due to mutations in the promoter region of the γ-tocopherol methytransferase-3 (γ-TMT3) gene. We surveyed a wild soybean germplasm collection and detected 16 accessions with stable seed α-tocopherol ratios of >20% under different growth conditions. The α-tocopherol ratios were greatly reduced when the plants were grown under cool temperatures during seed maturation, but increased to varying degrees at higher temperatures. Sequence analysis of the γ-TMT3 promoter of 11 of the accessions identified four haplotypes, one of which corresponded to that of cultivars with higher contents. These wild accessions can thus serve as novel donors for breeding cultivars with high α-tocopherol ratios and for better understanding the genetic basis of α-tocopherol synthesis in soybean.

Tocochromanols in wood: a potential new tool for dendrometabolomics

Tocochromanols are the most abundant lipid-soluble antioxidants in plants. Among them, α-tocopherol (α-Toc) shows a particularly high sensitivity to environmental stressors and its content is used as a stress biomarker even in non-photosynthetic tissues. Nevertheless, the presence of tocochromanols has not been described yet in the xylem of woody plants, even when their functions regarding cell membrane protection and the transport of photoassimilates may be crucial in this tissue and despite its potential utility in dendrometabolomics. Considering all these, we aimed to determine the presence and distribution of tocochromanols in the xylem of woody plants, to examine their responsiveness to high temperature and to evaluate their potential as environmental bioindicators. The analysis of 29 phyllogenetically diverse species showed that α-Toc is the most abundant and frequent tocochromanol in the xylem and is ubiquitously present in all the studied species, with a concentration ranging from 0.5 to 39.3 μg g(-1) of dry weight. α-Tocopherol appeared to be mainly located in the parenchyma rays and was found in both the sapwood and the heartwood, suggesting that it is present even in dead parenchyma cells. The levels of α-Toc in the xylem did not change in response to locally induced xylem heating, but responded positively to the 3-year moving average of annual precipitation. The present findings suggest that α-Toc may be linked to changes in climatic stress. This should enhance further research on the environmental controls of α-Toc variation in the xylem as a first step towards a deeper understanding of dendrometabolomics.

Crop management, genotypes, and environmental factors affect soyasaponin B concentration in soybean

Soybean [Glycine max (L.) Merr.] seeds contain soyasaponin B, which has putative health benefits. Studies were conducted in multiple environments in Quebec, Canada to determine the effects of genotypes, environments, and seeding dates on soyasaponin B concentration in mature seeds. A growth chamber study was also conducted to determine the impact of high air temperature imposed at specific growth development stages on soyasaponin B in soybeans. Concentrations of individual and total soyasaponin B were determined using high-performance liquid chromatography. Genotype and environment main effects were the main determinants of soyasaponin B concentration in soybean, genotype × environment interactions accounting for less than 5% of the variation for all soyasaponin. Ranking of 20 early maturing soybean genotypes was thus relatively consistent across four environments, total concentration varying between 2.31 and 6.59 μmol g(-1). Seeding date consistently impacted soyasaponin B concentrations, early seeding date resulting in the highest concentrations. There was an 11% difference in total soyasaponin B concentration of soybeans seeded in mid-May compared to that in late-June. The response to high air temperature was complex and cultivar specific. High temperature stress restricted to the seed filling stages increased total soyasaponin B concentration in one cultivar by 28% when compared to that in control nonstressed plants; however, in another cultivar high temperature applied during all growth stages reduced total concentration by 27%. Results from the present study thus demonstrate that environmental factors and crop management both impact soyasaponin B concentration in soybeans.

Vitamin E levels in soybean (Glycine max (L.) Merr.) expressing a p-hydroxyphenylpyruvate gene from oat (Avena sativa L.)

The enzyme p-hydroxyphenylpyruvate dioxygenase (HPPD) is ubiquitous in plants and functions in the tyrosine catabolic pathway, resulting in the formation of homogentisate. Homogentisate is the aromatic precursor of all plastoquinones and tocochromanols, including tocopherols and tocotrienols. Soybean (Glycine max (L.) Merr.) has been genetically modified to express the gene avhppd-03 that encodes the protein AvHPPD-03 derived from oat (Avena sativa L.). The AvHPPD-03 isozyme has an inherent reduced binding affinity for mesotrione, a herbicide that inhibits the wild-type soybean HPPD enzyme. Expression of avhppd-03 in soybean plants confers a mesotrione-tolerant phenotype. Seeds from three different avhppd-03-expressing soybean events were quantitatively assessed for content of eight vitamin E isoforms. Although increased levels of two tocopherol isoforms were identified for each of the three soybean events, they were within, or not substantially different from, the ranges of these isoforms found in nontransgenic soybean varieties. The increases of these tocopherols in the avhppd-03-expressing soybean events may have a slight benefit with regard to vitamin E nutrition but, given the commercial processing of soybeans, are unlikely to have a material impact on human nutrition with regard to vitamin E concentrations in soybean oil.

Investigation of biochemical diversity in a soybean lineage representing 35 years of breeding

This paper reports an assessment of seed biochemical and metabolite variability and diversity in a series of nine soybean varieties; all lines share the same genetic lineage but represent ∼35 years of breeding (launch years 1972-2008) and differing yield potentials. These varieties, including six conventional and three glyphosate-tolerant lines, were grown concurrently at two replicated field sites in the United States during the 2011 growing season, and seeds were harvested at maturity. A compositional assessment included measurement of proximates, amino acids, fatty acids, tocopherols, isoflavones, saccharides, organic acids, and selected phytohormones. Statistical analysis included application of principal variance component analysis (PVCA) to investigate the interrelationships among compositional components from these soybean varieties and the impacts of location (environment) and pedigree on variability of these components. Results demonstrated that (i) some biochemical analytes showed trends (either increased or decreased) with launch year and/or yield, (ii) some analytes varied according to variety but showed no trend with launch year and/or yield, and (iii) almost all analytes showed extensive variation within and across sites. In summary, varietal development of high-yielding soybean, as represented in this study, has been accompanied by compositional changes but these are typically modest relative to environmental factors.

Safety of GM crops: compositional analysis

The compositional analysis of genetically modified (GM) crops has continued to be an important part of the overall evaluation in the safety assessment program for these materials. The variety and complexity of genetically engineered traits and modes of action that will be used in GM crops in the near future, as well as our expanded knowledge of compositional variability and factors that can affect composition, raise questions about compositional analysis and how it should be applied to evaluate the safety of traits. The International Life Sciences Institute (ILSI), a nonprofit foundation whose mission is to provide science that improves public health and well-being by fostering collaboration among experts from academia, government, and industry, convened a workshop in September 2012 to examine these and related questions, and a series of papers has been assembled to describe the outcomes of that meeting.

Analysis of isoflavone, phenolic, soyasapogenol, and tocopherol compounds in soybean [ Glycine max (L.) Merrill] germplasms of different seed weights and origins

This study investigated the functional compounds, including isoflavones, phenolics, soyasapogenols, and tocopherols, that were detected in 204 soybean [ Glycine max (L.) Merrill] germplasms. The soybean samples were divided into three groups according to origin: America, China, and Korea. The soybean samples were also classified into three groups on the basis of 100-seed weight: small (<13 g), medium (13-24 g), and large (>24 g). Among the soybean germplasms, CSRV121 (Bosukkong) had the highest level of isoflavone content (4778.1 μg g(-1)), whereas CS01316 had the lowest isoflavone content (682.4 μg g(-1)). Of the soybeans from the three different countries of origin, those from Korea showed the highest average concentration of total isoflavones (2252.6 μg g(-1)). The small seeds had the highest average total isoflavone concentration (2520.0 μg g(-1)) of the three different seed sizes. Among the 204 soybean germplasms, CS01405 had the highest content of total phenolics (5219.6 μg g(-1)), and CSRV017 (Hwangkeumkong) had the lowest phenolic content (654.6 μg g(-1)). The mean concentrations of total phenolic compounds were 2729.1 μg g(-1) in American soybean seeds, 1680.4 μg g(-1) in Chinese soybean seeds, and 1977.6 μg g(-1) in Korean soybean seeds. Of the soybean seeds from the three different countries of origin, American soybean seeds had the highest average concentration of total phenolic compounds, and Korean varieties showed the second highest value. Small soybean seeds had the highest average content of total phenolic compounds (2241.7 μg g(-1)), whereas medium-sized (1926.8 μg g(-1)) and large (1949.9 μg g(-1)) soybeans had lower concentrations of phenolic compounds. In whole soybean germplasms, the level of total soyasapogenols was higher in CS01173 (1802.3 μg g(-1)) and CS01346 (1736.8 μg g(-1)) than in the other types of soybeans. The mean concentrations of total soyasapogenol were 1234.0 μg g(-1) in American, 1294.5 μg g(-1) in Chinese, and 1241.5 μg g(-1) in Korean soybean varieties. Chinese soybean varieties showed the highest mean concentration of total soyasapogenol, and Korean soybean seeds showed the second highest level. The medium-seed group had the highest soyasapogenol content (1269.3 μg g(-1)) of the seeds that were grouped by size. A larger amount of soyasapogenol B than soyasapogenol A was detected. In whole soybeans, CS01202 showed the highest level of total tocopherols (330.5 μg g(-1)), whereas CSRV056 (Pungsannamulkong) had the lowest content (153.3 μg g(-1)). Chinese soybeans had the highest average concentration of total tocopherols (255.1 μg g(-1)). By comparison, the medium-sized Chinese soybean group had the highest (256.1 μg g(-1)) average total tocopherol content.