Genome-wide transcriptome analyses of developing seeds from low and normal phytic acid soybean lines

Network Inference of Transcriptional Regulation in Germinating Low Phytic Acid Soybean Seeds

The low phytic acid (lpa) trait in soybeans can be conferred by loss-of-function mutations in genes encoding myo-inositol phosphate synthase and two epistatically interacting genes encoding multidrug-resistance protein ATP-binding cassette (ABC) transporters. However, perturbations in phytic acid biosynthesis are associated with poor seed vigor. Since the benefits of the lpa trait, in terms of end-use quality and sustainability, far outweigh the negatives associated with poor seed performance, a fuller understanding of the molecular basis behind the negatives will assist crop breeders and engineers in producing variates with lpa and better germination rate. The gene regulatory network (GRN) for developing low and normal phytic acid soybean seeds was previously constructed, with genes modulating a variety of processes pertinent to phytic acid metabolism and seed viability being identified. In this study, a comparative time series analysis of low and normal phytic acid soybeans was carried out to investigate the transcriptional regulatory elements governing the transitional dynamics from dry seed to germinated seed. GRNs were reverse engineered from time series transcriptomic data of three distinct genotypic subsets composed of lpa soybean lines and their normal phytic acid sibling lines. Using a robust unsupervised network inference scheme, putative regulatory interactions were inferred for each subset of genotypes. These interactions were further validated by published regulatory interactions found in Arabidopsis thaliana and motif sequence analysis. Results indicate that lpa seeds have increased sensitivity to stress, which could be due to changes in phytic acid levels, disrupted inositol phosphate signaling, disrupted phosphate ion (Pi) homeostasis, and altered myo-inositol metabolism. Putative regulatory interactions were identified for the latter two processes. Changes in abscisic acid (ABA) signaling candidate transcription factors (TFs) putatively regulating genes in this process were identified as well. Analysis of the GRNs reveal altered regulation in processes that may be affecting the germination of lpa soybean seeds. Therefore, this work contributes to the ongoing effort to elucidate molecular mechanisms underlying altered seed viability, germination and field emergence of lpa crops, understanding of which is necessary in order to mitigate these problems.

Transcriptome analysis identifies differentially expressed genes in the progenies of a cross between two low phytic acid soybean mutants

Phytic acid (PA) is a major antinutrient that cannot be digested by monogastric animals, but it can decrease the bioavailability of micronutrients (e.g., Zn and Fe). Lowering the PA content of crop seeds will lead to enhanced nutritional traits. Low-PA mutant crop lines carrying more than one mutated gene (lpa) have lower PA contents than mutants with a single lpa mutant gene. However, little is known about the link between PA pathway intermediates and downstream regulatory activities following the mutation of these genes in soybean. Consequently, we performed a comparative transcriptome analysis using an advanced generation recombinant inbred line with low PA levels [2mlpa (mips1/ipk1)] and a sibling line with homozygous non-mutant alleles and normal PA contents [2MWT (MIPS1/IPK1)]. An RNA sequencing analysis of five seed developmental stages revealed 7945 differentially expressed genes (DEGs) between the 2mlpa and 2MWT seeds. Moreover, 3316 DEGs were associated with 128 metabolic and signal transduction pathways and 4980 DEGs were annotated with 345 Gene Ontology terms related to biological processes. Genes associated with PA metabolism, photosynthesis, starch and sucrose metabolism, and defense mechanisms were among the DEGs in 2mlpa. Of these genes, 36 contributed to PA metabolism, including 22 genes possibly mediating the low-PA phenotype of 2mlpa. The expression of most of the genes associated with photosynthesis (81 of 117) was down-regulated in 2mlpa at the late seed developmental stage. In contrast, the expression of three genes involved in sucrose metabolism was up-regulated at the late seed developmental stage, which might explain the high sucrose content of 2mlpa soybeans. Furthermore, 604 genes related to defense mechanisms were differentially expressed between 2mlpa and 2MWT. In this study, we detected a low PA content as well as changes to multiple metabolites in the 2mlpa mutant. These results may help elucidate the regulation of metabolic events in 2mlpa. Many genes involved in PA metabolism may contribute to the substantial decrease in the PA content and the moderate accumulation of InsP3-InsP5 in the 2mlpa mutant. The other regulated genes related to photosynthesis, starch and sucrose metabolism, and defense mechanisms may provide additional insights into the nutritional and agronomic performance of 2mlpa seeds.

Developmental and genomic architecture of plant embryogenesis: from model plant to crops

Embryonic development represents an important reproductive phase of sexually reproducing plant species. The fusion of egg and sperm produces the plant zygote, a totipotent cell that, through cell division and cell identity specification in early embryogenesis, establishes the major cell lineages and tissues of the adult plant. The subsequent morphogenesis phase produces the full-sized embryo, while the late embryogenesis maturation process prepares the seed for dormancy and subsequent germination, ensuring continuation of the plant life cycle. In this review on embryogenesis, we compare the model eudicot Arabidopsis thaliana with monocot crops, focusing on genome activation, paternal and maternal regulation of early zygote development, and key organizers of patterning, such as auxin and WOX transcription factors. While the early stages of embryo development are apparently conserved among plant species, embryo maturation programs have diversified between eudicots and monocots. This diversification in crop species reflects the likely effects of domestication on seed quality traits that are determined during embryo maturation, and also assures seed germination in different environmental conditions. This review describes the most important features of embryonic development in plants, and the scope and applications of genomics in plant embryo studies.

Diverse role of phytic acid in plants and approaches to develop low-phytate grains to enhance bioavailability of micronutrients

Natural or synthetic compounds that interfere with the bioavailability of nutrients are called antinutrients. Phytic acid (PA) is one of the major antinutrients present in the grains and acts as a chelator of micronutrients. The presence of six reactive phosphate groups in PA hinders the absorption of micronutrients in the gut of non-ruminants. Consumption of PA-rich diet leads to deficiency of minerals such as iron and zinc among human population. On the contrary, PA is a natural antioxidant, and PA-derived molecules function in various signal transduction pathways. Therefore, optimal concentration of PA needs to be maintained in plants to avoid adverse pleiotropic effects, as well as to ensure micronutrient bioavailability in the diets. Given this, the chapter enumerates the structure, biosynthesis, and accumulation of PA in food grains followed by their roles in growth, development, and stress responses. Further, the chapter elaborates on the antinutritional properties of PA and explains the conventional breeding and transgene-based approaches deployed to develop low-PA varieties. Studies have shown that conventional breeding methods could develop low-PA lines; however, the pleiotropic effects of these methods viz. reduced yield, embryo abnormalities, and poor seed quality hinder the use of breeding strategies. Overexpression of phytase in the endosperm and RNAi-mediated silencing of genes involved in myo-inositol biosynthesis overcome these constraints. Next-generation genome editing approaches, including CRISPR-Cas9 enable the manipulation of more than one gene involved in PA biosynthesis pathway through multiplex editing, and scope exists to deploy such tools in developing varieties with optimal PA levels.

Genetic interactions regulating seed phytate and oligosaccharides in soybean (Glycine max L.)

Two low-phytate soybean (Glycine max (L.) Merr.) mutant lines- V99-5089 (mips mutation on chromosome 11) and CX-1834 (mrp-l and mrp-n mutations on chromosomes 19 and 3, respectively) have proven to be valuable resources for breeding of low-phytate, high-sucrose, and low-raffinosaccharide soybeans, traits that are highly desirable from a nutritional and environmental standpoint. A recombinant inbred population derived from the cross CX1834 x V99-5089 provides an opportunity to study the effect of different combinations of these three mutations on soybean phytate and oligosaccharides levels. Of the 173 recombinant inbred lines tested, 163 lines were homozygous for various combinations of MIPS and two MRP loci alleles. These individuals were grouped into eight genotypic classes based on the combination of SNP alleles at the three mutant loci. The two genotypic classes that were homozygous mrp-l/mrp-n and either homozygous wild-type or mutant at the mips locus (MIPS/mrp-l/mrp-n or mips/mrp-l/mrp-n) displayed relatively similar ~55% reductions in seed phytate, 6.94 mg g -1 and 6.70 mg g-1 respectively, as compared with 15.2 mg g-1 in the wild-type MIPS/MRP-L/MRP-N seed. Therefore, in the presence of the double mutant mrp-l/mrp-n, the mips mutation did not cause a substantially greater decrease in seed phytate level. However, the nutritionally-desirable high-sucrose/low-stachyose/low-raffinose seed phenotype originally observed in soybeans homozygous for the mips allele was reversed in the presence of mrp-l/mrp-n mutations: homozygous mips/mrp-l/mrp-n seed displayed low-sucrose (7.70%), high-stachyose (4.18%), and the highest observed raffinose (0.94%) contents per gram of dry seed. Perhaps the block in phytic acid transport from its cytoplasmic synthesis site to its storage site, conditioned by mrp-l/mrp-n, alters myo-inositol flux in mips seeds in a way that restores to wild-type levels the mips conditioned reductions in raffinosaccharides. Overall this study determined the combinatorial effects of three low phytic acid causing mutations on regulation of seed phytate and oligosaccharides in soybean.

Transcriptomic Insights into Mechanisms of Early Seed Maturation in the Garden Pea (Pisum sativum L.)

The garden pea (Pisum sativum L.) is a legume crop of immense economic value. Extensive breeding has led to the emergence of numerous pea varieties, of which some are distinguished by accelerated development in various stages of ontogenesis. One such trait is rapid seed maturation, which, despite novel insights into the genetic control of seed development in legumes, remains poorly studied. This article presents an attempt to dissect mechanisms of early maturation in the pea line Sprint-2 by means of whole transcriptome RNA sequencing in two developmental stages. By using a de novo assembly approach, we have obtained a reference transcriptome of 25,756 non-redundant entries expressed in pea seeds at either 10 or 20 days after pollination. Differential expression in Sprint-2 seeds has affected 13,056 transcripts. A comparison of the two pea lines with a common maturation rate demonstrates that while at 10 days after pollination, Sprint-2 seeds show development retardation linked to intensive photosynthesis, morphogenesis, and cell division, and those at 20 days show a rapid onset of desiccation marked by the cessation of translation and cell anabolism and accumulation of dehydration-protective and -storage moieties. Further inspection of certain transcript functional categories, including the chromatin constituent, transcription regulation, protein turnover, and hormonal regulation, has revealed transcriptomic trends unique to specific stages and cultivars. Among other remarkable features, Sprint-2 demonstrated an enhanced expression of transposable element-associated open reading frames and an altered expression of major maturation regulators and DNA methyltransferase genes. To the best of our knowledge, this is the first comparative transcriptomic study in which the issue of the seed maturation rate is addressed.

Expression regulation of myo-inositol 3-phosphate synthase 1 (INO1) in determination of phytic acid accumulation in rice grain

Phytic acid (PA) is the primary phosphorus (P) storage compound in the seeds of cereals and legumes. Low PA crops, which are considered an effective way to improve grain nutrient availability and combat environmental issues relating to seed P have been developed using mutational and reverse genetics approaches. Here, we identify molecular mechanism regulating PA content among natural rice variants. First, we performed genome-wide association (GWA) mapping of world rice core collection (WRC) accessions to understand the genetic determinants underlying PA trait in rice. Further, a comparative study was undertaken to identify the differences in PA accumulation, protein profiles, and gene expression in low (WRC 5) and high PA (WRC 6) accessions. GWA results identified myo-inositol 3-phosphate synthase 1 (INO1) as being closely localized to a significant single nucleotide polymorphism. We found high rates of PA accumulation 10 days after flowering, and our results indicate that INO1 expression was significantly higher in WRC 6 than in WRC 5. Seed proteome assays found that the expression of INO1 was significantly higher in WRC 6. These results suggest that not only the gene itself but regulation of INO1 gene expression at early developmental stages is important in determining PA content in rice.

RNAi-Mediated Downregulation of Inositol Pentakisphosphate Kinase (IPK1) in Wheat Grains Decreases Phytic Acid Levels and Increases Fe and Zn Accumulation

Enhancement of micronutrient bioavailability is crucial to address the malnutrition in the developing countries. Various approaches employed to address the micronutrient bioavailability are showing promising signs, especially in cereal crops. Phytic acid (PA) is considered as a major antinutrient due to its ability to chelate important micronutrients and thereby restricting their bioavailability. Therefore, manipulating PA biosynthesis pathway has largely been explored to overcome the pleiotropic effect in different crop species. Recently, we reported that functional wheat inositol pentakisphosphate kinase (TaIPK1) is involved in PA biosynthesis, however, the functional roles of the IPK1 gene in wheat remains elusive. In this study, RNAi-mediated gene silencing was performed for IPK1 transcripts in hexaploid wheat. Four non-segregating RNAi lines of wheat were selected for detailed study (S3-D-6-1; S6-K-3-3; S6-K-6-10 and S16-D-9-5). Homozygous transgenic RNAi lines at T₄ seeds with a decreased transcript of TaIPK1 showed 28-56% reduction of the PA. Silencing of IPK1 also resulted in increased free phosphate in mature grains. Although, no phenotypic changes in the spike was observed but, lowering of grain PA resulted in the reduced number of seeds per spikelet. The lowering of grain PA was also accompanied by a significant increase in iron (Fe) and zinc (Zn) content, thereby enhancing their molar ratios (Zn:PA and Fe:PA). Overall, this work suggests that IPK1 is a promising candidate for employing genome editing tools to address the mineral accumulation in wheat grains.

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/.

Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with seed development and seed size/weight determination in chickpea

Seed development is an intricate process regulated via a complex transcriptional regulatory network. To understand the molecular mechanisms governing seed development and seed size/weight in chickpea, we performed a comprehensive analysis of transcriptome dynamics during seed development in two cultivars with contrasting seed size/weight (small-seeded, Himchana 1 and large-seeded, JGK 3). Our analysis identified stage-specific expression for a significant proportion (>13%) of the genes in each cultivar. About one half of the total genes exhibited significant differential expression in JGK 3 as compared with Himchana 1. We found that different seed development stages can be delineated by modules of coexpressed genes. A comparative analysis revealed differential developmental stage specificity of some modules between the two cultivars. Furthermore, we constructed transcriptional regulatory networks and identified key components determining seed size/weight. The results suggested that extended period of cell division during embryogenesis and higher level of endoreduplication along with more accumulation of storage compounds during maturation determine large seed size/weight. Further, we identified quantitative trait loci-associated candidate genes harboring single nucleotide polymorphisms in the promoter sequences that differentiate small- and large-seeded chickpea cultivars. The results provide a valuable resource to dissect the role of candidate genes governing seed development and seed size/weight in chickpea.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Tissue specific transcript profiling of wheat phosphate transporter genes and its association with phosphate allocation in grains

Approaches enabling efficient phosphorus utilization in crops are of great importance. In cereal crop like wheat, utilization of inorganic phosphate (Pi) is high and mature grains are the major sink for Pi utilization and storage. Research that addresses the importance of the Pi homeostasis in developing grains is limited. In an attempt to understand the Pi homeostasis in developing wheat grains, we identified twelve new phosphate transporters (PHT), these are phyologentically well distributed along with the members reported from Arabidopsis and rice. Enhanced expression of PHT1-subfamily genes was observed in roots subjected to the Pi starvation suggesting their active role in Pi homeostasis. Differential expression patterns of all the PHT genes during grain filling stages suggested their importance in the filial tissues. Additionally, high accumulation of Pi and total P in aleurone correlates well with the expression of TaPHTs and other phosphate starvation related genes. Tissue specific transcript accumulation of TaPHT1.1, TaPHT1.2, TaPHT1.4 in aleurone; TaPHT3.1 in embryo and TaPHT4.2 in the endosperm was observed. Furthermore, their transcript abundance was affected in low phytate wheat grains. Altogether, this study helps in expanding the knowledge and prioritize the candidate wheat Pi-transporters to modulate the Pi homeostasis in cereal grains.