MucoRice-CTB line 19A, a new marker-free transgenic rice-based cholera vaccine produced in an LED-based hydroponic system.

We previously established the selection-marker-free rice-based oral cholera vaccine (MucoRice-CTB) line 51A for human use by Agrobacterium-mediated co-transformation and conducted a double-blind, randomized, placebo-controlled phase I trial in Japan and the United States. Although MucoRice-CTB 51A was acceptably safe and well tolerated by healthy Japanese and U.S. subjects and induced CTB-specific antibodies neutralizing cholera toxin secreted by Vibrio cholerae, we were limited to a 6-g cohort in the U.S. trial because of insufficient production of MucoRice-CTB. Since MucoRice-CTB 51A did not grow in sunlight, we re-examined the previously established marker-free lines and selected MucoRice-CTB line 19A. Southern blot analysis of line 19A showed a single copy of the CTB gene. We resequenced the whole genome and detected the transgene in an intergenic region in chromosome 1. After establishing a master seed bank of MucoRice-CTB line 19A, we established a hydroponic production facility with LED lighting to reduce electricity consumption and to increase production capacity for clinical trials. Shotgun MS/MS proteomics analysis of MucoRice-CTB 19A showed low levels of α-amylase/trypsin inhibitor-like proteins (major rice allergens), which was consistent with the data for line 51A. We also demonstrated that MucoRice-CTB 19A had high oral immunogenicity and induced protective immunity against cholera toxin challenge in mice. These results indicate that MucoRice-CTB 19A is a suitable oral cholera vaccine candidate for Phase I and II clinical trials in humans, including a V. cholerae challenge study.

Raman Multi-Omic Snapshots of Koshihikari Rice Kernels Reveal Important Molecular Diversities with Potential Benefits in Healthcare.

This study exploits quantitative algorithms of Raman spectroscopy to assess, at the molecular scale, the nutritional quality of individual kernels of the Japanese short-grain rice cultivar Koshihikari in terms of amylose-to-amylopectin ratio, fractions of phenylalanine and tryptophan aromatic amino acid residues, protein-to-carbohydrate ratio, and fractions of protein secondary structures. Statistical assessments on a large number of rice kernels reveal wide distributions of the above nutritional parameters over nominally homogeneous kernel batches. This demonstrates that genetic classifications cannot catch omic fluctuations, which are strongly influenced by a number of extrinsic factors, including the location of individual grass plants within the same rice field and the level of kernel maturation. The possibility of collecting nearly real-time Raman "multi-omic snapshots" of individual rice kernels allows for the automatic (low-cost) differentiation of groups of kernels with restricted nutritional characteristics that could be used in the formulation of functional foods for specific diseases and in positively modulating the intestinal microbiota for protection against bacterial infection and cancer prevention.

In planta evidence that the HAK transporter OsHAK2 is involved in Na+ transport in rice.

HAK family transporters primarily function as K+ transporters and play major roles in K+ uptake and translocation in plants, whereas several HAK transporters exhibit Na+ transport activity. OsHAK2, a rice HAK transporter, was shown to mediate Na+ transport in Escherichia coli in a previous study. In this study, we investigated whether OsHAK2 is involved in Na+ transport in the rice plant. Overexpression of OsHAK2 increased Na+ translocation from the roots to the shoots of transgenic rice. It also increased both root and whole-plant Na+ content, and enhanced shoot length under low Na+ and K+ conditions. Meanwhile, OsHAK2 overexpression increased salt sensitivity under a long-term salt stress condition, indicating that OsHAK2 is not involved in salt tolerance, unlike in the case of ZmHAK4 in maize. These results suggest that OsHAK2 is permeable to Na+ and contributes to shoot growth in rice plants under low Na+ and K+ conditions.

Alleviation of high light stress in shade-treated tea leaves by acclimation to light before shade removal.

Shade cultivation of tea plants (Camellia sinensis L.) is employed for the production of high-quality green tea which increases the content of chlorophylls and free amino acids, including theanine. However, shaded tea plants suffer from photooxidative stress caused by sudden exposure to high light (HL) when the shade is removed. In this study, we tried to acclimatize shaded tea plants to light prior to shade removal to alleviate HL-induced stress. Acclimated tea plants showed milder photoinhibition in response to HL exposure than the shaded plants without acclimation. Moreover, there were no large differences in the total chlorophylls and free amino acids (including theanine) content between acclimated and non-acclimated plants. These results indicate that acclimation of shaded tea plants can alleviate subsequent HL stress without causing large changes in the content of chemical components associated with tea quality.

Laser microdissection transcriptome data derived gene regulatory networks of developing rice endosperm revealed tissue- and stage-specific regulators modulating starch metabolism.

KEY MESSAGE: Laser microdissection applied on the developing rice endosperm revealed tissue- and stage-specific regulators modulating programmed cell death and desiccation tolerance mechanisms in the central starchy endosperm following starch metabolism. Rice (Oryza sativa L.) filial seed tissues are heterozygous in its function, which accumulate distinct storage compounds spatially in starchy endosperm and aleurone. In this study, we identified the 18 tissue- and stage-specific gene co-regulons in the developing endosperm by isolating four fine tissues dorsal aleurone layer (AL), central starchy endosperm (CSE), dorsal starchy endosperm (DSE), and lateral starchy endosperm (LSE) at two developmental stages (7 days after flowering, DAF and 12DAF) using laser microdissection (LM) coupled with gene expression analysis of a 44 K microarray. The derived co-expression regulatory networks depict that distinct set of starch biosynthesis genes expressed preferentially at first in CSE at 7 DAF and extend its spatial expression to LSE and DSE by 12 DAF. Interestingly, along with the peak of starch metabolism we noticed accumulation of transcripts related to phospholipid and glycolipid metabolism in CSE during 12 DAF. The spatial distribution of starch accumulation in distinct zones of starchy endosperm contains specific transcriptional factors and hormonal-regulated genes. Genes related to programmed cell death (PCD) were specifically expressed in CSE at 12DAF, when starch accumulation was already completed in that tissue. The aleurone layer present in the outermost endosperm accumulates transcripts of lipid, tricarboxylic acid metabolism, several transporters, while starch metabolism and PCD is not pronounced. These regulatory cascades are likely to play a critical role in determining the positional fate of cells and offer novel insights into the molecular physiological mechanisms of endosperm development from early to middle storage phase.

Raman Fingerprints of Rice Nutritional Quality: A Comparison between Japanese Koshihikari and Internationally Renowned Cultivars.

Raman spectroscopy was applied to characterize at the molecular scale the nutritional quality of the Japanese Koshihikari rice cultivar in comparison with other renowned rice cultivars including Carnaroli from Italy, Calrose from the USA, Jasmine rice from Thailand, and Basmati from both India and Pakistan. For comparison, two glutinous (mochigome) cultivars were also investigated. Calibrated and validated Raman analytical algorithms allowed quantitative determinations of: (i) amylopectin and amylose concentrations, (ii) fractions of aromatic amino acids, and (iii) protein content and secondary structure. The Raman assessments non-destructively linked the molecular composition of grains to key nutritional parameters and revealed a complex intertwine of chemical properties. The Koshihikari cultivar was rich in proteins (but with low statistical relevance as compared to other investigated cultivars) and aromatic amino acids. However, it also induced a clearly higher glycemic impact as compared to long-grain cultivars from Asian countries. Complementary to genomics and wet-chemistry analyses, Raman spectroscopy makes non-destructively available factual and data-driven information on rice nutritional characteristics, thus providing customers, dietitian nutritionists, and producers with a solid science-consolidated platform.

Raman Molecular Fingerprints of Rice Nutritional Quality and the Concept of Raman Barcode.

The nutritional quality of rice is contingent on a wide spectrum of biochemical characteristics, which essentially depend on rice genome, but are also greatly affected by growing/environmental conditions and aging during storage. The genetic basis and related identification of genes have widely been studied and rationally linked to accumulation of micronutrients in grains. However, genetic classifications cannot catch quality fluctuations arising from interannual, environmental, and storage conditions. Here, we propose a quantitative spectroscopic approach to analyze rice nutritional quality based on Raman spectroscopy, and disclose analytical algorithms for the determination of: (i) amylopectin and amylose concentrations, (ii) aromatic amino acids, (iii) protein content and structure, and (iv) chemical residues. The proposed Raman algorithms directly link to the molecular composition of grains and allow fast/non-destructive determination of key nutritional parameters with minimal sample preparation. Building upon spectroscopic information at the molecular level, we newly propose to represent the nutritional quality of labeled rice products with a barcode specially tailored on the Raman spectrum. The Raman barcode, which can be stored in databases promptly consultable with barcode scanners, could be linked to diet applications (apps) to enable a rapid, factual, and unequivocal product identification based on direct molecular screening.

Raman spectroscopic analysis of polysaccharides in popular Japanese rice cultivars.

Analytical algorithms based on Raman spectroscopy are proposed for the determination of amylopectin and amylose concentrations in polished white rice, and applied to characterize and compare linear and branched polysaccharide structures in nine different types of Japanese rice. A selected algorithm used symmetric bending vibrations of the COC glycosidic linkage from a relatively narrow spectral zone between 830 and 895 cm-1. It specifically compared the intensity of Raman signals from two types of bending common to both starch components (C1-O-C5 and C1-O-C4 at 868 and 855 cm-1, respectively) and that at the branch point peculiar to amylopectin (C1-O-C6 at 844 cm-1). Raman data were confronted with data collected by conventional amylose-iodine colorimetry method. Consistency was found between Raman and colorimetric methods over the entire series of tested rice cultivars, thus validating the newly proposed spectroscopic algorithm. The amylose content of the tested rice species broadly varied between 1.2 and 20.4%. The proposed Raman algorithm allows fast and nondestructive determination of amylose content in rice with minimal sample preparation. These characteristics might be key in the development of portable Raman devices capable to promptly screen polysaccharides in different rice cultivars with respect to their interannual and plantation-related fluctuations.

Rapidly evolving phosphoenolpyruvate carboxylase Gmppc1 and Gmppc7 are highly expressed in the external seed coat of immature soybean seeds.

Phosphoenolpyruvate carboxylase (PEPC) is a carbon fixation enzyme which probably plays crucial roles in seed development. A greater number of PEPC isoforms are encoded in the soybean genome, while most of the PEPC isoforms are functionally unknown. In this study, we investigated on soybean PEPC expressed in the external layer of seed coat (ELSC) during seed formation. PEPC activity in ELSC ranged from 0.24 to 1.0 U/g F.W., which could be comparable to those in whole seeds at U per dry matter. Public RNA-Seq data in separated soybean seed tissues revealed that six plant-type PEPC isogenes were substantially expressed in ELSC, and Gmppc1 and Gmppc7 were highly expressed in hourglass cells of ELSC. Gene Ontology enrichment of co-expressed genes with Gmppc1 and Gmppc7 implicated a role of these isogenes in assisting energy production and cellulose biosynthesis. Comparison of PEPC sequences from 16 leguminous species hypothesized adaptive evolution of the Gmppc1 and Gmppc7 lineage after divergence from the other plant-type PEPC lineages. Molecular diversification of these plant-type PEPC was possibly accomplished by adaptation to the functions of the soybean seed tissues. This study indicates that energy demand in immature seeds may be a driving force for the molecular evolution of PEPC.

Determination of free fatty acids in plasma by gas chromatography.

A method was developed for determination of free fatty acids (FFAs) in plasma by gas chromatography. Plasma was extracted with 3 vol of methanol. Most cholesterol esters and triacylglycerols did not dissolve in the aqueous methanol. FFAs in the crude lipid solution were directly and selectively methylated with (trimethylsilyl)diazomethane at room temperature. Fatty acid methyl esters (FAMEs) formed were extracted with hexane, and nonreactive phospholipids were washed out with 95% methanol. The partially purified FAME preparation was analyzed by gas chromatography. The composition and amount of plasma FFAs closely approximated those obtained using two different methods.

Stress Responses of Shade-Treated Tea Leaves to High Light Exposure after Removal of Shading.

High-quality green tea is produced from buds and young leaves grown by the covering-culture method, which employs shading treatment for tea plants (Camellia sinensis L.). Shading treatment improves the quality of tea, but shaded tea plants undergo sudden exposures to high light (HL) at the end of the treatment by shade removal. In this study, the stress response of shaded tea plants to HL illumination was examined in field condition. Chl a/b ratio was lower in shaded plants than nonshaded control, but it increased due to exposure to HL after 14 days. Rapid decline in Fv/Fm values and increases in carbonylated protein level were induced by HL illumination in the shaded leaves on the first day, and they recovered thereafter between a period of one and two weeks. These results revealed that shaded tea plants temporarily suffered from oxidative damages caused by HL exposure, but they could also recover from these damages in 2 weeks. The activities of antioxidant enzymes, total ascorbate level, and ascorbate/dehydroascorbate ratio were decreased and increased in response to low light and HL conditions, respectively, suggesting that the upregulation of antioxidant defense systems plays a role in the protection of the shaded tea plants from HL stress.

The plant-type phosphoenolpyruvate carboxylase Gmppc2 is developmentally induced in immature soy seeds at the late maturation stage: a potential protein biomarker for seed chemical composition.

Phosphoenolpyruvate carboxylase (PEPC) is a carbon-fixing enzyme with critical roles in seed development. Previously we observed a positive correlation between PEPC activity and protein content in mature seeds among soybean cultivars and varietal differences of PEPC activity in immature seeds, which is concordant with seed protein accumulation. Here, we report a PEPC isoform (Gmppc2) which is preferentially expressed in immature soybean seeds at the late maturation stage. Gmppc2 was co-expressed with enzyme genes involved in starch degradation: α-amylase, hexokinase, and α-glucan phosphorylase. Gmppc2 was developmentally induced in the external seed coats, internal seed coats, hypocotyls, and cotyledons at the late maturation stage. The expression of Gmppc2 protein was negatively regulated by the application of a nitrogen fertilizer, which suppressed nodule formation. These results imply that Gmppc2 is involved in the metabolism of nitrogen originated from nodules into seeds, and Gmppc2 might be applicable as a biomarker of seed protein content.Abbreviations: PEP: phosphoenolpyruvate; PEPC: phosphoenolpyruvate carboxylase; RNA-Seq: RNA sequencing; PCA: principal component analysis; SE: standard error.

Essential roles of autophagy in metabolic regulation in endosperm development during rice seed maturation.

Autophagy plays crucial roles in the recycling of metabolites, and is involved in many developmental processes. Rice mutants defective in autophagy are male sterile due to immature pollens, indicating its critical role in pollen development. However, physiological roles of autophagy during seed maturation had remained unknown. We here found that seeds of the rice autophagy-deficient mutant Osatg7-1, that produces seeds at a very low frequency in paddy fields, are smaller and show chalky appearance and lower starch content in the endosperm at the mature stage under normal growth condition. We comprehensively analyzed the effects of disruption of autophagy on biochemical properties, proteome and seed quality, and found an abnormal activation of starch degradation pathways including accumulation of α-amylases in the endosperm during seed maturation in Osatg7-1. These results indicate critical involvement of autophagy in metabolic regulation in the endosperm of rice, and provide insights into novel autophagy-mediated regulation of starch metabolism during seed maturation.

Pattern analysis suggests that phosphoenolpyruvate carboxylase in maturing soybean seeds promotes the accumulation of protein.

The protein and oil contents in soybean seeds are major factors in seed quality. Seed proteins and oils are synthesized from sucrose and nitrogenous compounds transported into maturing seeds. In this study, we compared changes in the activity of phosphoenolpyruvate carboxylase (PEPC) and the accumulation profiles of protein and oil in maturing seeds of two soybean cultivars, which exhibit different protein and oil contents in seeds, to determine the interrelationships of them. A principal component analysis indicated a concordance of seed PEPC activity with the protein content, but did not with the oil content. PEPC activity per seed was highest in the late maturation stage, when the physiological status of the vegetative organs drastically changed. The high-protein cultivar had higher PEPC activity compared to the low-protein cultivar. These results highlight the biological role of PEPC in the synthesis of protein, therefore it was implied that PEPC could be a biomarker in soybean breeding.Abbreviations: ANOVA: analysis of variance; DS: developmental stage; DW: dry weight; FW: fresh weight; NIR: near infrared; PEP(C): phosphoenolpyruvate (carboxylase); PC(A): principal component (analysis); S.E.: standard error; WC: water content.

Laser Microdissection-Based Tissue-Specific Transcriptome Analysis Reveals a Novel Regulatory Network of Genes Involved in Heat-Induced Grain Chalk in Rice Endosperm.

Heat stress occurrence during seed filling leads to the formation of a chalky portion in the limited zone of the starchy endosperm of rice grains. In this study, isolation of aleurone, dorsal, central and lateral tissues of developing endosperm by laser-microdissection (LM) coupled with gene expression analysis of a 44 K microarray was performed to identify key regulatory genes involved in the formation of milky-white (MW) and white-back (WB) grains during heat stress. Gene regulatory network analysis classified the genes changed under heat stress into five modules. The most distinct expression pattern was observed in modules where most of the small heat shock proteins and cellular organization genes were changed under heat stress in dorsal aleurone cells and dorsal starchy endosperm zones. The histological observation supported the significant increase in cell number and size of dorsal aleurone cells in WB grains. With regard to the central starchy endosperm zone, preferential down-regulation of high molecular weight heat shock proteins (HMW HSPs), including a prominent member encoding endoplasmic reticulum (ER) chaperones, by heat stress was observed, while changes in expression of starch biosynthesis genes were minimal. Characterization of transgenic plants suppressing endosperm lumenal binding protein gene (BiP1), an ER chaperone preferentially down-regulated at the MW zone under heat stress, showed evidence of forming the chalky grains without disturbing the expression of starch biosynthesis genes. The present LM-based comprehensive expression analysis provides novel inferences that HMW HSPs play an important role in controlling redox, nitrogen and amino acid metabolism in endosperm leading to the formation of MW and WB chalky grains under heat stress.

The localization of rice prolamin species in protein body type I is determined by the temporal control of gene expression of the respective prolamin promoters.

Rice prolamin species form a layered structure in the protein body type I (PB-I) storage organelle. Rice prolamins are classified as 10 kDa, 13a-1, 13a-2, 13b-1, 13b-2 and 16 kDa prolamin. Prolamin species form layer structure in PB-I in order of 10 kDa core, 13b-1 layer, 13a (13a-1 and 13a-2) and 16 kDa middle layer and 13b-2 outer-most layer. In a previous study, we showed that the fusion proteins in 13b-2 prolamin-GFP, 13a-1 prolamin-GFP and 10 kDa prolamin-GFP were localized in the same layer of PB-I as the native prolamin, when they were expressed by their respective native prolamin promoters. Our preliminary study suggested that the temporal control of the native prolamin promoters was responsible for the localization of the respective prolamins. The aim of this study was to determine whether the use of a prolamin promoter other than the native prolamin promoter would change the localization of prolamin-GFP fusion proteins. For this purpose, we generated transgenic lines expressing 13b-2 prolamin-GFP and 13a-1 prolamin-GFP fusion proteins driven by each prolamin promoter other than the native prolamin promoter. As a result, the localization of the fusion protein in PB-I was changed. Based on our results, foreign protein localization in PB-I can be achieved by the temporal control of the different prolamin promoters.

Accumulation of foreign polypeptides to rice seed protein body type I using prolamin portion sequences.

KEY MESSAGE: Rice prolamins are accumulated in endoplasmic reticulum (ER)-derived proteins bodies, although conserved sequences retained in ER are not confirmed. We investigated portion sequences of prolamins that must accumulate in PB-Is. Rice seed prolamins are accumulated in endoplasmic reticulum (ER)-derived protein body type I (PB-I), but ER retention sequences in rice prolamin polypeptides have not been confirmed. Here we investigated the lengths of the prolamin portion sequences required for accumulation in PB-Is. Of the rice prolamins, we compared 13a and 13b prolamins because the amino acid sequences of these prolamins are quite similar except for the presence or absence of Cys-residues. We also generated and analyzed transgenic rice expressing several prolamin portion sequence-GFP fusion proteins. We observed that in 13a prolamin, when the portion sequences were extended more than the 68th amino acid residue from the initiating methionine, the prolamin portion sequence-GFP fusion proteins were accumulated in PB-Is. In 13b prolamin, when the portion sequences were extended by more than the 82nd amino acid residue from the initiating methionine, the prolamin portion sequence-GFP fusion proteins were accumulated in PB-Is. When those fusion proteins were extracted under non-reduced or reduced conditions, the 13a prolamin portion sequence-GFP fusion proteins in PB-Is were soluble under only the reduced condition. In contrast, 13b prolamin portion sequence-GFP fusion proteins were soluble under both non-reduced and reduced conditions. These results suggest that the accumulation of 13a prolamin in PB-Is is associated with the formation of disulfide bonds and/or hydrophobicity in 13a prolamin polypeptide, whereas the accumulation of 13b prolamin in PB-Is was less involved in the formation of disulfide bonds.

Control of foreign polypeptide localization in specific layers of protein body type I in rice seed.

KEY MESSAGE: Prolamin-GFP fusion proteins, expressed under the control of native prolamin promoters, were localized in specific layers of PB-Is. Prolamin-GFP fusion proteins were gradually digested from outside by pepsin digestion. In rice seed endosperm, protein body type I (PB-I) has a layered structure consisting of prolamin species and is the resistant to digestive juices in the intestinal tract. We propose the utilization of PB-Is as an oral vaccine carrier to induce mucosal immune response effectively. If vaccine antigens are localized in a specific layer within PB-Is, they could be protected from gastric juice and be delivered intact to the small intestine. We observed the localization of GFP fluorescence in transgenic rice endosperm expressing prolamin-GFP fusion proteins with native prolamin promoters, and we confirmed that the foreign proteins were located in specific layers of PB-Is artificially. Each prolamin-GFP fusion protein was localized in specific layers of PB-Is, such as the outer-most layer, middle layer, and core region. Furthermore, to investigate the resistance of prolamin-GFP fusion proteins against pepsin digestion, we performed in vitro pepsin treatment. Prolamin-GFP fusion proteins were gradually digested from the peripheral region and the contours of PB-Is were made rough by in vitro pepsin treatment. These findings suggested that prolamin-GFP fusion proteins accumulating specific layers of PB-Is were gradually digested and exposed from the outside by pepsin digestion.

Expression of a rice glutaredoxin in aleurone layers of developing and mature seeds: subcellular localization and possible functions in antioxidant defense.

MAIN CONCLUSION: A rice glutaredoxin isoform (OsGrxC2;2) with antioxidant capacity is expressed abundantly in seed tissues and is localized to storage vacuoles in aleurone layers in developing and mature seeds. Seed tissues undergo drastic water loss at the late stage of seed development, and thus need to tolerate oxidative injuries associated with desiccation. We previously found a rice glutaredoxin isoform, OsGrxC2;2, as a gene expressed abundantly in developing seeds. Since glutaredoxin is involved in antioxidant defense, in the present study we investigated the subcellular localization and expression profile of OsGrxC2;2 and whether OsGrxC2;2 has a role in the defense against reactive oxygen species. Western blotting and immunohistochemistry revealed that the OsGrxC2;2 protein accumulated at a high level in the embryo and aleurone layers of developing and mature seeds. The OsGrxC2;2 in developing seeds was particularly localized to aleurone grains, which are storage organelles derived from vacuoles. Overexpression of OsGrxC2;2 resulted in an enhanced tolerance to menadione in yeast and methyl viologen in green leaves of transgenic rice plants. These results suggest that OsGrxC2;2 participates in the defense against oxidative stress in developing and mature seeds.

Laser microdissection-based gene expression analysis in the aleurone layer and starchy endosperm of developing rice caryopses in the early storage phase.

BACKGROUND: Rice endosperm is composed of aleurone cells in the outermost layers and starchy endosperm cells in the inner part. The aleurone layer accumulates lipids, whereas starchy endosperm mainly accumulates starch. During the ripening stage, the starch accumulation rate is known to be asynchronous, depending on the position of the starchy endosperm. Different physiological and molecular mechanisms are hypothesized to underlie the qualitative and quantitative differences in storage products among developing rice endosperm tissues. RESULTS: Target cells in aleurone layers and starchy endosperm were isolated by laser microdissection (LM), and RNAs were extracted from each endosperm tissue in the early storage phase. Genes important for carbohydrate metabolism in developing endosperm were analyzed using qRT-PCR, and some of the genes showed specific localization in either tissue of the endosperm. Aleurone layer-specific gene expression of a sucrose transporter, OsSUT1, suggested that the gene functions in sucrose uptake into aleurone cells. The expression levels of ADP-glucose pyrophosphorylase (AGPL2 and AGPS2b) in each endosperm tissue spatially corresponded to the distribution of starch granules differentially observed among endosperm tissues. By contrast, expressions of genes for sucrose cleavage-hexokinase, UDP-glucose pyrophosphorylase, and phosphoglucomutase-were observed in all endosperm tissues tested. Aleurone cells predominantly expressed mRNAs for the TCA cycle and oxidative phosphorylation. This finding was supported by the presence of oxygen (8 % concentration) and large numbers of mitochondria in the aleurone layers. In contrast, oxygen was absent and only a few mitochondria were observed in the starchy endosperm. Genes for carbon fixation and the GS/GOGAT cycle were expressed highly in aleurone cells compared to starchy endosperm. CONCLUSIONS: The transcript level of AGPL2 and AGPS2b encoding ADP-glucose pyrophosphorylase appears to regulate the asynchronous development of starch granules in developing caryopses. Aleurone cells appear to generate, at least partially, ATP via aerobic respiration as observed from specific expression of identified genes and large numbers of mitochondria. The LM-based expression analysis and physiological experiments provide insight into the molecular basis of the spatial and nutritional differences between rice aleurone cells and starchy endosperm cells.