Search

Root meristem activity is the most critical process influencing root development. Although several factors that regulate meristem activity have been identified in rice, studies on the enhancement of meristem activity in roots are limited. We identified a T-DNA activation tagging line of a zinc-finger homeobox gene, OsZHD2, which has longer seminal and lateral roots due to increased meristem activity. The phenotypes were confirmed in transgenic plants overexpressing OsZHD2. In addition, the overexpressing plants showed enhanced grain yield under low nutrient and paddy field conditions. OsZHD2 was preferentially expressed in the shoot apical meristem and root tips. Transcriptome analyses and quantitative real-time PCR experiments on roots from the activation tagging line and the wild type showed that genes for ethylene biosynthesis were up-regulated in the activation line. Ethylene levels were higher in the activation lines compared with the wild type. ChIP assay results suggested that OsZHD2 induces ethylene biosynthesis by controlling ACS5 directly. Treatment with ACC (1-aminocyclopropane-1-carboxylic acid), an ethylene precursor, induced the expression of the DR5 reporter at the root tip and stele, whereas treatment with an ethylene biosynthesis inhibitor, AVG (aminoethoxyvinylglycine), decreased that expression in both the wild type and the OsZHD2 overexpression line. These observations suggest that OsZHD2 enhances root meristem activity by influencing ethylene biosynthesis and, in turn, auxin.

Transcription factors (TFs) are an important class of regulatory molecules. Despite their importance, only a small number of genes encoding TFs have been characterized in Oryza sativa (rice), often because gene duplication and functional redundancy complicate their analysis. To address this challenge, we developed a web-based tool called the Rice Transcription Factor Phylogenomics Database (RTFDB) and demonstrate its application for predicting TF function. The RTFDB hosts transcriptome and co-expression analyses. Sources include high-throughput data from oligonucleotide microarray (Affymetrix and Agilent) as well as RNA-Seq-based expression profiles. We used the RTFDB to identify tissue-specific and stress-related gene expression. Subsequently, 273 genes preferentially expressed in specific tissues or organs, 455 genes showing a differential expression pattern in response to 4 abiotic stresses, 179 genes responsive to infection of various pathogens and 512 genes showing differential accumulation in response to various hormone treatments were identified through the meta-expression analysis. Pairwise Pearson correlation coefficient analysis between paralogous genes in a phylogenetic tree was used to assess their expression collinearity and thereby provides a hint on their genetic redundancy. Integrating transcriptome with the gene evolutionary information reveals the possible functional redundancy or dominance played by paralog genes in a highly duplicated genome such as rice. With this method, we estimated a predominant role for 83.3% (65/78) of the TF or transcriptional regulator genes that had been characterized via loss-of-function studies. In this regard, the proposed method is applicable for functional studies of other plant species with annotated genome.

Flooding stress caused by climate change is a serious threat to crop productivity. To enhance our understanding of flooding stress in soybean, we analyzed the transcriptome of the roots of soybean plants after waterlogging treatment for 10 days at the V2 growth stage. Through RNA sequencing analysis, 870 upregulated and 1129 downregulated differentially expressed genes (DEGs) were identified and characterized using Gene Ontology (GO) and MapMan software (version 3.6.0RC1). In the functional classification analysis, "alcohol biosynthetic process" was the most significantly enriched GO term in downregulated DEGs, and phytohormone-related genes such as ABA, cytokinin, and gibberellin were upregulated. Among the transcription factors (TFs) in DEGs, AP2/ERFs were the most abundant. Furthermore, our DEGs encompassed eight soybean orthologs from Arabidopsis and rice, such as 1-aminocyclopropane-1-carboxylate oxidase. Along with a co-functional network consisting of the TF and orthologs, the expression changes of those genes were tested in a waterlogging-resistant cultivar, PI567343. These findings contribute to the identification of candidate genes for waterlogging tolerance in soybean, which can enhance our understanding of waterlogging tolerance. [ABSTRACT FROM AUTHOR]

BackgroundProtein kinases catalyze the transfer of a phosphate moiety from a phosphate donor to the substrate molecule, thus playing critical roles in cell signaling and metabolism. Although plant genomes contain more than 1000 genes that encode kinases, knowledge is limited about the function of each of these kinases. A major obstacle that hinders progress towards kinase characterization is functional redundancy. To address this challenge, we previously developed the rice kinase database (RKD) that integrated omics-scale data within a phylogenetics context.ResultsAn updated version of rice kinase database (RKD) that contains metadata derived from NCBI GEO expression datasets has been developed. RKD 2.0 facilitates in-depth transcriptomic analyses of kinase-encoding genes in diverse rice tissues and in response to biotic and abiotic stresses and hormone treatments. We identified 261 kinases specifically expressed in particular tissues, 130 that are significantly up- regulated in response to biotic stress, 296 in response to abiotic stress, and 260 in response to hormones. Based on this update and Pearson correlation coefficient (PCC) analysis, we estimated that 19 out of 26 genes characterized through loss-of-function studies confer dominant functions. These were selected because they either had paralogous members with PCC values of

The completion of whole genome sequence of rice (Oryza sativa) has significantly accelerated functional genomics studies. Prior to the release of the sequence, only a few genes were assigned a function each year. Since sequencing was completed in 2005, the rate has exponentially increased. As of 2014, 1,021 genes have been described and added to the collection at The Overview of functionally characterized Genes in Rice online database (OGRO). Despite this progress, that number is still very low compared with the total number of genes estimated in the rice genome. One limitation to progress is the presence of functional redundancy among members of the same rice gene family, which covers 51.6 % of all non-transposable element-encoding genes. There remain a significant portion or rice genes that are not functionally redundant, as reflected in the recovery of loss-of-function mutants. To more accurately analyze functional redundancy in the rice genome, we have developed a phylogenomics databases for six large gene families in rice, including those for glycosyltransferases, glycoside hydrolases, kinases, transcription factors, transporters, and cytochrome P450 monooxygenases. In this review, we introduce key features and applications of these databases. We expect that they will serve as a very useful guide in the post-genomics era of research.

Cross-species translation of genomic information may play a pivotal role in applying biological knowledge gained from relatively simple model system to other less studied, but related, genomes. The information of abiotic stress (ABS)-responsive genes in Arabidopsis was identified and translated into the legume model system, Medicago truncatula. Various data resources, such as TAIR/AtGI DB, expression profiles and literatures, were used to build a genome-wide list of ABS genes. tBlastX/BlastP similarity search tools and manual inspection of alignments were used to identify orthologous genes between the two genomes. A total of 1,377 genes were finally collected and classified into 18 functional criteria of gene ontology (GO). The data analysis according to the expression cues showed that there was substantial level of interaction among three major types (i.e., drought, salinity and cold stress) of abiotic stresses. In an attempt to translate the ABS genes between these two species, genomic locations for each gene were mapped using an in-house-developed comparative analysis platform. The comparative analysis revealed that fragmental colinearity, represented by only 37 synteny blocks, existed between Arabidopsis and M. truncatula. Based on the combination of E-value and alignment remarks, estimated translation rate was 60.2% for this cross-family translation. As a prelude of the functional comparative genomic approaches, in-silico gene network/interactome analyses were conducted to predict key components in the ABS responses, and one of the sub-networks was integrated with corresponding comparative map. The results demonstrated that core members of the sub-network were well aligned with previously reported ABS regulatory networks. Taken together, the results indicate that network-based integrative approaches of comparative and functional genomics are important to interpret and translate genomic information for complex traits such as abiotic stresses.

Inorganic phosphate (Pi) homeostasis plays an important role in plant growth and abiotic stress tolerance. Several MYB-CC transcription factors involved in Pi homeostasis have been identified in rice (Oryza sativa). PHOSPHATE STARVATION RESPONSE-LIKE 7 (PHL7) is a class II MYC-CC protein, in which the MYC-CC domain is located at the N terminus. In this study, we established that OsPHL7 is localized to the nucleus and that the encoding gene is induced by Pi deficiency. The Pi-responsive genes and Pi transporter genes are positively regulated by OsPHL7. The overexpression of OsPHL7 enhanced the tolerance of rice plants to Pi starvation, whereas the RNA interference-based knockdown of this gene resulted in increased sensitivity to Pi deficiency. Transgenic rice plants overexpressing OsPHL7 produced more roots than wild-type plants under both Pi-sufficient and Pi-deficient conditions and accumulated more Pi in the shoots and roots. In addition, the overexpression of OsPHL7 enhanced rice tolerance to salt stress. Together, these results demonstrate that OsPHL7 is involved in the maintenance of Pi homeostasis and enhances tolerance to Pi deficiency and salt stress in rice. [ABSTRACT FROM AUTHOR]

Rice is an important cereal crop worldwide, the growth of which is affected by rice blast disease, caused by the fungal pathogen Magnaporthe oryzae. As climate change increases the diversity of pathogens, the disease resistance genes (R genes) in plants must be identified. The major blast-resistance genes have been identified in indica rice varieties; therefore, japonica rice varieties with R genes now need to be identified. Because leucine-rich repeat (LRR) domain proteins possess R-gene properties, we used bioinformatics analysis to identify the rice candidate LRR domain receptor-like proteins (OsLRR-RLPs). OsLRR-RLP2, which contains six LRR domains, showed differences in the DNA sequence, containing 43 single-nucleotide polymorphisms (SNPs) in indica and japonica subpopulations. The results of the M. oryzae inoculation analysis indicated that indica varieties with partial deletion of OsLRR-RLP2 showed susceptibility, whereas japonica varieties with intact OsLRR-RLP2 showed resistance. The oslrr-rlp2 mutant, generated using clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), showed increased pathogen susceptibility, whereas plants overexpressing this gene showed pathogen resistance. These results indicate that OsLRR-RLP2 confers resistance to rice, and OsLRR-RLP2 may be useful for breeding resistant cultivars. [ABSTRACT FROM AUTHOR]

Phosphate (Pi) starvation is a critical factor limiting crop growth, development, and productivity. Rice (Oryza sativa) R2R3-MYB transcription factors function in the transcriptional regulation of plant responses to various abiotic stresses and micronutrient deprivation, but little is known about their roles in Pi starvation signaling and Pi homeostasis. Here, we identified the R2R3-MYB transcription factor gene OsMYB58, which shares high sequence similarity with AtMYB58. OsMYB58 expression was induced more strongly by Pi starvation than by other micronutrient deficiencies. Overexpressing OsMYB58 in Arabidopsis thaliana and rice inhibited plant growth and development under Pi-deficient conditions. In addition, the overexpression of OsMYB58 in plants exposed to Pi deficiency strongly affected root development, including seminal root, lateral root, and root hair formation. Overexpressing OsMYB58 strongly decreased the expression of the rice microRNAs OsmiR399a and OsmiR399j. By contrast, overexpressing OsMYB58 strongly increased the expression of rice PHOSPHATE 2 (OsPHO2), whose expression is repressed by miR399 during Pi starvation signaling. OsMYB58 functions as a transcriptional repressor of the expression of its target genes, as determined by a transcriptional activity assay. These results demonstrate that OsMYB58 negatively regulates OsmiR399-dependent Pi starvation signaling by enhancing OsmiR399s expression. [ABSTRACT FROM AUTHOR]

3D-printed bone scaffolds hold great promise for the individualized treatment of critical-size bone defects. Among the resorbable polymers available for use as 3D-printable scaffold materials, poly(ε-caprolactone) (PCL) has many benefits. However, its relatively low stiffness and lack of bioactivity limit its use in load-bearing bone scaffolds. This study tests the hypothesis that surface-oxidized cellulose nanocrystals (SO-CNCs), decorated with carboxyl groups, can act as multi-functional scaffold additives that (1) improve the mechanical properties of PCL and (2) induce biomineral formation upon PCL resorption. To this end, an in vitro biomineralization study was performed to assess the ability of SO-CNCs to induce the formation of calcium phosphate minerals. In addition, PCL nanocomposites containing different amounts of SO-CNCs (1, 2, 3, 5, and 10 wt%) were prepared using melt compounding extrusion and characterized in terms of Young's modulus, ultimate tensile strength, crystallinity, thermal transitions, and water contact angle. Neither sulfuric acid-hydrolyzed CNCs (SH-CNCs) nor SO-CNCs were toxic to MC3T3 preosteoblasts during a 24 h exposure at concentrations ranging from 0.25 to 3.0 mg/mL. SO-CNCs were more effective at inducing mineral formation than SH-CNCs in simulated body fluid (1x). An SO-CNC content of 10 wt% in the PCL matrix caused a more than 2-fold increase in Young's modulus (stiffness) and a more than 60% increase in ultimate tensile strength. The matrix glass transition and melting temperatures were not affected by the SO-CNCs but the crystallization temperature increased by about 5.5°C upon addition of 10 wt% SO-CNCs, the matrix crystallinity decreased from about 43 to about 40%, and the water contact angle decreased from 87 to 82.6°. The abilities of SO-CNCs to induce calcium phosphate mineral formation and increase the Young's modulus of PCL render them attractive for applications as multi-functional nanoscale additives in PCL-based bone scaffolds.

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.

Microarray technologies facilitate high-throughput gene expression analysis. However, the diversity of platforms for rice gene expression analysis hinders efficient analysis. Tools to broadly integrate microarray data from different platforms are needed.In this study, we developed the Rice Oligonucleotide Array Database (ROAD, http://www.ricearray.org) to explore gene expression across 1,867 publicly available rice microarray hybridizations. The ROAD's user-friendly web interface and variety of visualization tools facilitate the extraction of gene expression profiles using gene and microarray element identifications. The ROAD supports meta-analysis of genes expressed in different tissues and at developmental stages. Co-expression analysis tool provides information on co-regulation between genes under general, abiotic and biotic stress conditions. Additionally, functional analysis tools, such as Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthology, are embedded in the ROAD. These tools facilitate the identification of meaningful biological patterns in a list of query genes.The Rice Oligonucleotide Array Database provides comprehensive gene expression profiles for all rice genes, and will be a useful resource for researchers of rice and other grass species.

Efficient and cost-effective conversion of plant biomass to usable forms of energy requires a thorough understanding of cell wall biosynthesis, modification and degradation. To elucidate these processes, we assessed the expression dynamics during enzymatic removal and regeneration of rice cell walls in suspension cells over time. In total, 928 genes exhibited significant up-regulation during cell wall removal, whereas, 79 genes were up-regulated during cell wall regeneration. Both gene sets are enriched for kinases, transcription factors and genes predicted to be involved in cell wall-related functions. Integration of the gene expression datasets with a catalog of known and/or predicted biochemical pathways from rice, revealed metabolic and hormonal pathways involved in cell wall degradation and regeneration. Rice lines carrying Tos17 mutations in genes up-regulated during cell wall removal exhibit dwarf phenotypes. Many of the genes up-regulated during cell wall development are also up-regulated in response to infection and environmental perturbations indicating a coordinated response to diverse types of stress.

Rice (Oryza sativa) is a staple food for more than half the world and a model for studies of monocotyledonous species, which include cereal crops and candidate bioenergy grasses. A major limitation of crop production is imposed by a suite of abiotic and biotic stresses resulting in 30%-60% yield losses globally each year. To elucidate stress response signaling networks, we constructed an interactome of 100 proteins by yeast two-hybrid (Y2H) assays around key regulators of the rice biotic and abiotic stress responses. We validated the interactome using protein-protein interaction (PPI) assays, co-expression of transcripts, and phenotypic analyses. Using this interactome-guided prediction and phenotype validation, we identified ten novel regulators of stress tolerance, including two from protein classes not previously known to function in stress responses. Several lines of evidence support cross-talk between biotic and abiotic stress responses. The combination of focused interactome and systems analyses described here represents significant progress toward elucidating the molecular basis of traits of agronomic importance.

Determination of gene function is particularly problematic when studying large-gene families because redundancy limits the ability to assess the contributions of individual genes experimentally. Phylogenomics is a phylogenetic approach used in comparative genomics to predict the biological functions of members of large gene-families by assessing the similarity among gene products. In this report, we describe the application of the Rice Kinase Database for elucidating functions of individual members of this gene family.

Bacterial leaf blight (BLB), caused by Xanthomonas oryzae pv. oryzae (Xoo), gives rise to devastating crop losses in rice. Disease resistant rice cultivars are the most economical way to combat the disease. The TP309 cultivar is susceptible to infection by Xoo strain PXO99. A transgenic variety, TP309_Xa21, expresses the pattern recognition receptor Xa21, and is resistant. PXO99 big up tri, openraxST, a strain lacking the raxST gene, is able to overcome Xa21-mediated immunity. We used a single extraction solvent to demonstrate comprehensive metabolomics and transcriptomics profiling under sample limited conditions, and analyze the molecular responses of two rice lines challenged with either PXO99 or PXO99 big up tri, openraxST. LC-TOF raw data file filtering resulted in better within group reproducibility of replicate samples for statistical analyses. Accurate mass match compound identification with molecular formula generation (MFG) ranking of 355 masses was achieved with the METLIN database. GC-TOF analysis yielded an additional 441 compounds after BinBase database processing, of which 154 were structurally identified by retention index/MS library matching. Multivariate statistics revealed that the susceptible and resistant genotypes possess distinct profiles. Although few mRNA and metabolite differences were detected in PXO99 challenged TP309 compared to mock, many differential changes occurred in the Xa21-mediated response to PXO99 and PXO99 big up tri, openraxST. Acetophenone, xanthophylls, fatty acids, alkaloids, glutathione, carbohydrate and lipid biosynthetic pathways were affected. Significant transcriptional induction of several pathogenesis related genes in Xa21 challenged strains, as well as differential changes to GAD, PAL, ICL1 and Glutathione-S-transferase transcripts indicated limited correlation with metabolite changes under single time point global profiling conditions. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-010-0218-7) contains supplementary material, which is available to authorized users.

Bacterial leaf blight (BLB), caused by Xanthomonas oryzae pv. oryzae (Xoo), gives rise to devastating crop losses in rice. Disease resistant rice cultivars are the most economical way to combat the disease. The TP309 cultivar is susceptible to infection by Xoo strain PXO99. A transgenic variety, TP309_Xa21, expresses the pattern recognition receptor Xa21, and is resistant. PXO99△raxST, a strain lacking the raxST gene, is able to overcome Xa21-mediated immunity. We used a single extraction solvent to demonstrate comprehensive metabolomics and transcriptomics profiling under sample limited conditions, and analyze the molecular responses of two rice lines challenged with either PXO99 or PXO99△raxST. LC–TOF raw data file filtering resulted in better within group reproducibility of replicate samples for statistical analyses. Accurate mass match compound identification with molecular formula generation (MFG) ranking of 355 masses was achieved with the METLIN database. GC–TOF analysis yielded an additional 441 compounds after BinBase database processing, of which 154 were structurally identified by retention index/MS library matching. Multivariate statistics revealed that the susceptible and resistant genotypes possess distinct profiles. Although few mRNA and metabolite differences were detected in PXO99 challenged TP309 compared to mock, many differential changes occurred in the Xa21-mediated response to PXO99 and PXO99△raxST. Acetophenone, xanthophylls, fatty acids, alkaloids, glutathione, carbohydrate and lipid biosynthetic pathways were affected. Significant transcriptional induction of several pathogenesis related genes in Xa21 challenged strains, as well as differential changes to GAD, PAL, ICL1 and Glutathione-S-transferase transcripts indicated limited correlation with metabolite changes under single time point global profiling conditions.

Biofuels derived from plant cell wall lignocellulose have the potential to serve as an alternative source of energy, relieving dependence on finite petroleum reserves and reducing production of climate-changing greenhouse gases. To better elucidate cell wall structure, the plant research community has developed databases to host the accumulated information on plant cell wall-related enzymes. The goal of this review is to provide a comprehensive catalog of these databases, as well as to describe the data and tools that they make available. A Directory of Databases for Plant Cell Wall-Related Enzymes (plantcellwalls.ucdavis.edu) was also developed to provide links to all the databases reviewed in this paper and relevant publications.

We previously characterized the rice (Oryza sativa) Submergence1 (Sub1) locus encoding three ethylene-responsive factor (ERF) transcriptional regulators. Genotypes carrying the Sub1A-1 allele are tolerant of prolonged submergence. To elucidate the mechanism of Sub1A-1-mediated tolerance, we performed transcriptome analyses comparing the temporal submergence response of Sub1A-1-containing tolerant M202(Sub1) with the intolerant isoline M202 lacking this gene. We identified 898 genes displaying Sub1A-1-dependent regulation. Integration of the expression data with publicly available metabolic pathway data identified submergence tolerance-associated pathways governing anaerobic respiration, hormone responses, and antioxidant systems. Of particular interest were a set of APETALA2 (AP2)/ERF family transcriptional regulators that are associated with the Sub1A-1-mediated response upon submergence. Visualization of expression patterns of the AP2/ERF superfamily members in a phylogenetic context resolved 12 submergence-regulated AP2/ERFs into three putative functional groups: (1) anaerobic respiration and cytokinin-mediated delay in senescence via ethylene accumulation during submergence (three ERFs); (2) negative regulation of ethylene-dependent gene expression (five ERFs); and (3) negative regulation of gibberellin-mediated shoot elongation (four ERFs). These results confirm that the presence of Sub1A-1 impacts multiple pathways of response to submergence.

Three subfamilies of grasses, the Erhardtoideae (rice), the Panicoideae (maize, sorghum, sugar cane and millet), and the Pooideae (wheat, barley and cool season forage grasses) provide the basis of human nutrition and are poised to become major sources of renewable energy. Here we describe the complete genome sequence of the wild grass Brachypodium distachyon (Brachypodium), the first member of the Pooideae subfamily to be completely sequenced. Comparison of the Brachypodium, rice and sorghum genomes reveals a precise sequence- based history of genome evolution across a broad diversity of the grass family and identifies nested insertions of whole chromosomes into centromeric regions as a predominant mechanism driving chromosome evolution in the grasses. The relatively compact genome of Brachypodium is maintained by a balance of retroelement replication and loss. The complete genome sequence of Brachypodium, coupled to its exceptional promise as a model system for grass research, will support the development of new energy and food crops

Studies of gene function are often hampered by gene-redundancy, especially in organisms with large genomes such as rice (Oryza sativa). We present an approach for using transcriptomics data to focus functional studies and address redundancy. To this end, we have constructed and validated an inexpensive and publicly available rice oligonucleotide near-whole genome array, called the rice NSF45K array. We generated expression profiles for light- vs. dark-grown rice leaf tissue and validated the biological significance of the data by analyzing sources of variation and confirming expression trends with reverse transcription polymerase chain reaction. We examined trends in the data by evaluating enrichment of gene ontology terms at multiple false discovery rate thresholds. To compare data generated with the NSF45K array with published results, we developed publicly available, web-based tools (www.ricearray.org). The Oligo and EST Anatomy Viewer enables visualization of EST-based expression profiling data for all genes on the array. The Rice Multi-platform Microarray Search Tool facilitates comparison of gene expression profiles across multiple rice microarray platforms. Finally, we incorporated gene expression and biochemical pathway data to reduce the number of candidate gene products putatively participating in the eight steps of the photorespiration pathway from 52 to 10, based on expression levels of putatively functionally redundant genes. We confirmed the efficacy of this method to cope with redundancy by correctly predicting participation in photorespiration of a gene with five paralogs. Applying these methods will accelerate rice functional genomics.

Glycosyltransferases (GTs; EC 2.4.x.y) constitute a large group of enzymes that form glycosidic bonds through transfer of sugars from activated donor molecules to acceptor molecules. GTs are critical to the biosynthesis of plant cell walls, among other diverse functions. Based on the Carbohydrate-Active enZymes (CAZy) database and sequence similarity searches, we have identified 609 potential GT genes (loci) corresponding to 769 transcripts (gene models) in rice (Oryza sativa), the reference monocotyledonous species. Using domain composition and sequence similarity, these rice GTs were classified into 40 CAZy families plus an additional unknown class. We found that two Pfam domains of unknown function, PF04577 and PF04646, are associated with GT families GT61 and GT31, respectively. To facilitate functional analysis of this important and large gene family, we created a phylogenomic Rice GT Database (http://ricephylogenomics.ucdavis.edu/cellwalls/gt/). Through the database, several classes of functional genomic data, including mutant lines and gene expression data, can be displayed for each rice GT in the context of a phylogenetic tree, allowing for comparative analysis both within and between GT families. Comprehensive digital expression analysis of public gene expression data revealed that most ( approximately 80%) rice GTs are expressed. Based on analysis with Inparanoid, we identified 282 'rice-diverged' GTs that lack orthologs in sequenced dicots (Arabidopsis thaliana, Populus tricocarpa, Medicago truncatula, and Ricinus communis). Combining these analyses, we identified 33 rice-diverged GT genes (45 gene models) that are highly expressed in above-ground, vegetative tissues. From the literature and this analysis, 21 of these loci are excellent targets for functional examination toward understanding and manipulating grass cell wall qualities. Study of the remainder may reveal aspects of hormone and protein metabolism that are critical for rice biology. This list of 33 genes and the Rice GT Database will facilitate the study of GTs and cell wall synthesis in rice and other plants.

Functional redundancy limits detailed analysis of genes in many organisms. Here, we report a method to efficiently overcome this obstacle by combining gene expression data with analysis of gene-indexed mutants. Using a rice NSF45K oligo-microarray to compare 2-week-old light- and dark-grown rice leaf tissue, we identified 365 genes that showed significant 8-fold or greater induction in the light relative to dark conditions. We then screened collections of rice T-DNA insertional mutants to identify rice lines with mutations in the strongly light-induced genes. From this analysis, we identified 74 different lines comprising two independent mutant lines for each of 37 light-induced genes. This list was further refined by mining gene expression data to exclude genes that had potential functional redundancy due to co-expressed family members (12 genes) and genes that had inconsistent light responses across other publicly available microarray datasets (five genes). We next characterized the phenotypes of rice lines carrying mutations in ten of the remaining candidate genes and then carried out co-expression analysis associated with these genes. This analysis effectively provided candidate functions for two genes of previously unknown function and for one gene not directly linked to the tested biochemical pathways. These data demonstrate the efficiency of combining gene family-based expression profiles with analyses of insertional mutants to identify novel genes and their functions, even among members of multi-gene families.

Xanthomonas oryzae pv. oryzae (Xoo) and X. oryzae pv. oryzicola (Xoc) are bacterial pathogens of the worldwide staple and grass model, rice. Xoo and Xoc are closely related but Xoo invades rice vascular tissue to cause bacterial leaf blight, a serious disease of rice in many parts of the world, and Xoc colonizes the mesophyll parenchyma to cause bacterial leaf streak, a disease of emerging importance. Both pathogens depend on hrp genes for type III secretion to infect their host. We constructed a 50-70 mer oligonucleotide microarray based on available genome data for Xoo and Xoc and compared gene expression in Xoo strains PXO99A and Xoc strain BLS256 grown in the rich medium PSB vs. XOM2, a minimal medium previously reported to induce hrp genes in Xoo strain T7174.Three biological replicates of the microarray experiment to compare global gene expression in representative strains of Xoo and Xoc grown in PSB vs. XOM2 were carried out. The non-specific error rate and the correlation coefficients across biological replicates and among duplicate spots revealed that the microarray data were robust. 247 genes of Xoo and 39 genes of Xoc were differentially expressed in the two media with a false discovery rate of 5% and with a minimum fold-change of 1.75. Semi-quantitative-RT-PCR assays confirmed differential expression of each of 16 genes each for Xoo and Xoc selected for validation. The differentially expressed genes represent 17 functional categories.We describe here the construction and validation of a two-genome microarray for the two pathovars of X. oryzae. Microarray analysis revealed that using representative strains, a greater number of Xoo genes than Xoc genes are differentially expressed in XOM2 relative to PSB, and that these include hrp genes and other genes important in interactions with rice. An exception was the rax genes, which are required for production of the host resistance elicitor AvrXa21, and which were expressed constitutively in both pathovars.

Rice, one of the most important food crops for humans, is the first crop plant to have its genome sequenced. Rice whole-genome microarrays, genome tiling arrays and genome-wide gene-indexed mutant collections have recently been generated. With the availability of these resources, discovering the function of the estimated 41,000 rice genes is now within reach. Such discoveries have broad practical implications for understanding the biological processes of rice and other economically important grasses such as cereals and bioenergy crops.