It was isozyme studies that first suggested widespread overlap of gametophytic and sporophytic gene expression. In various species, 60% to 72% of isozymes analyzed were expressed in both gametophytic and sporophytic tissues (Tanksley et al., 1981; Sari-Gorla et al., 1986; Pedersen et al., 1987). Studies of the kinetics of [³H]cDNA to poly(A⁺) RNA hybridization confirmed this overlap and quantified the complexity of haploid gene expression. These estimates suggested that Tradescantia paludosa and maize (Zea mays) pollen contain approximately 20,000 to 24,000 different mRNA sequences (Willing and Mascarenhas, 1984; Willing et al., 1988). The complexity of pollen mRNA populations were significantly lower when compared with roots, in both T. paludosa and maize, where it was estimated that 30,000 genes were expressed. The overlap between pollen and root mRNAs in these species was estimated to be 65%.

Hybridization studies of pollen cDNA libraries also suggested that the majority of pollen-expressed mRNAs showed gametophytic-sporophytic overlap. In fact, only 10% of pollen-expressed mRNAs were considered to be pollen specific (Stinson et al., 1987; Mascarenhas, 1990). However, these estimations did not account for cross-hybridization between closely related gene family members. Subsequent studies have shown that a number of pollen-specific genes have closely related counterparts expressed in sporophytic tissues (Belostotsky and Meagher, 1993, 1996; Brander and Kuhlemeier, 1995; Lopez et al., 1996), but the extent to which this characterizes the male transcriptome is unknown.

The assembly of gene-specific data concerning the male gametophytic transcriptome has revealed approximately 150 different genes, assigned to 16 distinct function groups, with strong evidence for pollen-specific expression in about 30 (Twell, 2002). Thus, compared with the estimated 20,000 pollen-expressed mRNAs and 2,000 to 7,000 pollen-specific mRNAs (Willing and Mascarenhas, 1984; Willing et al., 1988), current knowledge of the pollen transcriptome is strongly deficient. Microarray technology now provides the opportunity to reveal this unknown fraction and to compile comprehensive data concerning the extent of overlap between male gametophytic and sporophytic gene expression.

In addition to its vital role in sexual reproduction, pollen provides a microcosm of cellular development, which is an attractive system in which to dissect the fundamental processes of cell growth and division, cellular differentiation, and intercellular communication (Bedinger, 1992; Twell, 1994, 2002). An essential step toward the detailed understanding of these processes is to define the transcriptome at the cellular level. In plants, the male gametophyte is a uniquely accessible cell type for such studies, enabling RNA isolation from a pure cell population and transcriptome analysis of what is essentially a single cell type poised for an explosive growth phase.

Here, we report the first genome-wide view of the pollen transcriptome for the model species Arabidopsis. We demonstrate the unique composition of the pollen transcriptome based on comparative analysis with sporophytic tissues. The identities of 992 pollen-expressed mRNAs were identified, nearly 40% of which were detected specifically in pollen. We also analyze the functional composition of the pollen transcriptome and demonstrate common patterns of mutually exclusive gene expression between pollen and sporophytic tissues for different gene family members.

Affymetrix Arabidopsis 8K GeneChip arrays were used to explore the pollen transcriptome in comparison with that of the developing sporophyte at various stages of plant development. Microarrays were hybridized with cRNA probes made from total RNA isolated from mature pollen. Microarray hybridization data for pollen RNA prepared from two independently grown plant populations were compared. Only genes with a positive hybridization signal accompanied by a detection call value of 1 in both experiments were scored as pollen expressed. Sporophytic Affymetrix 8K GeneChip data were downloaded from the GARNet Web site (http://www.Arabidopsis.info). This provided transcriptome data for Arabidopsis plants at open cotyledon stage (stage 1.00; Boyes et al., 2001), four rosette leaf stage (1.08), first visible flower bud stage (5.10), and first open flower stage (6.00). Sporophytically expressed genes were identified using the same algorithm used for pollen microarray data.

The Affymetrix 8K GeneChip arrays contain oligonucleotide probes representing 7,792 annotated genes with Arabidopsis Genome Initiative (AGI) numbers. This represents approximately 28% of the total of 27,117 annotated protein-coding genes in Arabidopsis (The Institute for Genomic Research, 2002). In pollen, 992 genes consistently gave a positive hybridization signal, representing nearly 13% of the unigene targets present on the micorarray (Fig. 1). Complete microarray data are publicly available at the Nottingham Arabidopsis Stock Centre (NASC) microarray database (http://Arabidopsis.info/prototype). The fully annotated list of pollen-expressed mRNAs is available as a supplemental information (Supplemental Fig. 1; supplemental data can be viewed at http://www.plantphysiol.org). Given that the contribution of sperm cells to the detectable pollen transcriptome is likely to be negligible (Xu et al., 2002), these data essentially represent the transcriptome of the vegetative cell present in mature pollen.

Figure 1. Quantitative evaluation of gametophytic and sporophytic gene expression. A, Relative proportion of stage-specific mRNAs in mature pollen and in four sporophytic stages: open cotyledon (S 1.00), four leaf rosette (S 1.08), first visible flower bud (S (more ...)

The number of diverse mRNAs expressed in pollen calculated from microarray analysis was significantly lower by 30% to 60% (Fig. 1A) than that of other sporophytic samples. This may not be surprising given that all sporophytic RNA samples represented more complex tissues. However, a striking feature was that 39% of pollen-expressed mRNAs were detected exclusively in pollen. In all other sporophytic samples, the percentage of stage-specific mRNAs was less than 8%, reflecting the broadly similar tissue compositions and cellular functions throughout sporophytic development. Although less pronounced than between sporophytic stages, there was a significant proportion (61%) of pollen-expressed mRNAs that were also expressed in one or more sporophytic stages.

Analysis of the distribution of mRNAs among three mRNA abundance classes; high (up to 10-fold less than the maximum signal), medium (10- to 100-fold less), and low (more than 100-fold less; Fig. 1C) revealed that the proportion of mRNAs forming the high-abundance class in pollen (17%) was significantly above that in sporophytic samples (5.2%–8.2%). There was a tendency toward reduction of the medium-abundance class and expansion of the low-abundance class during sporophytic development. Fifty-five percent of pollen-expressed mRNAs forming the highly abundant class were only detected in pollen. Moreover, the proportion of pollen-specific genes (55%) was much higher than in the two lower abundance classes (36%). The high proportion of the pollen transcriptome contributed by pollen-specific mRNAs highlights the potential significance and regulatory specialization of this abundant subset of the pollen transcriptome.

Figure 2. Verification of 8K GeneChip expression data by RT-PCR. Expression of four genes considered pollen-specific according to 8K GeneChip data in roots (R), stems (S), leaves (L), flowers (F), and pollen (P) were verified by RT-PCR. These are compared with (more ...)

The purity of pollen populations and isolated RNA was confirmed by the lack of expression of the Rubisco small subunit 1b gene At5g38430 by RT-PCR (Fig. 2). Moreover, all six genes encoding light-harvesting chlorophyll a/b-binding proteins (LHCB2, At3g27690; LHCB3, At5g54270; LHCB4, At5g01530; LHCB4.3, At2g40100; LHCB5, At4g10340; and LHCB6, At1g15820) present on the microarray were abundantly expressed in the sporophyte and were among the 50 most abundantly expressed genes. However, none of these genes gave positive signals in pollen (data not shown).

The recent compilation of approximately 150 pollen-expressed genes from various species included 23 genes found in Arabidopsis (Twell, 2002). A positive microarray hybridization signal in pollen was verified for 10 of 11 of these genes that were present on the 8K GeneChip. The only exception was the lack of a microarray signal for AtSUC1 gene (At1g71880) in mature pollen. However, AtSUC1 expression was not examined in mature dehiscent pollen (Stadler et al., 1999). With regard to Arabidopsis homologs of published pollen-expressed or -specific genes, two Arabidopsis homologs (At1g55560 and At1g55570) of the abundantly expressed tobacco (Nicotiana tabacum) pollen-specific gene ntp303 (Weterings et al., 1992) were detected exclusively in pollen and were among the 25 most abundantly expressed genes. The same was true for At1g14420 (At59 in Kulikauskas and McCormick, 1997), a homolog of the tomato (Lycopersicon esculentum) late anther-specific gene lat59 encoding a pectate lyase protein (Wing et al., 1989), whose expression ranked slightly lower at 31 of 992. Similarly, several members of the polygalacturonase gene family (At2g23900, At3g07820, At3g07850, and At4g33440) were detected exclusively in pollen, in accord with data previously published for maize and tobacco (Allen and Lonsdale, 1993). Moreover, four of five actin genes expressed in reproductive tissues (see McDowell et al., 1996) were present on the 8K GeneChip. All four belonged into the high-abundance class in pollen including two that were pollen-specific (ACT4 [At5g59370] rank 8 and ACT12 [At3g46520] rank 172), and the remaining two were highly pollen-enriched (ACT1 [At2g37620] rank 53 and ACT11 [At3g12110] rank 74). Of two profilin genes present on the 8K GeneChip, PRF4 (At4g29340) was expressed exclusively in pollen and was ranked 84, whereas PRF1 (At2g19760) showed low-level constitutive expression. Again, these data are in accord with previously published data (Huang et al., 1996). Data presented in this paragraph can be found in Supplemental Figure 1.

Figure 3. Scatter plots comparing gene expression levels in pollen and at different sporophytic developmental stages. The expression levels of individual genes were normalized using a scale of 0 to 100. Genes co-expressed in pairs of transcriptome data sets (more ...)

We also performed quantitative analysis of the extent of gametophytic-sporophytic overlap that resolved all 31 possible intersections between the five independent transcriptome data sets. The outcome is displayed using a Venn diagram (Fig. 4; Ruskey, 1997), directly outlining numbers of genes comprising each of the 31 categories. For quantitative analyses, gene expression levels were not taken into account, each gene was treated only as expressed or not. Of 3,282 genes expressed in at least one sample, only 287 (8.7%) were shared by all five. Pollen expressed the highest percentage (39%) of specific genes compared with all sporophytic samples. In contrast, 830 genes not expressed in pollen, representing 29% of the total number of sporophytically expressed genes, were co-expressed in all four sporophytic data sets. The percentage of expressed genes shared by the three post-seedling data sets reached 65% of the total number of sporophytically expressed genes. Pollen and seedlings co-expressed 312 genes, representing the minimum overlap between all data sets. Moreover, 92% of these represented constitutive genes expressed in all five data sets. Taken together, Venn and scatter plot analyses highlight the unique composition of the pollen transcriptome in comparison with the developing sporophyte.

Figure 4. Quantification of gametophytic-sporophytic and sporophytic-sporophytic overlaps of transcriptome profiles. Venn diagram illustrating all 31 possible intersections between five independent transcriptome data sets including the number of genes composing (more ...)

Figure 5. Distribution of expressed mRNAs among gene function categories. Pie charts display each gene function category according to the number of expressed genes (B, D, F, and H) or their relative signal intensity (C, E, G, and I). The number of expressed (more ...)

Sporophytic tissues, represented here by the four rosette leaf stage 1.08 (Fig. 5, B and C), displayed transcriptome profiles that were broadly similar to each other, but distinct from pollen. Compared with the distribution of genes on 8K GeneChip, the most significant increase was observed in the number of genes involved in basic and energy metabolism, transport, and protein synthesis. Unlike pollen, genes involved in energy metabolism, transport, protein synthesis, and stress were preferentially expressed. The dramatic increase in expression of genes involved in energy metabolism from 2% to 11% was due to the extraordinary expression of genes involved in the light phase of photosynthesis. Apart from cell wall proteins, these categories account for the major differences between gametophytic and sporophytic transcriptomes.

The preferential expression of genes belonging to particular functional categories became even more apparent when only genes detected specifically in pollen were considered (Fig. 5, F and G). Narrowing of both transcription and protein synthesis categories was accompanied by a dramatic overrepresentation of genes involved in cell wall metabolism, signaling, and the cytoskeleton. In particular, the sum of the hybridization signals for genes in these last three classes reached 45% of the total signal intensity for all pollen-specific genes. On the contrary, the complementary set of pollen-expressed genes showing gametophytic-sporophytic overlap resembled the distribution of functional categories in the sporophyte (Fig. 5, H and I). Apparent under-representation was found for only two categories, energy metabolism and protein synthesis, reflecting their significant under-representation in the complete mature pollen transcriptome. Conversely, mRNAs encoding cytoskeletonrelated proteins were over-represented among subsets showing pollen-specific or gametophytic-sporophytic overlap.

Table I. Gene families selected for the hierarchical clustering analysis (EPCLUST)

Figure 6. Expression profiles of the glycoside hydrolase gene superfamily. Hierarchical cluster analysis of 51 members of the glycoside hydrolase superfamily in mature pollen (1) and in four sporophytic data sets: S 1.00 (2), S 1.08 (3), S 5.10 (4), and S 6.00 (more ...)

Figure 7. Expression profiles of the RLK gene family. Hierarchical cluster analysis of 62 members of RLK gene family in mature pollen (1) and in four sporophytic data sets (2–5) according to legend of Figure 6. For each expression profile, the AGI number (more ...)

Four selected families and superfamilies contained proteins involved in cell wall metabolism. Within three superfamilies of carbohydrate-modifying enzymes and the family of expansins, large and unique subsets of genes were expressed in pollen (Fig. 6; Supplemental Figs. 2, 3, 4). Approximately 40% of genes composing glycoside hydrolase families were present on 8K GeneChip, one-third of which were expressed at least at one sample. Four members of glycoside hydrolase family GHF28 (polygalacturonases) were only detected in pollen. These polygalacturonase genes were highly expressed within high (At3g07820, At3g07850, and At2g23900) and medium (At4g33440) mRNA abundance classes. Within five pollen-expressed glycoside hydrolase families, GHF1 (β-glycosidases), GHF16 (xyloglucan endotransglycosylases), GHF28 (polygalacturonases), GHF35 (β-galactosidases), and GHF9 (glucanases/cellulases), there was no overlap between gene family members expressed in pollen and the sporophyte.

Among glycosyltransferase gene families, we identified 11 pollen-expressed genes, six of which were highly expressed and pollen specific (Table I; Supplemental Fig. 2). Four of these belonged to the GTF2 family (glucosyltransferases), including two encoding cellulose synthases (At2g33100 and At4g38190). Therefore pollen expresses specific subsets of enzymes involved in cellulose synthesis as well as degradation (At4g11050 from GHF9). We observed the same general pattern for a gene family encoding other abundant cell wall proteins, pectinesterases (Supplemental Fig. 3), but not for expansins (Supplemental Fig. 4).

The RLK superfamily was well represented on the 8K GeneChip, with 226 RLK genes or 37% of the total number of RLKs in the Arabidopsis genome. Sixty-two RLK genes distributed among 25 subfamilies were expressed in at least one sample (Fig. 7). The large majority (21 of 23) of pollen-expressed RLKs were detected exclusively in pollen. Furthermore, specificity for pollen versus sporophytic expression was observed not only for individual RLK family members, but also at the level of RLK subfamilies. Of a total of nine pollen-expressed RLK subfamilies, four did not show sporophytic expression. These were Pro-rich extensin-like receptor kinases, the receptor-like cytoplasmic kinase subfamily IX, Crinkly4-like, and the Leu-rich repeat subfamily VI. Five pollen-expressed RLK subfamilies were coexpressed in sporophyte, and the remaining 16 subfamilies were expressed only in the sporophyte.

Genes encoding translation initiation factors more frequently exhibited overlap between pollen and sporophytic expression compared with the other gene families analyzed (Supplemental Fig. 5). This overlap was relatively high (two-thirds of pollen-expressed genes) with the exception of poly(A)binding proteins (PAB). Two PABs were pollen specific (PAB-like [At2g36660] and PAB6 [At3g16380]), one gene was expressed in the sporophyte, but not in pollen (PAB1 [At2g23350]), and one was not detectably expressed (PAB3 [At1g22760]).

To identify, on a genome-wide scale, genes involved in pollen functions and those with overlapping expression between the male gametophyte and the sporophyte, we analyzed the transcriptome of mature Arabidopsis pollen using microarrays. This represents the first attempt to characterize the transcriptome of gametophytic cell types on a genomewide scale. Gene-by-gene approaches have previously identified only 20 different genes expressed in Arabidopsis pollen (Twell, 2002), such that the data sets generated here provide a 50-fold increase in knowledge of the number, identity, and relative expression levels of pollen-expressed genes in Arabidopsis. In the sporophyte, transcriptome profiling at the single-cell level has recently been reported for leaf epidermal and mesophyll cells in Arabidopsis using expressed sequence tag filter arrays (Brandt et al., 2002). Of 16,000 expressed sequence tags, 680 showed expression in one or both cells types, with 3% to 14% of sequences uniquely expressed in epidermis or mesophyll, respectively.

The pollen transcriptome was analyzed using Affymetrix Arabidopsis 8K GeneChip arrays harboring probe sets for 7,792 annotated genes (Institute for Genomic Research, 2002). Of these, 992 genes gave positive hybridization signals in pollen. Taking into account that the proportion of genes embedded on the 8K microarray is approximately 28%, we estimate the total number of pollen-expressed genes in Arabidopsis to be over 3,500. According to the criteria given, 39% of pollen-expressed genes were considered pollen specific. Although there were minor differences in the growth conditions of plants used for gametophytic and sporophytic samples, given the wide variety of developmental stages and the timing of sample collection, their comparison with pollen is unlikely to account for much of the variation of gene expression between pollen and sporophytic tissues. Therefore, we estimate the total number of genes preferentially or specifically expressed in Arabidopsis to be more than 1,400. On the contrary, a significantly lower percentage (3% and 14%, respectively) of mRNAs uniquely expressed in epidermal and mesophyll cells was reported (Brandt et al., 2002). This degree of cell-specific expression is likely to be severely overestimated because it does not take into account a number of other cell types that make up the sporophyte.

The fact that 1 in 20 Arabidopsis genes is preferentially or specifically active in pollen highlights the evolution and potential importance of transcriptional regulatory mechanisms that operate to specify gametophytic gene expression. Previous functional studies have defined multiple cis-regulatory elements involved in such transcriptional regulation in several species (for review, see Twell, 2002). Microarray data in a single species now provide the opportunity to use bioinformatic approaches to accelerate the identification of conserved regulatory elements that regulate pollen-specific expression.

We propose that 1,400 pollen-specific genes represent an ample pool for selection for pollen fitness at microgametophytic level as suggested by Mulcahy et al. (1996). On the contrary, the observed 61% gametophytic-sporophytic overlap provides the potential to improve the fitness of the sporophytic generation through gametophytic competition and selection. Our results strongly support the extent of gametophytic-sporophytic overlap previously estimated from isozyme and hybridization kinetic studies. There was 72% overlap for isozyme profiles in maize (Sari-Gorla et al., 1986) and 60% in barley (Pedersen et al., 1987) and tomato (Tanksley et al., 1981). Similarly, hybridization kinetic studies in T. paludosa (Willing and Mascarenhas, 1984) and maize (Willing et al., 1988) estimated that 65% of pollen-expressed transcripts were also present in roots. Here, the identity and relative expression levels of 605 genes that show gametophytic-sporophytic overlap are reported.

There were significant differences between estimates for the number of pollen-expressed genes in Arabidopsis based on microarray data and previous estimates based on mRNA reassociation kinetics in T. paludosa and maize (Willing and Mascarenhas, 1984). Although our data are likely to be minimum estimates that would lead to the underestimation of very low-abundance transcripts, we suggest that 20,000 pollen-expressed genes may be a significant overestimation. Despite a large difference in genome size, the total number of genes in rice is estimated to be only 1.2- to 1.9-fold larger than in Arabidopsis (Goff et al., 2002). We suggest that it is unlikely that T. paludosa and maize pollen will express 5- to 6-fold more genes than Arabidopsis.

Our results demonstrated reduced complexity of the pollen transcriptome compared with the sporophyte that was associated with the expression of a greater number of abundant mRNAs. This is in accord with previous findings in T. paludosa and maize (Willing and Mascarenhas, 1984). These observations are consistent with the limited, but highly specialized activities of mature pollen grains that are prepared for rapid germination and pollen tube growth (see Mascarenhas, 1990; Taylor and Hepler, 1997). Mature pollen is charged with mRNA and proteins in preparation for extensive cell wall synthesis and the establishment of a prominent actin-rich cytoskeleton. Although comparable transcriptome data are not available for other differentiated plant cell types, there could be similar requirements for abundant mRNAs sets in rapidly growing cells such as root hairs.

We demonstrated the pollen-specific expression of discrete subsets of genes within several gene families. The pollen-specific expression of specific members of gene families, such as those encoding cytoskeletal proteins and proteins involved in carbohydrate metabolism, could influence pollen fitness through regulatory and/or functional specialization.

For the RLK gene family, microarray data demonstrated pollen-specific expression not only of individual genes, but also for several RLK subfamilies. In total, 21 of 23 pollen-expressed RLK genes were only detected in pollen, suggesting multiple and specialized roles in pollen functions. Our description of 23 pollen-expressed receptor-like protein kinases extends data recently published on this gene family (Kim et al., 2002) and confirms the reported pollen-specific expression of AtPRKc (At07040) and AtPRKd (At5g35390). For two other pollen RLKs named in this article, AtPRKg (At2g26730) was not detected and AtPRKh (At1g48480) was detected only in sporophytic tissues. The significance of RLKs for pollen function was recently highlighted in a study showing interaction of the essential tomato pollen-specific LAT52 protein with the extracellular domain of the pollen-specific receptor kinase LePRK2 (Tang et al., 2002).

The majority of pollen-expressed genes showing gametophytic-sporophytic overlap belonged to the two lowest mRNA abundance classes. In contrast, 55% of genes forming the highly abundant pollen-expressed class were specifically expressed in pollen. This division of abundantly expressed pollen-specific transcripts and poorly expressed transcripts showing gametophytic-sporophytic overlap may reflect different evolutionary constraints in gametophyte and sporophyte. We speculate that the evolution of highly expressed pollen-specific genes depends on initial gene redundancy caused by gene duplication, which allows rapid gametophytic selection accompanied by mutations leading to the loss of the sporophytic component of expression.

The singularity of the pollen transcriptome became evident after allocation of pollen-expressed genes into function groups in comparison with their sporophytic equivalents. Pollen-tube walls are the product of gametophytic expression and possess a unique composition distinct from that of most sporophytic cell types (Li et al., 1999; Hepler et al., 2001). Moreover, the specialized tip growth mechanism of the pollen tube is associated with extensive cell wall synthesis, dynamic changes in cell wall structure, and continuous cell-cell interactions in the pistil (Hepler et al., 2001). Our results revealed a significant investment in the gametophytic expression of genes involved in cell wall metabolism, cytoskeleton, and signaling pathways, consistent with such functional requirements. With regard to the gametophytic synthesis of cell wall polysaccharides, there is limited information concerning the pollen phenotypes of functional knockouts of polysaccharide synthetic enzymes. A cellulose-synthase-like glycosyltransferase, CslA7, has recently been shown to be required for efficient pollen tube growth in Arabidopsis and may synthesize a β-linked polymer (Goubet et al., 2003). The corresponding gene, At2g35650, was expressed in pollen and sporophytic tissues according to RTPCR (Goubet et al., 2003), which was confirmed in our microarray analysis. Pollen microarray data will help to target further knockouts to discover the significance of pollen-expressed polysaccharide synthetic enzymes.

Genes involved in transcription and protein synthesis were generally expressed at lower levels in pollen than in the sporophyte. These data are consistent with the fact that pollen germination and early tube growth in many species are largely independent of transcription, but vitally dependent on translation (Hoekstra and Bruinsma, 1979). mRNA accumulates during pollen maturation and is stored for use in protein synthesis during pollen germination (Schrauwen et al., 1990; Honys et al., 2000). Moreover, in physiologically advanced species, such as T. paludosa, rRNA gene transcription is repressed during pollen maturation and preformed rRNA is stored for use in protein synthesis during pollen germination (see Mascarenhas, 1989). Our data support the concept that Arabidopsis pollen is charged with stored mRNA and preformed translational apparatus enabling rapid activation upon hydration and germination.

Genome-wide analysis of the pollen transcriptome unequivocally demonstrates the unique state of differentiation that distinguishes the mature male gametophyte from the sporophyte. This study revealed the identities and relative expression levels of genes comprising approximately 30% of the pollen transcriptome, including genes expressed specifically in pollen and those showing gametophytic-sporophytic overlap. The functional specialization of the mature male gametophyte for recognition of target tissues and rapid directional growth was highlighted by the overrepresentation of genes expected to be vital for fulfilling these tasks. These included genes involved in cell wall metabolism, the cytoskeleton, and cell signaling. Our data also highlight the diminished role for transcription and the important role of mRNA storage in pollen function. The main impact of this work is that it provides the first genomic and complex view of the male gametophyte. This knowledge now provides a significant opportunity to construct new genome-targeted questions concerning plant cellular functions and the regulation of male gametophyte development and evolution.

For the isolation of mature pollen, inflorescences from over 500 plants were harvested in a large Erlenmeyer flask, 300 mL of ice-cold 0.3 m mannitol was added, and the flask was vigorously shaken for 1 min. The pollen suspension was sequentially filtered through 100- and 53-μm nylon mesh. Pollen grains were concentrated by repeated centrifugation steps (50-mL Falcon tubes, 450 g, 5 min, 4°C), and the final compact pollen pellet was stored at –80°C. The purity of isolated pollen was determined by light microscopy and 4′,6-diamino-phenylindole-staining according to Park et al. (1998). Arabidopsis ecotype Landsberg erecta plants used in sporophytic microarray analyses at NASC were grown in controlled-environment cabinets at 21°C to 22°C under constant light.

Biotinylated target RNA was prepared from 20 μg of total RNA as described in the Affymetrix GeneChip Technical Analysis Manual (Affymetrix, Santa Clara, CA). Double-stranded cDNA was synthesized using SuperScript Choice System (Invitrogen, Carlsbad, CA) with oligo(dT)₂₄ primer fused to T7 RNA polymerase promoter. Biotin-labeled target complementary RNA (cRNA) was prepared by cDNA in vitro transcription using the BioArray High-Yield RNA Transcript Labeling kit (Enzo Biochem, New York) in the presence of biotinylated UTP and CTP.

The Arabidopsis 8K GeneChip array was hybridized with 15 μg of labeled target cRNA for 16 h at 45°C as described in the Affymetrix Technical Analysis Manual. GeneChips were stained with Streptavidin-Phycoerythrin solution and scanned with an Agilent 2500A GeneArray Scanner (Agilent Technologies).

For PCR amplification, 1 μL of 50× diluted RT mix was used. The PCR reaction was carried out in 25 μL with 0.5 unit of BioTaq DNA polymerase (Bioline, London), 1.2 mm MgCl₂, and 20 pmol of each primer. The PCR program was as follows: 2 min at 95°C, 33 cycles of 15 s at 94°C, 15 s at the optimal annealing temperature (63°C to 67°C), and 30 s at 72°C, followed by 10 min at 72°C. As a reverse primer, NESTED primer (5′-AAGCAGTGGTAACAACGCAGAGT-3′) overlapping the 3′-RACE primer was used to eliminate genomic DNA amplification. The following gene-specific forward primers were designed using Primer3 software (http://www-genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi): At3g57690, 5′-TCCCTTGTCTCCCTCTTCAGCTACT-3′; At5g14380, 5′-CTAAGTTCTCTGTTGTCGGCACAGTC-3′; At1g02790, 5′-AACTTAGTTTCCCTATGTGTCCCAAA-3′; At4g27110, 5′-TGCAGTCGGTTATCTCGTTTAAGAAG-3′; At4g19440, 5′-TTGATGGATATGGTAAATTGGGTCAG-3′; At3g20770, 5′-ACATCGAGACAGTCGCTTACCGTAT-3′; and Rubisco, 5′-ATCTTACCTCCCTGACCTTACTGACGT-3′.

Microsoft Excel (Microsoft, Redmont, WA) was used to manage and filter the microarray data. For annotation of genes present on the 8K GeneChip, the Arabidopsis Genome Annotation Release 3.0 published by Institute for Genomic Research (2002) was used. Genes were sorted into functional categories created according to data mined from MAtDB (http://mips.gsf.de/proj/thal/db/index.html), KEGG (http://www.genome.ad.jp/kegg/), and TAIR (http://www.Arabidopsis.org/) databases. Hierarchical clustering of gene families was performed using expression profile data clustering and analysis software EPCLUST (http://ep.ebi.ac.uk/EP/EPCLUS), with correlation measure-based distance and average linkage clustering methods. SOTA analysis was performed using SOTA software with default parameters (Herrero et al., 2001; http://ep.ebi.ac.uk/EP/SOTA/).

We gratefully acknowledge the GARNet transcriptomic center at NASC for performing pollen microarray hybridizations and for providing public baseline 8K GeneChip data for sporophytic samples.