Endosperm-preferred expression of maize genes as revealed by transcriptome-wide analysis of expressed sequence tags

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.

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.

A PCR-based forward genetics screening, using expression domain-specific markers, identifies mutants in endosperm transfer cell development

Mutant collections are an invaluable source of material on which forward genetic approaches allow the identification of genes affecting a wide variety of biological processes. However, some particular developmental stages and morphological structures may resist analysis due to their physical inaccessibility or to deleterious effects associated to their modification. Furthermore, lethal mutations acting early in development may escape detection. We have approached the characterization of 101 maize seed mutants, selected from a collection of 27,500 visually screened Mu-insertion lines, using a molecular marker approach based on a set of genes previously ascribed to different tissue compartments within the early developing kernel. A streamlined combination of qRT-PCR assays has allowed us to preliminary pinpoint the affected compartment, establish developmental comparisons to WT siblings and select mutant lines with alterations in the different compartments. Furthermore, clusters of markers co-affected by the underlying mutation were identified. We have analyzed more extensively a set of lines presenting significant variation in transfer cell-associated expression markers, and have performed morphological observations, and immunolocalization experiments to confirm the results, validating this approach as an efficient mutant description tool.

Transcriptomic analysis of rice (Oryza sativa) endosperm using the RNA-Seq technique

The endosperm plays an important role in seed formation and germination, especially in rice (Oryza sativa). We used a high-throughput sequencing technique (RNA-Seq) to reveal the molecular mechanisms involved in rice endosperm development. Three cDNA libraries were taken from rice endosperm at 3, 6 and 10 days after pollination (DAP), which resulted in the detection of 21,596, 20,910 and 19,459 expressed gens, respectively. By ERANGE, we identified 10,371 differentially expressed genes (log(2)Ratio ≥1, FDR ≤0.001). The results were compared against three public databases (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and MapMan) in order to annotate the gene descriptions, associate them with gene ontology terms and to assign each to pathways. A large number of genes related to ribosomes, the spliceosome and oxidative phosphorylation were found to be expressed in the early and middle stages. Plant hormone, galactose metabolism and carbon fixation related genes showed a significant increase in expression at the middle stage, whereas genes for defense against disease or response to stress as well as genes for starch/sucrose metabolism were strongly expressed during the later stages of endosperm development. Interestingly, most metabolic pathways were down-regulated between 3 and 10 DAP except for those involved in the accumulation of material, such as starch/sucrose and protein metabolism. We also identified the expression of 1,118 putative transcription factor genes in endosperm development. The RNA-Seq results provide further systematic understanding of rice endosperm development at a fine scale and a foundation for future studies.

An EST-based analysis identifies new genes and reveals distinctive gene expression features of Coffea arabica and Coffea canephora

BACKGROUND

Coffee is one of the world's most important crops; it is consumed worldwide and plays a significant role in the economy of producing countries. Coffea arabica and C. canephora are responsible for 70 and 30% of commercial production, respectively. C. arabica is an allotetraploid from a recent hybridization of the diploid species, C. canephora and C. eugenioides. C. arabica has lower genetic diversity and results in a higher quality beverage than C. canephora. Research initiatives have been launched to produce genomic and transcriptomic data about Coffea spp. as a strategy to improve breeding efficiency.

RESULTS

Assembling the expressed sequence tags (ESTs) of C. arabica and C. canephora produced by the Brazilian Coffee Genome Project and the Nestlé-Cornell Consortium revealed 32,007 clusters of C. arabica and 16,665 clusters of C. canephora. We detected different GC3 profiles between these species that are related to their genome structure and mating system. BLAST analysis revealed similarities between coffee and grape (Vitis vinifera) genes. Using KA/KS analysis, we identified coffee genes under purifying and positive selection. Protein domain and gene ontology analyses suggested differences between Coffea spp. data, mainly in relation to complex sugar synthases and nucleotide binding proteins. OrthoMCL was used to identify specific and prevalent coffee protein families when compared to five other plant species. Among the interesting families annotated are new cystatins, glycine-rich proteins and RALF-like peptides. Hierarchical clustering was used to independently group C. arabica and C. canephora expression clusters according to expression data extracted from EST libraries, resulting in the identification of differentially expressed genes. Based on these results, we emphasize gene annotation and discuss plant defenses, abiotic stress and cup quality-related functional categories.

CONCLUSION

We present the first comprehensive genome-wide transcript profile study of C. arabica and C. canephora, which can be freely assessed by the scientific community at http://www.lge.ibi.unicamp.br/coffea. Our data reveal the presence of species-specific/prevalent genes in coffee that may help to explain particular characteristics of these two crops. The identification of differentially expressed transcripts offers a starting point for the correlation between gene expression profiles and Coffea spp. developmental traits, providing valuable insights for coffee breeding and biotechnology, especially concerning sugar metabolism and stress tolerance.

The Zea mays mutants opaque-2 and opaque-7 disclose extensive changes in endosperm metabolism as revealed by protein, amino acid, and transcriptome-wide analyses

Background

The changes in storage reserve accumulation during maize (Zea mays L.) grain maturation are well established. However, the key molecular determinants controlling carbon flux to the grain and the partitioning of carbon to starch and protein are more elusive. The Opaque-2 (O2) gene, one of the best-characterized plant transcription factors, is a good example of the integration of carbohydrate, amino acid and storage protein metabolisms in maize endosperm development. Evidence also indicates that the Opaque-7 (O7) gene plays a role in affecting endosperm metabolism. The focus of this study was to assess the changes induced by the o2 and o7 mutations on maize endosperm metabolism by evaluating protein and amino acid composition and by transcriptome profiling, in order to investigate the functional interplay between these two genes in single and double mutants.

Results

We show that the overall amino acid composition of the mutants analyzed appeared similar. Each mutant had a high Lys and reduced Glx and Leu content with respect to wild type. Gene expression profiling, based on a unigene set composed of 7,250 ESTs, allowed us to identify a series of mutant-related down (17.1%) and up-regulated (3.2%) transcripts. Several differentially expressed ESTs homologous to genes encoding enzymes involved in amino acid synthesis, carbon metabolism (TCA cycle and glycolysis), in storage protein and starch metabolism, in gene transcription and translation processes, in signal transduction, and in protein, fatty acid, and lipid synthesis were identified. Our analyses demonstrate that the mutants investigated are pleiotropic and play a critical role in several endosperm-related metabolic processes. Pleiotropic effects were less evident in the o7 mutant, but severe in the o2 and o2o7 backgrounds, with large changes in gene expression patterns, affecting a broad range of kernel-expressed genes.

Conclusion

Although, by necessity, this paper is descriptive and more work is required to define gene functions and dissect the complex regulation of gene expression, the genes isolated and characterized to date give us an intriguing insight into the mechanisms underlying endosperm metabolism.

Monitoring the expression of maize genes in developing kernels under drought stress using oligo-microarray

Preharvest aflatoxin contamination of grain grown on the US southeastern Coast Plain is provoked and aggravated by abiotic stress. The primary abiotic stress is drought along with high temperatures. The objectives of the present study were to monitor gene expression in developing kernels in response to drought stress and to identify drought-responsive genes for possible use in germplasm assessment. The maize breeding line Tex6 was used, and gene expression profiles were analyzed in developing kernels under drought stress verses well-watered conditions at the stages of 25, 30, 35, 40, 45 d after pollination (DAP) using the 70 mer maize oligo-arrays. A total of 9 573 positive array spots were detected with unique gene IDs, and 7 988 were common in both stressed and well-watered samples. Expression patterns of some genes in several stress response-associated pathways, including abscisic acid, jasmonic acid and phenylalanine ammonia-lyase, were examined, and these specific genes were responsive to drought stress positively. Real-time quantitative polymerase chain reaction validated microarray expression data. The comparison between Tex6 and B73 revealed that there were significant differences in specific gene expression, patterns and levels. Several defense-related genes had been downregulated, even though some defense-related or drought responsive genes were upregulated at the later stages.

Identification of differentially expressed genes at two key endosperm development stages using two maize inbreds with large and small grain and integration with detected QTL for grain weight

Maize endosperm accounts for more than 80% of the grain weight. Cell division and grain filling are the two key stages for endosperm development. Previous studies showed that gene expression during differential stages in endosperm development is greatly different. However, information on systematic identification and characterization of the differentially expressed genes between the two stages are limited. In this study, suppression subtractive hybridization (SSH) was used to generate four subtracted cDNA libraries for the two stages using two maize inbreds with large and small grain. Totally, 4,784 differentially expressed sequence tags (ESTs) were sequenced and 902 were non-redundant, which consisted of 344 unique ESTs. Among them 192 had high sequence similarity to the GenBank entries and represent diverse of functional categories, such as metabolism, cell growth/division, transcription, signal transduction, protein destination/storage, protein synthesis and others. The expression patterns of 75.7% SSH-derived cDNAs were confirmed by reverse Northern blot and semi-quantitative reverse transcription polymerase chain reaction, and exhibited the similar results (75.0%). Genes differentially expressed between two key stages for the two inbreds were involved in diverse physiological process pathway, which might be responsible for the formation of grain weight. 43.8% (70 of the 160 unique ESTs) of the identified ESTs were assigned to 39 chromosome bins distributed over all ten maize chromosomes. Eleven ESTs were found to co-localize with previous detected QTLs for grain weight, which might be considered as the candidate genes of grain weight for further study.

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.

Epistatic interactions between Opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize

Association genetics is a powerful method to track gene polymorphisms responsible for phenotypic variation, since it takes advantage of existing collections and historical recombination to study the correlation between large genetic diversity and phenotypic variation. We used a collection of 375 maize (Zea mays ssp. mays) inbred lines representative of tropical, American, and European diversity, previously characterized for genome-wide neutral markers and population structure, to investigate the roles of two functionally related candidate genes, Opaque2 and CyPPDK1, on kernel quality traits. Opaque2 encodes a basic leucine zipper transcriptional activator specifically expressed during endosperm development that controls the transcription of many target genes, including CyPPDK1, which encodes a cytosolic pyruvate orthophosphate dikinase. Using statistical models that correct for population structure and individual kinship, Opaque2 polymorphism was found to be strongly associated with variation of the essential amino acid lysine. This effect could be due to the direct role of Opaque2 on either zein transcription, zeins being major storage proteins devoid of lysine, or lysine degradation through the activation of lysine ketoglutarate reductase. Moreover, we found that a polymorphism in the Opaque2 coding sequence and several polymorphisms in the CyPPDK1 promoter nonadditively interact to modify both lysine content and the protein-versus-starch balance, thus revealing the role in quantitative variation in plants of epistatic interactions between a transcriptional activator and one of its target genes.

Structural and kinetic characterization of a maize aldose reductase

The aldo-keto reductases (AKRs) are classified as oxidoreductases and are found in organisms from prokaryotes to eukaryotes. The AKR superfamily consists of more than 120 proteins that are distributed throughout 14 families. Very few plant AKRs have been characterized and their biological functions remain largely unknown. Previous work suggests that AKRs may participate in stress tolerance by detoxifying reactive aldehyde species. In maize endosperm, the presence of an aldose reductase (AR; EC 1.1.1.21) enzyme has also been hypothesized based on the extensive metabolism of sorbitol. This manuscript identifies and characterizes an AKR from maize (Zea mays L.) with features of an AR. The cDNA clone, classified as AKR4C7, was expressed as a recombinant His-tag fusion protein in Escherichia coli. The product was purified by immobilized metal affinity chromatography followed by anion exchange chromatography. Circular dichroism spectrometry and SAXS analysis indicated that the AKR4C7 protein was stable, remained folded throughout the purification process, and formed monomers of a globular shape, with a molecular envelope similar to human AR. Maize AKR4C7 could utilize dl-glyceraldehyde and some pentoses as substrates. Although the maize AKR4C7 was able to convert sorbitol to glucose, the low affinity for this substrate indicated that AKR4C7 was probably a minimal contributor to sorbitol metabolism in maize seeds. Polyclonal antisera raised against AKR4C7 recognized at least three AR-like polypeptides in maize kernels, consistent with the presence of a small gene family. Diverse functions may have evolved for maize AKRs in association with specific physiological requirements of kernel development.

Dynamic proteomic analysis reveals a switch between central carbon metabolism and alcoholic fermentation in rice filling grains

Accumulation of reserve materials in filling grains involves the coordination of different metabolic and cellular processes, and understanding the molecular mechanisms underlying the interconnections remains a major challenge for proteomics. Rice (Oryza sativa) is an excellent model for studying grain filling because of its importance as a staple food and the available genome sequence database. Our observations showed that embryo differentiation and endosperm cellularization in developing rice seeds were completed approximately 6 d after flowering (DAF); thereafter, the immature seeds mainly underwent cell enlargement and reached the size of mature seeds at 12 DAF. Grain filling began at 6 DAF and lasted until 20 DAF. Dynamic proteomic analyses revealed 396 protein spots differentially expressed throughout eight sequential developmental stages from 6 to 20 DAF and determined 345 identities. These proteins were involved in different cellular and metabolic processes with a prominently functional skew toward metabolism (45%) and protein synthesis/destination (20%). Expression analyses of protein groups associated with different functional categories/subcategories showed that substantially up-regulated proteins were involved in starch synthesis and alcoholic fermentation, whereas the down-regulated proteins in the process were involved in central carbon metabolism and most of the other functional categories/subcategories such as cell growth/division, protein synthesis, proteolysis, and signal transduction. The coordinated changes were consistent with the transition from cell growth and differentiation to starch synthesis and clearly indicated that a switch from central carbon metabolism to alcoholic fermentation may be important for starch synthesis and accumulation in the developmental process.

Sugar levels modulate sorbitol dehydrogenase expression in maize

The first step in sucrose use by maize kernels produces fructose, regardless of whether the initial reaction is catalyzed by an invertase or the reversible sucrose synthase. This fructose can enter subsequent metabolism via hexokinase, or in maize kernels, by a sorbitol dehydrogenase that reversibly converts fructose + NADH to sorbitol + NAD. High levels of SDH activity suggest that kernels synthesize considerable amounts of sorbitol, but the molecular mechanism and functional role for this process have remained equivocal. To gain insights on the role of sorbitol synthesis in maize endosperm we cloned and characterized the transcriptional control of the maize sorbitol dehydrogenase (Sdh1) gene. Data indicated that Sdh1 was essentially kernel- and endosperm-specific, with maximal expression at both the mRNA and enzyme activity levels during early kernel development. Expression was elevated in high-sugar mutants (sugary1, shrunken2), also by sugar injections, and was more pronounced when transfected tissues were incubated at low oxygen concentrations. Control of Sdh1 expression in our transient assays was largely dependent on the first intron of Sdh1. We speculate that SDH activity may represent an adaptation to the high-sugar/low-oxygen environment of the endosperm. Under these conditions, the NADH-dependent reduction of fructose to sorbitol would regenerate NAD[+], thus contributing to the maintenance of the redox and energy status of the cell.

Organization of the prolamin gene family provides insight into the evolution of the maize genome and gene duplications in grass species

Zea mays, commonly known as corn, is perhaps the most greatly produced crop in terms of tonnage and a major food, feed, and biofuel resource. Here we analyzed its prolamin gene family, encoding the major seed storage proteins, as a model for gene evolution by syntenic alignments with sorghum and rice, two genomes that have been sequenced recently. Because a high-density gene map has been constructed for maize inbred B73, all prolamin gene copies can be identified in their chromosomal context. Alignment of respective chromosomal regions of these species via conserved genes allow us to identify the pedigree of prolamin gene copies in space and time. Its youngest and largest gene family, the alpha prolamins, arose about 22-26 million years ago (Mya) after the split of the Panicoideae (including maize, sorghum, and millet) from the Pooideae (including wheat, barley, and oats) and Oryzoideae (rice). The first dispersal of alpha prolamin gene copies occurred before the split of the progenitors of maize and sorghum about 11.9 Mya. One of the two progenitors of maize gained a new alpha zein locus, absent in the other lineage, to form a nonduplicated locus in maize after allotetraplodization about 4.8 Mya. But dispersed copies gave rise to tandem duplications through uneven expansion and gene silencing of this gene family in maize and sorghum, possibly because of maize's greater recombination and mutation rates resulting from its diploidization process. Interestingly, new gene loci in maize represent junctions of ancestral chromosome fragments and sites of new centromeres in sorghum and rice.

Nonadditive expression and parent-of-origin effects identified by microarray and allele-specific expression profiling of maize endosperm

Plant endosperm cells have a nuclear ratio of two maternal genomes to one paternal genome. This 2 to 1 dosage relationship provides a unique system for studying the additivity of gene expression levels in reciprocal hybrids. A combination of microarray profiling and allele-specific expression analysis was performed using RNA isolated from endosperm tissues of maize (Zea mays) inbred lines B73 and Mo17 and their reciprocal hybrids at two developmental stages, 13 and 19 d after pollination. The majority of genes exhibited additive expression in reciprocal hybrids based on microarray analyses. However, a substantial number of genes exhibited nonadditive expression patterns, including maternal like, paternal like, high parent like, low parent like, and expression patterns outside the range of the parental inbreds. The frequency of hybrid expression patterns outside of the parental range in maize endosperm tissue is much higher than that observed for vegetative tissues. For a set of 90 genes, allele-specific expression assays were employed to monitor allelic bias and regulatory variation. Eight of these genes exhibited evidence for maternally or paternally biased expression at multiple stages of endosperm development and are potential examples of differential imprinting. Our data indicate that parental effects on gene expression are much stronger in endosperm than in vegetative tissues.

Developmental analysis of maize endosperm proteome suggests a pivotal role for pyruvate orthophosphate dikinase

Although the morphological steps of maize (Zea mays) endosperm development are well described, very little is known concerning the coordinated accumulation of the numerous proteins involved. Here, we present a proteomic study of maize endosperm development. The accumulation pattern of 409 proteins at seven developmental stages was examined. Hierarchical clustering analysis allowed four main developmental profiles to be recognized. Comprehensive investigation of the functions associated with clusters resulted in a consistent picture of the developmental coordination of cellular processes. Early stages, devoted to cellularization, cell division, and cell wall deposition, corresponded to maximal expression of actin, tubulins, and cell organization proteins, of respiration metabolism (glycolysis and tricarboxylic acid cycle), and of protection against reactive oxygen species. An important protein turnover, which is likely associated with the switch from growth and differentiation to storage, was also suggested from the high amount of proteases. A relative increase of abundance of the glycolytic enzymes compared to tricarboxylic acid enzymes is consistent with the recent demonstration of anoxic conditions during starch accumulation in the endosperm. The specific late-stage accumulation of the pyruvate orthophosphate dikinase may suggest a critical role of this enzyme in the starch-protein balance through inorganic pyrophosphate-dependent restriction of ADP-glucose synthesis in addition to its usually reported influence on the alanine-aromatic amino acid synthesis balance.