Soybean Seed Development: Fatty Acid and Phytohormone Metabolism and Their Interactions

A spatially resolved multiomic single-cell atlas of soybean development

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

Identification of the domestication gene GmCYP82C4 underlying the major quantitative trait locus for the seed weight in soybean

KEY MESSAGE

A major quantitative trait locus (QTL) for the hundred-seed weight (HSW) was identified and confirmed in the two distinct soybean populations, and the target gene GmCYP82C4 underlying this locus was identified that significantly associated with soybean seed weight, and it was selected during the soybean domestication and improvement process. Soybean is a major oil crop for human beings and the seed weight is a crucial goal of soybean breeding. However, only a limited number of target genes underlying the quantitative trait loci (QTLs) controlling seed weight in soybean are known so far. In the present study, six loci associated with hundred-seed weight (HSW) were detected in the first population of 573 soybean breeding lines by genome-wide association study (GWAS), and 64 gene models were predicted in these candidate QTL regions. The QTL qHSW_1 exhibits continuous association signals on chromosome four and was also validated by region association study (RAS) in the second soybean population (409 accessions) with wild, landrace, and cultivar soybean accessions. There were seven genes in qHSW_1 candidate region by linkage disequilibrium (LD) block analysis, and only Glyma.04G035500 (GmCYP82C4) showed specifically higher expression in flowers, pods, and seeds, indicating its crucial role in the soybean seed development. Significant differences in HSW trait were detected when the association panels are genotyped by single-nucleotide polymorphisms (SNPs) in putative GmCYP82C4 promoter region. Eight haplotypes were generated by six SNPs in GmCYP82C4 in the second soybean population, and two superior haplotypes (Hap2 and Hap4) of GmCYP82C4 were detected with average HSW of 18.27 g and 18.38 g, respectively. The genetic diversity of GmCYP82C4 was analyzed in the second soybean population, and GmCYP82C4 was most likely selected during the soybean domestication and improvement process, leading to the highest proportion of Hap2 of GmCYP82C4 both in landrace and cultivar subpopulations. The QTLs and GmCYP82C4 identified in this study provide novel genetic resources for soybean seed weight trait, and the GmCYP82C4 could be used for soybean molecular breeding to develop desirable seed weight in the future.

Regulation of seed traits in soybean

Soybean (Glycine max) is an essential economic crop that provides vegetative oil and protein for humans, worldwide. Increasing soybean yield as well as improving seed quality is of great importance. Seed weight/size, oil and protein content are the three major traits determining seed quality, and seed weight also influences soybean yield. In recent years, the availability of soybean omics data and the development of related techniques have paved the way for better research on soybean functional genomics, providing a comprehensive understanding of gene functions. This review summarizes the regulatory genes that influence seed size/weight, oil content and protein content in soybean. We also provided a general overview of the pleiotropic effect for the genes in controlling seed traits and environmental stresses. Ultimately, it is expected that this review will be beneficial in breeding improved traits in soybean.

Cytokinin Translocation to, and Biosynthesis and Metabolism within, Cereal and Legume Seeds: Looking Back to Inform the Future

Early in the history of cytokinins, it was clear that Zea mays seeds contained not just trans-zeatin, but its nucleosides and nucleotides. Subsequently, both pods and seeds of legumes and cereal grains have been shown to contain a complex of cytokinin forms. Relative to the very high quantities of cytokinin detected in developing seeds, only a limited amount appears to have been translocated from the parent plant. Translocation experiments, and the detection of high levels of endogenous cytokinin in the maternal seed coat tissues of legumes, indicates that cytokinin does not readily cross the maternal/filial boundary, indicating that the filial tissues are autonomous for cytokinin biosynthesis. Within the seed, trans-zeatin plays a key role in sink establishment and it may also contribute to sink strength. The roles, if any, of the other biologically active forms of cytokinin (cis-zeatin, dihydrozeatin and isopentenyladenine) remain to be elucidated. The recent identification of genes coding for the enzyme that leads to the biosynthesis of trans-zeatin in rice (OsCYP735A3 and 4), and the identification of a gene coding for an enzyme (CPN1) that converts trans-zeatin riboside to trans-zeatin in the apoplast, further cements the key role played by trans-zeatin in plants.

Differential metabolic reprogramming in developing soybean embryos in response to nutritional conditions and abscisic acid

Seed storage compound deposition is influenced by both maternal and filial tissues. Within this framework, we analyzed strategies that operate during the development and filling of soybean embryos, using in vitro culture systems combined with metabolomics and proteomics approaches. The carbon:nitrogen ratio (C:N) of the maternal supply and the hormone abscisic acid (ABA) are specific and interacting signals inducing differential metabolic reprogrammings linked to changes in the accumulation of storage macromolecules like proteins or oils. Differences in the abundance of sugars, amino acids, enzymes, transporters, transcription factors, and proteins involved in signaling were detected. Embryos adapted to the nutritional status by enhancing the metabolism of both carbon and nitrogen under lower C:N ratio condition or only carbon under higher C:N ratio condition. ABA turned off multiple pathways especially in high availability of amino acids, prioritizing the storage compounds biosynthesis. Common responses induced by ABA involved increased sucrose uptake (to increase the sink force) and oleosin (oil body structural component) accumulation. In turn, ABA differentially promoted protein degradation under lower nitrogen supply in order to sustain the metabolic demands. Further, the operation of a citrate shuttle was suggested by transcript quantification and enzymatic activity measurements. The results obtained are useful to help define biotechnological tools and technological approaches to improve oil and protein yields, with direct impact on human and animal nutrition as well as in green chemistry.

Kinase CIPK9 integrates glucose and abscisic acid signaling to regulate seed oil metabolism in rapeseed

Rapeseed (Brassica napus), an important oil crop worldwide, provides large amounts of lipids for human requirements. Calcineurin B-like (CBL)-interacting protein kinase 9 (CIPK9) was reported to regulate seed oil content in the plant. Here, we generated gene-silenced lines through RNA interference biotechnology and loss-of-function mutant bnacipk9 using CRISPR/Cas9 to further study BnaCIPK9 functions in the seed oil metabolism of rapeseeds. We discovered that compared with wild-type (WT) lines, gene-silenced and bnacipk9 lines had substantially different oil contents and fatty acid compositions: seed oil content was improved by 3%-5% and 1%-6% in bnacipk9 lines and gene-silenced lines, respectively; both lines were with increased levels of monounsaturated fatty acids and decreased levels of polyunsaturated fatty acids. Additionally, hormone and glucose content analyses revealed that compared with WT lines the bnacipk9 lines showed significant differences: in bnacipk9 seeds, indoleacetic acid and abscisic acid (ABA) levels were higher; glucose and sucrose contents were higher with a higher hexose-to-sucrose ratio in bnacipk9 mid-to-late maturation development seeds. Furthermore, the bnacipk9 was less sensitive to glucose and ABA than the WT according to stomatal aperture regulation assays and the expression levels of genes involved in glucose and ABA regulating pathways in rapeseeds. Notably, in Arabidopsis (Arabidopsis thaliana), exogenous ABA and glucose imposed on developing seeds revealed the effects of ABA and glucose signaling on seed oil accumulation. Altogether, our results strongly suggest a role of CIPK9 in mediating the interaction between glucose flux and ABA hormone signaling to regulate seed oil metabolism in rapeseed.

Bioengineering of Soybean Oil and Its Impact on Agronomic Traits

Soybean is a major oil crop and is also a dominant source of nutritional protein. The 20% seed oil content (SOC) of soybean is much lower than that in most oil crops and the fatty acid composition of its native oil cannot meet the specifications for some applications in the food and industrial sectors. Considerable effort has been expended on soybean bioengineering to tailor fatty acid profiles and improve SOC. Although significant advancements have been made, such as the creation of high-oleic acid soybean oil and high-SOC soybean, those genetic modifications have some negative impacts on soybean production, for instance, impaired germination or low protein content. In this review, we focus on recent advances in the bioengineering of soybean oil and its effects on agronomic traits.

Genome-wide identification of the soybean cytokinin oxidase/dehydrogenase gene family and its diverse roles in response to multiple abiotic stress

Cytokinin oxidase/dehydrogenase (CKX) irreversibly degrades cytokinin, regulates growth and development, and helps plants to respond to environmental stress. Although the CKX gene has been well characterized in various plants, its role in soybean remains elusive. Therefore, in this study, the evolutionary relationship, chromosomal location, gene structure, motifs, cis-regulatory elements, collinearity, and gene expression patterns of GmCKXs were analyzed using RNA-seq, quantitative real-time PCR (qRT-PCR), and bioinformatics. We identified 18 GmCKX genes from the soybean genome and grouped them into five clades, each comprising members with similar gene structures and motifs. Cis-acting elements involved in hormones, resistance, and physiological metabolism were detected in the promoter regions of GmCKXs. Synteny analysis indicated that segmental duplication events contributed to the expansion of the soybean CKX family. The expression profiling of the GmCKXs genes using qRT-PCR showed tissue-specific expression patterns. The RNA-seq analysis also indicated that GmCKXs play an important role in response to salt and drought stresses at the seedling stage. The responses of the genes to salt, drought, synthetic cytokinin 6-benzyl aminopurine (6-BA), and the auxin indole-3-acetic acid (IAA) at the germination stage were further evaluated by qRT-PCR. Specifically, the GmCKX14 gene was downregulated in the roots and the radicles at the germination stage. The hormones 6-BA and IAA repressed the expression levels of GmCKX1, GmCKX6, and GmCKX9 genes but upregulated the expression levels of GmCKX10 and GmCKX18 genes. The three abiotic stresses also decreased the zeatin content in soybean radicle but enhanced the activity of the CKX enzymes. Conversely, the 6-BA and IAA treatments enhanced the CKX enzymes' activity but reduced the zeatin content in the radicles. This study, therefore, provides a reference for the functional analysis of GmCKXs in soybean in response to abiotic stresses.

Role of cytokinins in seed development in pulses and oilseed crops: Current status and future perspective

Cytokinins constitutes a vital group of plant hormones regulating several developmental processes, including growth and cell division, and have a strong influence on grain yield. Chemically, they are the derivatives of adenine and are the most complex and diverse group of hormones affecting plant physiology. In this review, we have provided a molecular understanding of the role of cytokinins in developing seeds, with special emphasis on pulses and oilseed crops. The importance of cytokinin-responsive genes including cytokinin oxidases and dehydrogenases (CKX), isopentenyl transferase (IPT), and cytokinin-mediated genetic regulation of seed size are described in detail. In addition, cytokinin expression in germinating seeds, its biosynthesis, source-sink dynamics, cytokinin signaling, and spatial expression of cytokinin family genes in oilseeds and pulses have been discussed in context to its impact on increasing economy yields. Recently, it has been shown that manipulation of the cytokinin-responsive genes by mutation, RNA interference, or genome editing has a significant effect on seed number and/or weight in several crops. Nevertheless, the usage of cytokinins in improving crop quality and yield remains significantly underutilized. This is primarily due to the multigene control of cytokinin expression. The information summarized in this review will help the researchers in innovating newer and more efficient ways of manipulating cytokinin expression including CKX genes with the aim to improve crop production, specifically of pulses and oilseed crops.

ABSCISIC ACID-INSENSITIVE 5-ω3 FATTY ACID DESATURASE3 module regulates unsaturated fatty acids biosynthesis in Paeonia ostii

Paeonia ostii is an authorized novel vegetable oil crop due to its seeds rich in unsaturated fatty acids (UFAs) especially α-linolenic acid (ALA), which overweight the current available edible oil. However, little is known on the regulation mechanism of UFAs biosynthesis during its seed development. Here, we used transcriptome and proteome data combining phytochemistry means to uncover the relationship between abscisic acid (ABA) signaling and UFAs biosynthesis during P. ostii seed development. Based on transcriptome and proteome analysis, two desaturases of omega-6 and omega-3 fatty acid, named as PoFAD2 and PoFAD3 responsible for ALA biosynthesis were identified. Then, an ABSCISIC ACID-INSENSITIVE 5 (ABI5) proteins was identified as an upstream transcriptional factor, which activated the expression of PoFAD3 instead of PoFAD2. Moreover, silencing of PoABI5 repressed the response of PoFAD3 to ABA. This study provides the first view on the connection between the function of ABA signaling factors and ALA biosynthesis in the P. ostii seed, which lays the foundation for studies on the regulatory mechanism of ABA signaling involved in the UFAs synthesis during seeds development, meanwhile, it will shed light on manipulation of ALA content for satisfying human demands on high quality of edible oil or healthy supplement.

Distribution analysis of jasmonic acid-related compounds in developing Glycine max L. (soybean) seeds using mass spectrometry imaging and liquid chromatography-mass spectrometry

INTRODUCTION

Jasmonic acid (JA) and its precursors are oxylipins derived from α-linolenic acid (αLA) and hexadecatrienoic acid, and regulate seed development. However, their spatial distribution in the developing Glycine max L. (soybean) seeds has not been elucidated.

OBJECTIVE

To investigate the distribution of JA-related compounds in the developing soybean seeds using desorption electrospray ionisation-mass spectrometry imaging (DESI-MSI) and liquid chromatography-electrospray ionisation-mass spectrometry (LC-ESI-MS) analyses.

METHODS

Cryosections of developing seeds were prepared using adhesive films, and subjected to DESI-MSI analysis. Verification of the DESI-MSI ion images were performed using DESI-tandem MSI (MS/MSI), LC-ESI-MS and tandem MS (MS/MS).

RESULTS

In the DESI-MSI mass spectrum, peaks matching the chemical formulae of αLA, 12-oxo-phytodienoic acid (OPDA), and 3-oxo-2-(2-(Z)-pentenyl)-cyclopentane-1-octanoic acid (OPC-8:0) were detected. These compounds were mainly distributed in the seed coat, especially near the hilum. This was consistent with the quantitative results obtained by LC-ESI-MS. While, DESI-MS/MSI and LC-ESI-MS/MS suggested the presence of isomers for OPDA and OPC-8:0. The effect of isomers on the DESI-MSI ion images was small for OPDA, and considerable for OPC-8:0.

CONCLUSION

These results demonstrated that free αLA, OPDA, and OPC-8:0 were the abundant JA-related compounds mainly distributed in the seed coat of the developing soybeans. OPDA and OPC-8:0 might exert a biological role in the seed coat. To the best of my knowledge, this is the first report on the accumulation of OPDA and OPC-8:0 in the seed coat. The combination of DESI-MSI and LC-ESI-MS is a useful tool for distribution analysis of JA-related compounds in the developing seeds.

The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis

Safflower (Carthamus tinctorius L.), a member of the Asteraceae, is a popular crop due to its high linoleic acid (LA) and flavonoid (such as hydroxysafflor yellow A) contents. Here, we report the first high-quality genome assembly (contig N50 of 21.23 Mb) for the 12 pseudochromosomes of safflower using single-molecule real-time sequencing, Hi-C mapping technologies and a genetic linkage map. Phyloge nomic analysis showed that safflower diverged from artichoke (Cynara cardunculus) and sunflower (Helianthus annuus) approximately 30.7 and 60.5 million years ago, respectively. Comparative genomic analyses revealed that uniquely expanded gene families in safflower were enriched for those predicted to be involved in lipid metabolism and transport and abscisic acid signalling. Notably, the fatty acid desaturase 2 (FAD2) and chalcone synthase (CHS) families, which function in the LA and flavonoid biosynthesis pathways, respectively, were expanded via tandem duplications in safflower. CarFAD2-12 was specifically expressed in seeds and was vital for high-LA content in seeds, while tandemly duplicated CarFAD2 genes were up-regulated in ovaries compared to CarFAD2-12, which indicates regulatory divergence of FAD2 in seeds and ovaries. CarCHS1, CarCHS4 and tandem-duplicated CarCHS5˜CarCHS6, which were up-regulated compared to other CarCHS members at early stages, contribute to the accumulation of major flavonoids in flowers. In addition, our data reveal multiple alternative splicing events in gene families related to fatty acid and flavonoid biosynthesis. Together, these results provide a high-quality reference genome and evolutionary insights into the molecular basis of fatty acid and flavonoid biosynthesis in safflower.

Inoculation with Bacillus amyloliquefaciens and mycorrhiza confers tolerance to drought stress and improve seed yield and quality of soybean plant

The present study aimed to evaluate the effect of Bacillus amyloliquefaciens and/or Arbuscular Mycorrhizal Fungi (AMF) as natural biofertilizers on biomass, yield, and seed nutritive quality of soybean (Giza 111). The conditions investigated include a well-watered (WW) control and irrigation withholding at the seed development stage (R5, after 90 days from sowing) (DS). Co-inoculation with B. amyloliquefaciens and AMF, resulted in the highest plant biomass and yield under WW and DS conditions. The nuclear DNA content analysis suggested that co-inoculation with B. amyloliquefaciens and AMF decreased the inhibition of drought stress on both the size and granularity of seed cells, which were comparable to the normal level. The single or co-inoculation with B. amyloliquefaciens and AMF increased the primary metabolites content and alleviated the drought-induced reduction in soluble sugars, lipids, protein and oil contents. Plant inoculation induced the expression of genes involved in lipid and protein biosynthesis, whereas an opposite trend was observed for genes involved in lipid and protein degradation, supporting the observed increase in lipid and protein content. Plant inoculated with B. amyloliquefaciens showed the highest α-amylase and β-amylase activities, indicating improved osmolyte (soluble sugar) synthesis, particularly under drought. Interestingly, single or co-inoculation further strengthen the positive effect of drought on the antioxidant and osmoprotectant levels, i.e. phenol, flavonoid, glycine betaine contents, and glutathione-S-transferase (GST) activity. As a result of stress release, there was a decrease in the level of stress hormones (abscisic acid, ABA) and an increase in gibberellin (GA), trans-zeatin-riboside (ZR), and indole acetic acid (IAA) in the seeds of inoculated plants. Additionally, the ATP content, hydrolytic activities of plasma membrane H+ -ATPase, Ca2+ -ATPase, and Mg2+ -ATPase were also increased by the inoculation.

Temporal changes in metabolism late in seed development affect biomass composition

The negative association between protein and oil production in soybean (Glycine max) seed is well-documented. However, this inverse relationship is based primarily on the composition of mature seed, which reflects the cumulative result of events over the course of soybean seed development and therefore does not convey information specific to metabolic fluctuations during developmental growth regimes. In this study, we assessed maternal nutrient supply via measurement of seed coat exudates and metabolite levels within the cotyledon throughout development to identify trends in the accumulation of central carbon and nitrogen metabolic intermediates. Active metabolic activity during late seed development was probed through transient labeling with 13C substrates. The results indicated: (1) a drop in lipid contents during seed maturation with a concomitant increase in carbohydrates, (2) a transition from seed filling to maturation phases characterized by quantitatively balanced changes in carbon use and CO2 release, (3) changes in measured carbon and nitrogen resources supplied maternally throughout development, (4) 13C metabolite production through gluconeogenic steps for sustained carbohydrate accumulation as the maternal nutrient supply diminishes, and (5) oligosaccharide biosynthesis within the seed coat during the maturation phase. These results highlight temporal engineering targets for altering final biomass composition to increase the value of soybeans and a path to breaking the inverse correlation between seed protein and oil content.

The soybean (Glycine max L.) cytokinin oxidase/dehydrogenase multigene family; Identification of natural variations for altered cytokinin content and seed yield

Cytokinins (CKs) play a fundamental role in regulating dynamics of organ source/sink relationships during plant development, including flowering and seed formation stages. As a result, CKs are key drivers of seed yield. The cytokinin oxidase/dehydrogenase (CKX) is one of the critical enzymes responsible for regulating plant CK levels by causing their irreversible degradation. Variation of CKX activity is significantly correlated with seed yield in many crop species while in soybean (Glycine max L.), the possible associations between CKX gene family members (GFMs) and yield parameters have not yet been assessed. In this study, 17 GmCKX GFMs were identified, and natural variations among GmCKX genes were probed among soybean cultivars with varying yield characteristics. The key CKX genes responsible for regulating CK content during seed filling stages of reproductive development were highlighted using comparative phylogenetics, gene expression analysis and CK metabolite profiling. Five of the seventeen identified GmCKX GFMs, showed natural variations in the form of single nucleotide polymorphisms (SNPs). The gene GmCKX7-1, with high expression during critical seed filling stages, was found to have a non-synonymous mutation (H105Q), on one of the active site residues, Histidine 105, previously reported to be essential for co-factor binding to maintain structural integrity of the enzyme. Soybean lines with this mutation had higher CK content and desired yield characteristics. The potential for marker-assisted selection based on the identified natural variation within GmCKX7-1, is discussed in the context of hormonal control that can result in higher soybean yield.

Cytokinin activity during early kernel development corresponds positively with yield potential and later stage ABA accumulation in field-grown wheat (Triticum aestivum L.)

Early cytokinin activity and late abscisic acid dynamics during wheat kernel development correspond to cultivars with higher yield potential. Cytokinins represent prime targets for marker development for wheat breeding programs. Two major phytohormone groups, abscisic acid (ABA) and cytokinins (CKs), are of crucial importance for seed development. Wheat (Triticum aestivum L.) yield is, to a high degree, determined during the milk and dough stages of kernel development. Therefore, understanding the hormonal regulation of these early growth stages is fundamental for crop-improvement programs of this important cereal. Here, we profiled ABA and 25 CK metabolites (including active forms, precursors and inactive conjugates) during kernel development in five field-grown wheat cultivars. The levels of ABA and profiles of CK forms varied greatly among the tested cultivars and kernel stages suggesting that several types of CK metabolites are involved in spatiotemporal regulation of kernel development. The seed yield potential was associated with the elevated levels of active CK levels (tZ, cZ). Interestingly, the increased kernel cZ levels were followed by higher ABA production, suggesting there is an interaction between these two phytohormones. Furthermore, we analyzed the expression patterns of representatives of the four main CK metabolic gene families. The unique transcriptional patterns of the IPT (biosynthesis) and ZOG (reversible inactivation) gene family members (GFMs) in the high and low yield cultivars additionally indicate that there is a significant association between CK metabolism and yield potential in wheat. Based on these results, we suggest that both CK metabolites and their associated genes, can serve as important, early markers of yield performance in modern wheat breeding programs.

Roles of metabolic regulation in developing Quercus variabilis acorns at contrasting geologically-derived phosphorus sites in subtropical China

BACKGROUND

Phosphorus (P) -rich soils develop in phosphorite residing areas while P-deficient soils are ubiquitous in subtropical regions. Little has been reported that how metabolites participate in the seed development and the processes involved in their coping with contrasting-nutrient environments.

RESULTS

Here we quantified the metabolites of Quercus variabilis acorns in the early (July), middle (August), late (September) development stages, and determined element (C, H, O, N, P, K, Ca, Mg, S, Fe, Al, Mn, Na, Zn, and Cu) concentrations of acorns in the late stage, at geologically-derived contrasting-P sites in subtropical China. The primary metabolic pathways included sugar metabolism, the TCA cycle, and amino acid metabolism. Most metabolites (especially C- and N-containing metabolites) increased and then decreased from July to September. Acorns between the two sites were significantly discriminated at the three stages, respectively, by metabolites (predominantly sugars and organic acids). Concentrations of P, orthophosphoric acid and most sugars were higher; erythrose was lower in late-stage acorns at P-rich sites than those at P-deficient sites. No significant differences existed in the size and dry mass of individual acorns between oak populations at the two sites.

CONCLUSIONS

Oak acorns at the two sites formed distinct metabolic phenotypes related to their distinct geologically-derived soil conditions, and the late-stage acorns tended to increase P-use-efficiency in the material synthesis process at P-deficient sites, relative to those at P-rich sites.

Comparative metabolic profiling of cultivated and wild black soybeans reveals distinct metabolic alterations associated with their domestication

Generally, cultivated black soybean (CBS) has been used as a major source of various nutrients for humans and animals. To assess the metabolic alterations induced by domestication in soybean, we performed a comprehensive metabolite profiling of 56 soybean varieties, including 28 CBS and 28 wild black soybean (WBS) varieties. A total of 48 metabolites were characterized, including 45 primary and 3 secondary metabolites, from CBS and WBS. The results of principal component analysis and hierarchical cluster analysis (HCA) revealed significant metabolic differences between CBS and WBS that were closely related to metabolic pathways. The results indicate that flavonoids correlated positively with phenylalanine, a precursor for phenylpropanoid biosynthesis; the contents of flavonoids and phenylpropanoids were higher in WBS. Pathway analysis revealed that CBS contained large amounts of TCA cycle intermediates, amino acids, and fatty acids as a result of increased energy metabolism, amino acid metabolism, and seed filling. The projection to latent structure method, using the partial least squares method, was applied to predict the flavonoid content in soybean seed, which indicated that sucrose, threonic acid, citric acid, and fatty acids are important in predicting the antioxidant content of samples. This work will provide important information for designing new soybean cultivars with enhanced nutritional and agricultural traits.

Exogenous Abscisic Acid Supplementation at Early Stationary Growth Phase Triggers Changes in the Regulation of Fatty Acid Biosynthesis in Chlorella vulgaris UMT-M1

Abscisic acid (ABA) has been known to exist in a microalgal system and serves as one of the chemical stimuli in various biological pathways. Nonetheless, the involvement of ABA in fatty acid biosynthesis, particularly at the transcription level in microalgae is poorly understood. The objective of this study was to determine the effects of exogenous ABA on growth, total oil content, fatty acid composition, and the expression level of beta ketoacyl-ACP synthase I (KAS I) and omega-3 fatty acid desaturase (ω-3 FAD) genes in Chlorella vulgaris UMT-M1. ABA was applied to early stationary C. vulgaris cultures at concentrations of 0, 10, 20, and 80 μM for 48 h. The results showed that ABA significantly increased biomass production and total oil content. The increment of palmitic (C16:0) and stearic (C18:0) acids was coupled by decrement in linoleic (C18:2) and α-linolenic (C18:3n3) acids. Both KAS I and ω-3 FAD gene expression were downregulated, which was negatively correlated to saturated fatty acid (SFAs), but positively correlated to polyunsaturated fatty acid (PUFA) accumulations. Further analysis of both KAS I and ω-3 FAD promoters revealed the presence of multiple ABA-responsive elements (ABREs) in addition to other phytohormone-responsive elements. However, the role of these phytohormone-responsive elements in regulating KAS I and ω-3 FAD gene expression still remains elusive. This revelation might suggest that phytohormone-responsive gene regulation in C. vulgaris and microalgae as a whole might diverge from higher plants which deserve further scientific research to elucidate its functional roles.

Alleviation of cadmium-induced genotoxicity and cytotoxicity by calcium chloride in faba bean (Vicia faba L. var. minor) roots

Alleviation of cadmium-induced root genotoxicity and cytotoxicity by calcium chloride (CaCl₂) in faba bean (Vicia faba L. var. minor) seedlings were studied. Faba bean seeds were treated with H₂O or 2% CaCl₂ for 6 h before germination. Seeds were then exposed to 0 and 50 µM CdCl₂ concentrations for 7 days. Genotoxic damaging effects of Cd was examined through the determination of the mitotic index (MI), chromosomal aberrations (CA) and micronucleus (MN) in the meristem cells of faba bean roots. Similarly, effects of Cd stress on metal accumulation, total membrane lipid contents, total fatty acid composition (TFA), lipid peroxidation as indicated by malondialdehyde production, soluble protein and non-protein thiols (NP-SH) contents, hydrogen peroxide production and the activities of superoxide dismutase (SOD), catalase (CAT) and guaiacol peroxidase (GPX) were evaluated after 7 days of Cd stress in the seedling roots. Cd stress resulted in the reduction of MI, in addition to MN formation and CA induction in the roots of non-primed seeds (treated with H₂O). Moreover, Cd induced lipid peroxidation, H₂O₂ overproduction and loss of membrane lipid amount and soluble protein content, and changes in the TFA composition in roots of faba bean seedlings. SOD activity declined, but CAT and GPX activities increased. However, seed pre-treatment with CaCl₂ attenuated the genotoxic and cytotoxic effects of Cd on Vicia faba roots. The results showed that CaCl₂ induced reduction of Cd accumulation, improved cell membrane stability and increased the antioxidant defence systems, thus reducing and alleviating Cd genotoxicity and oxidative damage.

Comparative transcriptome analysis of flower heterosis in two soybean F1 hybrids by RNA-seq

Heterosis has been widely exploited as an approach to enhance crop traits during breeding. However, its underlying molecular genetic mechanisms remain unclear. Recent advances in RNA sequencing technology (RNA-seq) have provided an opportunity to conduct transcriptome profiling for heterosis studies. We used RNA-seq to analyze the flower transcriptomes of two F1 hybrid soybeans (HYBSOY-1 and HYBSOY-5) and their parents. More than 385 million high-quality reads were generated and aligned against the soybean reference genome. A total of 681 and 899 genes were identified as being differentially expressed between HYBSOY-1 and HYBSOY-5 and their parents, respectively. These differentially expressed genes (DEGs) were categorized into four major expression categories with 12 expression patterns. Furthermore, gene ontology (GO) term analysis showed that the DEGs were enriched in the categories metabolic process and catalytic activity, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that metabolic pathway and biosynthesis of secondary metabolites were enriched in the two F1 hybrids. Comparing the DEGs of the two F1 hybrids by GO term and KEGG pathway analyses identified 26 common DEGs that showed transgressive up-regulation, and which could be considered potential candidate genes for heterosis in soybean F1 hybrids. This identification of an extensive transcriptome dataset gives a comprehensive overview of the flower transcriptomes in two F1 hybrids, and provides useful information for soybean hybrid breeding. These findings lay the foundation for future studies on molecular mechanisms underlying soybean heterosis.

Mutant Transcriptome Sequencing Provides Insights into Pod Development in Peanut (Arachis hypogaea L.)

Pod size is the major yield component and a key target trait that is selected for in peanut breeding. However, although numerous quantitative trait loci (QTLs) for peanut pod size have been described, the molecular mechanisms underlying the development of this characteristic remain elusive. A peanut mutant with a narrower pod was developed in this study using ethyl methanesulfonate (EMS) mutagenesis and designated as the "pod width" mutant line (pw). The fresh pod weight of pw was only about 40% of that seen in the wild-type (WT) Zhonghua16, while the hull and seed filling of the mutant both also developed at earlier stages. Pods from both pw and WT lines were sampled 20, 40, and 60 days after flowering (DAF) and used for RNA-Seq analysis; the results revealed highly differentially expressed lignin metabolic pathway genes at all three stages, but especially at DAF 20 and DAF 40. At the same time, expression of genes related to auxin signal transduction was found to be significantly repressed during the pw early pod developmental stage. A genome-wide comparative analysis of expression profiles revealed 260 differentially expressed genes (DEGs) across all three stages, and two candidate genes, c26901_g1 (CAD) and c37339_g1 (ACS), responsible for pod width were identified by integrating expression patterns and function annotation of the common DEGs within the three stages. Taken together, the information provided in this study illuminates the processes underlying peanut pod development, and will facilitate further identification of causal genes and the development of improved peanut varieties with higher yields.