Using genomics to study legume seed development

Domestication has altered gene expression and secondary metabolites in pea seed coat

The mature seed in legumes consists of an embryo and seed coat. In contrast to knowledge about the embryo, we know relatively little about the seed coat. We analyzed the gene expression during seed development using a panel of cultivated and wild pea genotypes. Gene co-expression analysis identified gene modules related to seed development, dormancy, and domestication. Oxidoreductase genes were found to be important components of developmental and domestication processes. Proteomic and metabolomic analysis revealed that domestication favored proteins involved in photosynthesis and protein metabolism at the expense of seed defense. Seed coats of wild peas were rich in cell wall-bound metabolites and the protective compounds predominated in their seed coats. Altogether, we have shown that domestication altered pea seed development and modified (mostly reduced) the transcripts along with the protein and metabolite composition of the seed coat, especially the content of the compounds involved in defense. We investigated dynamic profiles of selected identified phenolic and flavonoid metabolites across seed development. These compounds usually deteriorated the palatability and processing of the seeds. Our findings further provide resources to study secondary metabolism and strategies for improving the quality of legume seeds which comprise an important part of the human protein diet.

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.

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Ancestral function but divergent epigenetic regulation of HAIKU2 reveals routes of seed developmental evolution

Evolution is driven by various mechanisms. A directional increase in the embryo to endosperm ratio is an evolutionary trend within the angiosperms. The endosperm constitutes a major portion of the seed volume in Poales and some dicots. However, in other dicots such as Arabidopsis and soybean, the endosperm proliferates early, followed by embryo growth to replace the endosperm. The Arabidopsis leucine-rich repeat receptor protein kinase AtHAIKU2 (AtIKU2) is a key regulator of early endosperm proliferation. In this study, we found that IKU2s from Brachypodium, rice, and soybean can complement the abnormal seed developmental phenotype of Atiku2, while AtIKU2 also rescues the defective endosperm proliferation in the Brachypodium BdIKU2 knockout mutant seeds. AtIKU2 and soybean GmIKU2 are actively expressed a few days after fertilization. Thereafter, expression of AtIKU2 is suppressed by the FIS-PRC2 complex-mediated H3K27me3. The soybean GmIKU2 locus is also enriched with H3K27me3 marks. The histone methyltransferase AtMEA is unique to Brassicaceae, but one GmSWN in soybean plays a similar role in seed development as AtMEA. By contrast, the BdIKU2 and rice OsIKU2 loci are continuously expressed and are devoid of H3K27me3 marks. Taken together, these results suggest that IKU2 genes retain an ancestral function, but the duration of their expression that is controlled by PRC2-mediated epigenetic silencing contributes to silenced or persistent endosperm proliferation in different species. Our study reveals an epigenetic mechanism that drives the development of vastly different seed ontogenies.

The Seed Coat’s Impact on Crop Performance in Pea (Pisum sativum L.)

Seed development in angiosperms produces three genetically and developmentally distinct sub-compartments: the embryo, endosperm, and seed coat. The maternally derived seed coat protects the embryo and interacts closely with the external environment especially during germination and seedling establishment. Seed coat is a key contributor to seed composition and an important determinant of nutritional value for humans and livestock. In this review, we examined pea crop productivity through the lens of the seed coat, its contribution to several valued nutritional traits of the pea crop, and its potential as a breeding target. Key discoveries made in advancing the knowledge base for sensing and transmission of external signals, the architecture and chemistry of the pea seed coat, and relevant insights from other important legumes were discussed. Furthermore, for selected seed coat traits, known mechanisms of genetic regulation and efforts to modulate these mechanisms to facilitate composition and productivity improvements in pea were discussed, alongside opportunities to support the continued development and improvement of this underutilized crop. This review describes the most important features of seed coat development in legumes and highlights the key roles played by the seed coat in pea seed development, with a focus on advances made in the genetic and molecular characterization of pea and other legumes and the potential of this key seed tissue for targeted improvement and crop optimization.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Temporal Control of Seed Development in Dicots: Molecular Bases, Ecological Impact and Possible Evolutionary Ramifications

In flowering plants, seeds serve as organs of both propagation and dispersal. The developing seed passes through several consecutive stages, following a conserved general outline. The overall time needed for a seed to develop, however, may vary both within and between plant species, and these temporal developmental properties remain poorly understood. In the present paper, we summarize the existing data for seed development alterations in dicot plants. For genetic mutations, the reported cases were grouped in respect of the key processes distorted in the mutant specimens. Similar phenotypes arising from the environmental influence, either biotic or abiotic, were also considered. Based on these data, we suggest several general trends of timing alterations and how respective mechanisms might add to the ecological plasticity of the families considered. We also propose that the developmental timing alterations may be perceived as an evolutionary substrate for heterochronic events. Given the current lack of plausible models describing timing control in plant seeds, the presented suggestions might provide certain insights for future studies in this field.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development

BACKGROUND

Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement.

RESULTS

In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules.

CONCLUSIONS

Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Comparative transcriptome analysis during seeds development between two soybean cultivars

Soybean is one of the important economic crops, which supplies a great deal of vegetable oil and proteins for human. The content of nutrients in different soybean seeds is different, which is related to the expression of multiple genes, but the mechanisms are complicated and still largely uncertain. In this study, to reveal the possible causes of the nutrients difference in soybeans A7 (containing low oil and high protein) and A35 (containing high oil and low protein), RNA-seq technology was performed to compare and identify the potential differential expressed genes (DEGs) at different seed developmental stages. The results showed that DEGs mainly presented at the early stages of seeds development and more DEGs were up-regulated at the early stage than the late stages. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed that the DEGs have diverged in A7 and A35. In A7, the DEGs were mainly involved in cell cycle and stresses, while in A35 were the fatty acids and sugar metabolism. Specifically, when the DEGs contributing to oil and protein metabolic pathways were analyzed, the differences between A7 and A35 mainly presented in fatty acids metabolism and seeds storage proteins (SSPs) synthesis. Furthermore, the enzymes, fatty acid dehydrogenase 2, 3-ketoacyl-CoA synthase and 9S-lipoxygenase, in the synthesis and elongation pathways of fatty acids, were revealed probably to be involved in the oil content difference between A7 and A35, the SSPs content might be due to the transcription factors: Leafy Cotyledon 2 and Abscisic acid-intensitive 3, while the sugar transporter, SWEET10a, might contribute to both oil and protein content differences. Finally, six DEGs were selected to analyze their expression using qRT-PCR, and the results were consistent with the RNA-seq results. Generally, the study provided a comprehensive and dynamic expression trends for the seed development processes, and uncovered the potential DEGs for the differences of oil in A7 and A35.

Genome-wide profiling of miRNAs during seed development reveals their functional relevance in seed size/weight determination in chickpea

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at transcriptional and post-transcriptional levels. The role of miRNAs in seed development and seed size/weight determination is poorly understood in legumes. In this study, we profiled miRNAs at seven successive stages of seed development in a small-seeded and a large-seeded chickpea cultivar via small RNA sequencing. In total, 113 known and 243 novel miRNAs were identified. Gene ontology analysis revealed the enrichment of seed/reproductive/post-embryonic development and signaling pathways processes among the miRNA target genes. A large fraction of the target genes exhibited antagonistic correlation with miRNA expression. The sets of co-expressed miRNAs showing differential expression between the two cultivars were recognized. Known transcription factor (TF) encoding genes involved in seed size/weight determination, including SPL, GRF, MYB, ARF, HAIKU1, SHB1, KLUH/CYP78A5, and E2Fb along with novel genes were found to be targeted by the predicted miRNAs. Differential expression analysis revealed higher transcript levels of members of SPL and REVOLUTA TF families and lower expression of their corresponding miRNAs in the large-seeded cultivar. At least 19 miRNAs known to be involved in seed development or differentially expressed between small-seeded and large-seeded cultivars at late-embryogenesis and/or mid-maturation stages were located within known quantitative trait loci (QTLs) associated with seed size/weight determination. Moreover, 41 target genes of these miRNAs were also located within these QTLs. Altogether, we revealed important roles of miRNAs in seed development and identified candidate miRNAs and their target genes that have functional relevance in determining seed size/weight in chickpea.

Gene co-expression analysis of tomato seed maturation reveals tissue-specific regulatory networks and hubs associated with the acquisition of desiccation tolerance and seed vigour

BACKGROUND

During maturation seeds acquire several physiological traits to enable them to survive drying and disseminate the species. Few studies have addressed the regulatory networks controlling acquisition of these traits at the tissue level particularly in endospermic seeds such as tomato, which matures in a fully hydrated environment and does not undergo maturation drying. Using temporal RNA-seq analyses of the different seed tissues during maturation, gene network and trait-based correlations were used to explore the transcriptome signatures associated with desiccation tolerance, longevity, germination under water stress and dormancy.

RESULTS

During maturation, 15,173 differentially expressed genes were detected, forming a gene network representing 21 expression modules, with 3 being specific to seed coat and embryo and 5 to the endosperm. A gene-trait significance measure identified a common gene module between endosperm and embryo associated with desiccation tolerance and conserved with non-endospermic seeds. In addition to genes involved in protection such LEA and HSP and ABA response, the module included antioxidant and repair genes. Dormancy was released concomitantly with the increase in longevity throughout fruit ripening until 14 days after the red fruit stage. This was paralleled by an increase in SlDOG1-2 and PROCERA transcripts. The progressive increase in seed vigour was captured by three gene modules, one in common between embryo and endosperm and two tissue-specific. The common module was enriched with genes associated with mRNA processing in chloroplast and mitochondria (including penta- and tetratricopeptide repeat-containing proteins) and post-transcriptional regulation, as well several flowering genes. The embryo-specific module contained homologues of ABI4 and CHOTTO1 as hub genes associated with seed vigour, whereas the endosperm-specific module revealed a diverse set of processes that were related to genome stability, defence against pathogens and ABA/GA response genes.

CONCLUSION

The spatio-temporal co-expression atlas of tomato seed maturation will serve as a valuable resource for the in-depth understanding of the dynamics of gene expression associated with the acquisition of seed vigour at the tissue level.

Genome-wide identification and expression profiling of glutathione S-transferase family under multiple abiotic and biotic stresses in Medicago truncatula L

Glutathione transferases (GSTs) constitute an ancient, ubiquitous, multi-functional antioxidant enzyme superfamily that has great importance on cellular detoxification against abiotic and biotic stresses as well as plant development and growth. The present study aimed to a comprehensive genome-wide identification and functional characterization of GST family in one of the economically important legume plants-Medicago truncatula. Here, we have identified a total of ninety-two putative MtGST genes that code for 120 proteins. All these members were classified into twelve classes based on their phylogenetic relationship and the presence of structural conserved domain/motif. Among them, 7 MtGST gene pairs were identified to have segmental duplication. Expression profiling of MtGST transcripts revealed their high level of organ/tissue-specific expression in most of the developmental stages and anatomical tissues. The transcripts of MtGSTU5, MtGSTU8, MtGSTU17, MtGSTU46, and MtGSTU47 showed significant up-regulation in response to various abiotic and biotic stresses. Moreover, transcripts of MtGSTU8, MtGSTU14, MtGSTU28, MtGSTU30, MtGSTU34, MtGSTU46 and MtGSTF8 were found to be highly upregulated in response to drought treatment for 24h and 48h. Among the highly stress-responsive MtGST members, MtGSTU17 showed strong affinity towards its conventional substrates reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB) with the lowest binding energy of-5.7 kcal/mol and -6.5 kcal/mol, respectively. Furthermore, the substrate-binding site residues of MtGSTU17 were found to be highly conserved. These findings will facilitate the further functional and evolutionary characterization of GST genes in Medicago.

Comparative analysis of embryo proper and suspensor transcriptomes in plant embryos with different morphologies

An important question is what genes govern the differentiation of plant embryos into suspensor and embryo proper regions following fertilization and division of the zygote. We compared embryo proper and suspensor transcriptomes of four plants that vary in embryo morphology within the suspensor region. We determined that genes encoding enzymes in several metabolic pathways leading to the formation of hormones, such as gibberellic acid, and other metabolites are up-regulated in giant scarlet runner bean and common bean suspensors. Genes involved in transport and Golgi body organization are up-regulated within the suspensors of these plants as well, strengthening the view that giant specialized suspensors serve as a hormone factory and a conduit for transferring substances to the developing embryo proper. By contrast, genes controlling transcriptional regulation, development, and cell division are up-regulated primarily within the embryo proper. Transcriptomes from less specialized soybean and Arabidopsis suspensors demonstrated that fewer genes encoding metabolic enzymes and hormones are up-regulated. Genes active in the embryo proper, however, are functionally similar to those active in scarlet runner bean and common bean embryo proper regions. We uncovered a set of suspensor- and embryo proper-specific transcription factors (TFs) that are shared by all embryos irrespective of morphology, suggesting that they are involved in early differentiation processes common to all plants. Chromatin immunoprecipitation sequencing (ChIP-Seq) experiments with scarlet runner bean and soybean WOX9, an up-regulated suspensor TF, gained entry into a regulatory network important for suspensor development irrespective of morphology.

In silico Analysis of Acyl-CoA-Binding Protein Expression in Soybean

Plant acyl-CoA-binding proteins (ACBPs) form a highly conserved protein family that binds to acyl-CoA esters as well as other lipid and protein interactors to function in developmental and stress responses. This protein family had been extensively studied in non-leguminous species such as Arabidopsis thaliana (thale cress), Oryza sativa (rice), and Brassica napus (oilseed rape). However, the characterization of soybean (Glycine max) ACBPs, designated GmACBPs, has remained unreported although this legume is a globally important crop cultivated for its high oil and protein content, and plays a significant role in the food and chemical industries. In this study, 11 members of the GmACBP family from four classes, comprising Class I (small), Class II (ankyrin repeats), Class III (large), and Class IV (kelch motif), were identified. For each class, more than one copy occurred and their domain architecture including the acyl-CoA-binding domain was compared with Arabidopsis and rice. The expression profile, tertiary structure and subcellular localization of each GmACBP were predicted, and the similarities and differences between GmACBPs and other plant ACBPs were deduced. A potential role for some Class III GmACBPs in nodulation, not previously encountered in non-leguminous ACBPs, has emerged. Interestingly, the sole member of Class III ACBP in each of non-leguminous Arabidopsis and rice had been previously identified in plant-pathogen interactions. As plant ACBPs are known to play important roles in development and responses to abiotic and biotic stresses, the in silico expression profiles on GmACBPs, gathered from data mining of RNA-sequencing and microarray analyses, will lay the foundation for future studies in their applications in biotechnology.

Analysis of Spatio-Temporal Transcriptome Profiles of Soybean (Glycine max) Tissues during Early Seed Development

Soybean (Glycine max) is an important crop providing oil and protein for both human and animal consumption. Knowing which biological processes take place in specific tissues in a temporal manner will enable directed breeding or synthetic approaches to improve seed quantity and quality. We analyzed a genome-wide transcriptome dataset from embryo, endosperm, endothelium, epidermis, hilum, outer and inner integument and suspensor at the global, heart and cotyledon stages of soybean seed development. The tissue specificity of gene expression was greater than stage specificity, and only three genes were differentially expressed in all seed tissues. Tissues had both unique and shared enriched functional categories of tissue-specifically expressed genes associated with them. Strong spatio-temporal correlation in gene expression was identified using weighted gene co-expression network analysis, with the most co-expression occurring in one seed tissue. Transcription factors with distinct spatiotemporal gene expression programs in each seed tissue were identified as candidate regulators of expression within those tissues. Gene ontology (GO) enrichment of orthogroup clusters revealed the conserved functions and unique roles of orthogroups with similar and contrasting expression patterns in transcript abundance between soybean and Arabidopsis during embryo proper and endosperm development. Key regulators in each seed tissue and hub genes connecting those networks were characterized by constructing gene regulatory networks. Our findings provide an important resource for describing the structure and function of individual soybean seed compartments during early seed development.

Global transcriptome analysis of subterranean pod and seed in peanut (Arachis hypogaea L.) unravels the complexity of fruit development under dark condition

Peanut pods develop underground, which is the most salient characteristic in peanut. However, its developmental transcriptome remains largely unknown. In the present study, we sequenced over one billion transcripts to explore the developmental transcriptome of peanut pod using Illumina sequencing. Moreover, we identified and quantified the abundances of 165,689 transcripts in seed and shell tissues along with a pod developmental gradient. The dynamic changes of differentially expressed transcripts (DETs) were described in seed and shell. Additionally, we found that photosynthetic genes were not only pronouncedly enriched in aerial pod, but also played roles in developing pod under dark condition. Genes functioning in photomorphogenesis showed distinct expression profiles along subterranean pod development. Clustering analysis unraveled a dynamic transcriptome, in which transcripts for DNA synthesis and cell division during pod expansion were transitioning to transcripts for cell expansion and storage activity during seed filling. Collectively, our study formed a transcriptional baseline for peanut fruit development under dark condition.

A D-lactate dehydrogenase from rice is involved in conferring tolerance to multiple abiotic stresses by maintaining cellular homeostasis

D-lactate dehydrogenase (D-LDH) converts D-lactate (the end product of glyoxalase system) to pyruvate and thereby completes the detoxification process of methylglyoxal. D-LDH detoxifies and diverts the stress induced toxic metabolites, MG and D-lactate, towards energy production and thus, protects the cell from their deteriorating effects. In this study, a D-LDH enzyme from rice (OsD-LDH2, encoded by Os07g08950.1) was characterized for its role in abiotic stress tolerance. For this, a combination of in silico, molecular, genetic and biochemical approaches was used. The kinetic analysis revealed OsD-LDH2 to be the most efficient D-LDH enzyme in comparison to D-LDHs from other plant species. Heterologous overexpression of OsD-LDH2 provides tolerance against multiple abiotic stresses in E. coli, yeast and plant system. The analysis of D-LDH mutant and OsD-LDH2 overexpressing transgenic plants uncovered the crucial role of D-LDH in mitigation of abiotic stresses. OsD-LDH2 overexpressing plants maintained lower level of ROS and other toxic metabolites along with better functioning of antioxidant system. This is the first report on correlation of D-LDH with multiple abiotic stress tolerance. Overall, OsD-LDH2 emerged as a promising candidate which can open a new direction for engineering stress tolerant crop varieties by maintaining their growth and yield in unfavorable conditions.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Using Giant Scarlet Runner Bean (Phaseolus coccineus) Embryos to Dissect the Early Events in Plant Embryogenesis

The giant embryo of the scarlet runner bean (Phaseolus coccineus) has been used historically to investigate the molecular and developmental processes that control the early events of plant embryo development. In more recent years, our laboratory has been using scarlet runner bean embryos to uncover the genes and regulatory events that control embryo proper and suspensor region differentiation shortly after fertilization. In this chapter we describe methods that we have developed to isolate scarlet runner bean embryos at the globular stage of development, and capture embryo proper and suspensor regions by either hand dissection or laser capture microdissection (LCM) for use in downstream genomic analysis. These methods are also applicable for use in investigating the early events of common bean (Phaseolus vulgaris) embryo development, a close relative of scarlet runner bean, which also has a giant embryo in addition to a sequenced genome.

Microarray-Based Quality Assessment as a Supporting Criterion for de novo Transcriptome Assembly Selection

RNA-Sequencing and de novo assembly have enabled the analysis of species with non-available reference transcriptomes, although intrinsic features (biological and technical) induce errors in the reconstruction. A strategy to resolve these errors consists of varying assembling process parameters to generate multiple reconstructions. However, the best assembly selection remains a challenge. Quantitative metrics for quality assessment have been inconsistent when compared with pertinent references. In this paper, a criterion for supporting assembly selection based on mapping DNA microarray hybridized probes to assembly sets is proposed. Mouse and fruit fly RNA-Seq datasets were assembled with standard de novo procedures. Quality assessment was estimated using quantitative metrics and the proposed criterion. The assembly that best mapped to the available reference transcriptomes of these model species provided the highest quality assembly. The hybridized probes identified the best assemblies, whereas quantitative metrics remained inconsistent. For example, subtle probe mapping difference of 0.25 percent, but statistically significant (ANOVA, p < 0.05), enabled the assembly selection that led to identify 3,719 more contigs and led to 1,049 further mapped contigs to the mouse reference transcriptome. The microarray data availability for non-model species makes the proposed criterion suitable for quality assessment of multiple de novo assembly strategies.

Identification and characterization of a major QTL underlying soybean isoflavone malonylglycitin content

Isoflavones in soybean seeds are responsible for plant-microbe interactions and defend against pathogens, and are also beneficial to human health. We used two biparental populations and mini core collection of soybean germplasm to identify and validate QTLs underlying the content of isoflavone components. We identified a major QTL, qMGly_11, which regulates the content of malonylglycitin, on chromosome Gm11, in populations bred from parents with high, low, and null glycitein contents. qMGly_11 explained 44.5% of phenotypic variance in a population derived from a cross between 'Aokimame' (high) and 'Fukuyutaka' (low) and 79.9% of that in a population between 'Kumaji-1' (null) and 'Fukuyutaka' (low). The effect was observed only in the hypocotyl. We further confirmed the effect of qMGly_11 in a mini-core collection, where it explained 57.1% of the genetic diversity of glycitin production and 56.5% of malonylglycitin production. qMGly_11 increased the contents of glycitin and malonylglycitin at the expense of daidzin and malonyldaidzin in all analyzed populations. We discuss the gene responsible for this QTL and the availability of the null allele for metabolic engineering of soybean seed isoflavones.

Mapping developmental QTL for plant height in soybean [Glycine max (L.) Merr.] using a four-way recombinant inbred line population

Plant height (PH) is an important trait in soybean, as taller plants may have higher yields but may also be at risk for lodging. Many genes act jointly to influence PH throughout development. To map the quantitative trait loci (QTL) controlling PH, we used the unconditional variable method (UVM) and conditional variable method (CVM) to analyze PH data for a four-way recombinant inbred line (FW-RIL) population derived from the cross of (Kenfeng14 × Kenfeng15) × (Heinong48 × Kenfeng19). We identified 7, 8, 16, 19, 15, 27, 17, 27, 22, and 24 QTL associated with PH at 10 developmental stages, respectively. These QTL mapped to 95 genomic regions. Among these QTL, 9 were detected using UVM and CVM, and 89 and 66 were only detected by UVM or CVM, respectively. In total, 36 QTL controlling PH were detected at multiple developmental stages and these made unequal contributions to genetic variation throughout development. Among 19 novel regions discovered in our study, 7 could explain over 10% of the phenotypic variation and contained only one single QTL. The unconditional and conditional QTL detected here could be used in molecular design breeding across the whole developmental procedure.

The Soybean Sugar Transporter GmSWEET15 Mediates Sucrose Export from Endosperm to Early Embryo

Soybean (Glycine max) seed is primarily composed of a mature embryo that provides a major source of protein and oil for humans and other animals. Early in development, the tiny embryos grow rapidly and acquire large quantities of sugars from the liquid endosperm of developing seeds. An insufficient supply of nutrients from the endosperm to the embryo results in severe seed abortion and yield reduction. Hence, an understanding of the molecular basis and regulation of assimilate partitioning involved in early embryo development is important for improving soybean seed yield and quality. Here, we used expression profiling analysis to show that two paralogous sugar transporter genes from the SWEET (Sugars Will Eventually be Exported Transporter) family, GmSWEET15a and GmSWEET15b, were highly expressed in developing soybean seeds. In situ hybridization and quantitative real-time PCR showed that both genes were mainly expressed in the endosperm at the cotyledon stage. GmSWEET15b showed both efflux and influx activities for sucrose in Xenopus oocytes. In Arabidopsis (Arabidopsis thaliana), knockout of three AtSWEET alleles is required to see a defective, but not lethal, embryo phenotype, whereas knockout of both GmSWEET15 genes in soybean caused retarded embryo development and endosperm persistence, resulting in severe seed abortion. In addition, the embryo sugar content of the soybean knockout mutants was greatly reduced. These results demonstrate that the plasma membrane sugar transporter, GmSWEET15, is essential for embryo development in soybean by mediating Suc export from the endosperm to the embryo early in seed development.

Central role of the LEAFY COTYLEDON1 transcription factor in seed development

Seed development is a complex period of the flowering plant life cycle. After fertilization, the three main regions of the seed, embryo, endosperm and seed coat, undergo a series of developmental processes that result in the production of a mature seed that is developmentally arrested, desiccated, and metabolically quiescent. These processes are highly coordinated, both temporally and spatially, to ensure the proper growth and development of the seed. The transcription factor, LEAFY COTYLEDON1 (LEC1), is a central regulator that controls several aspects of embryo and endosperm development, including embryo morphogenesis, photosynthesis, and storage reserve accumulation. Thus, LEC1 regulates distinct sets of genes at different stages of seed development. Despite its critical importance for seed development, an understanding of the mechanisms underlying LEC1's multifunctionality is only beginning to be obtained. Recent studies describe the roles of specific transcription factors and the hormones, gibberellic acid and abscisic acid, in controlling the activity and transcriptional specificity of LEC1 across seed development. Moreover, studies indicate that LEC1 acts as a pioneer transcription factor to promote epigenetic reprogramming during embryogenesis. In this review, we discuss the mechanisms that enable LEC1 to serve as a central regulator of seed development.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean’s agronomic traits.

Genetic Mapping of Milling Quality Traits in Lentil (Lens culinaris Medik.)

Milling qualities are key traits for the red lentil ( Medik.) industry as price is largely determined by recovery yield. milling involves removal of the seed coat and splitting of the cotyledon to produce either splits or footballs (cotyledons still attached). The objectives of the study were to determine the heritability of the milling traits dehulling efficiency (DE), milling recovery (MR), and football recovery (FR) and to identify the genomic regions controlling them. We used a lentil recombinant inbred population from the cross 'CDC Robin' × '946a-46', which have contrasting seed characteristics. The mapping population consists of 127 F-derived lentil recombinant inbred lines that were phenotyped for milling quality parameters from four site-years in Saskatchewan, Canada. A total of 534 single nucleotide polymorphism markers, seven simple sequence repeat markers, and four morphological markers were used for quantitative trait locus (QTL) mapping. The broad-sense heritability was moderate for DE and MR and relatively low for FR. Milling quality traits were significantly correlated with seed shape (seed diameter and seed plumpness). Multiple QTLs for milling traits were detected in six of seven linkage groups (LGs). The most stable QTLs governing DE and MR were clustered on LGs 1, 2, 3, and 7, whereas FR QTLs were clustered on LGs 4, 5, 6, and 7. The molecular markers identified for these traits could be used for improving milling quality in lentil breeding programs.

A shared cis-regulatory module activates transcription in the suspensor of plant embryos

The mechanisms controlling the transcription of gene sets in specific regions of a plant embryo shortly after fertilization remain unknown. Previously, we showed that G564 mRNA, encoding a protein of unknown function, accumulates to high levels in the giant suspensor of both Scarlet Runner Bean (SRB) and Common Bean embryos, and a cis-regulatory module containing three unique DNA sequences, designated as the 10-bp, Region 2, and Fifth motifs, is required for G564 suspensor-specific transcription [Henry KF, et al. (2015) Plant Mol Biol 88:207-217; Kawashima T, et al. (2009) Proc Natl Acad Sci USA 106:3627-3632]. We tested the hypothesis that these motifs are also required for transcription of the SRB GA 20-oxidase gene, which encodes a gibberellic acid hormone biosynthesis enzyme and is coexpressed with G564 at a high level in giant bean suspensors. We used deletion and gain-of-function experiments in transgenic tobacco embryos to show that two GA 20-oxidase DNA regions are required for suspensor-specific transcription, one in the 5' UTR (+119 to +205) and another in the 5' upstream region (-341 to -316). Mutagenesis of sequences in these two regions determined that the cis-regulatory motifs required for G564 suspensor transcription are also required for GA 20-oxidase transcription within the suspensor, although the motif arrangement differs. Our results demonstrate the flexibility of motif positioning within a cis-regulatory module that activates gene transcription within giant bean suspensors and suggest that G564 and GA 20-oxidase comprise part of a suspensor gene regulatory network.

The stage of seed development influences iron bioavailability in pea (Pisum sativum L.)

Pea seeds are widely consumed in their immature form, known as garden peas and petit pois, mostly after preservation by freezing or canning. Mature dry peas are rich in iron in the form of ferritin, but little is known about the content, form or bioavailability of iron in immature stages of seed development. Using specific antibodies and in-gel iron staining, we show that ferritin loaded with iron accumulated gradually during seed development. Immunolocalization and high-resolution secondary ion mass spectrometry (NanoSIMS) revealed that iron-loaded ferritin was located at the surface of starch-containing plastids. Standard cooking procedures destabilized monomeric ferritin and the iron-loaded form. Iron uptake studies using Caco-2 cells showed that the iron in microwaved immature peas was more bioavailable than in boiled mature peas, despite similar levels of soluble iron in the digestates. By manipulating the levels of phytic acid in the digestates we demonstrate that phytic acid is the main inhibitor of iron uptake from mature peas in vitro. Taken together, our data show that immature peas and mature dry peas contain similar levels of ferritin-iron, which is destabilized during cooking. However, iron from immature peas is more bioavailable because of lower phytic acid levels compared to mature peas.

Ribosome profiling reveals changes in translational status of soybean transcripts during immature cotyledon development

To understand translational capacity on a genome-wide scale across three developmental stages of immature soybean seed cotyledons, ribosome profiling was performed in combination with RNA sequencing and cluster analysis. Transcripts representing 216 unique genes demonstrated a higher level of translational activity in at least one stage by exhibiting higher translational efficiencies (TEs) in which there were relatively more ribosome footprint sequence reads mapping to the transcript than were present in the control total RNA sample. The majority of these transcripts were more translationally active at the early stage of seed development and included 12 unique serine or cysteine proteases and 16 2S albumin and low molecular weight cysteine-rich proteins that may serve as substrates for turnover and mobilization early in seed development. It would appear that the serine proteases and 2S albumins play a vital role in the early stages. In contrast, our investigation of profiles of 19 genes encoding high abundance seed storage proteins, such as glycinins, beta-conglycinins, lectin, and Kunitz trypsin inhibitors, showed that they all had similar patterns in which the TE values started at low levels and increased approximately 2 to 6-fold during development. The highest levels of these seed protein transcripts were found at the mid-developmental stage, whereas the highest ribosome footprint levels of only up to 1.6 TE were found at the late developmental stage. These experimental findings suggest that the major seed storage protein coding genes are primarily regulated at the transcriptional level during normal soybean cotyledon development. Finally, our analyses also identified a total of 370 unique gene models that showed very low TE values including over 48 genes encoding ribosomal family proteins and 95 gene models that are related to energy and photosynthetic functions, many of which have homology to the chloroplast genome. Additionally, we showed that genes of the chloroplast were relatively translationally inactive during seed development.

Single Nucleotide Polymorphism Markers Associated with Seed Quality Characteristics of Cultivated Lentil

The dimensions of lentil ( Medik.) seeds are important quality parameters that are major determinants of market preference, cooking time, and post-harvest milling quality. Knowledge of the genetic control of traits related to seed dimensions would be useful for crop improvement. The principal aim of this study was to identify single nucleotide polymorphism (SNP) markers linked to genes that control seed diameter, seed thickness, and seed plumpness. Association mapping analysis with SNP markers was used to study the seed dimensions of 138 diverse cultivated lentil accessions grown at two locations in Saskatchewan, Canada, in 2011 and 2012. Six marker-trait associations were shown to be significant for the studied seed dimension characteristics. Two SNP markers closely associated with seed diameter across locations and years identified in previous work were validated in this study. Three additional marker-seed thickness associations were identified. Using the association mapping strategy, we confirmed the presence of two genomic regions controlling seed diameter and plumpness. This information can be used worldwide as a resource for lentil seed quality improvement programs.

Physiological and transcriptomic responses in the seed coat of field-grown soybean (Glycine max L. Merr.) to abiotic stress

BACKGROUND

Understanding how intensification of abiotic stress due to global climate change affects crop yields is important for continued agricultural productivity. Coupling genomic technologies with physiological crop responses in a dynamic field environment is an effective approach to dissect the mechanisms underpinning crop responses to abiotic stress. Soybean (Glycine max L. Merr. cv. Pioneer 93B15) was grown in natural production environments with projected changes to environmental conditions predicted for the end of the century, including decreased precipitation, increased tropospheric ozone concentrations ([O3]), or increased temperature.

RESULTS

All three environmental stresses significantly decreased leaf-level photosynthesis and stomatal conductance, leading to significant losses in seed yield. This was driven by a significant decrease in the number of pods per node for all abiotic stress treatments. To understand the underlying transcriptomic response involved in the yield response to environmental stress, RNA-Sequencing analysis was performed on the soybean seed coat, a tissue that plays an essential role in regulating carbon and nitrogen transport to developing seeds. Gene expression analysis revealed 49, 148 and 1,576 differentially expressed genes in the soybean seed coat in response to drought, elevated [O3] and elevated temperature, respectively.

CONCLUSIONS

Elevated [O3] and drought did not elicit substantive transcriptional changes in the soybean seed coat. However, this may be due to the timing of sampling and does not preclude impacts of those stresses on different tissues or different stages in seed coat development. Expression of genes involved in DNA replication and metabolic processes were enriched in the seed coat under high temperate stress, suggesting that the timing of events that are important for cell division and proper seed development were altered in a stressful growth environment.

Inference of Transcription Regulatory Network in Low Phytic Acid Soybean Seeds

A dominant loss of function mutation in myo-inositol phosphate synthase (MIPS) gene and recessive loss of function mutations in two multidrug resistant protein type-ABC transporter genes not only reduce the seed phytic acid levels in soybean, but also affect the pathways associated with seed development, ultimately resulting in low emergence. To understand the regulatory mechanisms and identify key genes that intervene in the seed development process in low phytic acid crops, we performed computational inference of gene regulatory networks in low and normal phytic acid soybeans using a time course transcriptomic data and multiple network inference algorithms. We identified a set of putative candidate transcription factors and their regulatory interactions with genes that have functions in myo-inositol biosynthesis, auxin-ABA signaling, and seed dormancy. We evaluated the performance of our unsupervised network inference method by comparing the predicted regulatory network with published regulatory interactions in Arabidopsis. Some contrasting regulatory interactions were observed in low phytic acid mutants compared to non-mutant lines. These findings provide important hypotheses on expression regulation of myo-inositol metabolism and phytohormone signaling in developing low phytic acid soybeans. The computational pipeline used for unsupervised network learning in this study is provided as open source software and is freely available at https://lilabatvt.github.io/LPANetwork/.

Similarity between soybean and Arabidopsis seed methylomes and loss of non-CG methylation does not affect seed development

We profiled soybean and Arabidopsis methylomes from the globular stage through dormancy and germination to understand the role of methylation in seed formation. CHH methylation increases significantly during development throughout the entire seed, targets primarily transposable elements (TEs), is maintained during endoreduplication, and drops precipitously within the germinating seedling. By contrast, no significant global changes in CG- and CHG-context methylation occur during the same developmental period. An Arabidopsis ddcc mutant lacking CHH and CHG methylation does not affect seed development, germination, or major patterns of gene expression, implying that CHH and CHG methylation does not play a significant role in seed development or in regulating seed gene activity. By contrast, over 100 TEs are transcriptionally de-repressed in ddcc seeds, suggesting that the increase in CHH-context methylation may be a failsafe mechanism to reinforce transposon silencing. Many genes encoding important classes of seed proteins, such as storage proteins, oil biosynthesis enzymes, and transcription factors, reside in genomic regions devoid of methylation at any stage of seed development. Many other genes in these classes have similar methylation patterns, whether the genes are active or repressed. Our results suggest that methylation does not play a significant role in regulating large numbers of genes important for programming seed development in both soybean and Arabidopsis. We conclude that understanding the mechanisms controlling seed development will require determining how cis-regulatory elements and their cognate transcription factors are organized in genetic regulatory networks.

Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development

Seed development is programmed by expression of many genes in plants. Seed maturation is an important developmental process to soybean seed quality and yield. DNA methylation is a major epigenetic modification regulating gene expression. However, little is known about the dynamic nature of DNA methylation and its effects on gene expression during plant development. Through whole-genome bisulfite sequencing, we showed that DNA methylation went through dynamic changes during seed maturation. An average of 66% CG, 45% CHG and 9% CHH contexts was methylated in cotyledons. CHH methylation levels in cotyledons changed greatly from 6% at the early stage to 11% at the late stage. Transcribed genes were approximately two-fold more likely to be differentially methylated than non-transcribed genes. We identified 40, 66 and 2136 genes containing differentially methylated regions (DMRs) with negative correlation between their expression and methylation in the CG, CHG and CHH contexts, respectively. The majority of the DMR genes in the CHH context were transcriptionally down-regulated as seeds mature: 99% of them during early maturation were down-regulated, and preferentially associated with DNA replication and cell division. The results provide novel insights into the dynamic nature of DNA methylation and its relationship with gene regulation in seed development.

Differential expression of a WRKY gene between wild and cultivated soybeans correlates to seed size

Soybean (Glycine max) probably originated from the wild soybean (Glycine soja). Glycine max has a significantly larger seed size, but the underlying genomic changes are largely unknown. Candidate regulatory genes were preliminarily proposed by data co-localizing RNA sequencing with the quantitative loci (QTLs) for seed size. The soybean gene locus SoyWRKY15a and its orthologous genes from G. max (GmWRKY15a) and G. soja (GsWRKY15a) were analyzed in detail. The coding sequences were nearly identical between the two orthologs, but GmWRKY15a was significantly more highly expressed than GsWRKY15a. Four haplotypes (H1-H4) were found and they varied in the size of a CT-core microsatellite locus in the 5'-untranslated region of this gene. H1 (with six CT-repeats) was the only allelic version found in G. max, while H3 (with five CT-repeats) was the dominant G. soja allele. Differential expression of this gene in soybean pods was correlated with CT-repeat variation, and manipulation of the CT copy number altered the reporter gene expression, suggesting a regulatory role for the simple sequence repeats. Seed weight of wild soybeans harboring H1 was significantly greater than that of soybeans having haplotypes H2, H3, or H4, and seed weight was correlated with gene expression, suggesting the influence of GsWRKY15a in controlling seed size. However, the seed size might be refractory to increased SoyWRKY15a expression in cultivated soybeans. The evolutionary significance of SoyWRKY15a variation in soybean seed domestication is discussed.

Integration of omics approaches to understand oil/protein content during seed development in oilseed crops

Oilseed crops, especially soybean (Glycine max) and canola/rapeseed (Brassica napus), produce seeds that are rich in both proteins and oils and that are major sources of energy and nutrition worldwide. Most of the nutritional content in the seed is accumulated in the embryo during the seed filling stages of seed development. Understanding the metabolic pathways that are active during seed filling and how they are regulated are essential prerequisites to crop improvement. In this review, we summarize various omics studies of soybean and canola/rapeseed during seed filling, with emphasis on oil and protein traits, to gain a systems-level understanding of seed development. Currently, most (80-85%) of the soybean and rapeseed reference genomes have been sequenced (950 and 850 megabases, respectively). Parallel to these efforts, extensive omics datasets from different seed filling stages have become available. Transcriptome and proteome studies have detected preponderance of starch metabolism and glycolysis enzymes to be the possible cause of higher oil in B. napus compared to other crops. Small RNAome studies performed during the seed filling stages have revealed miRNA-mediated regulation of transcription factors, with the suggestion that this interaction could be responsible for transitioning the seeds from embryogenesis to maturation. In addition, progress made in dissecting the regulation of de novo fatty acid synthesis and protein storage pathways is described. Advances in high-throughput omics and comprehensive tissue-specific analyses make this an exciting time to attempt knowledge-driven investigation of complex regulatory pathways.

Genome-Wide Analysis of Soybean LATERAL ORGAN BOUNDARIES Domain-Containing Genes: A Functional Investigation of GmLBD12

Plant-specific () genes play critical roles in various plant growth and development processes. However, the number and characteristics of genes in soybean [ (L.) Merr.] remain unknown. Here, we identified 90 homologous genes in the soybean genome that phylogenetically clustered into two classes (I and II). The majority of the genes were evenly distributed across all 20 soybean chromosomes, and 77 (81.11%) of them were detected in segmental duplicated regions. Furthermore, the exon-intron organization and motif composition for each were analyzed. A close phylogenetic relationship was identified between the soybean genes and 41 previously reported genes of different plants in the same group, providing insights into their putative functions. Expression analysis indicated that more than half of the genes were expressed, with the two gene classes showing differential tissue expression characteristics; in addition, they were differentially induced by biotic and abiotic stresses. To further explore the functions of genes in soybean, was selected for functional characterization. GmLBD12 was mainly localized to the nucleus and showed high expression in root and seed tissues. Overexpressing in (L.) Heynh resulted in increases in lateral root (LR) number and plant height. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis demonstrated that was induced by drought, salt, cold, indole acetic acid (IAA), abscisic acid (ABA), and salicylic acid SA treatments. This study provides the first comprehensive analysis of the soybean gene family and a valuable foundation for future functional studies of genes.

A Vernonia Diacylglycerol Acyltransferase Can Increase Renewable Oil Production

Increasing the production of plant oils such as soybean oil as a renewable resource for food and fuel is valuable. Successful breeding for higher oil levels in soybean, however, usually results in reduced protein, a second valuable seed component. This study shows that by manipulating a highly active acyl-CoA:diacylglycerol acyltransferase (DGAT) the hydrocarbon flux to oil in oilseeds can be increased without reducing the protein component. Compared to other plant DGATs, a DGAT from Vernonia galamensis (VgDGAT1A) produces much higher oil synthesis and accumulation activity in yeast, insect cells, and soybean. Soybean lines expressing VgDGAT1A show a 4% increase in oil content without reductions in seed protein contents or yield per unit land area. Incorporation of this trait into 50% of soybeans worldwide could result in an increase of 850 million kg oil/year without new land use or inputs and be worth ∼U.S.$1 billion/year at 2012 production and market prices.

Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes

Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Towards a better monitoring of seed ageing under ex situ seed conservation

Long-term conservation of 7.4 million ex situ seed accessions held in agricultural genebanks and botanic gardens worldwide is a challenging mission for human food security and ecosystem services. Recent advances in seed biology and genomics may have opened new opportunities for effective management of seed germplasm under long-term storage. Here, we review the current development of tools for assessing seed ageing and research advances in seed biology and genomics, with a focus on exploring their potential as better tools for monitoring of seed ageing. Seed ageing is found to be associated with the changes reflected in reactive oxygen species and mitochondria-triggered programmed cell deaths, expression of antioxidative genes and DNA and protein repair genes, chromosome telomere lengths, epigenetic regulation of related genes (microRNA and methylation) and altered organelle and nuclear genomes. Among these changes, the signals from mitochondrial and nuclear genomes may show the most promise for use in the development of tools to predict seed ageing. Non-destructive and non-invasive analyses of stored seeds through calorimetry or imaging techniques are also promising. It is clear that research into developing advanced tools for monitoring seed ageing to supplement traditional germination tests will be fruitful for effective conservation of ex situ seed germplasm.

Down-Regulating the Expression of 53 Soybean Transcription Factor Genes Uncovers a Role for SPEECHLESS in Initiating Stomatal Cell Lineages during Embryo Development

We used an RNA interference screen to assay the function of 53 transcription factor messenger RNAs (mRNAs) that accumulate specifically within soybean (Glycine max) seed regions, subregions, and tissues during development. We show that basic helix-loop-helix (bHLH) transcription factor genes represented by Glyma04g41710 and its paralogs are required for the formation of stoma in leaves and stomatal precursor complexes in mature embryo cotyledons. Phylogenetic analysis indicates that these bHLH transcription factor genes are orthologous to Arabidopsis (Arabidopsis thaliana) SPEECHLESS (SPCH) that initiate asymmetric cell divisions in the leaf protoderm layer and establish stomatal cell lineages. Soybean SPCH (GmSPCH) mRNAs accumulate primarily in embryo, seedling, and leaf epidermal layers. Expression of Glyma04g41710 under the control of the SPCH promoter rescues the Arabidopsis spch mutant, indicating that Glyma04g41710 is a functional ortholog of SPCH. Developing soybean embryos do not form mature stoma, and stomatal differentiation is arrested at the guard mother cell stage. We analyzed the accumulation of GmSPCH mRNAs during soybean seed development and mRNAs orthologous to MUTE, FAMA, and inducer of C-repeat/dehydration responsive element-binding factor expression1/scream2 that are required for stoma formation in Arabidopsis. The mRNA accumulation patterns provide a potential explanation for guard mother cell dormancy in soybean embryos. Our results suggest that variation in the timing of bHLH transcription factor gene expression can explain the diversity of stomatal forms observed during plant development.

A cis-regulatory module activating transcription in the suspensor contains five cis-regulatory elements

Little is known about the molecular mechanisms by which the embryo proper and suspensor of plant embryos activate specific gene sets shortly after fertilization. We analyzed the upstream region of the Scarlet Runner Bean (Phaseolus coccineus) G564 gene in order to understand how genes are activated specifically in the suspensor during early embryo development. Previously, we showed that a 54-bp fragment of the G564 upstream region is sufficient for suspensor transcription and contains at least three required cis-regulatory sequences, including the 10-bp motif (5'-GAAAAGCGAA-3'), the 10 bp-like motif (5'-GAAAAACGAA-3'), and Region 2 motif (partial sequence 5'-TTGGT-3'). Here, we use site-directed mutagenesis experiments in transgenic tobacco globular-stage embryos to identify two additional cis-regulatory elements within the 54-bp cis-regulatory module that are required for G564 suspensor transcription: the Fifth motif (5'-GAGTTA-3') and a third 10-bp-related sequence (5'-GAAAACCACA-3'). Further deletion of the 54-bp fragment revealed that a 47-bp fragment containing the five motifs (the 10-bp, 10-bp-like, 10-bp-related, Region 2 and Fifth motifs) is sufficient for suspensor transcription, and represents a cis-regulatory module. A consensus sequence for each type of motif was determined by comparing motif sequences shown to activate suspensor transcription. Phylogenetic analyses suggest that the regulation of G564 is evolutionarily conserved. A homologous cis-regulatory module was found upstream of the G564 ortholog in the Common Bean (Phaseolus vulgaris), indicating that the regulation of G564 is evolutionarily conserved in closely related bean species.

Proteomic insights into synthesis of isoflavonoids in soybean seeds

Soybean seeds are the major human dietary source of isoflavonoids, a class of plant natural products almost entirely exclusive to legumes. Isoflavonoids reduce the risk of a number of chronic human illnesses. Biosynthesis and accumulation of this class of compounds is a multigenic and complex trait, with a great deal of variability among soybean cultivars and with respect to the environment. There is a wealth of genomic, transcriptomic, and metabolomics data regarding isoflavonoid biosynthesis, but the connection between multigene families and their cognate proteins is a missing link that could provide us with a great deal of functional information. The changing proteome of the developing seed can shed light on the correlative increase in isoflavonoids, while the maternal seed coat proteome can provide the link with inherited metabolic and signaling machinery. In this effort, 'seed-filling' proteomics has revealed key secondary metabolite enzymes that quantitatively vary throughout seed development. Seed coat proteomics has revealed the existence of metabolic apparatus specific to isoflavonoid biosynthesis (isoflavonoid reductase) that could potentially influence the chemical content of this organ. The future of proteomic analysis of isoflavonoid biosynthesis should be centered on the development of quantitative, tissue-specific proteomes that emphasize low-abundance metabolic proteins to extract the whole suite of factors involved.

Using giant scarlet runner bean embryos to uncover regulatory networks controlling suspensor gene activity

One of the major unsolved issues in plant development is understanding the regulatory networks that control the differential gene activity that is required for the specification and development of the two major embryonic regions, the embryo proper and suspensor. Historically, the giant embryo of scarlet runner bean (SRB), Phaseolus coccineus, has been used as a model system to investigate the physiological events that occur early in embryogenesis-focusing on the question of what role the suspensor region plays. A major feature distinguishing SRB embryos from those of other plants is a highly enlarged suspensor containing at least 200 cells that synthesize growth regulators required for subsequent embryonic development. Recent studies have exploited the giant size of the SRB embryo to micro-dissect the embryo proper and suspensor regions in order to use genomics-based approaches to identify regulatory genes that may be involved in controlling suspensor and embryo proper differentiation, as well as the cellular processes that may be unique to each embryonic region. Here we review the current genomics resources that make SRB embryos a compelling model system for studying the early events required to program embryo development.

Comparative Transcriptomic Analyses of Vegetable and Grain Pea (Pisum sativum L.) Seed Development

Understanding the molecular mechanisms regulating pea seed developmental process is extremely important for pea breeding. In this study, we used high-throughput RNA-Seq and bioinformatics analyses to examine the changes in gene expression during seed development in vegetable pea and grain pea, and compare the gene expression profiles of these two pea types. RNA-Seq generated 18.7 G of raw data, which were then de novo assembled into 77,273 unigenes with a mean length of 930 bp. Our results illustrate that transcriptional control during pea seed development is a highly coordinated process. There were 459 and 801 genes differentially expressed at early and late seed maturation stages between vegetable pea and grain pea, respectively. Soluble sugar and starch metabolism related genes were significantly activated during the development of pea seeds coinciding with the onset of accumulation of sugar and starch in the seeds. A comparative analysis of genes involved in sugar and starch biosynthesis in vegetable pea (high seed soluble sugar and low starch) and grain pea (high seed starch and low soluble sugar) revealed that differential expression of related genes at late development stages results in a negative correlation between soluble sugar and starch biosynthetic flux in vegetable and grain pea seeds. RNA-Seq data was validated by using real-time quantitative RT-PCR analysis for 30 randomly selected genes. To our knowledge, this work represents the first report of seed development transcriptomics in pea. The obtained results provide a foundation to support future efforts to unravel the underlying mechanisms that control the developmental biology of pea seeds, and serve as a valuable resource for improving pea breeding.

An RNA-Seq based gene expression atlas of the common bean

BACKGROUND

Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.

RESULTS

Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.

CONCLUSIONS

These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

Screening and identification of soybean seed-specific genes by using integrated bioinformatics of digital differential display, microarray, and RNA-seq data

Soybean is one of the most economically important crops in the world. Soybean seeds have abundant protein and lipid content and very high economic value. In this study, a total of 184 seed-specific genes were obtained using online microarray databases, DDD, and RNA-seq data. The reported seed-specific genes in soybean and the 184 seed-specific genes analyzed in this paper were compared. Of the screened genes, 26 were common to both previous reports and the current screening. Meanwhile, 90 of the 184 genes have homologous counterparts in Arabidopsis, among which 24 have seed-specific expression, as indicated by microarray data for Arabidopsis. Furthermore, promoter analysis showed that almost all seed-specific genes contain at least one seed specific-related element. Seed-specific element Skn-1 motif exists in most, if not all, of the seed-specific genes screened. Five genes were randomly selected from 184 soybean seed specific gene pool and their expressions were quantified using quantitative real time polymerase chain reaction (qRT-PCR) to further confirm the specificity of the screened genes. The results indicated that all five genes showed seed-specific expression. Moreover, the identification of genes with seed-specific expression screened in this study provides information valuable to the in-depth study of soybean.

Legume genomics: understanding biology through DNA and RNA sequencing

BACKGROUND

The legume family (Leguminosae) consists of approx. 17 000 species. A few of these species, including, but not limited to, Phaseolus vulgaris, Cicer arietinum and Cajanus cajan, are important dietary components, providing protein for approx. 300 million people worldwide. Additional species, including soybean (Glycine max) and alfalfa (Medicago sativa), are important crops utilized mainly in animal feed. In addition, legumes are important contributors to biological nitrogen, forming symbiotic relationships with rhizobia to fix atmospheric N2 and providing up to 30 % of available nitrogen for the next season of crops. The application of high-throughput genomic technologies including genome sequencing projects, genome re-sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) by the legume research community has provided major insights into genome evolution, genomic architecture and domestication.

SCOPE AND CONCLUSIONS

This review presents an overview of the current state of legume genomics and explores the role that next-generation sequencing technologies play in advancing legume genomics. The adoption of next-generation sequencing and implementation of associated bioinformatic tools has allowed researchers to turn each species of interest into their own model organism. To illustrate the power of next-generation sequencing, an in-depth overview of the transcriptomes of both soybean and white lupin (Lupinus albus) is provided. The soybean transcriptome focuses on analysing seed development in two near-isogenic lines, examining the role of transporters, oil biosynthesis and nitrogen utilization. The white lupin transcriptome analysis examines how phosphate deficiency alters gene expression patterns, inducing the formation of cluster roots. Such studies illustrate the power of next-generation sequencing and bioinformatic analyses in elucidating the gene networks underlying biological processes.