A joint transcriptomic, proteomic and metabolic analysis of maize endosperm development and starch filling

Auxin and abscisic acid play important roles in promoting glucose metabolism of reactivated young kernels of maize (Zea mays L.)

In maize, young kernels that are less competitive and have poor sink activity often abort. Studies have indicated that such poor competitiveness depends, in part, on the regulation by auxin (IAA) and abscisic acid (ABA). However, the mechanisms for such effects remain unclear. We used pollination-blocking and hand-pollination treatments accompanied by multi-omics and physiological tests, to identify underlying mechanism by which IAA and ABA, along with sugar signaling affect kernel development. Results showed that preventing pollination of the primary ears reactivated kernels in the secondary ears and altered both sugar metabolism and hormone signaling pathways. This was accompanied by increased enzyme activities in carbon metabolism and concentrations of glucose and starch, as well as increased levels of IAA and decreased levels of ABA in the reactivated kernels. Positive and negative correlations were observed between IAA, ABA contents and cell wall invertase (CWIN) activity, and glucose contents, respectively. In vitro culture revealed that the expression of genes involved in glucose utilization was upregulated by IAA, but downregulated by ABA. IAA could promote the expression of ABA signaling genes ZmPP2C9 and ZmPP2C13 but downregulated the expression of Zmnced5, an ABA biosynthesis gene, and ZmSnRK2.10, which is involved in ABA signal transduction. However, these genes showed opposite trends when IAA transport was inhibited. To summarize, we suggest a regulatory model for how IAA inhibits ABA metabolism by promoting the smooth utilization of glucose in reactivated young kernels. Our findings highlight the importance of IAA in ABA signaling by regulating glucose production and transport in maize.

ZmDRR206 Regulates Nutrient Accumulation in Endosperm through Its Role in Cell Wall Biogenesis during Maize Kernel Development

Dirigent proteins (DIRs) contribute to plant fitness by dynamically reorganizing the cell wall and/or by generating defense compounds during plant growth, development, and interactions with environmental stresses. ZmDRR206 is a maize DIR, it plays a role in maintaining cell wall integrity during seedling growth and defense response in maize, but its role in regulating maize kernel development is unclear. Association analysis of candidate genes indicated that the natural variations of ZmDRR206 were significantly associated with maize hundred-kernel weight (HKW). ZmDRR206 plays a dominant role in storage nutrient accumulation in endosperm during maize kernel development, ZmDRR206 overexpression resulted in small and shrunken maize kernel with significantly reduced starch content and significantly decreased HKW. Cytological characterization of the developing maize kernels revealed that ZmDRR206 overexpression induced dysfunctional basal endosperm transfer layer (BETL) cells, which were shorter with less wall ingrowth, and defense response was constitutively activated in developing maize kernel at 15 and 18 DAP by ZmDRR206 overexpression. The BETL-development-related genes and auxin signal-related genes were down-regulated, while cell wall biogenesis-related genes were up-regulated in developing BETL of the ZmDRR206-overexpressing kernel. Moreover, the developing ZmDRR206-overexpressing kernel had significantly reduced contents of the cell wall components such as cellulose and acid soluble lignin. These results suggest that ZmDRR206 may play a regulatory role in coordinating cell development, storage nutrient metabolism, and stress responses during maize kernel development through its role in cell wall biogenesis and defense response, and provides new insights into understanding the mechanisms of kernel development in maize.

Genetic Engineering of Starch Biosynthesis in Maize Seeds for Efficient Enzymatic Digestion of Starch during Bioethanol Production

Maize accumulates large amounts of starch in seeds which have been used as food for human and animals. Maize starch is an importantly industrial raw material for bioethanol production. One critical step in bioethanol production is degrading starch to oligosaccharides and glucose by α-amylase and glucoamylase. This step usually requires high temperature and additional equipment, leading to an increased production cost. Currently, there remains a lack of specially designed maize cultivars with optimized starch (amylose and amylopectin) compositions for bioethanol production. We discussed the features of starch granules suitable for efficient enzymatic digestion. Thus far, great advances have been made in molecular characterization of the key proteins involved in starch metabolism in maize seeds. The review explores how these proteins affect starch metabolism pathway, especially in controlling the composition, size and features of starch. We highlight the roles of key enzymes in controlling amylose/amylopectin ratio and granules architecture. Based on current technological process of bioethanol production using maize starch, we propose that several key enzymes can be modified in abundance or activities via genetic engineering to synthesize easily degraded starch granules in maize seeds. The review provides a clue for developing special maize cultivars as raw material in the bioethanol industry.

Multiomics approach reveals a role of translational machinery in shaping maize kernel amino acid composition

Maize (Zea mays) seeds are a good source of protein, despite being deficient in several essential amino acids. However, eliminating the highly abundant but poorly balanced seed storage proteins has revealed that the regulation of seed amino acids is complex and does not rely on only a handful of proteins. In this study, we used two complementary omics-based approaches to shed light on the genes and biological processes that underlie the regulation of seed amino acid composition. We first conducted a genome-wide association study to identify candidate genes involved in the natural variation of seed protein-bound amino acids. We then used weighted gene correlation network analysis to associate protein expression with seed amino acid composition dynamics during kernel development and maturation. We found that almost half of the proteome was significantly reduced during kernel development and maturation, including several translational machinery components such as ribosomal proteins, which strongly suggests translational reprogramming. The reduction was significantly associated with a decrease in several amino acids, including lysine and methionine, pointing to their role in shaping the seed amino acid composition. When we compared the candidate gene lists generated from both approaches, we found a nonrandom overlap of 80 genes. A functional analysis of these genes showed a tight interconnected cluster dominated by translational machinery genes, especially ribosomal proteins, further supporting the role of translation dynamics in shaping seed amino acid composition. These findings strongly suggest that seed biofortification strategies that target the translation machinery dynamics should be considered and explored further.

Transcriptome profiling provides insights into the molecular mechanisms of maize kernel and silk development

Background

Maize kernel filling, which is closely related to the process of double fertilization and is sensitive to a variety of environmental conditions, is an important component of maize yield determination. Silk is an important tissue of maize ears that can discriminate pollen and conduct pollination. Therefore, investigating the molecular mechanisms of kernel development and silk senescence will provide important information for improving the pollination rate to obtain high maize yields.

Results

In this study, transcript profiles were determined in an elite maize inbred line (KA105) to investigate the molecular mechanisms functioning in self-pollinated and unpollinated maize kernels and silks. A total of 5285 and 3225 differentially expressed transcripts (DETs) were identified between self-pollinated and unpollinated maize in a kernel group and a silk group, respectively. We found that a large number of genes involved in key steps in the biosynthesis of endosperm storage compounds were upregulated after pollination in kernels, and that abnormal development and senescence appeared in unpollinated kernels (KUP). We also identified several genes with functions in the maintenance of silk structure that were highly expressed in silk. Further investigation suggested that the expression of autophagy-related genes and senescence-related genes is prevalent in maize kernels and silks. In addition, pollination significantly altered the expression levels of senescence-related and autophagy-related genes in maize kernels and silks. Notably, we identified some specific genes and transcription factors (TFs) that are highly expressed in single tissues.

Conclusions

Our results provide novel insights into the potential regulatory mechanisms of self-pollinated and unpollinated maize kernels and silks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00981-4.

The comparative gene expression concern to the seed pigmentation in maize (Zea mays L.)

Maize seed pigmentation is one of the important issue to develop maize seed breeding. The differently gene expression was characterized and compared for three inbred lines, such as the pigment accumulated seed (CM22) and non-pigmented seed (CM5 and CM19) at 10 days after pollination. We obtained a total of 63,870, 82,496, and 54,555 contigs by de novo assembly to identify gene expression in the CM22, CM5, and CM19, respectably. In differentially expressed gene analysis, it was revealed that 7,044 genes were differentially expressed by at least two-fold, with 4,067 upregulated in colored maize inbred lines and 2,977 upregulated in colorless maize inbred lines. Of them,18 genes were included to the anthocyanin biosynthesis pathways, while 15 genes were upregulated in both CM22/5 and CM22/19. Additionally, 37 genes were detected in the metabolic pathway concern to the seed pigmentation by BINs analysis using MAPMAN software. Finally, these differently expressed genes may aid in the research on seed pigmentation in maize breeding programs.

The defective effect of starch branching enzyme IIb from weak to strong induces the formation of biphasic starch granules in amylose-extender maize endosperm

Biphasic starch granules in maize ae mutant underwent the weak to strong SBEIIb-defective effect during endosperm development, leading to no birefringence in their exterior due to extended long branch-chains of amylopectin. Biphasic starch granules are usually detected regionally in cereal endosperm lacking starch branching enzyme (SBE). However, their molecular structure, formation mechanism, and regional distribution are unclear. In this research, biphasic starch granules were observed in the inner region of crown endosperm of maize ae mutant, and had poorly oriented structure with comb-like profiles in their exterior. The inner endosperm (IE) rich in biphasic starch granules and outer endosperm (OE) without biphasic starch granules were investigated. The starch had lower amylose content and higher proportion of long branch-chains of amylopectin in IE than in OE, and the exterior of biphasic starch granules had less amylose and more long branch-chains of amylopectin than the interior. Compared with OE, the expression pattern of starch synthesis related enzymes changed significantly in IE. The granule-bound starch synthase I activity within biphasic starch granules decreased slightly. The IE experienced more severe hypoxic stress than OE, and the up-regulated anaerobic respiration pathway indicated an increase in carbon consumption. The starch in IE underwent the SBEIIb-defective effect from weak to strong due to the lack of sufficient carbon inflow, leading to the formation of biphasic starch granules and their regional distribution in endosperm. The results provided information for understanding the biphasic starch granules.

Label-Free Comparative Proteomic Analysis Combined with Laser-Capture Microdissection Suggests Important Roles of Stress Responses in the Black Layer of Maize Kernels

The black layer (BL) is traditionally used as an indicator for kernel harvesting in maize, as it turns visibly dark when the kernel reaches physiological maturity. However, the molecular roles of BL in kernel development have not been fully elucidated. In this work, microscopy images showed that BL began to appear at a growth stage earlier than 10 days after pollination (DAP), and its color gradually deepened to become dark as the development period progressed. Scanning electron microscopy observations revealed that BL is a tissue structure composed of several layers of cells that are gradually squeezed and compressed during kernel development. Laser-capture microdissection (LCM) was used to sample BL and its neighboring inner tissue, basal endosperm transfer layer (BETL), and outer tissue, inner epidermis (IEP), from 20 DAP of kernels. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry profiling (MALDI-TOF MS profiling) detected 41, 104, and 120 proteins from LCM-sampled BL, BETL, and IEP, respectively. Gene ontology (GO) analysis indicated that the 41 BL proteins were primarily involved in the response to stress and stimuli. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis found that the BL proteins were enriched in several defense pathways, such as the ascorbate and aldarate metabolic pathways. Among the 41 BL proteins, six were BL-specific proteins that were only detected from BL. Annotations of five BL-specific proteins were related to stress responses. During kernel development, transcriptional expression of most BL proteins showed an increase, followed by a decrease, and reached a maximum zero to 20 DAP. These results suggest a role for BL in stress responses for protecting filial tissue against threats from maternal sides, which helps to elucidate the biological functions of BL.

Investigation and development of maize fused network analysis with multi-omics

Maize is a critically important staple crop in the whole world, which has contributed to both economic security and food in planting areas. The main target for researchers and breeding is the improvement of maize quality and yield. The use of computational biology methods combined with multi-omics for selecting biomolecules of interest for maize breeding has been receiving more attention. Moreover, the rapid growth of high-throughput sequencing data provides the opportunity to explore biomolecules of interest at the molecular level in maize. Furthermore, we constructed weighted networks for each of the omics and then integrated them into a final fused weighted network based on a nonlinear combination method. We also analyzed the final fused network and mined the orphan nodes, some of which were shown to be transcription factors that played a key role in maize development. This study could help to improve maize production via insights at the multi-omics level and provide a new perspective for maize researchers. All related data have been released at http://lab.malab.cn/∼jj/maize.htm.

Global Dynamics in Protein Disorder during Maize Seed Development

Intrinsic protein disorder is a physicochemical attribute of some proteins lacking tridimensional structure and is collectively known as intrinsically disordered proteins (IDPs). Interestingly, several IDPs have been associated with protective functions in plants and with their response to external stimuli. To correlate the modulation of the IDPs content with the developmental progression in seed, we describe the expression of transcripts according to the disorder content of the proteins that they codify during seed development, from the early embryogenesis to the beginning of the desiccation tolerance acquisition stage. We found that the total expression profile of transcripts encoding for structured proteins is highly increased during middle phase. However, the relative content of protein disorder is increased as seed development progresses. We identified several intrinsically disordered transcription factors that seem to play important roles throughout seed development. On the other hand, we detected a gene cluster encoding for IDPs at the end of the late phase, which coincides with the beginning of the acquisition of desiccation tolerance. In conclusion, the expression pattern of IDPs is highly dependent on the developmental stage, and there is a general reduction in the expression of transcripts encoding for structured proteins as seed development progresses. We proposed maize seeds as a model to study the regulation of protein disorder in plant development and its involvement in the acquisition of desiccation tolerance in plants.

Mitochondrial Zea mays Brittle1-1 Is a Major Determinant of the Metabolic Fate of Incoming Sucrose and Mitochondrial Function in Developing Maize Endosperms

Zea mays Brittle1-1 (ZmBT1-1) is an essential component of the starch biosynthetic machinery in maize endosperms, enabling ADPglucose transport from cytosol to amyloplast in exchange for AMP or ADP. Although ZmBT1-1 has been long considered to be an amyloplast-specific marker, evidence has been provided that ZmBT1-1 is dually localized to plastids and mitochondria (Bahaji et al., 2011b). The mitochondrial localization of ZmBT1-1 suggested that this protein may have as-yet unidentified function(s). To understand the mitochondrial ZmBT1-1 function(s), we produced and characterized transgenic Zmbt1-1 plants expressing ZmBT1-1 delivered specifically to mitochondria. Metabolic and differential proteomic analyses showed down-regulation of sucrose synthase (SuSy)-mediated channeling of sucrose into starch metabolism, and up-regulation of the conversion of sucrose breakdown products generated by cell wall invertase (CWI) into ethanol and alanine, in Zmbt1-1 endosperms compared to wild-type. Electron microscopic analyses of Zmbt1-1 endosperm cells showed gross alterations in the mitochondrial ultrastructure. Notably, the protein expression pattern, metabolic profile, and aberrant mitochondrial ultrastructure of Zmbt1-1 endosperms were rescued by delivering ZmBT1-1 specifically to mitochondria. Results presented here provide evidence that the reduced starch content in Zmbt1-1 endosperms is at least partly due to (i) mitochondrial dysfunction, (ii) enhanced CWI-mediated channeling of sucrose into ethanol and alanine metabolism, and (iii) reduced SuSy-mediated channeling of sucrose into starch metabolism due to the lack of mitochondrial ZmBT1-1. Our results also strongly indicate that (a) mitochondrial ZmBT1-1 is an important determinant of the metabolic fate of sucrose entering the endosperm cells, and (b) plastidic ZmBT1-1 is not the sole ADPglucose transporter in maize endosperm amyloplasts. The possible involvement of mitochondrial ZmBT1-1 in exchange between intramitochondrial AMP and cytosolic ADP is discussed.

Identification of genes alternatively spliced in developing maize endosperm

The process of alternative splicing is critical for the regulation of growth and development of plants. Thus far, little is known about the role of alternative splicing in the regulation of maize (Zea mays L.) endosperm development. RNA sequencing (RNA-seq) data of endosperms from two maize inbred lines, Mo17 and Ji419, at 15 and 25 days after pollination (DAP), respectively, were used to identify genes that were alternatively spliced during endosperm development. Intron retention (IR) in GRMZM2G005887 was further validated using PCR and re-sequencing technologies. In total, 49,000 alternatively spliced events and ca. 20,000 alternatively spliced genes were identified in the two maize inbred lines. Of these, 30 genes involved in amino acid biosynthesis and starch biosynthesis were identified, with IR occurring only in a specific sample, and were significantly co-expressed with ten well-known genes related to maize endosperm development. Moreover, IR in GRMZM2G005887, which encodes a cysteine synthase, was confirmed to occur only in the endosperm of Mo17 at 15 DAP, resulting in the retention of a 121-bp fragment in its 5' untranslated region. Two cis-acting regulatory elements, CAAT-box and TATA-box were observed in the retained fragment in Mo17 at 15 DAP; this could regulate the expression of this gene and influence endosperm development. The results suggest that the 30 genes with IR identified herein might be associated with maize endosperm development, and are likely to play important roles in the developing maize endosperm.

Functions of maize genes encoding pyruvate phosphate dikinase in developing endosperm

Maize opaque2 (o2) mutations are beneficial for endosperm nutritional quality but cause negative pleiotropic effects for reasons that are not fully understood. Direct targets of the bZIP transcriptional regulator encoded by o2 include pdk1 and pdk2 that specify pyruvate phosphate dikinase (PPDK). This enzyme reversibly converts AMP, pyrophosphate, and phosphoenolpyruvate to ATP, orthophosphate, and pyruvate and provides diverse functions in plants. This study addressed PPDK function in maize starchy endosperm where it is highly abundant during grain fill. pdk1 and pdk2 were inactivated individually by transposon insertions, and both genes were simultaneously targeted by endosperm-specific RNAi. pdk2 accounts for the large majority of endosperm PPDK, whereas pdk1 specifies the abundant mesophyll form. The pdk1- mutation is seedling-lethal, indicating that C4 photosynthesis is essential in maize. RNAi expression in transgenic endosperm eliminated detectable PPDK protein and enzyme activity. Transgenic kernels weighed the same on average as nontransgenic siblings, with normal endosperm starch and total N contents, indicating that PPDK is not required for net storage compound synthesis. An opaque phenotype resulted from complete PPDK knockout, including loss of vitreous endosperm character similar to the phenotype conditioned by o2-. Concentrations of multiple glycolytic intermediates were elevated in transgenic endosperm, energy charge was altered, and starch granules were more numerous but smaller on average than normal. The data indicate that PPDK modulates endosperm metabolism, potentially through reversible adjustments to energy charge, and reveal that o2- mutations can affect the opaque phenotype through regulation of PPDK in addition to their previously demonstrated effects on storage protein gene expression.

iTRAQ-Based Proteomics Analysis and Network Integration for Kernel Tissue Development in Maize

Grain weight is one of the most important yield components and a developmentally complex structure comprised of two major compartments (endosperm and pericarp) in maize (Zea mays L.), however, very little is known concerning the coordinated accumulation of the numerous proteins involved. Herein, we used isobaric tags for relative and absolute quantitation (iTRAQ)-based comparative proteomic method to analyze the characteristics of dynamic proteomics for endosperm and pericarp during grain development. Totally, 9539 proteins were identified for both components at four development stages, among which 1401 proteins were non-redundant, 232 proteins were specific in pericarp and 153 proteins were specific in endosperm. A functional annotation of the identified proteins revealed the importance of metabolic and cellular processes, and binding and catalytic activities for the tissue development. Three and 76 proteins involved in 49 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were integrated for the specific endosperm and pericarp proteins, respectively, reflecting their complex metabolic interactions. In addition, four proteins with important functions and different expression levels were chosen for gene cloning and expression analysis. Different concordance between mRNA level and the protein abundance was observed across different proteins, stages, and tissues as in previous research. These results could provide useful message for understanding the developmental mechanisms in grain development in maize.

Proteomics analysis of maize (Zea mays L.) grain based on iTRAQ reveals molecular mechanisms of poor grain filling in inferior grains

In maize, inferior grains (IG) located on the upper part of the ear have poor grain filling process compared to superior grains (SG) located on the middle and lower parts of the ear. This difference limits satisfactory yield and quality; however, the underlying molecular mechanisms remain unknown. Here, using the isobaric tag for relative and absolute quantification (iTRAQ) technology, the proteomes of IG and SG during early and middle grain filling stages were investigated. In total, 4720 proteins were identified in maize grain and 305 differentially accumulated proteins (DiAPs) were detected between IG and SG. These DiAPs were involved in diverse cellular and metabolic processes with preferred distribution in protein synthesis/destination and metabolism. Compared to SG, DiAPs related to cell growth/division and starch synthesis were lag-accumulated and down-regulated in IG, respectively, resulting in smaller sink sizes and lower sink activities in IG. Meanwhile, impediment of the glycolysis pathway in IG may lead to reduce energy supply and building materials for substance synthesis. Additionally, reactive oxygen species (ROS) homeostasis and the defense system were disturbed in IG, which might lead to reduce protection against various environmental stresses. The present study provides new information on the proteomic differences between IG and SG, and explains possible molecular mechanisms for poor grain filling in IG.

Responses to Hypoxia and Endoplasmic Reticulum Stress Discriminate the Development of Vitreous and Floury Endosperms of Conventional Maize (Zea mays) Inbred Lines

Major nutritional and agronomical issues relating to maize (Zea mays) grains depend on the vitreousness/hardness of its endosperm. To identify the corresponding molecular and cellular mechanisms, most studies have been conducted on opaque/floury mutants, and recently on Quality Protein Maize, a reversion of an opaque2 mutation by modifier genes. These mutant lines are far from conventional maize crops. Therefore, a dent and a flint inbred line were chosen for analysis of the transcriptome, amino acid, and sugar metabolites of developing central and peripheral endosperm that is, the forthcoming floury and vitreous regions of mature seeds, respectively. The results suggested that the formation of endosperm vitreousness is clearly associated with significant differences in the responses of the endosperm to hypoxia and endoplasmic reticulum stress. This occurs through a coordinated regulation of energy metabolism and storage protein (i.e., zein) biosynthesis during the grain-filling period. Indeed, genes involved in the glycolysis and tricarboxylic acid cycle are up-regulated in the periphery, while genes involved in alanine, sorbitol, and fermentative metabolisms are up-regulated in the endosperm center. This spatial metabolic regulation allows the production of ATP needed for the significant zein synthesis that occurs at the endosperm periphery; this finding agrees with the zein-decreasing gradient previously observed from the sub-aleurone layer to the endosperm center. The massive synthesis of proteins transiting through endoplasmic reticulum elicits the unfolded protein responses, as indicated by the splicing of bZip60 transcription factor. This splicing is relatively higher at the center of the endosperm than at its periphery. The biological responses associated with this developmental stress, which control the starch/protein balance, leading ultimately to the formation of the vitreous and floury regions of mature endosperm, are discussed.

Proteomic analysis of maize grain development using iTRAQ reveals temporal programs of diverse metabolic processes

BACKGROUND

Grain development in maize is an essential process in the plant's life cycle and is vital for use of the plant as a crop for animals and humans. However, little is known regarding the protein regulatory networks that control grain development. Here, isobaric tag for relative and absolute quantification (iTRAQ) technology was used to analyze temporal changes in protein expression during maize grain development.

RESULTS

Maize grain proteins and changes in protein expression at eight developmental stages from 3 to 50 d after pollination (DAP) were performed using iTRAQ-based proteomics. Overall, 4751 proteins were identified; 2639 of these were quantified and 1235 showed at least 1.5-fold changes in expression levels at different developmental stages and were identified as differentially expressed proteins (DEPs). The DEPs were involved in different cellular and metabolic processes with a preferential distribution to protein synthesis/destination and metabolism categories. A K-means clustering analysis revealed coordinated protein expression associated with different functional categories/subcategories at different development stages.

CONCLUSIONS

Our results revealed developing maize grain display different proteomic characteristics at distinct stages, such as numerous DEPs for cell growth/division were highly expressed during early stages, whereas those for starch biosynthesis and defense/stress accumulated in middle and late stages, respectively. We also observed coordinated expression of multiple proteins of the antioxidant system, which are essential for the maintenance of reactive oxygen species (ROS) homeostasis during grain development. Particularly, some DEPs, such as zinc metallothionein class II, pyruvate orthophosphate dikinase (PPDK) and 14-3-3 proteins, undergo major changes in expression at specific developmental stages, suggesting their roles in maize grain development. These results provide a valuable resource for analyzing protein function on a global scale and also provide new insights into the potential protein regulatory networks that control grain yield and quality.

Characterization of factors underlying the metabolic shifts in developing kernels of colored maize

Elucidation of the metabolic pathways determining pigmentation and their underlying regulatory mechanisms in maize kernels is of high importance in attempts to improve the nutritional composition of our food. In this study, we compared dynamics in the transcriptome and metabolome between colored SW93 and white SW48 by integrating RNA-Seq and non-targeted metabolomics. Our data revealed that expression of enzyme coding genes and levels of primary metabolites decreased gradually from 11 to 21 DAP, corresponding well with the physiological change of developing maize kernels from differentiation through reserve accumulation to maturation, which was cultivar independent. A remarkable up-regulation of anthocyanin and phlobaphene pathway distinguished SW93 from SW48, in which anthocyanin regulating transcriptional factors (R1 and C1), enzyme encoding genes involved in both pathways and corresponding metabolic intermediates were up-regulated concurrently in SW93 but not in SW48. The shift from the shikimate pathway of primary metabolism to the flavonoid pathway of secondary metabolism, however, appears to be under posttranscriptional regulation. This study revealed the link between primary metabolism and kernel coloration, which facilitate further study to explore fundamental questions regarding the evolution of seed metabolic capabilities as well as their potential applications in maize improvement regarding both staple and functional foods.

Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel

Starch is the major component in maize kernels, providing a stable carbohydrate source for humans and livestock as well as raw material for the biofuel industry. Increasing maize kernel starch content will help meet industry demands and has the potential to increase overall yields. We developed a pair of maize near-isogenic lines (NILs) with different alleles for a starch quantitative trait locus on chromosome 3 (qHS3), resulting in different kernel starch content. To investigate the candidate genes for qHS3 and elucidate their effects on starch metabolism, RNA-Seq was performed for the developing kernels of the NILs at 14 and 21 d after pollination (DAP). Analysis of genomic and transcriptomic data identified 76 genes with nonsynonymous single nucleotide polymorphisms and 384 differentially expressed genes (DEGs) in the introgressed fragment, including a hexokinase gene, ZmHXK3a, which catalyzes the conversion of glucose to glucose-6-phosphate and may play a key role in starch metabolism. The expression pattern of all DEGs in starch metabolism shows that altered expression of the candidate genes for qHS3 promoted starch synthesis, with positive consequences for kernel starch content. These results expand the current understanding of starch biosynthesis and accumulation in maize kernels and provide potential candidate genes to increase starch content.

Dynamic Proteomic Characteristics and Network Integration Revealing Key Proteins for Two Kernel Tissue Developments in Popcorn

The formation and development of maize kernel is a complex dynamic physiological and biochemical process that involves the temporal and spatial expression of many proteins and the regulation of metabolic pathways. In this study, the protein profiles of the endosperm and pericarp at three important developmental stages were analyzed by isobaric tags for relative and absolute quantification (iTRAQ) labeling coupled with LC-MS/MS in popcorn inbred N04. Comparative quantitative proteomic analyses among developmental stages and between tissues were performed, and the protein networks were integrated. A total of 6,876 proteins were identified, of which 1,396 were nonredundant. Specific proteins and different expression patterns were observed across developmental stages and tissues. The functional annotation of the identified proteins revealed the importance of metabolic and cellular processes, and binding and catalytic activities for the development of the tissues. The whole, endosperm-specific and pericarp-specific protein networks integrated 125, 9 and 77 proteins, respectively, which were involved in 54 KEGG pathways and reflected their complex metabolic interactions. Confirmation for the iTRAQ endosperm proteins by two-dimensional gel electrophoresis showed that 44.44% proteins were commonly found. However, the concordance between mRNA level and the protein abundance varied across different proteins, stages, tissues and inbred lines, according to the gene cloning and expression analyses of four relevant proteins with important functions and different expression levels. But the result by western blot showed their same expression tendency for the four proteins as by iTRAQ. These results could provide new insights into the developmental mechanisms of endosperm and pericarp, and grain formation in maize.

Proteomic analysis of the enzymes involved in the starch biosynthesis of maize with different endosperm type and characterization of the starch

BACKGROUND

The characterization of starch maize with different endosperm type and the proteomic analysis of its biosynthetic enzymes at 20 and 50 days after pollination (DAP) was carried out.

RESULTS

There were differences between both endosperm types at 20 DAP, mainly in starch accumulation, amylose content, granule size and crystallinity percentage, whereas at 50 DAP the differences found were not relevant in the case of starch content, granule size, chain length distribution and thermal properties. SSSI, SBEIIb and GBSSI enzymes were identified; however, SBEIIb was only identified in two samples: floury endosperm at 20 DAP and vitreous at 50 DAP.

CONCLUSION

Starch did not show differences in its morphological or structural characteristics in either endosperm on reaching maturity. Starch biosynthetic enzymes identified by matrix-assisted laser desorption/ionization-time of flight did not show a relationship to starch structure.

The differential transcription network between embryo and endosperm in the early developing maize seed

Transcriptome analysis of early-developing maize (Zea mays) seed was conducted using Illumina sequencing. We mapped 11,074,508 and 11,495,788 paired-end reads from endosperm and embryo, respectively, at 9 d after pollination to define gene structure and alternative splicing events as well as transcriptional regulators of gene expression to quantify transcript abundance in both embryo and endosperm. We identified a large number of novel transcribed regions that did not fall within maize annotated regions, and many of the novel transcribed regions were tissue-specifically expressed. We found that 50.7% (8,556 of 16,878) of multiexonic genes were alternatively spliced, and some transcript isoforms were specifically expressed either in endosperm or in embryo. In addition, a total of 46 trans-splicing events, with nine intrachromosomal events and 37 interchromosomal events, were found in our data set. Many metabolic activities were specifically assigned to endosperm and embryo, such as starch biosynthesis in endosperm and lipid biosynthesis in embryo. Finally, a number of transcription factors and imprinting genes were found to be specifically expressed in embryo or endosperm. This data set will aid in understanding how embryo/endosperm development in maize is differentially regulated.

Proteomic identification of genes associated with maize grain-filling rate

Grain filling during the linear phase contributes most of the dry matter accumulated in the maize kernel, which in turn determines the final grain yield. Endosperms and embryos of three elite maize hybrids (Zhengdan 958, Nongda 108, and Pioneer 335) were sampled 17, 22, 25, and 28 days after pollination, during the linear phase of grain filling, for proteomic analysis to explore the regulatory factors critical for grain filling rate. In total, 39 and 43 protein spots that showed more than 2-fold changes in abundance at P<0.01 between any two sampling stages in the endosperm and embryo were analyzed by protein mass spectrometry. The changing patterns in expression index of these proteins in the endosperm were evenly distributed, whereas up-regulation patterns predominated (74%) in the embryo. Functional analysis revealed that metabolism was the largest category, represented by nine proteins in the endosperm and 12 proteins in the embryo, of the proteins that significantly changed in abundance. Glycolysis, a critical process both for glucose conversion into pyruvate and for release of free energy and reducing power, and proteins related to redox homeostasis were emphasized in the endosperm. Additionally, lipid, nitrogen, and inositol metabolism related to fatty acid biosynthesis and late embryogenesis abundant proteins were emphasized in the embryo. One protein related to cellular redox equilibrium, which showed a more than 50-fold change in abundance and was co-localized with a quantitative trait locus for grain yield on chromosome 1, was further investigated by transcriptional profile implying consistent expression pattern with protein accumulation. The present results provide a first step towards elucidation of the gene network responsible for regulation of grain filling in maize.

Deep sequencing of RNA from ancient maize kernels

The characterization of biomolecules from ancient samples can shed otherwise unobtainable insights into the past. Despite the fundamental role of transcriptomal change in evolution, the potential of ancient RNA remains unexploited – perhaps due to dogma associated with the fragility of RNA. We hypothesize that seeds offer a plausible refuge for long-term RNA survival, due to the fundamental role of RNA during seed germination. Using RNA-Seq on cDNA synthesized from nucleic acid extracts, we validate this hypothesis through demonstration of partial transcriptomal recovery from two sources of ancient maize kernels. The results suggest that ancient seed transcriptomics may offer a powerful new tool with which to study plant domestication.

Cofactome analyses reveal enhanced flux of carbon into oil for potential biofuel production

To identify the underlying molecular basis of carbon partitioning between starch and oil we conducted 454 pyrosequencing, followed by custom microarrays to profile gene expression throughout endosperm development, of two closely related oat cultivars that differ in oil content at the expense of starch as determined by several approaches including non-invasive magnetic resonance imaging. Comparative transcriptome analysis in conjunction with metabolic profiling displays a close coordination between energy metabolism and carbon partitioning pathways, with increased demands for energy and reducing equivalents in kernels with a higher oil content. These studies further expand the repertoire of networks regulating carbon partitioning to those involved in metabolism of cofactors, suggesting that an elevated supply of cofactors, here called cofactomes, contribute to the allocation of higher carbon pools for production of oils and storage proteins. These data highlight a close association between cofactomes and carbon partitioning, thereby providing a biotechnological target for conversion of starch to oil.

Analysis of amino acid metabolism in the ear of maize mutants deficient in two cytosolic glutamine synthetase isoenzymes highlights the importance of asparagine for nitrogen translocation within sink organs

Nitrogen (N) metabolism was characterized in the developing ear of glutamine synthetase deficient mutants (gln1-3, gln1-4 and gln1-3/gln1-4) of maize exhibiting a reduction in kernel yield. During the grain-filling period, the metabolite contents, enzyme activities and steady-state levels of transcripts for marker genes of amino acid synthesis and interconversion were monitored in the cob and kernels. The ear of gln1-3 and gln1-3/gln1-4 had a higher free amino acid content and a lower C/N ratio, when compared to the wild type. The free ammonium concentrations were also much higher in gln1-3/gln1-4, and Asn accumulation was higher in gln1-3 and gln1-3/gln1-4. The level of transcripts of ZmAS3 and ZmAS4, two genes encoding asparagine synthetase, increased in the 'aborted kernels' of gln1-3 and gln1-3/gln1-4. The results show that N metabolism is clearly different in developing and 'aborted kernels'. The data support the hypothesis that N accumulated in 'aborted kernels' is remobilized via the cob to developing kernels using Asn as a transport molecule. The two genes ZmAS3 and ZmAS4 are likely to play an important role during this process.

Computational identification of root hair-specific genes in Arabidopsis

Activated cortical domains (ACDs) are regions of the plant cell cortex performing localized membrane turnover, delimited by concerted action of the cortical cytoskeleton and endomembrane compartments. Arabidopsis thaliana rhizodermis consists of two cell types differing by a single ACD (trichoblasts, carrying tip-growing root hairs, and hairless atrichoblasts), providing a model for the study of ACD determination. We compiled a set of genes specifically upregulated in root hairs from published transcriptome data, and compared it with a "virtual Arabidopsis root hair proteome", i.e. a list of computationally identified homologs of proteins from the published soybean root hair proteome. Both data sets were enriched in genes and proteins associated with root hairs in functional studies, but there was little overlap between the transcriptome and the proteome: the former captured gene products specific to root hairs, while the latter selected those abundant in root hairs but not necessarily specific to them. Decisive steps in ACD specification may be performed by signaling proteins of high expression specifity and low abundance. Nevertheless, 73 genes specifically transcribed in Arabidopsis trichoblasts or root hairs encode homologs of abundant root hair proteins from soybean. Most of them encode "housekeeping" proteins required for rapid tip growth. However, among the "candidates" is also a generative actin isoform, ACT11. Preliminary characterization of an act11 mutant allele indeed suggests a hitherto unexpected role for this gene in root and root hair development.

Integrated proteomic and cytological study of rice endosperms at the storage phase

The endosperm at the storage phase undergoes a series of coordinated cellular and metabolic events, including starchy endosperm cell death, starch synthesis, and starch granule packaging, which leads to efficient accumulation of starch. However, the mechanism underlying the interconnections remains unknown. We used integrated proteomic and cytological approaches to probe the interconnections in rice (Oryza sativa) endosperm at the storage phase from 12 to 18 days after flowering (DAF). Starch granule packaging was completed first in the inner part of endosperm at 15 DAF and spread to almost the entire endosperm at 18 DAF. Programmed starchy endosperm cell death occurred after the starch granule packaging. Endogenous H(2)O(2) was detectable in the inner part of endosperm at 12 DAF and the region beyond the inner part at 15 DAF, with an H(2)O(2) burst at 15 DAF. Proteomics analysis with 2-D fluorescent difference gel electrophoresis and matrix-assisted laser-desorption ionization time-of-flight/time-of-flight mass spectrometry revealed 317 proteins, including almost all known antioxidants, differentially expressed throughout the 3 stages of the developmental phase. More than two-thirds of the 317 proteins were potential thioredoxin targets, with a preferential skew toward central carbon metabolism, alcoholic fermentation, starch metabolism, amino acid metabolism, and protein synthesis or folding. These proteins implicated in starch synthesis and gluconeogenesis were upregulated, whereas those involved in anabolism of biomacromolecules such as proteins, lipids, and cell wall components were downregulated, with upregulated expression of proteins involved in catabolism of these biomacromolecules, which suggests remobilization of nutrients for starch synthesis. These data suggested important roles of the H(2)O(2)-antioxidant interface in coordinating starch accumulation, programmed cell death of starchy endosperm, and remobilization of nutrients during the cell death.

An expression analysis of 57 transcription factors derived from ESTs of developing seeds in Maize (Zea mays)

Maize seeds are an important source of food, animal feed, and industrial raw materials. To understand global gene expression and regulation during maize seed development, a normalized cDNA library, covering most of the developmental stages of maize seeds, was constructed. Sequencing analysis of 10,848 randomly selected clones identified 6,630 unique ESTs. Among them, 57 putative transcription factors (TFs) were identified. The TFs belong to seven different super-families, specifically 17 Zinc-finger, 13 bZIP, 8 bHLH, 6 MADS, 7 MYB, 3 Homedomain, and 3 AP2/EREBP. The spatial and temporal expression of the TFs was analyzed by semi-quantitative RT-PCR with representative tissue types and seeds at different developmental stages, revealing their diverse expression patterns and expression levels. One-third (19) of the maize TFs was found their putative orthologs in Arabidopsis. Similar expression patterns were observed in both maize and Arabidopsis for the majority of orthologous pairs (15 out of 19), suggesting their conserved functions during seed development. In conclusion, the systematic analysis of maize seed TFs has provided valuable insight into transcriptional regulation during maize seed development.

Systems analysis of seed filling in Arabidopsis: using general linear modeling to assess concordance of transcript and protein expression

Previous systems analyses in plants have focused on a single developmental stage or time point, although it is often important to additionally consider time-index changes. During seed development a cascade of events occurs within a relatively brief time scale. We have collected protein and transcript expression data from five sequential stages of Arabidopsis (Arabidopsis thaliana) seed development encompassing the period of reserve polymer accumulation. Protein expression profiling employed two-dimensional gel electrophoresis coupled with tandem mass spectrometry, while transcript profiling used oligonucleotide microarrays. Analyses in biological triplicate yielded robust expression information for 523 proteins and 22,746 genes across the five developmental stages, and established 319 protein/transcript pairs for subsequent pattern analysis. General linear modeling was used to evaluate the protein/transcript expression patterns. Overall, application of this statistical assessment technique showed concurrence for a slight majority (56%) of expression pairs. Many specific examples of discordant protein/transcript expression patterns were detected, suggesting that this approach will be useful in revealing examples of posttranscriptional regulation.

Barley grain development toward an integrative view

Seeds are complex structures composed of several maternal and filial tissues which undergo rapid changes during development. In this review, the barley grain is taken as a cereal seed model. Following a brief description of the developing grain, recent progress in grain development modeling is described. 3-D/4-D models based on histological sections or nondestructive NMR measurements can be used to integrate a variety of datasets. Extensive transcriptome data are taken as a frame to augment our understanding of various molecular-physiological processes. Discussed are maternal influences on grain development and the role of different tissues (pericarp, nucellus, nucellar projection, endosperm, endosperm transfer cells). Programmed cell death (PCD) is taken to pinpoint tissue specificities and the importance of remobilization processes for grain development. Transcriptome data have also been used to derive transcriptional networks underlying differentiation and maturation in endosperm and embryo. They suggest that the "maturation hormone" ABA is important also in early grain development. Massive storage product synthesis during maturation is dependent on sufficient energy, which can only be provided by specific metabolic adaptations due to severe oxygen deficiencies within the seed. To integrate the great variety of data from different research areas in complex, predictive computational modeling as part of a systems biology approach is an important challenge of the future. First attempts of modeling barley grain metabolism are summarized.