Gibberellins (GAs) play important roles in regulating reproductive development, especially anther development. Our previous studies revealed that the MYB transcriptional factor GAMYB, an important component of GA signaling in cereal aleurone cells, is also important for anther development. Here, we examined the physiological functions of GA during anther development through phenotypic analyses of rice (Oryza sativa) GA-deficient, GA-insensitive, and gamyb mutants. The mutants exhibited common defects in programmed cell death (PCD) of tapetal cells and formation of exine and Ubisch bodies. Microarray analysis using anther RNAs of these mutants revealed that rice GAMYB is involved in almost all instances of GA-regulated gene expression in anthers. Among the GA-regulated genes, we focused on two lipid metabolic genes, a cytochrome P450 hydroxylase CYP703A3 and b-ketoacyl reductase, both of which might be involved in providing a substrate for exine and Ubisch body. GAMYB specifically interacted with GAMYB binding motifs in the promoter regions in vitro, and mutation of these motifs in promoter-b-glucuronidase (GUS) transformants caused reduced GUS expression in anthers. Furthermore, a knockout mutant for CYP703A3 showed gamyb-like defects in exine and Ubisch body formation. Together, these results suggest that GA regulates exine formation and the PCD of tapetal cells and that direct activation of CYP703A3 by GAMYB is key to exine formation.

Plant male reproductive development is highly organized and sensitive to various environmental stressors, including high temperature. We have established an experimental procedure to evaluate high temperature injury in japonica rice plants. High temperature treatment (39 degrees C/30 degrees C) starting at the microspore stage repeatedly reduced spikelet fertility in our system. Morphological observations revealed that pollen viability in plants exposed to high temperatures was lower than that in control plants. Most pollen grains in high temperature-treated plants displayed a normal round shape and stained reddish purple with Alexander's reagent; however, the pollen grains were very poorly attached and displayed limited germination on the stigma. To investigate gene regulatory mechanisms in the anther in high temperature environments, DNA microarray analysis was performed by comparing non-treated samples with samples treated with 2-4 d of high heat. Genes responsive to high temperatures were identified from clustering of microarray data. Among these, at least 13 were designated as high temperature-repressed genes in the anther. Expression analyses revealed that these genes were expressed specifically in the immature anther mainly in the tapetum at the microspore stage and down-regulated after 1 d of high temperature. The expression levels of Osc6, OsRAFTIN and TDR, which are tapetum-specific genes, were unaffected by high temperatures. These results suggest that not all tapetal genes are inhibited by increased temperatures and the tapetum itself is not degraded in such an environment. However, high temperatures may disrupt some of the tapetum functions required for pollen adhesion and germination on the stigma.

Gene expression throughout the reproductive process in rice (Oryza sativa) beginning with primordia development through pollination/fertilization to zygote formation was analyzed. We analyzed 25 stages/organs of rice reproductive development including early microsporogenesis stages with 57,381 probe sets, and identified around 26,000 expressed probe sets in each stage. Fine dissection of 25 reproductive stages/organs combined with detailed microarray profiling revealed dramatic, coordinated and finely tuned changes in gene expression. A decrease in expressed genes in the pollen maturation process was observed in a similar way with Arabidopsis and maize. An almost equal number of ab initio predicted genes and cloned genes which appeared or disappeared coordinated with developmental stage progression. A large number of organ-/stage-specific genes were identified; notably 2,593 probe sets for developing anther, including 932 probe sets corresponding to ab initio predicted genes. Analysis of cell cycle-related genes revealed that several cyclin-dependent kinases (CDKs), cyclins and components of SCF E3 ubiquitin ligase complexes were expressed specifically in reproductive organs. Cell wall biosynthesis or degradation protein genes and transcription factor genes expressed specifically in reproductive stages were also newly identified. Rice genes homologous to reproduction-related genes in other plants showed expression profiles both consistent and inconsistent with their predicted functions. The rice reproductive expression atlas is likely to be the most extensive and most comprehensive data set available, indispensable for unraveling functions of many specific genes in plant reproductive processes that have not yet been thoroughly analyzed.

The male gametophyte and tapetum play different roles during anther development although they are differentiated from the same cell lineage, the L2 layer. Until now, it has not been possible to delineate their transcriptomes due to technical difficulties in separating the two cell types. In the present study, we characterized the separated transcriptomes of the rice microspore/pollen and tapetum using laser microdissection (LM)-mediated microarray. Spatiotemporal expression patterns of 28,141 anther-expressed genes were classified into 20 clusters, which contained 3,468 (12.3%) anther-enriched genes. In some clusters, synchronous gene expression in the microspore and tapetum at the same developmental stage was observed as a novel characteristic of the anther transcriptome. Noteworthy expression patterns are discussed in connection with gene ontology (GO) categories and gene annotations, which are related to important biological events in anther development, such as pollen maturation, pollen germination, pollen tube elongation and pollen wall formation.

Similarity of gene expression profiles provides important clues for understanding the biological functions of genes, biological processes and metabolic pathways related to genes. A gene expression network (GEN) is an ideal choice to grasp such expression profile similarities among genes simultaneously. For GEN construction, the Pearson correlation coefficient (PCC) has been widely used as an index to evaluate the similarities of expression profiles for gene pairs. However, calculation of PCCs for all gene pairs requires large amounts of both time and computer resources. Based on correspondence analysis, we developed a new method for GEN construction, which takes minimal time even for large-scale expression data with general computational circumstances. Moreover, our method requires no prior parameters to remove sample redundancies in the data set. Using the new method, we constructed rice GENs from large-scale microarray data stored in a public database. We then collected and integrated various principal rice omics annotations in public and distinct databases. The integrated information contains annotations of genome, transcriptome and metabolic pathways. We thus developed the integrated database OryzaExpress for browsing GENs with an interactive and graphical viewer and principal omics annotations (http://riceball.lab.nig.ac.jp/oryzaexpress/). With integration of Arabidopsis GEN data from ATTED-II, OryzaExpress also allows us to compare GENs between rice and Arabidopsis. Thus, OryzaExpress is a comprehensive rice database that exploits powerful omics approaches from all perspectives in plant science and leads to systems biology.

Cool temperature conditions are known to lead to pollen sterility in rice. Pollen sterility is an agriculturally important phenomenon because it imparts a large influence directly on rice yield. However, cool temperature stress tolerance varies among rice cultivars and avoidance of cool temperature stress is difficult by practical method of agriculture. In this study using two rice cultivars, Hitomebore (high tolerance) and Sasanishiki (low tolerance), we analyzed morphological features and gene expression profiles, under cool temperature stress, in anther development of rice. Hitomebore was given cool temperature stress (19°C) at flowering stage, and showed 87.3% seed fertility. Meanwhile, the seed fertility decreased to 41.7% in the case of Sasanishiki. A transverse section of Hitomebore anther revealed that the degradation of the tapetum started at the uninucleate microspore stage, and the tapetum had completely vanished at mature stage. The tapetum provides nutrients for pollen development, and its degradation occurs at a late stage in pollen development. In contrast, degradation of the tapetum did not occur at the uninucleate microspore stage in Sasanishiki, and the tapetum was clearly intact at mature stage, suggesting that tapetum degradation is critical for accurate pollen development and cool temperature tolerance correlated with the degree of tapetum degeneration. In gene expression analysis of anther, 356 genes that showed different expression levels between two cultivars at cool temperatures were found. These genes will lead to understanding the mechanism of cool temperature stress response in rice pollen development and the identification of genes involved in accurate tapetum degradation.

Many plant species are susceptible to high-temperature (HT) injury during reproductive development. We recently demonstrated that HT represses the expression of YUCCA auxin biosynthesis genes and reduces endogenous auxin in the developing anthers of barley and Arabidopsis. Here, we show that DNA proliferation in mitochondria, chloroplasts and nuclei of developing panicles is inhibited with increasing temperatures in barley. Following DNA proliferation suppression, terminal abnormalities were observed in the organelles of anther wall cells, including mitochondrial swelling and overdevelopment of chloroplasts. Comprehensive transcriptome analyses using both reproductive organs and vegetative tissues showed high and positive pairwise correlations between the expression profiles of auxin-induced genes, DNA replication-related genes and mitochondrial-related genes. In contrast, the expression profiles of auxin-repressed protein genes and photosynthesis-/chloroplast-related genes were negatively correlated with those of the previously mentioned genes. Under HT conditions, the former was repressed and the latter was up-regulated in the developing panicles. Furthermore, application of exogenous auxin promoted the expression of DNA replication-related genes under HT conditions, inducing anther cell proliferation. These suggest that compromised auxin biosynthesis/IAA level under HT condition results in nuclear and organellar DNA proliferation arrest due to co-transcriptional alterations.

Correlations of gene-to-gene co-expression and metabolite-to-metabolite co-accumulation calculated from large amounts of transcriptome and metabolome data are useful for uncovering unknown functions of genes, functional diversities of gene family members and regulatory mechanisms of metabolic pathway flows. Many databases and tools are available to interpret quantitative transcriptome and metabolome data, but there are only limited ones that connect correlation data to biological knowledge and can be utilized to find biological significance of it. We report here a new metabolic pathway database, KaPPA-View4 (http://kpv.kazusa.or.jp/kpv4/), which is able to overlay gene-to-gene and/or metabolite-to-metabolite relationships as curves on a metabolic pathway map, or on a combination of up to four maps. This representation would help to discover, for example, novel functions of a transcription factor that regulates genes on a metabolic pathway. Pathway maps of the Kyoto Encyclopedia of Genes and Genomes (KEGG) and maps generated from their gene classifications are available at KaPPA-View4 KEGG version (http://kpv.kazusa.or.jp/kpv4-kegg/). At present, gene co-expression data from the databases ATTED-II, COXPRESdb, CoP and MiBASE for human, mouse, rat, Arabidopsis, rice, tomato and other plants are available.