Transcriptome analysis of grain-filling caryopses reveals involvement of multiple regulatory pathways in chalky grain formation in rice
© Liu et al; licensee BioMed Central Ltd. 2010
Received: 10 February 2010
Accepted: 30 December 2010
Published: 30 December 2010
Grain endosperm chalkiness of rice is a varietal characteristic that negatively affects not only the appearance and milling properties but also the cooking texture and palatability of cooked rice. However, grain chalkiness is a complex quantitative genetic trait and the molecular mechanisms underlying its formation are poorly understood.
A near-isogenic line CSSL50-1 with high chalkiness was compared with its normal parental line Asominori for grain endosperm chalkiness. Physico-biochemical analyses of ripened grains showed that, compared with Asominori, CSSL50-1 contains higher levels of amylose and 8 DP (degree of polymerization) short-chain amylopectin, but lower medium length 12 DP amylopectin. Transcriptome analysis of 15 DAF (day after flowering) caryopses of the isogenic lines identified 623 differential expressed genes (P < 0.01), among which 324 genes are up-regulated and 299 down-regulated. These genes were classified into 18 major categories, with 65.3% of them belong to six major functional groups: signal transduction, cell rescue/defense, transcription, protein degradation, carbohydrate metabolism and redox homeostasis. Detailed pathway dissection demonstrated that genes involved in sucrose and starch synthesis are up-regulated, whereas those involved in non-starch polysaccharides are down regulated. Several genes involved in oxidoreductive homeostasis were found to have higher expression levels in CSSL50-1 as well, suggesting potential roles of ROS in grain chalkiness formation.
Extensive gene expression changes were detected during rice grain chalkiness formation. Over half of these differentially expressed genes are implicated in several important categories of genes, including signal transduction, transcription, carbohydrate metabolism and redox homeostasis, suggesting that chalkiness formation involves multiple metabolic and regulatory pathways.
Endosperm chalkiness is a varietal characteristic that negatively affects not only the appearance and milling properties but also the cooking texture and palatability of cooked rice . Chalky grains have a lower density of starch granules compared to vitreous ones, and are therefore more prone to breakage during milling . In many rice-producing areas, high chalkiness is a major concern that decreases grain quality. In China, many early-season indica and japonica varieties are of high grain endosperm chalkiness and their market values are seriously affected . Therefore, one of the goals in rice breeding is to reduce chalkiness in rice varieties.
In rice grains, starch is the predominant storage substance that account for over 80% of the total dry mass. Starch in rice endosperm is composed of relatively unbranched amylose (linear α-1, 4-polyglucans) and highly branched amylopectin (α-1, 6-branched polyglucans). Recent work showed that multiple factors contribute to the formation of grain chalkiness, including starch synthesis, starch granule structure and arrangement [3–5]. For example, mutations in the Wx (waxy) gene (encoding grandule-bound sucrose synthase, or GBSS) and its regulator DULL cause low amylose content (≤ 2%) and hence whole opaque endosperm [6, 7]. The amylose-extender mutant has reduced activity of branching enzyme II (BEIIb), causing alteration in the fine structure of grain amylopectin . The flo-2 floury endosperm mutant harbors mutations affecting rice branching enzyme I (RBEI) activity . The floury endosperm-4 mutant and the sugary-1 mutant are defective in pyruvate orthophosphate dikinase (PPDK) and debranching enzymes (DBE) activity respectively [4, 10].
The formation of grain chalkiness can also be influenced by various external stresses during the grain-filling stage. Temperatures higher than 26°C, for example, could easily cause chalky appearance and a reduction in grain weight . Microscopic observation showed that, compared with the translucent portion of rice endosperm that ripened under normal temperature which were filled with densely packed and polygonal granules, the chalky portion of high temperature-ripened grains were loosely packed with elliptical-shaped starch granules containing air spaces which caused random light reflection and hence chalky appearance [5, 12–14]. These observations demonstrated that environmental stresses represent another major cause for grain chalkiness in rice. Furthermore, imaging on endosperm amyloplast development of various japonica and indica rice lines indicated that starch synthesis in the rice grain may involve complicated genetic networks [4, 15]. Previous studies have detected many major quantitative trait loci (QTLs) that may underlie chalkiness in rice [16–19], however; only few QTLs have been isolated and functionally analyzed . Thus, the molecular mechanisms underlying the formation of rice grain endosperm chalkiness still remain poorly understood.
In this study, we performed a comparative transcriptome analysis of the caryopses of a near-isogenic line CSSL50-1 (with high chalkiness) and its low chalkiness parental line Asominori. Corroborated with the phenotypic and physico-biochemical observations, our genome-wide transcription analysis supports the notion that rice grain endosperm development is controlled by a delicate, but complex genetic network. Notably, several pathways related to signal transduction, cell rescue/defense, transcription, protein degradation, carbohydrate metabolism and redox homeostasis were found to be predominant among the differentially expressed genes.
Phenotypic and physiochemical properties of Asominori and CSSL50-1 grains
CSSL50-1 is derived from the near-isogenic line CSSL50 with a small substituted segment of chromosome 8 from the original donor IR24 in the largely Asominori background . CSSL50-1 displays high chalkiness under normal field conditions whereas its parental line Asominori has normal grains. Therefore, CSSL50-1 represents an ideal genetic material with relatively stabilized genetic background suitable for exploring the molecular mechanism of chalkiness formation.
Increased grain-filling rate and enhanced activities of starch enzymes in CSSL50-1
Since grains of CSSL50-1 contain higher starch, amylose and sucrose contents compared with Asominori, we speculated that enzymes involved in starch synthesis might be more robust in CSSL50-1 than in Asominori. To confirm this, the enzymatic activities of four major enzymes involved in grain starch synthesis were measured during the first 30 days after flowering. Similar patterns were observed for SuSy in CSSL50-1 and Asominori, namely, its highest activity was detected at 10 DAF, and its lowest activity detected at 30 DAF (Figure 2C). However, SuSy activity in CSSL50-1 was higher than that in Asominori at the 15 and 20 DAF. Similarly, at 15 DAF, AGPase, SBE and DBE activities were significantly higher in CSSL50-1 compared to those in Asominori (Figure 2D-F). Additionally, enzyme activities of SuSy at 30 DAF and DBE at 5 DAF were found to be lower in CSSL50-1 than those in Asominori (Figure 2C). These results indicate that 15 DAF is a critical time point for grain filling when many enzymes involved in starch synthesis exhibit maximum activities. We therefore used RNAs extracted from 15 DAF endosperms for subsequent microarray analysis.
Transcriptome analysis of 15 DAF caryopses of CSSL50-1 and Asominori
Functional classification of differentially expressed genes associated with enriched GO terms1
A. Biological process
response to stress
carbohydrate metabolic process
regulation of transcription, DNA-dependent
regulation of RNA metabolic process
RNA biosynthetic process
primary metabolic process
B. Molecular function
hydrolase activity, hydrolyzing O-glycosyl compounds
Enhanced sucrose and starch synthesis vs. disrupted cellulose, hemicellulose and pectin metabolism in CSSL50-1
Physio-biochemical analysis of chalky rice endosperm indicated that the change in starch composition is a major difference between chalky and non-chalky rice grains (Figure 1 & Figure 2). GO analysis also showed that genes associated with carbohydrate metabolism are significantly represented among the differentially expressed transcripts. As shown in Additional file 4 (Table S2), more than 50 genes are annotated to be associated with carbohydrate metabolism. Of particular interest were several key genes that are known to be directly involved in the synthesis of starch and cell-wall related polysaccharides (eg. sucrose phosphatase, sucrose phosphate synthase, glucose-6-phosphate isomerase, glycosyltransferae family 5, starch debranching enzyme, α-amylase, cellulose synthesis and α-D-xylosidase).
In CSSL50-1, the enhancement of sucrose and starch (Figure 1I) seems to be accompanied by a metabolic disorder of cell wall related polysaccharides. First, two cellulose synthase genes were down-regulated which may reduce cellulose synthesis; Second, up-regulation of α-L-arabinofuranosidase and α-D-xylosidase and down-regulation of an xylanase inhibitor protein may promote hydrolysis of hemicellulose . These observations seem to suggest that the enhancement of sucrose and starch synthesis is at the cost of cell wall related non-storage polysaccharides in CSSL50-1. Such carbohydrate metabolism disorders may significantly contribute to the endosperm chalkiness during grain ripening.
Increased expression of redox genes and a higher level ROS homeostasis in CSSL50-1
Previous studies showed that rice grains develop chalkiness under adverse environmental conditions such as high temperatures . GO analysis also indicated significant enrichment in oxidoreductase activity in Molecular Function (Table 1 & Additional file 3, Table S1). About 40 genes fell in the category of redox homeostasis in our manual classification of differentially expressed genes (Additional file 4 Table S2). Since reactive oxygen species (ROS) are well known to be involved in various stress responses, we first measured the concentration of H2O2, a common ROS, in CSSL50-1 and Asominori. The results showed that the 15 DAF grains of CSSL50-1 contained much higher H2O2 concentration (617 umol/g fresh weight) than that in Asominori (451.3 umol/g fresh weight) (P < 0.01). Such an imbalance in ROS concentrations and its consequence may contribute to the development of chalkiness in grain endosperm at later developmental stages.
A delicate but complicated gene network may underlie chalky rice grain formation
CSSL50-1 is an ideal material for exploring the molecular basis of rice grain chalkiness
Grain endosperm chalkiness is a complex quantitative genetic trait and is controlled by multiple factors . Previous studies showed that there are as many as 42 QTLs that may contribute to the percentages of grains with chalkiness (PGWC) and degrees of endosperm chalkiness. These genes spread among 10 rice chromosomes as being located using seven different genetic populations [16, 17, 19, 27–30]. In addition to genetic factors, rice grains are also sensitive to various environmental stresses and will easily form chalky endosperms under adverse conditions. Changes in temperature and light, for example, are most effective at the milky stage of grain filling [5, 11, 31]. Meanwhile, conditions of water and nutrient in the field also play important roles on grain chalkiness [32, 33].
Despite its economic importance, only few genes have been functionally identified to be associated with endosperm chalkiness. This analysis, together with a previous transcriptome analysis on chalkiness formation under higher temperature , has identified a set of differentially expressed genes that may contribute to endosperm chalkiness. The identified candidate genes may serve as excellent starting materials for dissecting the pathway controlling rice endosperm development at the molecular level. Notably, most of these genes identified in this study belong to six major categories including cell rescue/defense, free radical clearance and redox homeostasis, signal transduction, hormone response, protein biosynthesis and degradation, and carbohydrate metabolism. Our data also showed that, similar to the effect of high temperature, expression of numerous genes was affected under this genetic background, but surprisingly, quite a few of these genes were oppositely regulated in CSSL50-1 when compared with the effect under the high temperature conditions. Thus, the use of a genetically stabilized line for endosperm chalkiness study is complementary to the previous physical stress-based studies and should provide novel information regarding the molecular mechanism for chalky endosperm formation in rice.
Enhanced starch synthesis causes imbalanced starch composition in CSSL50-1
Previously, chalky grains were found to have lower starch content . For CSSL50-1, its grains contain higher percentage of sucrose, amylose, starch, and even protein content when compared with the normal rice cultivar Asominori. Enzymes, such as SuSy, AGPase, SBE and DBE, exhibit higher activities at 15 DAF in CSSL50-1, correlating with the higher expression levels of corresponding genes that were detected in our microarray data. Since the shape and the arrangement of starch granule are closely related to endosperm chalkiness in rice [3, 4, 34], it is reasonable to conclude that a coordinated and balanced action of starch synthetic enzymes are critical to the prevention of chalky endosperm formation.
Endosperm development is a process of proper starch composition and accumulation [35, 36]. Gradual and smooth grain filling pace is required to form normal, translucent grains as seen in Asominori. CSSL50-1 has a higher grain-filling rate which may give insufficient time for long chain amylopectin to be synthesized, resulting in a relative higher percentage of short chain amylopectin (8-9 DP) when compared with Asominori. Consistent with our results, the decrease of 10-14 DP amylopectin was also observed when rice grains ripened under high temperatures . This is in line with our findings that both gene expression and the enzymatic activity of DBE are increased in CSSL50-1. Interestingly, the synthesis of non-starch polysaccharides appears to be significantly retarded because of the down regulation of related genes such as cellulose synthase and the up-regulation of genes for degrading hemicellulose and pectin. Therefore, the synthesis of cell wall related sugar might be sacrificed in CSSL50-1. In light of the importance of normal synthesis of starch and related polysaccharides, it is very likely that disorders in the enzyme activity and the expression of genes responsible for these events are among the major causes for endosperm chalkiness in CSSL50-1.
Potential roles of ROS in rice grain chalkiness formation
Reactive oxygen species (ROS) are partially reduced forms of atmospheric oxygen (O2). They typically result from the excitation of O2 to form singlet oxygen (O-) or from the transfer of one, two or three electrons to O2 to form, respectively, a superoxide radical (O2-), hydrogen peroxide (H2O2) or a hydroxyl radical (HO-) . Among them, H2O2 is one of the most stable ROS . With both reducing and oxidizing properties, H2O2 has effects on almost all organisms, and can influence the life of every single cell. On one hand, H2O2 is highly reactive and toxic, and can lead to oxidative destruction of cells; on the other hand, it acts as a signaling molecule in regulating cell growth and development, cell proliferation, cell stress response, and signal transduction . When accumulated at high enough concentrations, H2O2 can directly or indirectly oxidize enough of the cellular ascorbic acid and glutathione pool to alter the overall redox state of the cells. Such high concentrations of H2O2 can also damage a large variety of biomolecules such as lipids, proteins and nucleic acids that are essential for the activity and integrity of the cells [37, 40]. As sessile organisms, plants have evolved a high degree of developmental plasticity to optimize their growth and reproduction in response to various biotic and abiotic stresses . Under these conditions, the excessive H2O2 is efficiently scavenged by various antioxidative defense mechanisms in plant cells. The major ROS-scavenging enzymes include ascorbate peroxidase (APX), catalase (CAT), superoxide dismutase (SOD), glutathione peroxidase (GPX), monodehydroascorbate reductase (MDAR), Glutaredoxin (GLR) and peroxiredoxin (PrxR). Together with the antioxidants ascorbic acid and glutathione, these enzymes provide plant cells with highly efficient machinery for detoxifying H2O2 and other ROS [37, 40, 42].
In the present study, the expression levels of five genes involved in reactive oxygen species production and hemeostasis including superoxide dismutase (SOD), ascorbate peroxidase (APX), glutathione peroxidase (GPX), monodehydroascorbate reductase (MDAR) and peroxiredoxin (PrxR) were found to be higher in CSSL50-1 than those in Asominori (Figure 5B & Additional file 4 Table S2), suggesting that the antioxidative network in CSSL50-1 is activated. This result is consistent with the higher concentration of H2O2 in CSSL50-1. In the 15th day of grain-filling period, the concentration of H2O2 in CSSL50-1 reached 617 umol/g whereas that in Asominori is only 451.3 umol/g. The high concentration of H2O2 can provoke the defense system responsive to ROS (mainly H2O2) stress in CSSL50-1 to induce the expression of these antioxidative genes (including SOD, APX, GPX, MDAR and PrxR). Unexpectedly, another one antioxidative gene (GLR) was found to be down-regulated in CSSL50-1 compared to Asominori (Figure 5B & Additional file 4, Table S2), possibly due to more GSH to be needed for enhancing the functions of GST and Glx genes (Figure 5A).
In this study, we also found that the expression levels of GST, Glx and Trx genes are significantly higher in CSSL50-1compared to those in Asominori (Figure 5A & Additional file 4, Table S2). In contrast, LOX gene is down-regulated in CSSL50-1 (Figure 5B & Additional file 4, Table S2). GST is an antioxidative protein together with glutathione (GSH) to reduce oxidized biological macromolecule, and its expression can be strongly enhanced by abiotic and biotic stresses . Glyoxalase I can convert toxic 2-oxoaldehydes into less reactive 2-hydroxyacids using GSH as a cofactor . Thioredoxin (Trx) can reduce the oxidized proteins and peroxidative lipids . However, lipoxygenase 5 (LOX5) is one member of a family of enzymes that deoxygenate unsaturated fatty acids, thus initiating lipoperoxidation of membranes . These results suggest that the antioxidative level in CSSL50-1 is higher than that in Asominori.
Grain chalkiness involves coordinated regulation of multiple pathways
It has been known for a long time that adverse environmental conditions can easily cause chalkiness in rice grain. High-temperatures, for example, have been shown to cause changes in the expression of genes involved in starch synthesis and directly correlated with the extent of gains chalkiness [5, 45]. Drought stress, as well as sulphur deficiency which also activates antioxidation-related enzymes, can cause increased sucrose synthase activity and finally lead to the emergence of chalkiness too [33, 46–49]. The reported studies show that exterior coercive conditions can break the oxidation-reduction balance, causing the change of carbohydrate metabolization in the rice plant, leading finally to the emergence of chalkiness. For example, 1) High-temperature stress can not only cause the change in the expression quantity of antioxidation-related genes, but also lead to the change in the expression quantity of carbon metabolism-related radical/protein, finally increasing the rice chalkiness [5, 45] (Lin et al 2005; Yamakawa et al 2007). 2) Drought and coldness coercive conditions can also change the expression levels of carbon metabolization-related genes in the rice plant [50, 51] (Gorantla et al 2007; Cheng et al 2007). Drought stress can also induce increase of sucrose synthase in the rice plant, finally leading to the emergence of chalkiness [33, 46, 47], which agrees with the result of this experiment; namely, the sucrose synthase activity and seed grain filling rate in high-chalkiness CSSL50-1 are remarkably higher than that in low-chalkiness Asominori. 3) For the rice plant with outside trauma treatment, salt stress, and ray irradiation, except the activation of cell-defense related genes, the expression quantity of carbon-metabolization related genes is also changed [52–54]. 4) Sulphur-deficiency can also bring about nutrition coercion of the rice plant, causing an increase in the activity of the rice's antioxidation-related enzymes, further leading to the emergence of rice chalkiness [48, 49].
Reactive oxygen species (ROS), as represented by their most stable form H2O2 , play important roles as signaling molecules in regulating plant growth and development including cell proliferation, cell stress response, and signal transduction . H2O2 is known to be involved in biotic and abiotic stress responses. The observed drastic increase in H2O2 levels in CSSL50-1 and the differential expression of several key regulatory genes involved on ROS production and scavenge collectively suggest that ROS may play a critical role in regulating rice endosperm chalkiness. Changes in H2O2 levels may affect multiple metabolism pathways in the rice endosperm, causing chalkiness phenotypic change. Further genetic and biochemical studies should further test such a possibility.
Consistent with previous studies on the effect of adverse environmental conditions in causing chalky rice grain, our comparative transcriptome analysis of the caryopses of a near-isogenic line CSSL50-1 (with high chalkiness) and its low chalkiness parental line Asominori supports the notion that rice grain endosperm development is controlled by delicate, but complex genetic networks. Notably, several pathways related to signal transduction, cell rescue/defense, transcription, protein degradation, carbohydrate metabolism and redox homeostasis were found to be predominant among the differentially expressed genes, suggesting that formation of rice endosperm chalkiness may involve coordinated regulation of multiple pathways. Further refining of CSSL50-1 as a useful genetic material will help eventual cloning and engineering the major genes underlying the formation of rice grain chalkiness.
Plant material and growth
A japonica cultivar, Asominori, and its chromosome segment substitution line (CSSL50-1) were used in this study. CSSL50-1 is a near-isogenic line of Asominori with a substituted segment from the donor IR24 (an indica cultivar). Seventy-one F7 RILs were derived from a cross between Asominori and IR24 by single-seed descent . To produce a series of CSSLs in a largely Asominori background, 19 selected RILs were crossed and then backcrossed with Asominori, without selection, until the BC3F1 generation. Sixty-six individuals were then selected at BC3F1 on the basis of a whole genome survey (116 RFLP loci) and were denoted as CSSL1-CSSL66 . CSSL50 was observed to have high grain chalkiness characteristics. To further reduce the introgressed segment, CSSL50 was backcrossed with Asominori followed by two generations of self-pollination, and the progeny were evaluated using marker-assisted selection strategy. The homozygous line CSSL50-1 (BC4F3, Asominori/IR24//4*Asominori) was found to have high chalky grains and contains a small segment of IR24 chromosome 8 in a largely Asominori genetic background.
Four batches of seeds were sowed for both Asominori and CSSL50-1. The distance between the auricle of the flag-leaf (last leaf) and that of the penultimate leaf, or so-called Auricle distance (AD) was used as a non-destructive measurement to gauge rice flowering stage. CSSL50-1 started to flower when its AD reached 17 cm which was marked as day zero after flowering or DAF, whereas the AD for Asominori is 17.5 cm. Grain endosperms of the batches that flowered simultaneously for CSSL50-1 and Asominori were collected at 5, 10, 15, 20, 25, 30 and 35 DAF. The samples were immediately frozen in liquid nitrogen and stored at -80°C. RNA samples from 15 DAF endosperms were used for the DNA microarray analysis.
Phenotype and physico-chemical properties of CSSL50-1 and Asominori grains
Seed starch granules was imaged as described in Kang et al. (2005) . Samples for scanning electron microscopy (SEM) were pre-fixed with 3% glutaraldehyde for 3 h at room temperature, rinsed three times (15 min each) with 0.1 M sodium phosphate buffer (pH 6.8), and fixed overnight with 2% OsO4 at 4°C. The fixed samples were then washed three times (15 min each) with 0.1 M sodium phosphate buffer, dehydrated through an ethanol series, and incubated in a 1:3 (v:v) ethanol-isoamyl acetate mixture for 1 h. These samples were dried to a critical point, mounted on SEM stubs, and coated with gold. The mounted specimens were observed under SEM with an accelerating voltage of 10-20 kV. Fine structure of amylopectin was determined according to Fujita et al. (2006) .
Determination of starch, amylase and sucrose content, chalkiness and RVA profiles
Starch content, amylose content and sucrose content of each accession were determined as Fujita et al. (2007) . Percentage of grains with chalkiness (PGWC), area of chalky endosperm (ACE) and degree of endosperm chalkiness (DEC) were measured according to the method of Wan et al. (2005) . To separate chalky from vitreous grains, 100 grains per entry were assessed on a chalkiness visualizer to calculate PGWC. Twenty chalky grains were then selected at random, and the ratio of the area of chalkiness to the area of the whole endosperm for each grain was evaluated by visual assessment on the chalkiness visualizer. The values were averaged and used as values for ACE. DEC was calculated as the product of PGWC × ACE. Rapid viscosity analyzer (RVA) profiles were characterized by six parameters such as peak paste viscosity (PKV), hot paste viscosity (HPV), cool paste viscosity (CPV), breakdown viscosity (BDV = PKV-HPV), consistency viscosity (CSV = CPV-HPV), and setback viscosity (SBV = CPV-PKV) as described in Brabender (1998) .
Determination of photosynthesis efficiency
Maximum quantum efficiency of PS II photochemistry (Fv/Fm) and noncyclic electron flow (ΦPSII) of rice leaves were measured using a PAM-2000 portable PAM fluorometer (Walz Effeltrich, Germany) with the software DA-2000 (Heinz Walz). For each sample, at least nine leaves were measured.
Measurement of sucrose synthase, AGPase, BE, and DBE
All enzymatic activity measurements were carried out in a 4°C cold chamber. In general, five immature rice grains without hull, pericarp, and embryo at the late-milking stage were homogenized in 1 mL of solution composed of 50 mM HEPES-NaOH (pH 7.4), 2 mM MgCl2, 50 mM 2-mercaptoethanol, and 12.5% (v/v) glycerol. The homogenate was centrifuged twice at 15, 000 g for 15 min. The supernatant was used as the crude enzyme extract [59, 60]. The activities of sucrose synthase, AGPase, and BE were assayed as described [8, 61]. Activity of DBE was measured using the methods of Nelson (1944)  and Somogyi (1952) .
Measurement of grain H2O2 levels
H2O2 concentrations in Asominori and CSSL50-1 endosperm were measured according to Wan and Liu (2008)  with minor modifications. Briefly, rice endosperm of 15 DAF (1.0 g) were ground with a mortar and pestle in liquid nitrogen to fine powders and added to a 10-ml cuvette containing 8 ml of double distilled H2O and 2 ml of 25 mM titanium sulfate and then incubated for 1 h at room temperature. Oxidation of titanium sulfate was recorded by reading A410. Readings were converted to corresponding concentrations using a standard calibration plot.
RNA extraction, GeneChip hybridization, and initial data analysis
RNA samples were processed according to Affymetrix manual. Total RNA was isolated using TRIzol reagent. RNA was then purified using an RNeasy spin column (Qiagen) and an on-column DNase treatment. Hybridization of Affrymetrix rice GeneChips and initial data collection were conducted at CapitalBio Corporation (Beijing, China). A total of 6 chips, with three biological replicates for each sample, were used in the assay. The hybridization data were analyzed using GeneChip Operating software (GCOS 1.4). A global scaling procedure was performed to normalize different arrays using dChip software, which incorporates a statistical model for expression array data at the probe level. The expression values were log2 transformed after calculating the expression index. Two-class unpaired method in the SAM (Significant Analysis of Microarrays) software (Tusher et al. 2001)  was used to identify the differentially expressed genes. One-way ANOVA was applied as an alternative statistic tool to further filter the differentially expressed genes.
The differentially expressed genes were classified using the Gene Ontology scheme http://www.geneontology.org/ and enriched GO terms were evaluated using Fisher's exact test. Rice annotation information was downloaded from AFFYMETRIX netaffx-annotation http://www.affymetrix.com/support/support_result.affx. Manual functional classification of differentially expressed genes was performed with the aid of the bioinformatics tool JAFA http://jafa.burnham.org. The microarray data has been deposited in the EMBL ArrayExpress database under the accession number E-MTAB-397.
Semi-quantitative RT-PCR analysis
Five micrograms of RNA (15 DAF endosperms) were used for reverse-transcription (SuperScript II; Invitrogen). An aliquot of the first-strand cDNA mixture corresponding to 6.25 ng of total RNA was used as a template with 0.5 units of Taq polymerase (ExTaq; TaKaRa) in 50 uL volume. In general, after initial 5 min at 94°C, 30 cycles of 94°C for 30 s, 55°C for 30 s, 72°C for 1 min were performed with a final extension at 72°C for 10 min. Sequences of primers are listed in Additional file 5 (Table S3). PCR products were separated by electrophoresis in 1.5% agarose gels, stained with ethidium bromide, and visualized using the BioDoc-It system (UVP).
This work was supported in part by grants from the National Program on the Development of Basic Research ("973", #2006CB101706) and the National High Technology Research and Development Program of China ("863", #2007AA10Z137).
- Zhang Q: Strategies for developing Green Super Rice. Proc Natl Acad Sci USA. 2007, 104: 16402-16409. 10.1073/pnas.0708013104.PubMed CentralPubMedView ArticleGoogle Scholar
- Del-Rosario AR, Briones VP, Vidal AJ, Juliano BO: Composition and endosperm structure of developing and mature rice kernel. Cereal Chem. 1968, 45: 225-235.Google Scholar
- Myers AM, Morell MK, James MG, Ball SG: Recent progress toward understanding biosynthesis of the amylopectin crystal. Plant Physiol. 2000, 122: 989-997. 10.1104/pp.122.4.989.PubMed CentralPubMedView ArticleGoogle Scholar
- Kang HG, Park S, Matsuoka M, An G: White-core endosperm floury endosperm-4 in rice is generated by knockout mutations in the C-type pyruvate orthophosphate dikinase gene (OsPPDKB). Plant J. 2005, 42: 901-911. 10.1111/j.1365-313X.2005.02423.x.PubMedView ArticleGoogle Scholar
- Yamakawa H, Hirose T, Kuroda M, Yamaguchi T: Comprehensive expression profiling of rice grain filling-related genes under high temperature using DNA microarray. Plant Physiol. 2007, 144: 258-277. 10.1104/pp.107.098665.PubMed CentralPubMedView ArticleGoogle Scholar
- Dung LV, Mikami I, Amano E, Sano Y: Study on the response of dull endosperm 2-2, du2-2, to two Wx alleles in rice. Breeding Sci. 2000, 50: 215-219.View ArticleGoogle Scholar
- Isshiki M, Nakajima M, Satoh H, Shimamoto K: dull: rice mutants with tissue-specific effects on the splicing of the waxy pre-mRNA. Plant J. 2000, 23: 451-460. 10.1046/j.1365-313x.2000.00803.x.PubMedView ArticleGoogle Scholar
- Nishi A, Nakamura Y, Tanaka N, Satoh H: Biochemical and genetic analysis of the effects of amylose-extender mutation in rice endosperm. Plant Physiol. 2001, 127: 459-472. 10.1104/pp.010127.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawasaki T, Mizuno K, Baba T: Coordinated regulation of the gene participating in starch biosynthesis by the rice floury-2 locus. Plant Physiol. 1996, 110: 459-472.Google Scholar
- Nakamura Y, Kubo AK, Shimamune T, Matsuda T, Harada K, Satoh H: Correlation between activities of starch debranching enzyme and α-polyglucan structure in endosperms of sugary-1 mutants of rice. Plant J. 1997, 12: 143-153. 10.1046/j.1365-313X.1997.12010143.x.View ArticleGoogle Scholar
- Tashiro T, Wardlaw IF: The effect of high temperature on kernel dimensions and the type and occurrence of kernel damage in rice. Aust J Agric Res. 1991, 42: 485-496. 10.1071/AR9910485.View ArticleGoogle Scholar
- Zakaria S, Matsuda T, Tajima S, Nitta Y: Effect of high temperature at ripening stage on the reserve accumulation in seed in some rice cultivars. Plant Prod Sci. 2002, 5: 160-168. 10.1626/pps.5.160.View ArticleGoogle Scholar
- Cheng FM, Hu DW, Ding YS: Dynamic change of chalkiness and observation of grain endosperm structure with scanning electron microscope under controlled temperature condition. Chinese J Rice Sci. 2000, 14: 83-87.Google Scholar
- Shen B: Observation on the starch grain development of early indica rice during chalkiness formation with scanning electronic microscope. Chinese J Rice Sci. 2000, 14: 225-228.Google Scholar
- Yang F, Song H, Cui XY, Gao W: Observation on amyloplast development in endosperm of different chalkiness Japonica rice with scanning electronic microscope. Acta Agronomica Sinica. 2004, 30: 406-408.Google Scholar
- He P, Li SG, Qian Q, Ma YQ, Li JZ, Wang WM, Chen Y, Zhu LH: Genetic analysis of rice grain quality. Theor Appl Genet. 1999, 98: 502-508. 10.1007/s001220051098.View ArticleGoogle Scholar
- Tan YF, Xing YZ, Li JX, Yu SB, Xu CG, Zhang Q: Genetic bases of appearance quality of rice grains in Shanyou 63, an elite rice hybrid. Theor Appl Genet. 2000, 10: 823-829. 10.1007/s001220051549.View ArticleGoogle Scholar
- Wan XY, Weng JF, Zhai HQ, Wan JM: Fine mapping of pgwc-8 gene affecting percentage of grains with chalkiness in rice (Oryza sativa. L). Rice Genetics Newsletter. 2004, 21: 54-56.Google Scholar
- Wan XY, Wan JM, Weng JF, Jiang L, Bi JC, Wang CM, Zhai HQ: Stability of QTLs for rice grain dimension and endosperm chalkiness characteristics across eight environments. Theor Appl Genet. 2005, 110: 1334-1346. 10.1007/s00122-005-1976-x.PubMedView ArticleGoogle Scholar
- Woo MO, Ham TH, Ji HS, Choi MS, Jiang W, Chu SH, Piao R, Chin JH, Kim JA, Park BS, Seo HS, Jwa NS, McCouch S, Koh HJ: Inactivation of the UGPase1 gene causes genic male sterility and endosperm chalkiness in rice (Oryza sativa L.). Plant J. 2008, 54: 190-204. 10.1111/j.1365-313X.2008.03405.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Tsunematsu H, Yoshimura A, Harushima Y, Nagamura Y, Kurata N, Yano M, Sasaki T, Iwata N: RFLP framework map using recombinant inbred lines in rice. Breeding Sci. 1996, 46: 279-284.Google Scholar
- Han XZ, Hamaker BR: Amylopectin fine structure and rice starch paste breakdown. J Cereal Sci. 2001, 34: 279-284. 10.1006/jcrs.2001.0374.View ArticleGoogle Scholar
- Boehlein SK, Shaw JR, Stewart JD, Hannah LC: Heat stability and allosteric properties of the maize endosperm ADP-glucose pyrophosphorylase are intimately intertwined. Plant Physiol. 2008, 146: 289-299. 10.1104/pp.107.109942.PubMed CentralPubMedView ArticleGoogle Scholar
- Saha BC: Hemicellulose bioconversion. J Ind Microbiol Biotechnol. 2003, 30: 279-291. 10.1007/s10295-003-0049-x.PubMedView ArticleGoogle Scholar
- Dixon DP, Lapthorn A, Edwards R: Plant glutathione transferases. Genome Biol. 2002, 3: REVIEWS3004-10.1186/gb-2002-3-3-reviews3004.PubMed CentralPubMedView ArticleGoogle Scholar
- Vieira Dos Santos C, Rey P: Plant thioredoxins are key actors in the oxidative stress response. Trends Plant Sci. 2006, 11: 329-334. 10.1016/j.tplants.2006.05.005.PubMedView ArticleGoogle Scholar
- Zen DL, Qian Q, Ruan LQ, Teng S, Y. K, Fujimoto H, Zhu LH: QTL analysis of chalkiness size in three dimensions. Chinese J Rice Sci. 2002, 16: 11-14.Google Scholar
- Li ZF, Wan JM, Xia JF, Zhai HQ: Mapping quantitative trait loci underlying appearance quality of rice grains (Oryza sativa L.). Yi Chuan Xue Bao. 2003, 30: 251-259.PubMedGoogle Scholar
- Septiningsih EM, Trijatmiko KR, Moeljopawiro S, McCouch SR: Identification of quantitative trait loci for grain quality in an advanced backcross population derived from the Oryza sativa variety IR64 and the wild relative O. rufipogon. Theor Appl Genet. 2003, 107: 1433-1441. 10.1007/s00122-003-1376-z.PubMedView ArticleGoogle Scholar
- Li J, Xiao J, Grandillo S, Jiang L, Wan Y, Deng Q, Yuan L, McCouch SR: QTL detection for rice grain quality traits using an interspecific backcross population derived from cultivated Asian (O. sativa L.) and African (O. glaberrima S.) rice. Genome. 2004, 47: 697-704. 10.1139/g04-029.PubMedView ArticleGoogle Scholar
- Shouichi Y: Effects of air temperature and light on grain filling of an indica and a japonica under controlled environmental conditions. Soil sci plant Nut. 1987, 23: 93-104.Google Scholar
- Guo YP, Mu P, Liu JF, Lu YX, Li ZC: Comparative Studies on Quality Characters of Rice under Water-and Dry-cultivation Conditions. Acta Agronomica Sinica. 2005, 17: 1443-1448.Google Scholar
- Cai Y, Wang W, Zhu Z, Zhang Z, Langm Y, Zhu Q: Effects of water stress during grain-filling period on rice grain yield and its quality under different nitrogen levels. Chinese Journal of Applied Ecology. 2006, 17: 1201-1206.PubMedGoogle Scholar
- Patindol J, Wang YJ: Fine structures and physicochemical properties of starches from chalky and translucent rice kernels. J Agric Food Chem. 2003, 51: 2777-2784. 10.1021/jf026101t.PubMedView ArticleGoogle Scholar
- Zhong LJ, Cheng FM: Varietal differences in amylose accumulation and activities of major enzymes associated with starch synthesis during grain filling in rice. Acta Agronomica Sinica. 2003, 29: 452-456.Google Scholar
- Yang JC, Peng SB, Gu SL, Visperas RM, Zhu QS: Changes in activities of three enzymes associated with starch synthesis in rice grains during grain filling. Acta Agronomica Sinica. 2007, 27: 157-164.Google Scholar
- Mittler R: Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002, 7: 405-410. 10.1016/S1360-1385(02)02312-9.PubMedView ArticleGoogle Scholar
- Bienert GP, Schjoerring JK, Jahn TP: Membrane transport of hydrogen peroxide. Biochim Biophys Acta. 2006, 1758: 994-1003. 10.1016/j.bbamem.2006.02.015.PubMedView ArticleGoogle Scholar
- Wan XY, Liu JY: Comparative proteomic analysis reveals an intimate protein network provoked by hydrogen peroxide stress in rice seedling leaves. Molecular and Cellular Proteomics. 2008, 7: 1469-1488. 10.1074/mcp.M700488-MCP200.PubMed CentralPubMedView ArticleGoogle Scholar
- Mittler R, Vanderauwera S, Gollery M, Van Breusegem F: Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9: 490-498. 10.1016/j.tplants.2004.08.009.PubMedView ArticleGoogle Scholar
- Yang J, Lin R, Sullivan J, Hoecker U, Liu B, Xu L, Deng XW, Wang H: Light regulates COP1-mediated degradation of HFR1, a transcription factor essential for light signaling in Arabidopsis. Plant Cell. 2005, 17: 804-821. 10.1105/tpc.104.030205.PubMed CentralPubMedView ArticleGoogle Scholar
- Apel K, Hirt H: Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004, 55: 373-399. 10.1146/annurev.arplant.55.031903.141701.PubMedView ArticleGoogle Scholar
- Espartero J, Sanchez-Aguayo I, Pardo JM: Molecular characterization of glyoxalase-I from a higher plant; upregulation by stress. Plant Mol Biol. 1995, 29: 1223-1233. 10.1007/BF00020464.PubMedView ArticleGoogle Scholar
- Maccarrone M, Melinol G, Finazzi-Agro A: Lipoxygenases and their involvement in programmed cell death. Cell Death and Differentiation. 2001, 8: 776-784. 10.1038/sj.cdd.4400908.PubMedView ArticleGoogle Scholar
- Lin SK, Chang MC, Tsai YG, Lur HS: Proteomic analysis of the expression of proteins related to rice quality during caryopsis development and the effect of high temperature on expression. Proteomics. 2005, 5: 2140-2156. 10.1002/pmic.200401105.PubMedView ArticleGoogle Scholar
- Yang J, Zhang J, Wang Z, Zhu Q: Activities of starch hydrolytic enzymes and sucrose-phosphate synthase in the stems of rice subjected to water stress during grain filling. J Exp Bot. 2001, 52: 2169-2179.PubMedGoogle Scholar
- Yang J, Zhang J, Wang Z, Zhu Q, Liu L: Abscisic acid and cytokinins in the root exudates and leaves and their relationship to senescence and remobilization of carbon reserves in rice subjected to water stress during grain filling. Planta. 2002, 215: 645-652. 10.1007/s00425-002-0789-2.PubMedView ArticleGoogle Scholar
- Schaeffer GW, Sharpe FT: Increased Lysine and Seed Storage Protein in Rice Plants Recovered from Calli Selected with Inhibitory Levels of Lysine plus Threonine and S-(2-Aminoethyl)cysteine. Plant Physiol. 1987, 84: 509-515. 10.1104/pp.84.2.509.PubMed CentralPubMedView ArticleGoogle Scholar
- Lunde C, Zygadlo A, Simonsen HT, Nielsen PL, Blennow A, Haldrup A: Sulfur starvation in rice: the effect on photosynthesis, carbohydrate metabolism, and oxidative stress protective pathways. Physiol Plant. 2008, 134: 508-521. 10.1111/j.1399-3054.2008.01159.x.PubMedView ArticleGoogle Scholar
- Gorantla M, Babu PR, Lachagari VB, Reddy AM, Wusirika R, Bennetzen JL, Reddy AR: Identification of stress-responsive genes in an indica rice (Oryza sativa L.) using ESTs generated from drought-stressed seedlings. J Exp Bot. 2007, 58: 253-265. 10.1093/jxb/erl213.PubMedView ArticleGoogle Scholar
- Cheng C, Yun KY, Ressom HW, Mohanty B, Bajic VB, Jia Y, Yun SJ, de los Reyes BG: An early response regulatory cluster induced by low temperature and hydrogen peroxide in seedlings of chilling-tolerant japonica rice. BMC Genomics. 2007, 8: 175-10.1186/1471-2164-8-175.PubMed CentralPubMedView ArticleGoogle Scholar
- Cho SK, Jeung JU, Kang KH, Shim KS, Jung KW, You MK, Ok SH, Chung YS, Hwang HG, Choi HC, Moon HP, Shin JS: Identification of genes induced in wound-treated wild rice (Oryza minuta). Mol Cells. 2004, 17: 230-236.PubMedGoogle Scholar
- Kawasaki S, Borchert C, Deyholos M, Wang H, Brazille S, Kawai K, Galbraith D, Bohnert HJ: Gene expression profiles during the initial phase of salt stress in rice. Plant Cell. 2001, 13: 889-905. 10.1105/tpc.13.4.889.PubMed CentralPubMedView ArticleGoogle Scholar
- Rakwal R, Kimura S, Shibato J, Nojima K, Kim YK, Nahm BH, Jwa NS, Endo S, Tanaka K, Iwahashi H: Growth retardation and death of rice plants irradiated with carbon ion beams is preceded by very early dose- and time-dependent gene expression changes. Mol Cells. 2008, 25: 272-278.PubMedGoogle Scholar
- Kubo T, Nakamura K, Yoshimura A: Development of a series of Indica chromosome segment substitution lines in Japonica background of rice. Rice Genet Newsl. 1999, 16: 104-106.Google Scholar
- Fujita N, Yoshida M, Asakura N, Ohdan T, Miyao A, Hirochika H, Nakamura Y: Function and characterization of starch synthase I using mutants in rice. Plant Physiol. 2006, 140: 1070-1084. 10.1104/pp.105.071845.PubMed CentralPubMedView ArticleGoogle Scholar
- Fujita N, Yoshida M, Kondo T, Saito K, Utsumi Y, Tokunaga T, Nishi A, Satoh H, Park JH, Jane JL, Miyao A, Hirochika H, Nakamura Y: Characterization of SSIIIa-deficient mutants of rice: the function of SSIIIa and pleiotropic effects by SSIIIa deficiency in the rice endosperm. Plant Physiol. 2007, 144: 2009-2023. 10.1104/pp.107.102533.PubMed CentralPubMedView ArticleGoogle Scholar
- Brabender M: The new MICRO-VISCO-AMYLO-GRAPH: comparison of some results with those of the Viscograph. Poster presentation at 1998 American Association of Cereal Chemists Annual Meeting Minneapoils. 1998Google Scholar
- Yamanouchi H, Nakamura Y: Organ specificity of isoforms of starch branching enzyme (Q-enzyme) in rice. Plant Cell Physiol. 1992, 33: 985-991.Google Scholar
- Umemoto T, Nakamura Y, Ishikura N: Activity of starch synthase and the amylose content in rice endosperm. Phytochemistry. 1995, 40: 1613-1616. 10.1016/0031-9422(95)00380-P.View ArticleGoogle Scholar
- Nakamura Y, Yuki K, Park SY, Ohya T: Carbohydrate metabolism in the developing endosperm of rice grain. Plant Cell Physiol. 1989, 30: 833-839.Google Scholar
- Nelson N: A photometric adaptation of Somogyi method for the determination of glucose. J Biol Chem. 1944, 153: 375-380.Google Scholar
- Somogyi M: Notes on sugar determination. J Biol Chem. 1952, 195: 19-23.Google Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.