Evolution of the Grain Dispersal System in Barley.

About 12,000 years ago in the Near East, humans began the transition from hunter-gathering to agriculture-based societies. Barley was a founder crop in this process, and the most important steps in its domestication were mutations in two adjacent, dominant, and complementary genes, through which grains were retained on the inflorescence at maturity, enabling effective harvesting. Independent recessive mutations in each of these genes caused cell wall thickening in a highly specific grain "disarticulation zone," converting the brittle floral axis (the rachis) of the wild-type into a tough, non-brittle form that promoted grain retention. By tracing the evolutionary history of allelic variation in both genes, we conclude that spatially and temporally independent selections of germplasm with a non-brittle rachis were made during the domestication of barley by farmers in the southern and northern regions of the Levant, actions that made a major contribution to the emergence of early agrarian societies.

Deciphering and prediction of transcriptome dynamics under fluctuating field conditions.

Determining the drivers of gene expression patterns is more straightforward in laboratory conditions than in the complex fluctuating environments where organisms typically live. We gathered transcriptome data from the leaves of rice plants in a paddy field along with the corresponding meteorological data and used them to develop statistical models for the endogenous and external influences on gene expression. Our results indicate that the transcriptome dynamics are predominantly governed by endogenous diurnal rhythms, ambient temperature, plant age, and solar radiation. The data revealed diurnal gates for environmental stimuli to influence transcription and pointed to relative influences exerted by circadian and environmental factors on different metabolic genes. The model also generated predictions for the influence of changing temperatures on transcriptome dynamics. We anticipate that our models will help translate the knowledge amassed in laboratories to problems in agriculture and that our approach to deciphering the transcriptome fluctuations in complex environments will be applicable to other organisms.

Genome-wide analysis of spatiotemporal gene expression patterns during early embryogenesis in rice.

Embryogenesis in rice is different from that of most dicotolydonous plants in that it shows a non-stereotypic cell division pattern, formation of dorsal-ventral polarity, and endogenous initiation of the radicle. To reveal the transcriptional features associated with developmental events during rice early embryogenesis, we used microarray analysis coupled with laser microdissection to obtain both spatial and temporal transcription profiles. Our results allowed us to determine spatial expression foci for each expressed gene in the globular embryo, which revealed the importance of phytohormone-related genes and a suite of transcription factors to early embryogenesis. Our analysis showed the polarized expression of a small number of genes along the apical-basal and dorsal-ventral axes in the globular embryo, which tended to fluctuate in later developmental stages. We also analyzed gene expression patterns in the early globular embryo and how this relates to expression in embryonic organs at later stages. We confirmed the accuracy of the expression patterns found by microarray analysis of embryo subdomains using in situ hybridization. Our study identified homologous genes from Arabidopsis thaliana with known functions in embryogenesis in addition to unique and uncharacterized genes that show polarized expression patterns during embryogenesis. The results of this study are presented in a database to provide a framework for spatiotemporal gene expression during rice embryogenesis, to serve as a resource for future functional analysis of genes, and as a basis for comparative studies of plant embryogenesis.

Serotonin attenuates biotic stress and leads to lesion browning caused by a hypersensitive response to Magnaporthe oryzae penetration in rice.

The hypersensitive response (HR) of plants is one of the earliest responses to prevent pathogen invasion. A brown dot lesion on a leaf is visual evidence of the HR against the blast fungus Magnaporthe oryzae in rice, but tracking the browning process has been difficult. In this study, we induced the HR in rice cultivars harboring the blast resistance gene Pit by inoculation of an incompatible M. oryzae strain, which generated a unique resistance lesion with a brown ring (halo) around the brown fungal penetration site. Inoculation analysis using a plant harboring Pit but lacking an enzyme that catalyzes tryptamine to serotonin showed that high accumulation of the oxidized form of serotonin was the cause of the browning at the halo and penetration site. Our analysis of the halo browning process in the rice leaf revealed that abscisic acid enhanced biosynthesis of serotonin under light conditions, and serotonin changed to the oxidized form via hydrogen peroxide produced by light. The dramatic increase in serotonin, which has a high antioxidant activity, suppressed leaf damage outside the halo, blocked expansion of the browning area and attenuated inhibition of plant growth. These results suggest that serotonin helps to reduce biotic stress in the plant by acting as a scavenger of oxygen radicals to protect uninfected tissues from oxidative damage caused by the HR. The deposition of its oxide at the HR lesion is observed as lesion browning.

Overexpression of TIFY genes promotes plant growth in rice through jasmonate signaling.

Because environmental stress can reduce crop growth and yield, the identification of genes that enhance agronomic traits is increasingly important. Previous screening of full-length cDNA overexpressing (FOX) rice lines revealed that OsTIFY11b, one of 20 TIFY proteins in rice, affects plant size, grain weight, and grain size. Therefore, we analyzed the effect of OsTIFY11b and nine other TIFY genes on the growth and yield of corresponding TIFY-FOX lines. Regardless of temperature, grain weight and culm length were enhanced in lines overexpressing TIFY11 subfamily genes, except OsTIFY11e. The TIFY-FOX plants exhibited increased floret number and reduced days to flowering, as well as reduced spikelet fertility, and OsTIFY10b, in particular, enhanced grain yield by minimizing decreases in fertility. We suggest that the enhanced growth of TIFY-transgenic rice is related to regulation of the jasmonate signaling pathway, as in Arabidopsis. Moreover, we discuss the potential application of TIFY overexpression for improving crop yield.

TAWAWA1, a regulator of rice inflorescence architecture, functions through the suppression of meristem phase transition.

Inflorescence structures result from the activities of meristems, which coordinate both the renewal of stem cells in the center and organ formation at the periphery. The fate of a meristem is specified at its initiation and changes as the plant develops. During rice inflorescence development, newly formed meristems acquire a branch meristem (BM) identity, and can generate further meristems or terminate as spikelets. Thus, the form of rice inflorescence is determined by a reiterative pattern of decisions made at the meristems. In the dominant gain-of-function mutant tawawa1-D, the activity of the inflorescence meristem (IM) is extended and spikelet specification is delayed, resulting in prolonged branch formation and increased numbers of spikelets. In contrast, reductions in TAWAWA1 (TAW1) activity cause precocious IM abortion and spikelet formation, resulting in the generation of small inflorescences. TAW1 encodes a nuclear protein of unknown function and shows high levels of expression in the shoot apical meristem, the IM, and the BMs. TAW1 expression disappears from incipient spikelet meristems (SMs). We also demonstrate that members of the SHORT VEGETATIVE PHASE subfamily of MADS-box genes function downstream of TAW1. We thus propose that TAW1 is a unique regulator of meristem activity in rice and regulates inflorescence development through the promotion of IM activity and suppression of the phase change to SM identity.

Rice MEL2, the RNA recognition motif (RRM) protein, binds in vitro to meiosis-expressed genes containing U-rich RNA consensus sequences in the 3'-UTR.

Post-transcriptional gene regulation by RNA recognition motif (RRM) proteins through binding to cis-elements in the 3'-untranslated region (3'-UTR) is widely used in eukaryotes to complete various biological processes. Rice MEIOSIS ARRESTED AT LEPTOTENE2 (MEL2) is the RRM protein that functions in the transition to meiosis in proper timing. The MEL2 RRM preferentially associated with the U-rich RNA consensus, UUAGUU[U/A][U/G][A/U/G]U, dependently on sequences and proportionally to MEL2 protein amounts in vitro. The consensus sequences were located in the putative looped structures of the RNA ligand. A genome-wide survey revealed a tendency of MEL2-binding consensus appearing in 3'-UTR of rice genes. Of 249 genes that conserved the consensus in their 3'-UTR, 13 genes spatiotemporally co-expressed with MEL2 in meiotic flowers, and included several genes whose function was supposed in meiosis; such as Replication protein A and OsMADS3. The proteome analysis revealed that the amounts of small ubiquitin-related modifier-like protein and eukaryotic translation initiation factor3-like protein were dramatically altered in mel2 mutant anthers. Taken together with transcriptome and gene ontology results, we propose that the rice MEL2 is involved in the translational regulation of key meiotic genes on 3'-UTRs to achieve the faithful transition of germ cells to meiosis.

Whole-genome analysis of herbicide-tolerant mutant rice generated by Agrobacterium-mediated gene targeting.

Gene targeting (GT) is a technique used to modify endogenous genes in target genomes precisely via homologous recombination (HR). Although GT plants are produced using genetic transformation techniques, if the difference between the endogenous and the modified gene is limited to point mutations, GT crops can be considered equivalent to non-genetically modified mutant crops generated by conventional mutagenesis techniques. However, it is difficult to guarantee the non-incorporation of DNA fragments from Agrobacterium in GT plants created by Agrobacterium-mediated GT despite screening with conventional Southern blot and/or PCR techniques. Here, we report a comprehensive analysis of herbicide-tolerant rice plants generated by inducing point mutations in the rice ALS gene via Agrobacterium-mediated GT. We performed genome comparative genomic hybridization (CGH) array analysis and whole-genome sequencing to evaluate the molecular composition of GT rice plants. Thus far, no integration of Agrobacterium-derived DNA fragments has been detected in GT rice plants. However, >1,000 single nucleotide polymorphisms (SNPs) and insertion/deletion (InDels) were found in GT plants. Among these mutations, 20-100 variants might have some effect on expression levels and/or protein function. Information about additive mutations should be useful in clearing out unwanted mutations by backcrossing.

Inflorescence meristem identity in rice is specified by overlapping functions of three AP1/FUL-like MADS box genes and PAP2, a SEPALLATA MADS box gene.

In plants, the transition to reproductive growth is of particular importance for successful seed production. Transformation of the shoot apical meristem (SAM) to the inflorescence meristem (IM) is the crucial first step in this transition. Using laser microdissection and microarrays, we found that expression of PANICLE PHYTOMER2 (PAP2) and three APETALA1 (AP1)/FRUITFULL (FUL)-like genes (MADS14, MADS15, and MADS18) is induced in the SAM during meristem phase transition in rice (Oryza sativa). PAP2 is a MADS box gene belonging to a grass-specific subclade of the SEPALLATA subfamily. Suppression of these three AP1/FUL-like genes by RNA interference caused a slight delay in reproductive transition. Further depletion of PAP2 function from these triple knockdown plants inhibited the transition of the meristem to the IM. In the quadruple knockdown lines, the meristem continued to generate leaves, rather than becoming an IM. Consequently, multiple shoots were formed instead of an inflorescence. PAP2 physically interacts with MAD14 and MADS15 in vivo. Furthermore, the precocious flowering phenotype caused by the overexpression of Hd3a, a rice florigen gene, was weakened in pap2-1 mutants. Based on these results, we propose that PAP2 and the three AP1/FUL-like genes coordinately act in the meristem to specify the identity of the IM downstream of the florigen signal.

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.

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.

NYC4, the rice ortholog of Arabidopsis THF1, is involved in the degradation of chlorophyll - protein complexes during leaf senescence.

Yellowing/chlorophyll breakdown is a prominent phenomenon in leaf senescence, and is associated with the degradation of chlorophyll - protein complexes. From a rice mutant population generated by ionizing radiation, we isolated nyc4-1, a stay-green mutant with a defect in chlorophyll breakdown during leaf senescence. Using gene mapping, nyc4-1 was found to be linked to two chromosomal regions. We extracted Os07g0558500 as a candidate for NYC4 via gene expression microarray analysis, and concluded from further evidence that disruption of the gene by a translocation-related event causes the nyc4 phenotype. Os07g0558500 is thought to be the ortholog of THF1 in Arabidopsis thaliana. The thf1 mutant leaves show variegation in a light intensity-dependent manner. Surprisingly, the Fv /Fm value remained high in nyc4-1 during the dark incubation, suggesting that photosystem II retained its function. Western blot analysis revealed that, in nyc4-1, the PSII core subunits D1 and D2 were significantly retained during leaf senescence in comparison with wild-type and other non-functional stay-green mutants, including sgr-2, a mutant of the key regulator of chlorophyll degradation SGR. The role of NYC4 in degradation of chlorophyll and chlorophyll - protein complexes during leaf senescence is discussed.

Strigolactone and cytokinin act antagonistically in regulating rice mesocotyl elongation in darkness.

Strigolactones (SLs) are a group of phytohormones that control plant growth and development including shoot branching. Previous studies of the phenotypes of SL-related rice (Oryza sativa) dwarf (d) mutants demonstrated that SLs inhibit mesocotyl elongation by controlling cell division. Here, we found that the expression of cytokinin (CK)-responsive type-A RESPONSE REGULATOR (RR) genes was higher in d10-1 and d14-1 mutants than in the wild type. However, CK levels in mesocotyls of the d mutants were not very different from those in the wild type. On the other hand, application of a synthetic CK (kinetin) enhanced mesocotyl elongation in the d mutants and the wild type. d10-1 and d14-1 mesocotyls were more sensitive to CK than wild-type mesocotyls, suggesting that the up-regulation of the CK-responsive type-A RR genes and the higher elongation of mesocotyls in the d mutants are mainly due to the increased sensitivity of the d mutants to CK. Co-treatment with kinetin and a synthetic SL (GR24) confirmed the antagonistic functions of SL and CK on mesocotyl elongation. OsTCP5, which encodes a transcription factor belonging to the cell division-regulating TCP family, was also regulated by SL and CK and its expression was negatively correlated with mesocotyl length. These findings suggest that OsTCP5 contributes to the SL- and CK-controlled mesocotyl elongation in darkness.

A guanine nucleotide exchange factor for Rab5 proteins is essential for intracellular transport of the proglutelin from the Golgi apparatus to the protein storage vacuole in rice endosperm.

Rice (Oryza sativa) glutelins are synthesized on the endoplasmic reticulum as a precursor, which are then transported via the Golgi to protein storage vacuoles (PSVs), where they are proteolytically processed into acidic and basic subunits. The glutelin precursor mutant6 (glup6) accumulates abnormally large amounts of proglutelin. Map-base cloning studies showed that glup6 was a loss-of-function mutant of guanine nucleotide exchange factor (GEF), which activates Rab GTPase, a key regulator of membrane trafficking. Immunofluorescence studies showed that the transport of proglutelins and α-globulins to PSV was disrupted in glup6 endosperm. Secreted granules of glutelin and α-globulin were readily observed in young glup6 endosperm, followed by the formation of large dilated paramural bodies (PMBs) containing both proteins as the endosperm matures. The PMBs also contained membrane biomarkers for the Golgi and prevacuolar compartment as well as the cell wall component, ß-glucan. Direct evidence was gathered showing that GLUP6/GEF activated in vitro GLUP4/Rab5 as well as several Arabidopsis (Arabidopsis thaliana) Rab5 isoforms to the GTP-bound form. Therefore, loss-of-function mutations in GEF or Rab5 disrupt the normal transport of proglutelin from the Golgi to PSVs, resulting in the initial extracellular secretion of these proteins followed, in turn, by the formation of PMBs. Overall, our results indicate that GLUP6/GEF is the activator of Rab5 GTPase and that the cycling of GTP- and GDP-bound forms of this regulatory protein is essential for the intracellular transport of proglutelin and α-globulin from the Golgi to PSVs and in the maintenance of the general structural organization of the endomembrane system in rice seeds.

Microarray analysis of laser-microdissected tissues indicates the biosynthesis of suberin in the outer part of roots during formation of a barrier to radial oxygen loss in rice (Oryza sativa).

Internal aeration is crucial for root growth in waterlogged soil. A barrier to radial oxygen loss (ROL) can enhance long-distance oxygen transport via the aerenchyma to the root tip; a higher oxygen concentration at the apex enables root growth into anoxic soil. The ROL barrier is formed within the outer part of roots (OPR). Suberin and/or lignin deposited in cell walls are thought to contribute to the barrier, but it is unclear which compound is the main constituent. This study describes gene expression profiles during ROL barrier formation in rice roots to determine the relative responses of suberin and/or lignin biosyntheses for the barrier. OPR tissues were isolated by laser microdissection and their transcripts were analysed by microarray. A total of 128 genes were significantly up- or downregulated in the OPR during the barrier formation. Genes associated with suberin biosynthesis were strongly upregulated, whereas genes associated with lignin biosynthesis were not. By an ab initio analysis of the promoters of the upregulated genes, the putative cis-elements that could be associated with transcription factors, WRKY, AP2/ERF, NAC, bZIP, MYB, CBT/DREB, and MADS, were elucidated. They were particularly associated with the expression of transcription factor genes containing WRKY, AP2, and MYB domains. A semiquantitative reverse-transcription PCR analysis of genes associated with suberin biosynthesis (WRKY, CYP, and GPAT) confirmed that they were highly expressed during ROL barrier formation. Overall, these results suggest that suberin is a major constituent of the ROL barrier in roots of rice.

Os-GIGANTEA confers robust diurnal rhythms on the global transcriptome of rice in the field.

The circadian clock controls physiological traits such as flowering time, photosynthesis, and growth in plants under laboratory conditions. Under natural field conditions, however, little is known about the significance of the circadian clock in plants. By time-course transcriptome analyses of rice (Oryza sativa) leaves, using a newly isolated rice circadian clock-related mutant carrying a null mutation in Os-GIGANTEA (Os-GI), we show here that Os-GI controlled 75% (false discovery rate = 0.05) of genes among 27,201 genes tested and was required for strong amplitudes and fine-tuning of the diurnal rhythm phases of global gene expression in the field. However, transcripts involved in primary metabolism were not greatly affected by osgi. Time-course metabolome analyses of leaves revealed no trends of change in primary metabolites in osgi plants, and net photosynthetic rates and grain yields were not affected. By contrast, some transcripts and metabolites in the phenylpropanoid metabolite pathway were consistently affected. Thus, net primary assimilation of rice was still robust in the face of such osgi mutation-related circadian clock defects in the field, unlike the case with defects caused by Arabidopsis thaliana toc1 and ztl mutations in the laboratory.

Low-affinity cation transporter (OsLCT1) regulates cadmium transport into rice grains.

Accumulation of cadmium (Cd) in rice (Oryza sativa L.) grains poses a potential health problem, especially in Asia. Most Cd in rice grains accumulates through phloem transport, but the molecular mechanism of this transport has not been revealed. In this study, we identified a rice Cd transporter, OsLCT1, involved in Cd transport to the grains. OsLCT1-GFP was localized at the plasma membrane in plant cells, and OsLCT1 showed Cd efflux activity in yeast. In rice plants, strong OsLCT1 expression was observed in leaf blades and nodes during the reproductive stage. In the uppermost node, OsLCT1 transcripts were detected around large vascular bundles and in diffuse vascular bundles. RNAi-mediated knockdown of OsLCT1 did not affect xylem-mediated Cd transport but reduced phloem-mediated Cd transport. The knockdown plants of OsLCT1 accumulated approximately half as much Cd in the grains as did the control plants. The content of other metals in rice grains and plant growth were not negatively affected by OsLCT1 suppression. These results suggest that OsLCT1 functions at the nodes in Cd transport into grains and that in a standard japonica cultivar, the regulation of OsLCT1 enables the generation of "low-Cd rice" without negative effects on agronomical traits. These findings identify a transporter gene for phloem Cd transport in plants.

Genome-wide transcriptome dissection of the rice root system: implications for developmental and physiological functions.

The root system is a crucial determinant of plant growth potential because of its important functions, e.g. uptake of water and nutrients, structural support and interaction with symbiotic organisms. Elucidating the molecular mechanism of root development and functions is therefore necessary for improving plant productivity, particularly for crop plants, including rice (Oryza sativa). As an initial step towards developing a comprehensive understanding of the root system, we performed a large-scale transcriptome analysis of the rice root via a combined laser microdissection and microarray approach. The crown root was divided into eight developmental stages along the longitudinal axis and three radial tissue types at two different developmental stages, namely: epidermis, exodermis and sclerenchyma; cortex; and endodermis, pericycle and stele. We analyzed a total of 38 microarray data and identified 22,297 genes corresponding to 17,010 loci that showed sufficient signal intensity as well as developmental- and tissue type-specific transcriptome signatures. Moreover, we clarified gene networks associated with root cap function and lateral root formation, and further revealed antagonistic and synergistic interactions of phytohormones such as auxin, cytokinin, brassinosteroids and ethylene, based on the expression pattern of genes related to phytohormone biosynthesis and signaling. Expression profiling of transporter genes defined not only major sites for uptake and transport of water and nutrients, but also distinct signatures of the radial transport system from the rhizosphere to the xylem vessel for each nutrient. All data can be accessed from our gene expression profile database, RiceXPro (http://ricexpro.dna.affrc.go.jp), thereby providing useful information for understanding the molecular mechanisms involved in root system development of crop plants.

Roles of pollen-specific boron efflux transporter, OsBOR4, in the rice fertilization process.

Arabidopsis thaliana BOR1 was the first boron (B) transporter identified in living systems. There are four AtBOR1-like genes, OsBOR1, 2, 3 and 4, present in the rice genome. We characterized the activity, expression and physiological function of OsBOR4. OsBOR4 is an active efflux transporter of B. Quantitative PCR analysis and OsBOR4 promoter-green fluorescent protein (GFP) fusion revealed that OsBOR4 was both highly and specifically expressed in pollen. We obtained five Tos17 insertion mutants of osbor4. The pollen grains were viable and development of floral organs was normal in the homozygous osbor4 mutants. We observed that in all Tos17 insertion lines tested, the frequency of osbor4 homozygous plants was lower than expected in the progeny of self-fertilized heterozygous plants. These results establish that OsBOR4 is essential for normal reproductive processes. Pollen from osbor4 homozygous plants elongated fewer tubes on wild-type stigmas, and tube elongation of mutant pollen was less efficient compared with the wild-type pollen, suggesting reduced competence of osbor4 mutant pollen. The reduced competence of mutant pollen was further supported by the crosses of independent Tos17-inserted alleles of OsBOR4. Our results suggest that OsBOR4, a boron efflux transporter, is required for normal pollen germination and/or tube elongation.

Identification of transcription factors involved in rice secondary cell wall formation.

Using co-expression network analysis, we identified 123 transcription factors (TFs) as candidate secondary cell wall regulators in rice. To validate whether these TFs are associated with secondary cell wall formation, six TF genes belonging to the MYB, NAC or homeodomain-containing TF families were overexpressed or downregulated in rice. With the exception of OsMYB58/63-RNAi plants, all transgenic plants showed phenotypes possibly related to secondary cell wall alteration, such as dwarfism, narrow and dark green leaves, and also altered rice cinnamyl alcohol dehydrogenase 2 (OsCAD2) gene expression and lignin content. These results suggest that many of the 123 candidate secondary cell wall-regulating TFs are likely to function in secondary cell wall formation in rice. Further analyses were performed for the OsMYB55/61 and OsBLH6 TFs, the former being a TF in which the Arabidopsis ortholog is known to participate in lignin biosynthesis (AtMYB61) and the latter being one for which no previous involvement in cell wall formation has been reported even in Arabidopsis (BLH6). OsMYB55/61 and OsBLH6-GFP fusion proteins localized to the nucleus of onion epidermal cells. Moreover, expression of a reporter gene driven by the OsCAD2 promoter was enhanced in rice calli when OsMYB55/61 or OsBLH6 was transiently expressed, demonstrating that they function in secondary cell wall formation. These results show the validity of identifying potential secondary cell wall TFs in rice by the use of rice co-expression network analysis.