Skip Navigation


Plant and Cell Physiology Advance Access originally published online on September 6, 2008
Plant and Cell Physiology 2008 49(10):1417-1428; doi:10.1093/pcp/pcn128
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
49/10/1417    most recent
pcn128v2
pcn128v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Related articles in PCP
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (4)
Google Scholar
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Agricola
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved.
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and the Japanese Society of Plant Physiologists are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact journals.permissions@oxfordjournals.org


Special Issue - Regular Paper

Various Spatiotemporal Expression Profiles of Anther-Expressed Genes in Rice

Tokunori Hobo1,11, Keita Suwabe2,11, Koichiro Aya1, Go Suzuki3, Kentaro Yano4, Takeshi Ishimizu5, Masahiro Fujita6, Shunsuke Kikuchi4, Kazuki Hamada4, Masumi Miyano2,7, Tomoaki Fujioka2, Fumi Kaneko2, Tomohiko Kazama2,7, Yoko Mizuta6, Hirokazu Takahashi8, Katsuhiro Shiono8, Mikio Nakazono8, Nobuhiro Tsutsumi8, Yoshiaki Nagamura9, Nori Kurata6, Masao Watanabe2,7,* and Makoto Matsuoka1,*

1Bioscience and Biotechnology Center, Nagoya University, Furocho, Chikusa, Nagoya, 464-8601 Japan
2Graduate School of Life Sciences, Tohoku University, Sendai, 980-8577 Japan
3Division of Natural Science, Osaka Kyoiku University, Kashiwara, 582-8582 Japan
4Faculty of Agriculture, Meiji University, Kawasaki, 214-8571 Japan
5Department of Chemistry, Graduate School of Science, Osaka University, Osaka, 560-0043 Japan
6Plant Genetics Laboratory, National Institute of Genetics, Mishima, 411-8540 Japan
7The 21st Century Center of Excellence Program, Iwate University, Morioka, 020-8550 Japan
8Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, 113-8657 Japan
9Division of Genome and Biodiversity Research, National Institute of Agrobiological Sciences, Tsukuba, 305-8602 Japan
10Faculty of Science, Tohoku University, Sendai, 980-8587 Japan

*Corresponding authors: Masao Watanabe, E-mail, nabe{at}ige.tohoku.ac.jp; Fax, +81-22-217-5683; Makoto Matsuoka, E-mail, makoto{at}nuagr1.agr.nagoya-u.ac.jp; Fax, +81-52-789-5226


    Abstract
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
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.

Keywords: Anther - Gene ontology - Laser microdissection - Microarray - Oryza sativa L - Synchronous gene expression

Abbreviations: BC, bicellular; GO, gene ontology; LM, laser microdissection; MEI, meiosis; miRNA, microRNA; TC, tri­cellular; TET, tetrad; UN, uninuclear; VSN, variance stabiliza­tion normalization


    Introduction
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
Anther tissues in flowering plants are important for microsporogenesis and dispersing the male gamete, pollen. Anther tissues consist of three layers, L1, L2 and L3. The epidermis and stomium are differentiated at the final devel­opmental stage in the L1 layer. The microspore/pollen, endothecium and outer tapetum are formed in the L2 layer, while the L3 layer differentiates into vascular bundle and inner tapetum (Goldberg et al. 1993Go). Although the sporo­genous cells and parietal cells originate from the same L2 layer, their final fates are different. In the case of sporogenous cells, following meiosis, microspores are released into the anther locule, and each microspore divides into vegetative and generative cells. In rice, the generative cell further divides into two sperm cells. Desiccated mature pollen grains are dispersed from anthers onto the stigma surface ready for fertilization. In contrast, parietal cells differentiate into endothecium and outer tapetum. Tapetum cells function as nurse cells for pollen development; the tapetum supports pollen development by secreting nutrients and other secondary metabolites. The tapetum constitutes the pollen outer surface, but degrades before anther dehiscence (Goldberg et al. 1993Go). Therefore, while microspore/pollen and outer tapetum originate from the same cell lineage, the L2 layer, their fetal functions are completely different.

In order to understand the molecular mechanisms of anther development, particularly microspore/pollen differentiation and tapetum development, transcriptome analyses have been conducted in anthers of rice and Lotus japonicus (Endo et al. 2002Go, Endo et al. 2004Go). Some anther-specific genes were further characterized by RNA in situ hybridization analysis, which is extremely laborious and not comprehensive. These studies demonstrated the existence of microspore/pollen-specific and tapetum-specific genes (Endo et al. 2004Go, Masuko et al. 2006Go). Briefly, in the mature pollen grain, coordinated expression of pollen-specific transcripts was detected. These transcripts are related to the cell cytoskeleton, cell wall re-organization, methionine metabolism, proton pumps and sugar transporters, and are presumably important for pollen germination and pollen tube elongation (Endo et al. 2002Go, Endo et al. 2004Go, Masuko et al. 2006Go). In contrast, tapetum-specific transcripts were detected relating to secondary metabolism, fatty acid biosynthesis, protein secretion and the gibberellin signaling cascade (Endo et al. 2004Go, Masuko et al. 2006Go). The different characteristics of the two transcriptomes suggested that gene functions in pollen and tapetum are distinct. However, further comprehensive analysis of microspore/pollen-specific and tapetum-specific genes is necessary for a complete functional understanding of these cell types.

Because laser microdissection (LM) is a useful technique to separate cell types within complex plant tissues (Asano et al. 2002Go, Kerk et al. 2003Go, Nakazono et al. 2003Go, Day et al. 2005Go, Nelson et al. 2006Go, Ohtsu et al. 2007Go), we performed a 44K microarray with RNAs derived from LM-separated microspore/pollen and tapetum cells in rice. The reliability of the data was confirmed by comparison with the published gene expression profile of 156 previously identified anther-specific genes in rice (refer to Suwabe et al. 2008Go in this special issue). Having validated the technique, the genome-wide expression profiles of microspore/pollen and tapetum in rice were determined in this study according to the developmental stages. From the cluster analysis, 28,141 genes were categorized into 20 clusters according to their characteristic expression patterns, and gene ontology (GO) was applied to each cluster. We discuss noteworthy characteristics of the microspore/pollen and tapetum transcriptomes in rice.


    Results and Discussion
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
Various spatiotemporal expression patterns of genes in rice anthers
We conducted a 44K LM-microarray analysis of microspore/pollen and tapetum cells of japonica rice, Oryza sativa cv. Nipponbare, at each of the five stages of anther development: meiosis (MEI), tetrad (TET), uninuclear (UN) microspore, bicellular (BC) pollen and tricellular (TC) pollen (Suwabe et al. 2008Go). Because tapetum degradation occurs at the BC stage, the tapetum transcriptome was only characterized at the MEI, TET and UN stages. After normalization of Z-scored 44K array data, a cluster analysis was performed with TIGR MeV software (Fig. 1, Supplementary Figs. S1–S20, Tables S1–S20). It should be noted that the microspore/pollen and tapetum transcriptomes were successfully separated without contamination in the present LM-microarray, as described by Suwabe et al. (2008Go) in this issue. This conclusion was further supported by the existence of the anther-specific and microspore-specific clades in the heat map (Fig. 2). In total, 28,141 anther-expressed genes were categorized into 20 clusters, and the average number of genes in each cluster was 1,407. Expression levels of genes in six clusters (Nos. 4, 5, 6, 8, 10 and 19) were higher in the microspore/pollen than in the tapetum, indicating that these clusters were microspore/pollen dominant. On the other hand, no obvious tapetum-dominant profiles were observed in the 20 clusters, mainly because only a small proportion of genes showed the strictly tapetum-dominant expression pattern, and these were distributed among several clusters in this analysis.


Figure 1
View larger version (82K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1 Gene expression patterns of 20 clusters in microspore/pollen and tapetum. In each cluster, Z-scored data with a standard error bar were plotted according to developmental stages. The number of genes in the cluster is indicated at the upper left side. The calculated percentage of anther-enriched genes in each cluster is shown at the lower right side. The high synchronous gene expression pattern between microspore/pollen and tapetum is shown by red circles in clusters 2 and 20. The low synchronous pattern is shown by blue circles in clusters 3, 7, 12, 13, 16 and 17.

 

Figure 2
View larger version (70K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2 Heat map view of the microspore/pollen-specific clade (A) and tapetum-specific clade (B). Microspore/pollen- and tapetum-specific clades were extracted from clusters 3 and 2, respectively. Red indicates higher, while green represents lower expression.

 
To discriminate genes preferentially expressed in anthers, we compared the LM-microarray data of microspore/pollen and tapetum with microarray data sets obtained for rice roots and leaves (Y. Nagamura, unpublished results). In the present study, anther-expressed genes with a vegetative expression signal score <50 were defined as ‘anther-specific’ genes. According to this definition, 3,468 anther-specific genes were identified, accounting for 12.3% of the total 28,141 anther-expressed genes. The percentage of anther-specific genes was then calculated in each cluster (Fig. 1), and relatively lower ratios of anther-specific genes (~2–4%) were observed for the clusters of constitutively tapetum-expressed genes (Nos. 13, 14, 15 and 18). In contrast, in the clusters of genes expressed in MEI microspore (clusters 12 and 17) and in late developmental stages of pollen (clusters 5 and 10), >20% of genes were anther specific. The increased proportion of anther-specific transcripts in clusters 5 and 10 is consistent with analysis of the Arabidopsis male gametophyte-specific transcriptome (Honys and Twell 2004Go), which showed that the transition from BC to TC pollen was accompanied by an increase in the proportion of male gamete-specific transcripts.

When we focused on the tapetum transcriptome, genes in clusters 12, 20 and 2 were preferentially expressed only at one stage (MEI, TET and UN stages, respectively), suggesting that the biological functions of tapetum genes differed in the three stages. In other clusters, genes in cluster 7 were predominantly expressed in both MEI and TET stages, while those in clusters 13, 14, 15 and 18 were expressed continuously in all the three stages.

In contrast, analysis of the microspore/pollen transcriptome revealed complex profiles. Expression patterns in the BC pollen were variable in some clusters, e.g. in cluster 1, expression was high in the UN microspore, gradually down-regulated in the late BC pollen, and low in the TC pollen. Such gradual down-regulation within the BC stage was also observed in clusters 15 and 16, whereas gradual up-regulation was seen in clusters 5, 6, 10, 11 and 19. In view of these patterns, it would be possible to divide the BC stage in the pollen transcriptome into early BC and late BC stages.

Throughout the developmental stages from MEI to TC, there are many expression patterns among the microspore/pollen genes. For example, gradual up-regulation of genes from MEI to TC was observed in clusters 5 and 6, while gradual down-regulation occurred in cluster 16. In the case of clusters 9 and 19, genes were gradually down-regulated from MEI to UN, but then gradually up-regulated again from BC to TC pollen. Another pattern was shown in cluster 3, with constant expression in the microspore/pollen at the UN to TC pollen stages. The variety of expression profiles in the microspore/pollen indicated that biological events in the male gametophyte are developmentally regulated by many anther-expressed genes with various expression patterns.

Interestingly, synchronous gene expression in both micro­spore and tapetum was observed in eight clusters: Nos. 2 and 20 with a high degree of synchrony (indicated by the red color in Fig. 1), and Nos. 3, 7, 12, 13, 16 and 17 with low synchronization (indicated by the blue color). These synchronous profiles cannot be attributed to contamination artifacts during the LM procedure, because microspore-specific and tapetum-specific clades could be identified (Fig. 2, Suwabe et al. 2008Go), and the expression patterns of microspore and tapetum were fundamentally independent in the other 12 clusters. The number of genes classified into the eight synchronous clusters was 10,810 (38.4%). Synchronous expression in the microspore/pollen and tapetum is one of the major features of the expression profiles in anthers, and argues against the conventional idea that tapetum cells function solely as nurse cells for feeding microspores in anther locules, with genes in different categories being expressed in the microspore and tapetum. The finding that synchronous transcription of some groups of genes occurs in both the micro­spore and tapetum suggests that these cells do share many developmental pathways as described in Scott et al. (2004Go), and is also consistent with the fact that both the microspore and tapetum are derived from the same cell lineage, L2.

Characterization of specific gene clusters based on GO
Genes with functional annotation were classified in each cluster based on GO by using Gene Ontology Slim terms (Fig. 3, http://www.tigr.org/tdb/e2k1/osa1/batch_downloads.html). We focused on five categories with likely functional significance in anther development: transcription, translation, lipid metabolism, secondary metabolism and the cell cycle (Fig. 4).


Figure 3
View larger version (89K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3 Gene classification based on gene ontology (GO) in 20 clusters. The genes in the 20 clusters were re-classified according to GO. For five GO terms (transcription, translation, lipid metabolic process, secondary metabolic process and the cell cycle), significant terms are highlighted on the circled graph. For the lipid metabolism process category, three additional clusters (Nos. 13, 14 and 15) are indicated by asterisks.

 

Figure 4
View larger version (35K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4 Representative results of the GO classification in each cluster. In each category, clusters with a high percentage coverage, shown in Fig. 3, are highlighted. For the lipid metabolism process category, we selected three additional clusters (Nos. 13, 14 and 15). They are indicated by asterisks in the three cells.

 
Genes related to transcription occupied 1.8–7.2% of genes in each cluster, and were frequently observed in clusters 7, 13, 14 and 20 (6.0–7.2% in each cluster, Fig. 4). These clusters were all characterized as tapetum expression, indicating that transcriptional regulation might be significantly active in the tapetum. In particular, clusters 13 and 14, showing similar expression profiles, included 312 genes (28.0%) out of the 1,114 transcription-related genes in all 20 clusters. Clusters 13 and 14 contained various transcription factor genes, with annotations including bHLH, MYB, GRAS, MADS, HD-Zip and Zn-finger-type transcription factors.


View this table:
[in this window]
[in a new window]

 
Table 1 Representative common motifs in the 140 male gamete- and tapetum-expressed genes

 
Translation-related genes were frequently observed in clusters 1 and 4 (6.3 and 7.3% in each cluster, respectively); the two clusters contained 149 genes (34.6%) out of the 431 translation-related genes in all the clusters. Common characteristics of the expression profiles in clusters 1 and 4 were higher gene expression in early BC pollen. The profiles in the clusters with a high percentage of transcription-related genes (Nos. 7, 13, 14 and 20) and translation-related genes (Nos. 1 and 4) were considerably different; in particular, the expression profiles of clusters 4 and 14 were almost opposite (Fig. 1). Thus, translation in the BC pollen might be achieved using transcripts produced in earlier stages of pollen (TET and UN microspores).

Lipid metabolism is important in pollen development, and failure of lipid synthesis frequently causes male sterility (Mariani and Wolters-Arts 2000Go, Ariizumi et al. 2003Go). In our GO analysis, lipid metabolism genes often appeared in cluster 2, whose profiles were characterized by high expression in the UN microspores and tapetum. This indicates that regulation of lipid metabolism-related genes might function at the UN stage. Several fatty acid synthesis genes, such as genes for β-ketoacyl-CoA synthase and β-keto acyl reductase, were also observed in clusters 13, 14 and 15 (described further in the following section), although the percentage coverage of lipid metabolism genes in these clusters was not so high.

Genes related to secondary metabolism were the most frequently observed in clusters 2 and 20 (3.9 and 3.2% in each cluster), including chalcone synthase, chalcone reductase and phytoene synthase. These two clusters were categorized as showing synchronous expression in the TET or UN microspore and tapetum.

A relatively small portion of cell cycle-related genes (0.0–1.4%) was categorized in the GO analysis. However, most of these (19.8%; 17 genes out of all the 86 cell cycle genes) were included in cluster 17. This might be due to the existence of genes related to meiosis, because notable gene expression was observed at the MEI stage in cluster 17. Genes annotated as cyclin, cyclin–F-box and rad51/recA were identified in cluster 17.

Functional relationships between other characteristic genes in each cluster and biological events in the microspore/pollen and tapetum
During the maturation of anther tissues, many characteristic biological events can be related to the microspore/pollen and tapetum transcriptomes. For example, important biological events for the male gamete are pollen maturation, pollen germination and pollen tube elongation. To date, many pollen-specific genes, related to carbohydrate metabolism (glycosyltransferase, glycoside hydrolase, sugar transporter, etc.), cytoskeleton re-organization (actin), Ca+ signaling (calmodulin) and tip growth (small GTPase), have been identified and characterized in terms of these biological events (Franklin-Tong 1999Go, Endo et al. 2000Go, Endo et al. 2002Go, Golovkin and Reddy 2003Go, Honys and Twell 2003Go, Endo et al.2004Go, Jiang et al. 2005Go, Huang and Yang 2006Go, Malhó et al. 2006Go, Shimamura et al. 2007Go, Cheung and Wu 2008Go). These genes, whose transcripts are present in mature pollen grains and are translated during pollen germination (Hoekstra and Bruinsma 1979Go, Mascarenhas et al. 1984Go, Schrauwen et al. 1990Go), were confirmed in our study to be predominantly expressed in the late BC and TC pollen (clusters 5, 6, 10, 11 and 19).

During pollen tube elongation, the secretion of products to the tube tip is important, but only a few genes involved in the transport of secretory vesicles during tube growth have been identified to date (Sanderfoot et al. 2001Go). In our array data, genes encoding secretory vesicle molecules (SNARE, SNAP, synaptobrevin, clathrin, etc.) were identified in the late BC and TC pollen (clusters 5, 6, 10, 11 and 19). This new finding should be important in understanding the molecular machinery for vesicle trafficking in the pollen tube.

Clusters 13, 14 and 15 showed similar predominant expression in the MEI, TET and UN tapetum, and contained known genes related to pollen wall formation and fatty acid synthesis. Rice WDA1 (Jung et al. 2006Go), which is involved in the biosynthesis of very long chain fatty acids in the anther and is necessary for pollen exine formation, was categorized in cluster 14. As potential orthologs of known fatty acid biosynthesis genes, genes similar to Arabidopsis CER6/CUT1 (Fiebig et al. 2000Go) and FIDDLEHEAD (Pruitt et al. 2000Go) were categorized in cluster 13, while genes homologous to maize glossy8 (Xu et al. 1997Go) and Arabidopsis ECR (Zheng et al. 2005Go) were in cluster 15. The tapetum cells have a central role in pollen wall formation, and transfer lipids onto the exine surface. Transcriptional regulation of the lipid-related genes observed in these clusters is undoubtedly important in male gametophyte development.

Clusters 13, 14 and 15 also contained ethylene signaling-related genes (ethylene receptor genes, EIN2, EIN3, EIL3, ETO1, ERF genes, etc.) and ubiquitin–proteasome system-related genes [genes for 26S proteasome, ubiquitin, E2 (E3)-ligase, F-box protein, SUMO, etc.]. Ethylene is an important phytohormone in senescence, fruit ripening and stress responses, and it has been reported that ethylene signaling is regulated by the proteolytic degradation of transcription factor EIN3 (Guo and Ecker 2004Go, Schwechheimer and Villalobos 2004Go). Taken together, involvement of these sets of genes in clusters 13, 14 and 15 suggests that the phenomenon of gradual tapetum shrinking from the MEI to UN stages may be regulated by ethylene-mediated senescence-like events in the tapetum.

The argonaute gene and micro RNAs (miRNAs) were observed in 13 clusters (clusters 3, 5, 7–13, 15, 16, 18 and 19). The expression profiles of these clusters were diverse, with active transcription patterns distributed through all the developmental stages of the microspore/pollen and tapetum. The argonaute proteins are integral players in all known small RNA-directed regulatory pathways (Vaucheret 2008Go), and miRNAs have emerged as key guide molecules in the regulation of various biological processes (Ramachandran and Chen 2008Go). Thus, our LM-microarray analysis indicated that transcriptional regulation by small RNAs may be important throughout the process of anther development in rice. It is noteworthy that cluster 19 contained miR166, consistent with a recent report in which miR166 was detected in mature rice anthers by RNA gel blot analysis (Fujioka et al. 2008Go).

Many other biological events still remain to be analyzed, including tapetum cell death, sporopollenin formation, pollen–stigma interaction and asymmetric mitosis. The separated transcriptomes of microspore/pollen and tapetum contain considerable useful information that will shed light on many events during anther development, and may potentially reveal novel interactions between the male gametophyte and tapetum.

Characterization of cis-elements in the male gamete- and tapetum-specific genes
Using LM-microarray technology, we found that 28,141 genes were ‘anther-expressed’. Of these, 3,468 ‘anther-specific’ genes were identified by comparison with gene expression in vegetative tissues, as described above. We were interested in identifying specific cis-elements regulating spatiotemporal gene expression in the promoter regions of anther-expressed genes. A representative sample of 140 genes, characterized as male gamete-specific, tapetum-specific, and genes expressed in both tissues in our LM-micro­array (Suwabe et al. 2008Go), were analyzed for cis-elements. A total of 164,188 cis-elements, classified into 256 cis-motifs, were found using the PLACE database. Among the 19 clusters, the frequency of the 11 common elements, as a percentage of the total cis-elements in that clusters, ranged from 51.2 to 86.3% (Table 1). Among these common elements, two were pollen-specific elements, GTGANTG10 and POLLEN1LELAT52. The GTGA motif in the promoter of the tobacco late pollen gene g10 was found to be an element responsible for gene expression in the late pollen stage (Rogers et al. 2001Go). The AGAAA motif is one of two co-dependent regulatory elements responsible for pollen-specific expression of the tomato LAT52 gene (Bate and Twell 1998Go). Because these two elements were originally characterized in tobacco and tomato, Solanaceae, and are also found as common cis-elements in pollen small cysteine-rich proteins in rice (Park et al. 2006Go), they can be regarded as common elements responsible for gene expression in pollen and widely conserved among plant species. By comparing cis-elements between the microspore/pollen-expressed clusters (Nos. 4, 5, 6, 8, 10 and 19) and the tapetum-expressed clusters (Nos 2, 7, 12, 13, 14, 15, 18 and 20), five male gamete- and 35 tapetum-specific cis-elements were found (Table 2). They included various kinds of motif elements from a wide variety of plant species, but we are far from deducing how male gamete- or tapetum-specific gene expression is controlled by these cis-elements. Interestingly, >130,000 elements were found in cluster 15, and the ATATT motif accounted for approximately 70% of elements in this cluster. In contrast, the ATATT motif was not found in any of the cis-elements identified by the PLACE database in clusters 7 and 9, even though synchronous and predominant gene expression in MEI and TET stages, and preferential gene expression at MEI and TC pollen stages, was observed in microspore and tapetum in clusters 7 and in 9, respectively. It is difficult to rule out experimental bias in these clusters, but the difference in cis-element frequencies might truly reflect biological characteristics. As discussed above, there were many kinds of cis-elements identified, indicating considerable complexity in the regulation of gene expression during male gamete and tapetum development.


View this table:
[in this window]
[in a new window]

 
Table 2 Specific motifs in the male gamete- or tapetum-expressed genes

 
Characterization of Tos17 insertion mutant lines in male gamete- and tapetum-specific genes
As a systematic rice mutant resource, Tos17 insertion lines are convenient for field cultivation, although there is preference for a Tos17 insertion sites within the rice genome (Miyao et al. 2003Go). By using the Tos17 mutant panel database (http://tos.nias.affrc.go.jp/), Tos17 mutant lines were surveyed among the 140 representative genes (Suwabe et al. 2008Go). A total of 114 independent lines were identified among 17 of these genes (~12.1%). From the mutant phenotype descriptions available on the database, several phenotypes were observed (e.g. sterile, low fertility, semi-dwarf, dense panicle and pale green leaf). In particular, nine lines (52.9%) showed sterile or low fertility phenotypes, indicating that these genes are likely to be important in pollen development (data not shown).

Conclusion
In this study, using cluster analysis, we found various spatiotemporal expression patterns in anther-expressed genes. Interestingly, many genes were synchronously expressed in microspore/pollen and tapetum. This synchronous expression in different cell types should shed light on the mechanisms of anther development.

These global gene expression data from the LM-array in microspore/pollen and tapetum will form an important basis for future research in plant reproduction, as well as other fields in plant science. In order to utilize this large body of data effectively, it will be necessary to combine several methodologies (e.g. in silico, morphological, genetic and biochemical analyses).


    Materials and Methods
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
Plant materials
Rice (O. sativa L. ssp. japonica cv. Nipponbare) plants were grown in a greenhouse under normal conditions. Developing anthers at each of five stages, MEI, TET, UN microspore, BC pollen and TC pollen, were collected by confirming their developmental stage by microscopy, using one of six anthers from each flower. The stage classification of anthers is described in detail in Suwabe et al. (2008Go) in this issue.

Laser microdissection
The fixation, dehydration and paraffin embedding of anthers were described in Suwabe et al. (2008Go) in this issue. Briefly, paraffin-embedded sections were cut to a thickness of 16 µm and mounted on PEN membrane glass slides (Molecular Devices, Ontario, Canada) for LM. To remove paraffin, slides were immersed in 100% xylene (twice), 50% xylene/50% ethanol (v/v) and 100% ethanol for 5 min for each step, and then air-dried completely at room temperature. Three or four individual flowers were used for each LM experiment. LM was performed using the Veritas Laser Microdissection System LCC1704 (Molecular Devices). The target cells that fused to the CapSure Macro LCM cap (Molecular Devices) were collected by removing the cap from the tissue section.

Microarray and cluster analysis
Microarray analysis was performed as described by Suwabe et al. (2008Go) in this issue. Briefly, total RNAs were extracted from LM cells with a PicoPureTM RNA isolation kit (Molecular Devices), and were quantified with a Quant-iTTM RiboGreen RNA reagent and kit (Invitrogen, San Diego, CA, USA). A rice 44K oligo microarray (Agilent Technologies, Palo Alto, CA, USA) with ~42,000 oligonucleotides based on the Rice Annotation Project (RAP) was used. Fluorescent probe labeling using the oligo-dT-T7 strand-specific amplification method and hybridization were performed according to the manufacturer's instructions (Agilent Technologies) with slight modifications. The slide images were scanned with a DNA microarray scanner (Agilent Technologies, http://www.home.agilent.com/) using the manufacturer's Feature Extraction software. Microarray experiments using each stage sample were replicated at least three times.

Microarray data were statistically analyzed by the vari­ance stabilization normalization (VSN) method using R software (http://www.r-project.org/). Normalization by VSN, Z-scoring and one-way ANOVA analysis of variance were performed by Bioconductor (http://www.bioconductor.org/). Advanced data analyses including cluster analysis (K-means) were performed using Microsoft Excel and TIGR MeV (Saeedet al. 2003Go; http://www.jcvi.org/cms/research/software/#c622).

Gene ontology analysis
GO Slim terms for genes (probes) were obtained from the TIGR database site (http://www.tigr.org/tdb/e2k1/osa1/batch_downloads.html). To obtain GO terms, locus identifiers (IDs) of TIGR were assigned for probes by the RAP-DB and our perl scripts. For each cluster, GOSlim terms for locus IDs of TIGR were searched by the TIGR web site and counted.

Cis-element analysis
For identifying the cis-elements of anther-specific genes, the promoter region of each gene was compared with the PLACE database (http://www.dna.affrc.go.jp/PLACE/). Identified cis-elements were classified by sequence similarity.

Identification of Tos17 insertion mutant lines
The Tos17 mutant panel database (http://tos.nias.affrc.go.jp/) was used for identification of Tos17 insertion mutant lines. The number of mutant lines was counted, and their mutant phenotypes were surveyed.


    Supplementary Material
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
Supplementary Material are available at PCP Online.


    Funding
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
The Ministry of Education, Culture, Sports, Science, and Technology of Japan (MEXT) Grants-in-Aid for Special Research on Priority Areas (Nos. 18075003 and 18075012 to M.W., 18075005 to N.T., 18075006 to M.M., 18075009 to N.K., 19043009 to T.I., 19043015 and 19651084 to K.Y.); the Japan Society for Promotion of Science (JSPS) Grants-in-Aid for the 21st Century Center of Excellence Program (to M.W.); the Ministry of Agriculture, Forestry and Fisheries of Japan (MAFF) grant for the Integrated Research Project for Plant, Insect, and Animal using Genome Technology (Nos. IPG-0003 to M.M., IPG-0012 to M.N., IPG-0019 to M.W.); a Research Fellowship for Postdoctoral Fellowships for Young Researchers of JSPS (to K.S.).


    Acknowledgments
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
The authors thank Drs. Jun Hidema, Yukihiro Ito and Jong-In Park (Tohoku University) for their valuable comments, and Hisae Kamakura (University of Tokyo), Ayako Chiba (Iwate University) Kosuke Matsumoto, Hiromi Shoji, Yuta Tsunaga and Ayumi Yamakawa (Tohoku University) for technical assistance.


    Footnotes
 
11These authors contributed equally to this work. Back


    References
 Top
 Abstract
 Introduction
 Results and Discussion
 Materials and Methods
 Supplementary Material
 Funding
 Acknowledgments
 References
 
Ariizumi T, Hatakeyama K, Hinata K, Sato S, Kato T, Tabata S, et al. A novel male-sterile mutant of Arabidopsis thaliana, faceless pollen-1, produces pollen with a smooth surface and an acetolysis-sensitive exine. Plant Mol. Biol (2003) 53:107–116.[CrossRef][Web of Science][Medline]

Asano T, Masumura T, Kusano H, Kikuchi S, Kurita A, Shimada H, et al. Construction of a specialized cDNA library from plant cells isolated by laser capture microdissection: toward compre-hen­sive analysis of the genes expressed in the rice phloem. Plant J (2002) 32:401–408.[CrossRef][Web of Science][Medline]

Bate N, Twell D. Functional architecture of a late pollen promoter: pollen-specific transcription is developmentally regulated by multiple stage-specific and co-dependent activator elements. Plant Mol. Biol (1998) 37:859–869.[CrossRef][Web of Science][Medline]

Cheung AY, Wu H-M. Structural and signaling networks for the polar cell growth machinery in pollen tubes. Annu. Rev. Plant Biol (2008) 59:547–572.[CrossRef][Medline]

Day RC, Grossniklaus U, Macknight RC. Be more specific! Laser-assisted microdissection of plant cells. Trends Plant Sci (2005) 10:397–406.[CrossRef][Web of Science][Medline]

Endo M, Kokubun T, Takahata Y, Higashitani A, Tabata S, Watanabe M. Analysis of expressed sequence tags of flower buds in Lotus japonicus. DNA Res (2000) 7:213–216.[Abstract]

Endo M, Matsubara H, Kokubun T, Masuko H, Takahata Y, Tsuchiya T, et al. The advantages of cDNA microarray as an effective tool for identification of reproductive organ-specific genes in a model legume, Lotus japonicus. FEBS Lett (2002) 514:229–237.[CrossRef][Web of Science][Medline]

Endo M, Tsuchiya T, Saito H, Matsubara H, Hakozaki H, et al. Identification and molecular characterization of novel anther-specific genes in japonica rice, Oryza sativa L. by using cDNA microarray. Genes Genet. Syst (2004) 79:213–226.[CrossRef][Web of Science][Medline]

Fiebig A, Mayfield JA, Miley NL, Chau S, Fischer RL, Preuss D. Alterations in CER6, a gene identical to CUT1, differentially affect long-chain lipid content on the surface of pollen and stem. Plant Cell (2000) 12:2001–2008.[Abstract/Free Full Text]

Franklin-Tong VE. Signaling in pollination. Curr. Opin. Plant Biol (1999) 2:490–495.[CrossRef][Web of Science][Medline]

Fujioka T, Kaneko F, Kazama T, Suwabe K, Suzuki G, Makino A. Identification of small RNA products in the late developmental stage of rice anthers. Genes Genet. Syst (2008) 83:281–284.[CrossRef][Web of Science][Medline]

Goldberg RB, Beals TP, Sanders PM. Anther development: basic principles and practical application. Plant Cell (1993) 5:1217–1229.[Free Full Text]

Golovkin M, Reddy ASN. A calmodulin-binding protein from Arabidopsis has an essential role in pollen germination. Proc. Natl Acad. Sci. USA (2003) 100:10558–10563.[Abstract/Free Full Text]

Guo H, Ecker JR. The ethylene signaling pathway: new insights. Curr. Opin. Plant Biol (2004) 7:40–49.[CrossRef][Web of Science][Medline]

Hoekstra FA, Bruinsma J. Protein synthesis of binucleate and trinucleate pollen and its relationship to tube emergence and growth. Planta (1979) 146:559–566.[CrossRef][Web of Science]

Honys D, Twell D. Comparative analysis of the Arabidopsis pollen transcriptome. Plant Physiol (2003) 132:640–652.[Abstract/Free Full Text]

Honys D, Twell D. Transcriptome analysis of haploid male gametophyte development in Arabidopsis. Genome Biol (2004) 5:R85.[CrossRef][Medline]

Huang J-U, Yang Z. Small GTPases and spatiotemporal regulation of pollen tube growth. In: In The Pollen Tube Growth.—Malh R, ed. (2006) Berlin: Springer-Verlag. 95–116.

Jiang SY, Cai M, Ramachandran S. The Oryza sativa no pollen (Osnop) gene plays a role in male gametophyte development and most likely encodes a C2-GRAM domain-containing protein. Plant Mol. Biol (2005) 57:835–853.[CrossRef][Web of Science][Medline]

Jung K-I, Han M-J, Lee D-Y, Lee Y-S, Schreiber L, et al. Wax-deficient anther1 is involved in cuticle and wax production in rice anther walls and is required for pollen development. Plant Cell (2006) 18:3015–3032.[Abstract/Free Full Text]

Kerk NM, Ceserani T, Tausta SL, Sussex IM, Nelson TM. Laser capture microdissection of cells from plant tissues. Plant Physiol (2003) 132:27–35.[Abstract/Free Full Text]

Malhó R, Liu Q, Rato C, Camacho L, Dinis A. Signalling pathways in pollen germination and tube growth. Protoplasma (2006) 228:21–30.[CrossRef][Web of Science][Medline]

Mariani C, Wolters-Arts M. Complex waxes. Plant Cell (2000) 12:1795–1798.[Free Full Text]

Mascarenhas NT, Bashe D, Eisenberg A, Willing RP, Xiao CM, Mascarenhas JP. Messenger RNAs in corn pollen and protein synthesis during germination and pollen tube growth. Theor. Appl. Genet (1984) 68:323–326.[CrossRef][Web of Science]

Masuko H, Endo M, Saito H, Hakozaki H, Park J-I, et al. Anther-specific genes, which are expressed through microsporogenesis, are temporally and spatially regulated in model legume, Lotus japonicus. Genes Genet. Syst (2006) 81:57–62.[CrossRef][Web of Science][Medline]

Miyao A, Tanaka K, Murata K, Sawaki H, Takeda S, Abe K, et al. Target site specificity of the Tos17 retrotransposon shows a preference for insertion within genes and against insertion in retro-transposon-rich regions of the genome. Plant Cell (2003) 15:1771–1780.[Abstract/Free Full Text]

Nakazono M, Qui F, Brsuk LA, Schnable PS. Laser-capture microdissection, a tool for the global analysis of gene expression in specific plant cell type: identification of genes expressed differentially in epidermal cell or vascular tissues of maize. Plant Cell (2003) 15:583–596.[Abstract/Free Full Text]

Nelson T, Tausta SL, Gandotra N, Liu T. Laser microdissection of plant tissue: what you see is what you get. Annu. Rev. Plant Biol (2006) 57:181–201.[CrossRef][Medline]

Ohtsu K, Takahashi H, Schnable PS, Nakazono M. Cell type-specific gene expression profiling in plants by using a combination of laser microdissection and high-throughput technologies. Plant Cell Physiol (2007) 48:3–7.[Abstract/Free Full Text]

Park JI, Hakozaki H, Endo M, Takada Y, Ito H, Uchida M, et al. Molecular characterization of mature pollen-specific genes encoding novel small cysteine-rich proteins in rice (Oryza sativa L.). Plant Cell Rep (2006) 25:466–474.[CrossRef][Web of Science][Medline]

Pruitt RE, Vielle-Calzada J-P, Ploense SE, Grossniklaus U, Lolle SJ. FIDDLEHEAD, a gene required to suppress epidermal cell interactions in Arabidopsis, encodes a putative lipid biosynthetic enzyme. Proc. Natl Acad. Sci. USA (2000) 97:1311–1316.[Abstract/Free Full Text]

Ramachandran V, Chen X. Small RNA metabolism in Arabidopsis. Trends Plant Sci (2008) 13:369–374.

Rogers HJ, Bate N, Combe J, Sullivan J, Sweetman J, Swan C, et al. Functional analysis of cis-regulatory elements within the promoter of the tobacco late pollen gene g10. Plant Mol. Biol (2001) 45:577–585.[CrossRef][Web of Science][Medline]

Saeed AI, Sharov V, White J, Li J, Liang W. TM4: a free, open-source system for microarray data management and analysis. Biotechniques (2003) 34:374–378.[Web of Science][Medline]

Sanderfoot AA, Pilgrim M, Adam L, Raikhel NV. Disruption of individual membranes of Arabidopsis syntaxin gene families indicates each has essential functions. Plant Cell (2001) 13:659–666.[Abstract/Free Full Text]

Schwechheimer C, Villalabos LIAC. Cullin-containing E3 ubiquitin ligases in plant development. Curr. Opin. Plant Biol (2004) 7:677–686.[CrossRef][Web of Science][Medline]

Schrauwen JAM, Groot PFM, Herpen MMA, Lee T, Reynen WH, Weterings KAP, et al. Stage-related expression of mRNAs during pollen development in lily and tobacco. Planta (1990) 182:298–304.[Web of Science]

Scott RJ, Spielman M, Dickinson HG. Stamen structure and function. Plant Cell (2004) 16:S46–S60.[Free Full Text]

Shimamura K, Ishimizu T, Nishimura K, Matsubara K, Kodama H, Watanabe H, et al. Analysis of expressed sequence tags from petunia flowers. Plant Sci (2007) 173:495–500.[CrossRef][Web of Science]

Suwabe K, Suzuki G, Takahashi H, Shiono K, Endo M. Separated transcriptomes of male gametophyte and tapetum in rice: validity of a laser microdissection (LM) microarray. Plant Cell Physiol (2008) 49:1407–1416.[Abstract/Free Full Text]

Vaucheret H. Plant ARGONAUTES. Trends Plant Sci (2008) 13:350–358.[CrossRef][Web of Science][Medline]

Xu X, Dietrich CR, Delledonne M, Xia Y, Wen T-J, Robertson DS, et al. Sequence analysis of the cloned glossy8 gene of maize suggests that it may code for a β-ketoacyl reductase required for the biosynthesis of cuticular wax. Plant Physiol (1997) 115:501–510.[Abstract]

Zheng H, Rowland O, Kunst L. Disruptions of the Arabidopsis enoyl-CoA reductase gene reveal an essential role for very-long-chain fatty acid synthesis in cell expansion during plant morphogenesis. Plant Cell (2005) 17:1467–1481.[Abstract/Free Full Text]


(Received July 20, 2008; Accepted August 30, 2008)
Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?

Related articles in PCP:

Towards a Comprehensive Understanding of Molecular Mechanisms of Sexual Reproduction in Higher Plants
Masao Watanabe
PCP 2008 49: 1404-1406. [Extract] [FREE Full Text]  



This article has been cited by other articles:


Home page
Plant Cell PhysiolHome page
M. Endo, T. Tsuchiya, K. Hamada, S. Kawamura, K. Yano, M. Ohshima, A. Higashitani, M. Watanabe, and M. Kawagishi-Kobayashi
High Temperatures Cause Male Sterility in Rice Plants with Transcriptional Alterations During Pollen Development
Plant Cell Physiol., November 1, 2009; 50(11): 1911 - 1922.
[Abstract] [Full Text] [PDF]


Home page
Plant Cell PhysiolHome page
G. Suzuki
Recent Progress in Plant Reproduction Research: The Story of the Male Gametophyte through to Successful Fertilization
Plant Cell Physiol., November 1, 2009; 50(11): 1857 - 1864.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
M. Mihara, T. Itoh, and T. Izawa
SALAD database: a motif-based database of protein annotations for plant comparative genomics
Nucleic Acids Res., October 23, 2009; (2009) gkp831v1.
[Abstract] [Full Text] [PDF]


Home page
J Exp BotHome page
Compiled by, F. Tooke, T. Chiurugwi, and N. Battey
Flowering Newsletter bibliography for 2008
J. Exp. Bot., June 23, 2009; (2009) erp154v1.
[Full Text] [PDF]


Home page
J Exp BotHome page
Z. A. Wilson and D.-B. Zhang
From Arabidopsis to rice: pathways in pollen development
J. Exp. Bot., April 1, 2009; 60(5): 1479 - 1492.
[Abstract] [Full Text] [PDF]


Home page
Plant Cell PhysiolHome page
K. Kurasawa, A. Matsui, R. Yokoyama, T. Kuriyama, T. Yoshizumi, M. Matsui, K. Suwabe, M. Watanabe, and K. Nishitani
The AtXTH28 Gene, a Xyloglucan Endotransglucosylase/Hydrolase, is Involved in Automatic Self-Pollination in Arabidopsis thaliana
Plant Cell Physiol., February 1, 2009; 50(2): 413 - 422.
[Abstract] [Full Text] [PDF]


Home page
Plant Cell PhysiolHome page
K. Suwabe, G. Suzuki, H. Takahashi, K. Shiono, M. Endo, K. Yano, M. Fujita, H. Masuko, H. Saito, T. Fujioka, et al.
Separated Transcriptomes of Male Gametophyte and Tapetum in Rice: Validity of a Laser Microdissection (LM) Microarray
Plant Cell Physiol., October 1, 2008; 49(10): 1407 - 1416.
[Abstract] [Full Text] [PDF]


Home page
Plant Cell PhysiolHome page
M. Mihara, T. Itoh, and T. Izawa
In Silico Identification of Short Nucleotide Sequences Associated with Gene Expression of Pollen Development in Rice
Plant Cell Physiol., October 1, 2008; 49(10): 1451 - 1464.
[Abstract] [Full Text] [PDF]


Home page
Plant Cell PhysiolHome page
M. Watanabe
Towards a Comprehensive Understanding of Molecular Mechanisms of Sexual Reproduction in Higher Plants
Plant Cell Physiol., October 1, 2008; 49(10): 1404 - 1406.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
49/10/1417    most recent
pcn128v2
pcn128v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Related articles in PCP
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (4)
Google Scholar
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Agricola
Right arrow Articles by Hobo, T.
Right arrow Articles by Matsuoka, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?