Stages Of Endosperm Development Research Articles

Variety-dependent accumulation of glucomannan in the starchy endosperm and aleurone cell walls of rice grains and its possible genetic basis.

Plant cell wall plays important roles in the regulation of plant growth/development and affects the quality of plant-derived food and industrial materials. On the other hand, genetic variability of cell wall structure within a plant species has not been well understood. Here we show that the endosperm cell walls, including both starchy endosperm and aleurone layer, of rice grains with various genetic backgrounds are clearly classified into two groups depending on the presence/absence of β-1,4-linked glucomannan. All-or-none distribution of the glucomannan accumulation among rice varieties is very different from the varietal differences of arabinoxylan content in wheat and barley, which showed continuous distributions. Immunoelectron microscopic observation suggested that the glucomannan was synthesized in the early stage of endosperm development, but the synthesis was down-regulated during the secondary thickening process associated with the differentiation of aleurone layer. Significant amount of glucomannan in the cell walls of the glucomannan-positive varieties, i.e., 10% or more of the starchy endosperm cell walls, and its close association with the cellulose microfibril suggested possible effects on the physicochemical/biochemical properties of these cell walls. Comparative genomic analysis indicated the presence of striking differences between OsCslA12 genes of glucomannan-positive and negative rice varieties, Kitaake and Nipponbare, which seems to explain the all-or-none glucomannan cell wall trait in the rice varieties. Identification of the gene responsible for the glucomannan accumulation could lead the way to clarify the effect of the accumulation of glucomannan on the agronomic traits of rice by using genetic approaches.

Influence of high-molecular-weight glutenin subunit deletions at the Glu-A1 and Glu-D1 loci on protein body development, protein components and dough properties of wheat (Triticum aestivum L.)

High-molecular-weight glutenin subunits (HMW-GSs) play a critical role in determining the viscoelastic properties of wheat. As the organelle where proteins are stored, the development of protein bodies (PBs) reflects the status of protein synthesis and also affects grain quality to a great extent. In this study, with special materials of four near-isogenic lines in a Yangmai 18 background we created, the effects of Glu-A1 and Glu-D1 loci deletions on the development and morphological properties of the protein body, protein components and dough properties were investigated. The results showed that the deletion of the HMW-GS subunit delayed the development process of the PBs, and slowed the increases of volume and area of PBs from 10 days after anthesis (DAA) onwards. In contrast, the areas of PBs at 25 DAA, the middle or late stage of endosperm development, showed no distinguishable differences among the four lines. Compared to the wild type and single null type in Glu-A1, the ratios of HMW-GSs to low-molecular-weight glutenin subunits (LMW-GSs), glutenin macropolymer (GMP) content, mixograph parameters as well as extension parameters decreased in the single null type in Glu-D1 and double null type in Glu-A1 and Glu-D1, while the ratios of gliadins (Gli)/glutenins (Glu) in those types increased. The absence of Glu-D1 subunits decreased both dough strength and extensibility significantly compared to the Glu-A1 deletion type. These results provide a detailed description of the effect of HMW-GS deletion on PBs, protein traits and dough properties, and contribute to the utilization of Glu-D1 deletion germplasm in weak gluten wheat improvement for use in cookies, cakes and southern steamed bread in China and liquor processing.

Changes in fine structure of amylopectin and internal structures of starch granules in developing endosperms and culms caused by starch branching enzyme mutations of japonica rice.

BEIIb plays a specific role in determining the structure of amylopectin in rice endosperm, whereas BEIIa plays the similar role in the culm where BEIIb is absent. Cereals have three types of starch branching enzymes (BEs), BEI, BEIIa, and BEIIb. It is widely known that BEIIb is specifically expressed in the endosperm and plays a distinct role in the structure of amylopectin because in its absence the amylopectin type changes to the amylose-extender-type (ae-type) or B-type from the wild-type or A-type and this causes the starch crystalline allomorph to the B-type from the wild-type A-type. This study aimed to clarify the role of BEIIa in the culm where BEIIb is not expressed, by using a be2a mutant in comparison with results with be2b and be1 mutants. The results showed that the amylopectin structure exhibited the B-type in the be2a culm compared with the A-type in the wild-type culm. The starch granules from the be2a culm also showed the B-type like allomorph when examined by X-ray diffraction analysis and optical sum frequency generation spectroscopy. Both amylopectin chain-length profile and starch crystalline properties were found to be the A-type at the very early stage of endosperm development at 4-6days after pollination (DAP) even in the be2b mutant. All these results support a view that in the culm as well as in the endosperm at 4-6 DAP, BEIIa can play the role of BEIIb which has been well documented in maturing endosperm. The possible mechanism as to how BEIIa can play its role is discussed.

CTP synthase is essential for early endosperm development by regulating nuclei spacing.

SummaryCereal grain endosperms are an important source of human nutrition. Nuclear division in early endosperm development plays a major role in determining seed size; however, this development is not well understood. We identified the rice mutant endospermless 2 (enl2), which shows defects in the early stages of endosperm development. These phenotypes arise from mutations in OsCTPS1 that encodes a cytidine triphosphate synthase (CTPS). Both wild‐type and mutant endosperms were normal at 8 h after pollination (HAP). In contrast, at 24 HAP, enl2 endosperm had approximately 10–16 clumped nuclei while wild‐type nuclei had increased in number and migrated to the endosperm periphery. Staining of microtubules in endosperm at 24 HAP revealed that wild‐type nuclei were evenly distributed by microtubules while the enl2‐2 nuclei were tightly packed due to their reduction in microtubule association. In addition, OsCTPS1 interacts with tubulins; thus, these observations suggest that OsCTPS1 may be involved in microtubule formation. OsCTPS1 transiently formed macromolecular structures in the endosperm during early developmental stages, further supporting the idea that OsCTPS1 may function as a structural component during endosperm development. Finally, overexpression of OsCTPS1 increased seed weight by promoting endosperm nuclear division, suggesting that this trait could be used to increase grain yield.

Starch content and activities of starch biosynthetic enzymes in wheat, rice and millets

The aim of the study was to correlate the accumulation of amylose and total starch with activities of starch biosynthetic enzymes (AGPase, SSS and GBSS) in cereals like rice (PB-2, PD-19), wheat (UP-262, PBW-343) and millets (Barnyard millet, Finger millet and Foxtail millet) at different stages of endosperm development. ADP-glucose pyrophosphorylase (AGPase) activity was found higher in all varieties during developmental stage followed by granule bound starch synthase (GBSS) and soluble starch synthase (SSS). Accumulation of amylose was positively influenced by activity of AGPase and GBSS. Total starch accumulation was positively correlated with activity of AGPase and SSS. In all these varieties rate of enzymes activity was found higher at S3 stage. Higher ratio of AGPase and GBSS contributes lesser resistant starch, amylose content and am/ap ratio. Ratio of AGPase and GBSS can be used as an index for amylose content when it was studied in large sample size of similar varieties.

MADS78 and MADS79 Are Essential Regulators of Early Seed Development in Rice.

MADS box transcription factors (TFs) are subdivided into type I and II based on phylogenetic analysis. The type II TFs regulate floral organ identity and flowering time, but type I TFs are relatively less characterized. Here, we report the functional characterization of two type I MADS box TFs in rice (Oryza sativa), MADS78 and MADS79 Transcript abundance of both these genes in developing seed peaked at 48 h after fertilization and was suppressed by 96 h after fertilization, corresponding to syncytial and cellularized stages of endosperm development, respectively. Seeds overexpressing MADS78 and MADS 79 exhibited delayed endosperm cellularization, while CRISPR-Cas9-mediated single knockout mutants showed precocious endosperm cellularization. MADS78 and MADS 79 were indispensable for seed development, as a double knockout mutant failed to make viable seeds. Both MADS78 and 79 interacted with MADS89, another type I MADS box, which enhances nuclear localization. The expression analysis of Fie1, a rice FERTILIZATION-INDEPENDENT SEED-POLYCOMB REPRESSOR COMPLEX2 component, in MADS78 and 79 mutants and vice versa established an antithetical relation, suggesting that Fie1 could be involved in negative regulation of MADS78 and MADS 79 Misregulation of MADS78 and MADS 79 perturbed auxin homeostasis and carbon metabolism, as evident by misregulation of genes involved in auxin transport and signaling as well as starch biosynthesis genes causing structural abnormalities in starch granules at maturity. Collectively, we show that MADS78 and MADS 79 are essential regulators of early seed developmental transition and impact both seed size and quality in rice.

An image processing method for investigating the morphology of cereal endosperm cells

ABSTRACTIn cereal seeds, the number, morphology and development of endosperm cells are closely related to grain quality, weight and yield. Endosperm cells differ morphologically in different regions of the seed. Nevertheless, it is important to be able to analyze the morphology of cereal endosperm cells. We established an image processing method to enhance the outlines of endosperm cells. The endosperm cell wall was traced precisely using the “pen tool” in Photoshop software (PS). The tracing was defined as the “work path” and was highlighted using the PS “brush tool.” Images of mature rice, maize and wheat endosperm sections stained with different methods were analyzed using this method. Combined with the whole sections of mature and developing cereal kernels, the processed image exhibited clearly the morphology of endosperm cells in any region of endosperm and at any stage of endosperm development. The processed image was more accurate and efficient for analyzing morphological characteristics than the unprocessed image.

Post pollination events in a self-sterile hexaploid monoembryonic lineage of Handroanthus serratifolius (Vahl) S. Grose (Bignoniaceae)

ABSTRACT In diploid late-acting self-incompatible (LSI) species of Bignoniaceae, self-pollinated pistils show a marked delay in ovule penetration by pollen tubes, followed by delayed, but otherwise normal, initial stages of endosperm development and subsequent pistil abscission. Most polyploid species of Bignoniaceae are apomictic and set selfed fruits with viable polyembryonic seeds carrying adventitious embryos. Handroanthus serratifolius is a polyploid species with sporophytic apomictic individuals that produce polyembryonic seeds after either self- or cross-pollination, and other LSI individuals that produce exclusively monoembryionic seeds. Our aim was to determine whether self-sterility events in the latter are similar to those observed in other diploid LSI species of Bignoniaceae. No conspicuous differences in the progress of ovule penetration, fertilization and early stages of endosperm development were observed between selfed and crossed pistils, which contrasts with the marked delay observed in these events after self-pollination in diploid LSI species of the family. Adventitious embryo precursor cells (AEPs) were observed in some ovules, but they apparently do not develop into embryos. We believe that the AEPs in these plants represent a ‘potential’ for sporophytic apomixis, which would explain the similar behaviour of post-pollination events in selfed vs. crossed pistils despite their inability to set fruits after selfing.

Detailed observation on expression dynamics of Polycomb group genes during rice early endosperm development in subspecies hybridization reveals their characteristics of parent-of-origin genes

BackgroundParent-of-origin gene expression and its role in seed development have drown a great attention in recent years. Genome-wide analysis has identified hundreds of candidate imprinted genes, a major type of parent-of-origin genes, in rice hybrid endosperms at the stage of 5 days after pollination (dap). However, the expression of these genes in early endosperm have been never confirmed due to technique limitations and the behavior of the imprinted genes in different rice hybridizations are still largely unknown.ResultsHere, based on our elaborate technique established previously, the expression patterns of PcG genes in the early stages of endosperm development (within 3 dap), were comprehensively analyzed. We revealed that the free nucleus stage of endosperm development is critical for parent-of-origin gene analysis. The expression of the imprinted genes are highly dynamic, likely corresponding to the critical developmental events during this period. Hybridizations between Oryza sativa japonica and indica showed that the expression patterns of the same imprinted gene could be varied by crossing with different parental cultivars, indicative of their parent-dependent character. There are strong alleles that often showed predominant expression over other alleles regardless of the parental origin, which provides a possible explanation for the cultivar-dependent predominant phenotype in crop hybridizations. In addition, we found that the transcripts of the same gene behave differently, with imprinting or non-imprinting patterns, suggesting the existence of not only imprinted and non-imprinted genes but also imprinted or non-imprinted transcripts, which reveals new aspects of the genomic imprinting.ConclusionsThese findings on the characters of parent-of-origin genes shed light on the understanding the real role of gene imprinting in endosperm development.

Transcriptome Dynamics of Double Recessive Mutant, o2o2o16o16, Reveals the Transcriptional Mechanisms in the Increase of Its Lysine and Tryptophan Content in Maize

In maize, pyramiding of o2 and o16 alleles can greatly improve the nutritional quality of grains. To dissect its molecular mechanism, we created a double recessive mutant line, o2o2o16o16, by introgression of the o2 and o16 alleles into the wild-type maize inbred line, by molecular marker-assisted backcross selection. The kernels (18 day after pollination (DAP), 28 DAP, and 38 DAP) of the o2o2o16o16 mutant and its parent lines were subject to RNA sequencing (RNA-Seq). The RNA-Seq analysis revealed that 59 differentially expressed genes (DEGs) were involved in lysine metabolism and 43 DEGs were involved in tryptophan metabolism. Among them, the genes encoding AK, ASADH, and Dap-F in the lysine synthesis pathway were upregulated at different stages of endosperm development, promoting the synthesis of lysine. Meanwhile, the genes encoding LKR/SDH and L-PO in the lysine degradation pathway were downregulated, inhibiting the degradation of lysine. Moreover, the genes encoding TAA and YUC in the tryptophan metabolic pathway were downregulated, restraining the degradation of tryptophan. Thus, pyramiding o2 and o16 alleles could increase the lysine and tryptophan content in maize. These above results would help to uncover the molecular mechanisms involved in the increase in lysine and the tryptophan content, through the introgression of o2 and o16 alleles into the wild-type maize.

Cytosolic calcium localization and dynamics during early endosperm development in the genus Agave (Asparagales, Asparagaceae).

Calcium is a secondary messenger that regulates and coordinates the cellular responses to environmental cues. Despite calcium being a key player during fertilization in plants, little is known about its role during the development of the endosperm. For this reason, the distribution, abundance, and dynamics of cytosolic calcium during the first stages of endosperm development of Agave tequilana and Agave salmiana were analyzed. Cytosolic calcium and actin filaments detected in the embryo sacs of Agave tequilana and A. salmiana revealed that they play an important role during the division and nuclear migration of the endosperm. After fertilization, a relatively high concentration of cytosolic calcium was located in the primary nucleus of the endosperm, as well as around migrating nuclei during the development of the endosperm. Cytosolic calcium participates actively during the first mitosis of the endosperm mother cell and interacts with the actin filaments that generate the motor forces during the migration of the nuclei through the large cytoplasm of the central cell.

A Connection between Lysine and Serotonin Metabolism in Rice Endosperm.

Cereal endosperms produce a vast array of metabolites, including the essential amino acid lysine (Lys). Enhanced accumulation of Lys has been achieved via metabolic engineering in cereals, but the potential connection between metabolic engineering and Lys fortification is unclear. In mature seeds of engineered High Free Lysine (HFL) rice (Oryza sativa), the endosperm takes on a characteristic dark-brown appearance. In this study, we use an integrated metabolomic and transcriptomic approach combined with functional validation to elucidate the key metabolites responsible for the dark-brown phenotype. Importantly, we found that serotonin biosynthesis was elevated dramatically and closely linked with dark-brown endosperm color in HFL rice. A functional connection between serotonin and endosperm color was confirmed via overexpression of TDC3, a key enzyme of serotonin biosynthesis. Furthermore, we show that both the jasmonate signaling pathway and TDC expression were strongly induced in the late stage of endosperm development of HFL rice, coinciding with serotonin accumulation and dark-brown pigmentation. We propose a model for the metabolic connection between Lys and serotonin metabolism in which elevated 2-aminoadipate from Lys catabolism may play a key role in the connection between the jasmonate signaling pathway, serotonin accumulation, and the brown phenotype in rice endosperm. Our data provide a deeper understanding of amino acid metabolism in rice. In addition, the finding that both Lys and serotonin accumulate in HFL rice grains should promote efforts to create a nutritionally favorable crop.

Laser-Capture Microdissection of Maize Kernel Compartments for RNA-Seq-Based Expression Analysis.

Laser-capture microdissection (LCM) enables isolation of single cells or groups of cells for a variety of downstream applications including transcriptome profiling. Recently, this methodology has found a more widespread use particularly with the advent of next-generation sequencing techniques that enable deep profiling of the limited amounts of RNA obtained from fixed or frozen sections. When used with fixed tissues, a major experimental challenge is to balance the tissue integrity needed for microscopic visualization of the cell types of interest with that of the RNA quality necessary for deep profiling. Complex biological structures such as seeds or kernels pose an especially difficult case in this context as in many instances the key internal structures such as the embryo and the endosperm are relatively inaccessible. Here, we present an optimized LCM protocol for maize kernel that has been developed specifically to enable profiling of the early stages of endosperm development using RNA-Seq.

The isolation of early nuclear endosperm of Oryza sativa to facilitate gene expression analysis and screening imprinted genes.

BackgroundSince the quality and yield of rice production depends on endosperm development, previous studies have focused on the molecular mechanism that regulates this developmental process. Recently, how this process is epigenetically regulated has become an important topic. However, the gene expression analysis and screening imprinted genes during early endosperm development remain challenging since the isolation of early endosperm has not been possible. Here, we report a procedure for the isolation of endosperm at 24 or 48 HAP (hours after pollination) during the free nuclear stage of endosperm development.ResultsThis technique allows for rapid and convenient collection of pure free nuclear endosperm. Early endosperm RNA can then be extracted from the isolated endosperm cells using dynabeads. Our results showed that the quality of RNA is satisfactory for gene expression analysis and screening the parental-of-origin specific genes in early endosperm.ConclusionsThus, we offer a reliable method to overcome one of the major obstacles in the investigation of the molecular mechanisms of early endosperm development. Our approach can be used for accurate gene expression analysis and screening of imprinted genes, and facilitates the confirmation of endosperm-specific gene expression at the very early stages of endosperm development. This method could also be used in other species to collect early free nuclear endosperm.

A maize ADP-ribosylation factor ZmArf2 increases organ and seed size by promoting cell expansion in Arabidopsis.

ADP-ribosylation factors (ARFs) are small GTP-binding proteins that regulate a wide variety of cell functions. Previously, we isolated a new ARF, ZmArf2, from maize (Zea mays). Sequence and expression characteristics indicated that ZmArf2 might play a critical role in the early stages of endosperm development. In this study, we investigated ZmArf2 function by analysis of its GTP-binding activity and subcellular localization. We also over-expressed ZmArf2 in Arabidopsis and measured organ and cell size and counted cell numbers. The expression levels of five organ size-associated genes were also determined in 35S::ZmArf2 transgenic and wild-type plants. Results showed that the recombinant ZmArf2 protein purified from Escherichia coli exhibited GTP-binding activity. Subcellular localization revealed that ZmArf2 was localized in the cytoplasm and plasma membrane. ZmArf2 over-expression in Arabidopsis showed that 35S::ZmArf2 transgenic plants were taller and had larger leaves and seeds compared to wild-type plants, which resulted from cell expansions, not an increase in cell numbers. In addition, three cell expansion-related genes, AtEXP3, AtEXP5 and AtEXP10, were upregulated in 35S::ZmArf2 transgenic lines, while the expression levels of AtGIF1 and AtGRF5, were unchanged. Collectively, our studies suggest that ZmArf2 has an active GTP-binding function, and plays a crucial role in growth and development in Arabidopsis through cell expansion mediated by cell expansion genes.

Expression, regulation and activity of a B2-type cyclin in mitotic and endoreduplicating maize endosperm.

Cyclin-dependent kinases, the master regulators of the eukaryotic cell cycle, are complexes comprised of a catalytic serine/threonine protein kinase and an essential regulatory cyclin. The maize genome encodes over 50 cyclins grouped in different types, but they have been little investigated. We characterized a type B2 cyclin (CYCB2;2) during maize endosperm development, which comprises a cell proliferation phase based on the standard mitotic cell cycle, followed by an endoreduplication phase in which DNA replication is reiterated in the absence of mitosis or cytokinesis. CYCB2;2 RNA was present throughout the period of endosperm development studied, but its level declined as the endosperm transitioned from a mitotic to an endoreduplication cell cycle. However, the level of CYCB2;2 protein remained relatively constant during both stages of endosperm development. CYCB2;2 was recalcitrant to degradation by the 26S proteasome in endoreduplicating endosperm extracts, which could explain its sustained accumulation during endosperm development. In addition, although CYCB2;2 was generally localized to the nucleus of endosperm cells, a lower molecular weight form of the protein accumulated specifically in the cytosol of endoreduplicating endosperm cells. In dividing cells, CYCB2;2 appeared to be localized to the phragmoplast and may be involved in cytokinesis and cell wall formation. Kinase activity was associated with CYCB2;2 in mitotic endosperm, but was absent or greatly reduced in immature ear and endoreduplicating endosperm. CYCB2;2-associated kinase phosphorylated maize E2F1 and the “pocket” domains of RBR1 and RBR3. CYCB2;2 interacted with both maize CDKA;1 and CDKA;3 in insect cells. These results suggest CYCB2;2 functions primarily during the mitotic cell cycle, and they are discussed in the context of the roles of cyclins, CDKs and proteasome activity in the regulation of the cell cycle during endosperm development.

Development of endosperm cells and starch granules in common wheat

The objective of the present study was to understand the developmental regularity of wheat endosperm cells at different Days After Pollination (DAP) using microscopic and histochemical methods. Resin semi-thin sections of the endosperm and the enzymatically dissociated Starchy Endosperm Cells (SECs) were observed under a light microscope. The results showed that: (1) SECs were irregular-shaped and had two types of starch granules: large oval-shaped A-type starch granules and small spherical B-type starch granules. (2) The growth shape of SECs was referred to as S-curve and the fastest cell growth period was at 16–24 DAP. (3) The largest increase and growth of A-type starch granules were mainly at 4–16 DAP. B-type starch granules increased rapidly after 16 DAP and made up over 90% of the total starch granules in SEC during the late stage of endosperm development. (4) The nuclei of SEC deformed and degenerated during the middle and late stages of endosperm development and eventually disappeared. However, starch granules still increased and grew after the cell nuclei had degenerated. The investigations showed the development regularity of starch endosperm cells and starch granules, thereby improving the understanding of wheat endosperm development.

RNA-Seq analysis of laser-capture microdissected cells of the developing central starchy endosperm of maize.

Endosperm is a product of double fertilization, and provides nutrients and signals to the embryo during seed development in flowering plants. Early stages of endosperm development are critical for the development of its storage capacity through synthesis and accumulation of starch and storage proteins. Here we report on the isolation and sequencing of mRNAs from the central portion of the starchy endosperm of Zea mays (maize) B73 at 6 days after pollination. We detected a high level of correlation among the four biological replicates of RNAs isolated using laser-capture microdissection of the cell type. Because the assayed developmental stage precedes the synthesis and accumulation of the major storage proteins and starch in the endosperm, our dataset likely include mRNAs for genes that are involved in control and establishment of these storage programs. The mRNA-Seq data has been deposited in Gene Expression Omnibus (accession number GSE58504).

Dissection of two major components of the post‐zygotic hybridization barrier in rice endosperm

A post-zygotic hybridization barrier is often observed in the endosperm of seeds produced by interspecific or interploidy crosses. In Arabidopsis thaliana, for example, hybrid endosperm from both types of cross shows altered timing of cellularization and an altered rate of nuclear divisions. Therefore, it has been proposed that interspecific and interploidy crosses share common molecular mechanisms for establishment of an effective species barrier. However, these two types of hybridization barrier may be initiated by different intrinsic cues: the interspecific cross barrier arises after hybridization of genomes with differences in DNA sequences, while the interploidy cross barrier arises after hybridization of genomes with the same DNA sequences but differences in ploidy levels. In this study, we performed interploidy crosses to identify components of the post-hybridization barrier in the endosperm of rice. We performed an intra-cultivar cross of autotetraploid (4n) × diploid (2n) rice, and found precocious cellularization and a decreased rate of nuclear division in the syncytial endosperm. By contrast, seeds from the reciprocal cross showed delayed cellularization and an increased rate of nuclear division. This differential effect on nuclear division rates contrasts with the outcome of rice interspecific crosses, which were previously shown to have altered timing of cellularization without any change in nuclear division rates. Thus, we propose that the post-zygotic hybridization barrier in rice endosperm has two separable components, namely control of the timing of cellularization and control of the nuclear division rates in the syncytial stage of endosperm development.

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 10days 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.