The Rice Oligonucleotide Array Database: an atlas of rice gene expression

The mitochondrion-located protein OsB12D1 enhances flooding tolerance during seed germination and early seedling growth in rice

B12D belongs to a function unknown subgroup of the Balem (Barley aleurone and embryo) proteins. In our previous work on rice seed germination, we identified a B12D-like protein encoded by LOC_Os7g41350 (named OsB12D1). OsB12D1 pertains to an ancient protein family with an amino acid sequence highly conserved from moss to angiosperms. Among the six OsB12Ds, OsB12D1 is one of the major transcripts and is primarily expressed in germinating seed and root. Bioinformatics analyses indicated that OsB12D1 is an anoxic or submergence resistance-related gene. RT-PCR results showed OsB12D1 is induced remarkably in the coleoptiles or roots by flooding during seed germination and early seedling growth. The OsB12D1-overexpressed rice seeds could protrude radicles in 8 cm deep water, further exhibiting significant flooding tolerance compared to the wild type. Moreover, this tolerance was not affected by the gibberellin biosynthesis inhibitor paclobutrazol. OsB12D1 was identified in the mitochondrion by subcellular localization analysis and possibly enhances electron transport through mediating Fe and oxygen availability under flooded conditions. This work indicated that OsB12D1 is a promising gene that can help to enhance rice seedling establishment in farming practices, especially for direct seeding.

Evaluation of rice promoters conferring pollen-specific expression in a heterologous system, Arabidopsis

Promoters can direct gene expression specifically to targeted tissues or cells. Effective with both crop species and model plant systems, these tools can help researchers overcome the practical obstacles associated with transgenic protocols. Here, we identified promoters that allow one to target the manipulation of gene expression during pollen development. Utilizing published transcriptomic databases for rice, we investigated the promoter activity of selected genes in Arabidopsis. From various microarray datasets, including those for anthers and pollen grains at different developmental stages, we selected nine candidate genes that showed high levels of expression in the late stages of rice pollen development. We named these Oryza sativa late pollen-specific genes. Their promoter regions contained various cis-acting elements that could be responsible for anther-/pollen-specific expression. Promoter::GUS-GFP reporters were constructed and introduced into Arabidopsis plants. Histochemical GUS staining revealed that six of the nine rice promoters conferred strong GUS expression that was restricted to the anthers in Arabidopsis. Further analysis showed that although the GUS signals were not detected at the unicellular stage, they strengthened in the bicellular or tricellular stages, peaking at the mature pollen stage. This paralleled their transcriptomic profiles in rice. Based on our results, we proposed that these six rice promoters, which are active in the late stages of pollen formation in the dicot Arabidopsis, can aid molecular breeders in generating new varieties of a monocot plant, rice.

Auxin and Cell Wall Invertase Related Signaling during Rice Grain Development

Indole-3-acetic acid (IAA) synthesis is required for grain-fill in maize and appears to be regulated by cell-wall invertase (CWIN) activity. OsYUC12 is one of three IAA biosynthesis genes we previously reported as expressed during early rice grain development, correlating with a large increase in IAA content of the grain. This work aimed to investigate further the role of OsYUC12 and its relationship to CWIN activity and invertase inhibitors (INVINH). The analysis shows a brief peak of OsYUC12 expression early in endosperm development. Meta-analysis of microarray data, confirmed by quantitative expression analysis, revealed that OsYUC12 is coexpressed with OsIAA29, which encodes an unusual AUX/IAA transcription factor previously reported as poorly expressed. Maximum expression of OsYUC12 and OsIAA29 coincided with maximum CWIN activity, but also with a peak in INVINH expression. Unlike ZmYUC1, OsYUC12 expression is not reduced in the rice CWIN mutant, gif1. Several reports have investigated CWIN expression in rice grains but none has reported on expression of INVINH in this species. We show that rice has 54 genes encoding putative invertase/pectin methylesterase inhibitors, seven of which are expressed exclusively during grain development. Our results suggest a more complex relationship between IAA, CWIN, and INVINH than previously proposed.

Uncovering legumain genes in rice

Legumains are Asn specific cysteine proteases physiologically related to the biosynthesis of vacuolar components, degradation of storage proteins and programmed cell death. The present work identifies and characterizes the genic family of legumains in rice (Oryza sativa), which comprises five different loci. Rice legumains (OsaLegs) were ubiquitously detected in all plant tissues analyzed. However, phylogenetic analyses and gene expression studies demonstrated greater association of OsaLeg2 and OsaLeg3 to seed-related legumains, whereas OsaLeg1, 4 and 5 would act as vegetative-related proteases. Additionally, OsaLeg1 mRNA is strongly induced in senescent leaves. All rice legumain genes respond in different ways to environmental conditions such as wounding, salt and abscisic acid treatments. Mainly, wounding is capable of inducing all the four expressed genes OsaLeg1, 2, 3 and 4. Alternative splicing isoforms, with potential to generate pre-activated OsaLeg1 and OsaLeg2 nonvacuolar enzymes under different environmental situations were also observed.

Quantitative trait locus mapping and candidate gene analysis for plant architecture traits using whole genome re-sequencing in rice

Plant breeders have focused on improving plant architecture as an effective means to increase crop yield. Here, we identify the main-effect quantitative trait loci (QTLs) for plant shape-related traits in rice (Oryza sativa) and find candidate genes by applying whole genome re-sequencing of two parental cultivars using next-generation sequencing. To identify QTLs influencing plant shape, we analyzed six traits: plant height, tiller number, panicle diameter, panicle length, flag leaf length, and flag leaf width. We performed QTL analysis with 178 F7 recombinant in-bred lines (RILs) from a cross of japonica rice line 'SNUSG1' and indica rice line 'Milyang23'. Using 131 molecular markers, including 28 insertion/deletion markers, we identified 11 main- and 16 minor-effect QTLs for the six traits with a threshold LOD value > 2.8. Our sequence analysis identified fifty-four candidate genes for the main-effect QTLs. By further comparison of coding sequences and meta-expression profiles between japonica and indica rice varieties, we finally chose 15 strong candidate genes for the 11 main-effect QTLs. Our study shows that the whole-genome sequence data substantially enhanced the efficiency of polymorphic marker development for QTL fine-mapping and the identification of possible candidate genes. This yields useful genetic resources for breeding high-yielding rice cultivars with improved plant architecture.

Transcriptomic analysis of UV-treated rice leaves reveals UV-induced phytoalexin biosynthetic pathways and their regulatory networks in rice

Rice produces diterpenoid and flavonoid phytoalexins for defense against pathogen attack. The production of phytoalexins in rice is also induced by UV-irradiation. To understand the metabolic networks involved in UV-induced phytoalexin biosynthesis and their regulation, phytochemical and transcriptomic analyses of UV-treated rice leaves were performed. In response to UV treatment, the accumulation of flavonoids was observed in rice leaves, which may serve as antioxidants against UV-induced oxidative stress. The phytochemical analysis confirmed sakuranetin accumulation and also demonstrated the induction of phenylamide synthesis in rice leaves by UV-irradiation. Transcriptomic analysis established that aromatic amino acid biosynthetic genes were immediately up-regulated after UV treatment. The genes involved in the phenylpropanoid pathway and flavonoid biosynthesis were also up-regulated. These findings suggest that the aromatic amino acid and flavonoid biosynthetic pathways are coordinately activated for the production of flavonoids and phenolic phytoalexins such as sakuranetin and phenylamides. An in silico analysis of UV-induced O-methyltransferase and acyltransferase genes suggested that these genes may be implicated in sakuranetin and phenylamide synthesis, respectively. The transcriptomic analysis also showed up-regulation of both methylerythritol phosphate pathway and the diterpenoid phytoalexin biosynthetic genes in response to UV treatment. A functional gene network analysis of phytoalexin biosynthetic and UV-induced genes for signaling components and transcription factors using RiceNet suggested that regulatory networks comprising signal perceiving receptor kinases, G-proteins, signal transducing mitogen-activated protein kinases and calcium signaling components, and transcription factors control flavonoid and phytoalexin biosynthesis in rice leaves under UV-C stress conditions.

RNA-sequencing reveals previously unannotated protein- and microRNA-coding genes expressed in aleurone cells of rice seeds

The rice genome annotation has been greatly improved in recent years, largely due to the availability of full length cDNA sequences derived from many tissues. Among those yet to be studied is the aleurone layer, which produces hydrolases for mobilization of seed storage reserves during seed germination and post germination growth. Herein, we report transcriptomes of aleurone cells treated with the hormones abscisic acid, gibberellic acid, or both. Using a comprehensive approach, we identified hundreds of novel genes. To minimize the number of false positives, only transcripts that did not overlap with existing annotations, had a high level of expression, and showed a high level of uniqueness within the rice genome were considered to be novel genes. This approach led to the identification of 553 novel genes that encode proteins and/or microRNAs. The transcriptome data reported here will help to further improve the annotation of the rice genome.

Molecular characterization and functional analysis of "fruit-weight 2.2-like" gene family in rice

Tomato fruit-weight 2.2 (FW2.2) was reported to control up to 30 % fruit weight. Recent studies demonstrated that FW2.2-like (FWL) genes also play important roles in plant growth and development. For instance, a maize homolog of FW2.2, named cell number regulator 1 (CNR1), negatively regulates plant and organ size. However, FWL genes in rice have not been characterized yet. In this study, eight FWL genes were identified in rice genome and designated as OsFWL1-8. The chromosome location, gene structure, protein motif, and phylogenetic relationship of OsFWL genes were analyzed. RT-PCR result and microarray data revealed that OsFWL genes exhibited diverse expression patterns and the detailed expression patterns of OsFWL5, 6, and 7 negatively correlated with leaf growth activity. Rice protoplast transient transformation experiment showed that most OsFWL proteins locate at cell membrane but OsFWL8 is present in the nucleus. In addition, the functions of OsFWL genes were investigated by analyzing two T-DNA insertion lines for OsFWL3 and 5. Compared with wild type, the grain weight of osfwl3 mutant and the plant height of osfwl5 mutant were increased by 5.3 and 12.5 %, respectively. We also found that the increase in grain length of osfwl3 mutant was due chiefly to incremental cell number, not cell size and the expression of OsFWL3 negatively correlated with glume growth activity. These results provide a comprehensive foundation for further study of OsFWL functions in rice.

A systematic view of rice heat shock transcription factor family using phylogenomic analysis

The heat shock transcription factor (Hsf) family includes key regulators of the physiological response to heat stress. Here, we present a systematic analysis of the Hsf family in rice using a phylogenomics-based approach that integrates multi-omics data into the context of a phylogenetic tree. For 25 previously identified Hsfs, we integrated anatomical meta-profiling data from 983 Affymetrix arrays into a phylogenetic tree, revealing a global view of the functional redundancy within this family. Interestingly, most of the Hsfs showed significant fluctuation in gene expression patterns, suggesting that they have condition- or stress-dependent roles. Therefore, we further analyzed the abiotic stress responses of the Hsfs using log(2-)fold change data in response to heat, cold, drought and salt stresses. Subsequently, we identified 19 Hsfs that are positively associated with heat stress, 11 with drought, 9 with salt, and 7 with cold stress, as indicated by at least a 2-fold change and coefficient of variation less than 1. The Hsf subfamily A2 was conserved in the heat stress response. The Hsf subfamily C showed a strong positive association with drought, salt and cold stresses. Downregulation of three members in the Hsf subfamily B in response to cold stress is characteristic. More interestingly, half of the Hsf subfamily B genes were upregulated by heat, drought and salt stresses, while one gene in the other half was downregulated by drought, salt, and cold stresses. Finally, we developed a hypothetical functional gene network mediated by OsHsfA2e/OsHsf-12 that is involved in thermotolerance as well as upregulated in response to heat. We expect that our data will help researchers design more efficient strategies to study the rice Hsf family with information about probable functional redundancy.

Construction of a rice glycoside hydrolase phylogenomic database and identification of targets for biofuel research

Glycoside hydrolases (GH) catalyze the hydrolysis of glycosidic bonds in cell wall polymers and can have major effects on cell wall architecture. Taking advantage of the massive datasets available in public databases, we have constructed a rice phylogenomic database of GHs (http://ricephylogenomics.ucdavis.edu/cellwalls/gh/). This database integrates multiple data types including the structural features, orthologous relationships, mutant availability, and gene expression patterns for each GH family in a phylogenomic context. The rice genome encodes 437 GH genes classified into 34 families. Based on pairwise comparison with eight dicot and four monocot genomes, we identified 138 GH genes that are highly diverged between monocots and dicots, 57 of which have diverged further in rice as compared with four monocot genomes scanned in this study. Chromosomal localization and expression analysis suggest a role for both whole-genome and localized gene duplications in expansion and diversification of GH families in rice. We examined the meta-profiles of expression patterns of GH genes in twenty different anatomical tissues of rice. Transcripts of 51 genes exhibit tissue or developmental stage-preferential expression, whereas, seventeen other genes preferentially accumulate in actively growing tissues. When queried in RiceNet, a probabilistic functional gene network that facilitates functional gene predictions, nine out of seventeen genes form a regulatory network with the well-characterized genes involved in biosynthesis of cell wall polymers including cellulose synthase and cellulose synthase-like genes of rice. Two-thirds of the GH genes in rice are up regulated in response to biotic and abiotic stress treatments indicating a role in stress adaptation. Our analyses identify potential GH targets for cell wall modification.

RiceFREND: a platform for retrieving coexpressed gene networks in rice

Similarity of gene expression across a wide range of biological conditions can be efficiently used in characterization of gene function. We have constructed a rice gene coexpression database, RiceFREND (http://ricefrend.dna.affrc.go.jp/), to identify gene modules with similar expression profiles and provide a platform for more accurate prediction of gene functions. Coexpression analysis of 27 201 genes was performed against 815 microarray data derived from expression profiling of various organs and tissues at different developmental stages, mature organs throughout the growth from transplanting until harvesting in the field and plant hormone treatment conditions, using a single microarray platform. The database is provided with two search options, namely, 'single guide gene search' and 'multiple guide gene search' to efficiently retrieve information on coexpressed genes. A user-friendly web interface facilitates visualization and interpretation of gene coexpression networks in HyperTree, Cytoscape Web and Graphviz formats. In addition, analysis tools for identification of enriched Gene Ontology terms and cis-elements provide clue for better prediction of biological functions associated with the coexpressed genes. These features allow users to clarify gene functions and gene regulatory networks that could lead to a more thorough understanding of many complex agronomic traits.

RiceXPro version 3.0: expanding the informatics resource for rice transcriptome

A wide range of resources on gene expression profiling enhance various strategies in plant molecular biology particularly in characterization of gene function. We have updated our gene expression profile database, RiceXPro (http://ricexpro.dna.affrc.go.jp/), to provide more comprehensive information on the transcriptome of rice encompassing the entire growth cycle and various experimental conditions. The gene expression profiles are currently grouped into three categories, namely, ‘field/development’ with 572 data corresponding to 12 data sets, ‘plant hormone’ with 143 data corresponding to 13 data sets and ‘cell- and tissue-type’ comprising of 38 microarray data. In addition to the interface for retrieving expression information of a gene/genes in each data set, we have incorporated an interface for a global approach in searching an overall view of the gene expression profiles from multiple data sets within each category. Furthermore, we have also added a BLAST search function that enables users to explore expression profile of a gene/genes with similarity to a query sequence. Therefore, the updated version of RiceXPro can be used more efficiently to survey the gene expression signature of rice in sufficient depth and may also provide clues on gene function of other cereal crops.