Transcriptome analysis of haploid male gametophyte development in Arabidopsis
© Honys and Twell; licensee BioMed Central Ltd. 2004
Received: 10 August 2004
Accepted: 28 September 2004
Published: 27 October 2004
The haploid male gametophyte generation of flowering plants consists of two- or three-celled pollen grains. This functional specialization is thought to be a key factor in the evolutionary success of flowering plants. Moreover, pollen ontogeny is also an attractive model in which to dissect cellular networks that control cell growth, asymmetric cell division and cellular differentiation. Our objective, and an essential step towards the detailed understanding of these processes, was to comprehensively define the male haploid transcriptome throughout development.
We have developed staged spore isolation procedures for Arabidopsis and used Affymetrix ATH1 genome arrays to identify a total of 13,977 male gametophyte-expressed mRNAs, 9.7% of which were male-gametophyte-specific. The transition from bicellular to tricellular pollen was accompanied by a decline in the number of diverse mRNA species and an increase in the proportion of male gametophyte-specific transcripts. Expression profiles of regulatory proteins and distinct clusters of coexpressed genes were identified that could correspond to components of gametophytic regulatory networks. Moreover, integration of transcriptome and experimental data revealed the early synthesis of translation factors and their requirement to support pollen tube growth.
The progression from proliferating microspores to terminally differentiated pollen is characterized by large-scale repression of early program genes and the activation of a unique late gene-expression program in maturing pollen. These data provide a quantum increase in knowledge concerning gametophytic transcription and lay the foundations for new genomic-led studies of the regulatory networks and cellular functions that operate to specify male gametophyte development.
Development of eukaryotic cells towards particular cell fates is regulated by complex and dynamic changes in gene expression. These changes, when monitored on a genome-wide scale, provide a detailed framework for the analysis and modeling of cellular development. To monitor patterns of gene expression it is important to be able to isolate cells at precise stages along a developmental pathway. Well-developed procedures for cell culture and single-cell PCR techniques have allowed genome-wide changes in gene expression to be monitored during animal cell differentiation [1–4].
Transcriptomic studies of single cell types in plants have focused on diploid sporophytic cell types, including undifferentiated cell suspensions [5, 6], leaf epidermal and mesophyll cells , stomatal guard cells  and cultured mesophyll cells [9, 10]. These studies have provided valuable information about gene expression in single cell types; however, their coverage of the transcriptome has been limited and/or hampered by low RNA yields from individual cells, requiring mRNA preamplification steps that can bias the complementary RNA (cRNA) [11, 12]. Moreover, such studies have not involved the use of the most comprehensive tools for monitoring gene expression that are now available for Arabidopsis - which include the Affymetrix ATH1 gene arrays. Recently, a significant advance in transcriptome analysis of plant cell types has been achieved through fluorescence-activated cell sorting of cell-type marked and protoplasted root cells using Affymetrix ATH1 micorarrays . This has provided a near-comprehensive transcriptomic view of cell-fate determination at three developmental stages in five different domains of the root apex.
In contrast to such enabling technologies and procedures developed for sporophytic cell types there have been no studies that provide a genome-wide perspective of cell fate determination and differentiation during haploid gametophyte development. The haploid male gametophyte generation of flowering plants has a simple and well-defined pathway of development and consists of two- or three-celled pollen grains that deliver two sperm cells via the pollen tube to the embryo sac at fertilization. The highly reduced cell lineage and functional specialization of the male gametophyte are thought to be key factors in the reproductive fitness and evolutionary success of flowering plants. Moreover, pollen ontogeny provides an attractive model of cellular development in which to dissect the regulation of cell growth and division, cellular differentiation and intercellular communication (for reviews see [14–17]).
Recent progress in understanding of molecular and cellular aspects of pollen development has emerged from genetic studies that have identified mutants in Arabidopsis that affect all phases of pollen development [18–29]. In parallel, cDNA libraries and databanks have been obtained for sperm cells in maize, lily, tobacco and Plumbago zelanica [30–33]. Despite such advances there is limited information about developmental changes in gene expression associated with particular phases of male gametophyte development. Our objective was to develop procedures to enable the isolation of populations of microspores and developing pollen grains at precise developmental stages in Arabidopsis and to analyze changes in gene expression from unicellular microspores to mature differentiated pollen grains. A particular advantage of the male gametophyte generation is that developing microspores and pollen grains are symplastically isolated. This facilitates access to viable cell populations at different stages of haploid development without contaminating sporophytic cells.
Some initial progress has been made towards the definition of the male gametophytic transcriptome of Arabidopsis using serial analysis of gene expression (SAGE)  and Affymetrix AG microarrays that harbor probes for approximately 8,000 different genes [35, 36]. These studies have provided valuable insight into the complexity of gene expression in mature pollen and the extent of overlap between male gametophytic and sporophytic gene expression (reviewed in ). However, these studies monitored the expression of only 30% of the annotated genes in Arabidopsis and analyzed mRNA populations only in mature differentiated pollen grains. These studies do not, therefore, provide a developmental perspective of gene expression during development and differentiation of the male gametophyte.
Here we describe spore isolation procedures for Arabidopsis and the use of Affymetrix ATH1 Genome Arrays to analyze transcript expression profiles throughout four successive stages of male gametophyte development in Arabidopsis. Isolated spore populations were large enough to enable RNA extraction for direct microarray hybridization without any preceding amplification step that could lead to bias in expression signals between stages or between genes within individual stages. Progression from proliferating microspores to terminally differentiated pollen was characterized by large-scale repression and the activation of a unique collection of late-program genes during pollen maturation. Putative male gametophyte-specific genes and distinct clusters of coexpressed genes are identified, including key groups of regulatory factors including cell cycle, transcription and translation factors. Bioinformatic and experimental data are used to address the importance of transcription and translation during pollen germination and tube growth
Isolation and characterization of developing spores
Transcriptome profiling throughout microgametogenesis in Arabidopsis required the introduction of a procedure for the isolation of homogeneous populations of viable spores at precisely defined stages of development. The method was based on centrifugation of isolated mixed spores in a Percoll step-gradient [38, 39]. Large homogeneous spore populations at three developmental stages were collected: uninucleate microspores (UNM), bicellular pollen (BCP) and immature tricellular pollen (TCP). In addition, a homogeneous mature pollen grain (MPG) population was isolated from open flowers according to Honys and Twell .
Developmental changes in the male gametophytic transcriptome
Arabidopsis ATH1 Genome Arrays were used to explore the dynamics of gene expression throughout male gametophyte development in comparison with sporophytic tissues. Microarrays were hybridized with cRNA probes made from total RNA purified from isolated spores. Hybridization data from two biological replicates derived from independently grown populations of plants were compared. Only genes with a positive hybridization signal and a detection call value of 1 in both experiments were scored as expressed. Microarray data from each pair of replicates were highly correlated, with correlation coefficients of 0.986 (UNM), 0.972 (BCP), 0.991 (TCP) and 0.971 (MPG). Complete microarray data are publicly available at the European Arabidopsis Stock Centre (NASC) microarray database . Sporophytic ATH1 Genome Array datasets were downloaded from the NASC website . This provided transcriptome data for seedlings at open cotyledon stage (COT, stage 0.7 ), leaves (LEF, stage 6.0), petiole (PET, stage 3.9), stems (STM, stage 6.1), roots (ROT), root hair zone (RHR, stage 1.02), and suspension cell cultures (SUS). Genes that were consistently expressed in replicate sporophytic datasets were identified using the same algorithm used for gametophytic data.
We have previously confirmed and validated the expression pattern of 15 putative pollen-specific genes identified using Affymetrix AG arrays by reverse transcription-PCR analysis . Similarly we validated the current ATH1 datasets by RT-PCR analysis in two separate experiments that included analysis of 41 genes encoding predicted glycosylphosphotidylinositol-anchored proteins (GAPs)  and 16 cation/proton exchanger proteins . In both experiments the expression patterns of all genes tested that were identified as pollen-expressed, or pollen-specific by ATH1 analysis were confirmed by RT-PCR.
The ATH1 Genome Array harbors oligonucleotide probes representing 22,591 genes based on the Arabidopsis Genome Initiative annotation. This represents 80.7% of the most recent estimate of 28,000 protein-coding genes in Arabidopsis . Of these, 13,977 genes gave a consistently positive expression signal in at least one stage of male gametophyte development, representing 61.9% of the unigene targets on the microarray. The majority of these were expressed in the two earliest developmental stages; 11,565 in microspores and 11,909 in bicellular pollen (Figure 1f). After pollen mitosis II, there was a sharp decline in the number of diverse transcripts to 8,788 in tricellular pollen and 7,235 in mature pollen.
To identify genes expressed preferentially or specifically in developing male gametophytes, hybridization data was compared with sporophytic ATH1 datasets (COT, LEF, PET, STM, ROT and RHR; see Additional data file 1). Transcripts with a consistent positive expression signal in at least one stage of male gametophyte development and a zero signal in any sporophytic dataset were considered male gametophyte-specific. In total, 1,355 specific transcripts were identified, representing 9.7% of the male gametophytic transcriptome. The number of male gametophyte-specific transcripts ranged from 857 (BCP) to 625 (MPG). Thus, in contrast to the decline in the total number of diverse transcripts expressed, the representation of male gametophyte-specific transcripts increased, from 6.9% and 7.2% at UNM and BCP-stages to 8.0% and 8.6% at TCP and MPG-stages respectively.
Analysis of the distribution of transcripts among three 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) (Figure 1f), revealed a decrease in the proportion of transcripts forming the high-abundance class during development from 20% to 12%. On the contrary, there was sharp increase in the proportion of mRNAs forming the low-abundance class after pollen mitosis II from 4% (UNM) to 14% (MPG). Moreover, 55% of low-abundance transcripts at MPG stage represented repressed mRNAs expressed more abundantly at earlier stages. Thus, the dramatic decrease in the number of transcripts expressed between bicellular and tricellular stages is paralleled by redistribution of mRNA from the high to the low abundance classes. These changes may be associated with reduced cellular activities and cell differentiation processes together with preferential expression of certain classes of genes during pollen maturation. This finding is in accord with the over-representation of cytoskeleton, cell-wall and signaling-related genes that comprise 26% of the high-abundance transcripts at MPG stage. In particular, the average expression signals of cytoskeleton, cell-wall and signaling-related transcripts were increased by 3.1, 3.7 and 2.3-fold, respectively, compared with the UNM stage.
The relationship between cell proliferation activities and transcriptome profiles was examined by comparison of early UNM and late MPG stages with a publicly available suspension cell culture dataset. These comparisons demonstrated that the microspore transcriptome was significantly more similar to that of cell suspensions (R = 0.474) than to mature pollen (R = 0.194). This is also in accord with the lack of correlation between transcriptome profiles of mature pollen and cell suspensions (R = 0.13).
Co-regulated clusters of gametophytic genes
The majority of male gametophyte-expressed genes (52%) were grouped into four clusters (25, 27, 29 and 35) comprising early expressed genes repressed after PMII. Several large gene clusters collectively containing 1,899 genes (13.6%) were associated with pollen maturation. These were activated or upregulated between BCP and TCP stages, forming clusters 5, 7, 11, 13, 18-24, 26, 28, 38 and 39. In contrast, a discrete set of 298 genes forming cluster 17 was upregulated only after TCP stage. In total, 3,342 late genes (24%), forming clusters 1-3, 6, 8 and in particular, cluster 17, encode proteins that are likely to function during post-pollination development.
Expression of regulatory genes throughout male gametophyte development
Core cell-cycle genes
Among 61 core cell-cycle genes, 55 genes were present on the ATH1 GeneChip and 45 (82%) of these were expressed in the male gametophyte (see Additional data file 1). Representative(s) of all families and subfamilies were expressed. The majority of gametophytic core cell-cycle genes showed similar expression profiles (Figure 4a), with a decline in mRNA abundance after UNM stage to zero (or low levels) at TCP and MPG stages. This pattern is consistent with the termination of proliferation of the microspore and generative cell before pollen maturation.
Putative transcription factors
We identified 1,594 genes encoding putative transcription factors that were divided into 34 gene families (see Additional data file 1). Their representation on the ATH1 GeneChip was 1,350 (85%). Of these, 608 (45%) were expressed in the male gametophyte, including 54 (15.7%) that were male gametophyte-specific. There were distinct differences in the representation of large transcription factor families (with over 25 members) in the gametophyte. Among those over-represented were the p-coumarate 3-hydroxylase (C3H) family (67% of family members present on the ATH1 GeneChip), the CCAAT family (64%), C2H2 zinc finger proteins (57%), the WRKY family (53%), the bZIP family (51%), the TCP family (50%) and the GRAS family (50%). In contrast, the AUX/IAA (20%), HSF (33%), bHLH (34%), NAC (34%), AP2-EREBP (35%), HB (36%), R2R3-MYB (37%), MADS (37%) and C2C2 zinc finger (37%) gene families were all under-represented.
The dominant expression pattern of transcription factor genes reflected the general repression of mRNA diversity between BCP and TCP stages (Figure 4b). Besides a limited number of constitutively expressed genes, two major transcription factor gene groups could be distinguished. One contained a major group of early-expressed genes and the second a smaller group of genes that were more abundantly expressed late during pollen maturation. The same general tendency was apparent when the profiles of individual transcription factor families were analyzed (exemplified by the C3H family, Figure 4d). Several gene families comprised predominantly early-expressed genes. These were the NAC, WRKY, TCP, ARF, Aux/IAA, HMG-box and Alfin-like gene families (Figure 4c-e, Additional data file 3). Complete lists of transcription factor gene families and their expression profiles are presented in Additional data files 1 and 3.
Core translation factors
Among 100 annotated core translation factor genes, 82 were present on the ATH1 GeneChip and 75 (91%) of these were expressed in the male gametophyte (see Additional data file 1). The vast majority of translation factor genes belonged to the early group and these were strongly expressed (Figure 4g). Reflecting the constitutive requirement for protein synthesis, only six genes showed male gametophyte-specific expression. These were: AtPAB3 (At1g22760), AtPAB6 (At3g16380), AtPAB7 (At2g36660), AteIF2-B3 (At3g07920), AteIF4G-like (At4g30680) and AteIF6-2 (At2g39820). There was a striking over-representation of poly(A)-binding (PAB) proteins among the male gametophyte-specific genes; seven out of eight PAB genes were male gametophyte-expressed, three of which were specific. Moreover, two of these gametophyte-specific PAB genes were among the few late pollen genes encoding translation initiation factors (Figure 4h).
Integrating transcriptomic and experimental data
To identify patterns of gene expression involved in Arabidopsis male gametophyte development, we compared the transcriptomes of isolated spores at four discrete developmental stages using ATH1 microarrays. ATH1 microarrays harbor probe sets for 22,591 annotated genes . Of these, 61.9% (13,977 genes) gave positive hybridization signals in at least one stage of male gametophyte development. A comparable proportion of active genes was reported for isolated root cells which expressed 10,492 genes (46%) on ATH1 microarrays . Moreover, in similar studies of animal cell development, 53% of 13,179 arrayed genes were found to be expressed during early murine adipocyte differentiation .
As the proportion of known genes embedded on the ATH1 array is 80.7%, we estimate the total number of genes expressed throughout Arabidopsis male gametophyte development to be more than 17,000. Similarly, the total number of genes expressed at individual developmental stages is estimated to be 14,300 at UNM stage, 14,800 at BCP stage, 10,900 at TCP stage and 9,000 at MPG stage. Previous gene-by-gene approaches identified only 21 different genes expressed during Arabidopsis male gametophyte development (for a review see ). Moreover, only three of these genes were shown to be expressed at microspore stage [53–55]. The data sets reported here include more than 11,000 microspore-expressed genes, representing a 3,600-fold increase in knowledge of gene expression in haploid microspores.
Two recent studies of the Arabidopsis mature pollen transcriptome using Affymetrix 8K AG arrays led to the identification of 992 and 1,584 pollen-expressed mRNAs, respectively [35, 36]. Results obtained with ATH1 and AG arrays are considered comparable and largely independent of the different probe sets used . However, there was a significant discrepancy in the number of incorrectly annotated genes between both arrays, with 6.3% of probe sets on the AG array being incorrectly annotated in comparison with only 0.4% on the ATH1 array . Therefore, results from ATH1 arrays are more accurate as well as more comprehensive. Accordingly, the use of the more complete ATH1 array and more accurate microarray normalization protocols led to an increase in the estimated total number of genes expressed in mature pollen from around 3,500  to around 9,000 (this study). The proportion of these genes that are considered male-gametophyte specific is strongly dependent on the choice of the set of reference sporophytic datasets. In the work reported here, the availability of more comprehensive sporophytic datasets and the application of more stringent criteria therefore led to a decrease in the estimated number of putative pollen-specific genes from around 1,400  to around 800 (this study). This number could be reduced further if cell-type-specific expression within an organ limits detection of overlap with pollen expression. Our data highlight the extensive overlap between sporophytic and gametophyte gene expression and reveal the subset of the transcriptome that is strongly enhanced or specifically expressed during male gametophyte development. Considering all stages of microsporogenesis the total number of putative male-gametophyte-specific genes was 1,355 with the proportion of specific genes increasing from 6.9% at UNM-stage to 8.6% at MPG-stage. Among the male-gametophyte-specific genes identified there was an increase in the collective proportion of cell-wall, cytoskeleton, signaling and transport-related genes from 22% at UNM stage to 34% in MPG stage. This reflects the increasing functional specialization of mature pollen in preparation for a dramatic change in the pattern of cell growth during pollen germination and pollen tube growth.
Developmental analysis of transcriptome data revealed two striking features, a sharp reduction in transcript diversity after BCP stage and a major shift in mRNA populations between BCP and TCP stages. The decline in mRNA diversity after BCP stage is associated with terminal differentiation as well as the documented phenomenon of protein storage in pollen (see , and this study). Moreover, this large-scale repression associated with termination of cell proliferation after PMII is accompanied by the selective activation of new groups of genes that are likely to function during pollen maturation and post-pollination development.
It is interesting that the expression profiles of UNM stage and BCP stages are similar despite the presence of two different cell types in pollen grains at BCP stage - the larger vegetative cell and the smaller generative cell. Given the limited volume of cytoplasm associated with the generative cell, developmental changes in gene expression in the gametic or male germline cells are likely to be masked by the predominant contribution of the vegetative cell cytoplasm. Therefore, our male gametophytic gene expression profiles largely reflect the passage of the microspore through cell division and changes in gene expression associated with the differentiation of the vegetative cell. Large-scale changes in gene expression occur between BCP and TCP stages, and therefore do not coincide with asymmetric division of the microspore. UNM expression patterns persist into the bicellular stage, which is consistent with experiments that demonstrate that vegetative cell fate is specified independently of cell division at pollen mitosis I .
In contrast, generative cell fate appears to be dependent on asymmetric division at pollen mitosis I [25, 58]. Sperm-cell cDNAs and databanks recently established in maize, lily, tobacco and Plumbago zelanica [30–33] provide valuable gametic gene-expression data in other species. Although our data do not provide direct information about gametic gene expression in Arabidopsis, further development of cell gamete isolation sorting  would allow genome-wide identification of generative- and sperm-cell-specific genes in comparison with the datasets generated here.
Hierarchical cluster analysis provided detailed evidence for the dramatic switch between early and late developmental programs. We identified 39 gene clusters that could correspond to co-regulated genes. These included early clusters, several clusters of late genes, those with constitutive expression profiles and clusters showing transient expression with peaks at BCP or TCP stages. The large size of cluster 29 (4,464 genes) documents the homogeneity in expression profiles of most early genes. In contrast, late gene clusters included a significant number of genes with similar profiles between BCP and TCP stages, followed by expression profiles that deviated between TCP and MPG stages. Cluster 1, and in particular cluster 17, contained genes strongly upregulated in TCP and MPG, with likely functions in post-pollination events. The differential fate of certain late gene clusters is likely to be a feature of their requirement during pollen maturation or post-pollination events.
Our analysis revealed completely different expression profiles of transcription factors when compared to core translation factors. The majority of core translation factors belonged to the early-group genes with few that were male gametophyte-specific. This may be expected, given that many genes are involved in general cellular activities. However, genes encoding PAB proteins did not follow the general trend. Seven out of eight Arabidopsis PAB mRNAs were gametophytically expressed. Three PAB genes (PAB3, PAB6 and PAB7) appeared to be male gametophyte-specific and PAB5 was preferentially expressed in pollen. Moreover, PAB3 and PAB5 are the most abundant early and constitutive PAB mRNAs and PAB6 and PAB7 belong among the few late core translation-factor genes. Although these data suggest-specific expression, our data do not rule out expression in other sporophytic tissues, particularly in flowers. Indeed, previously published expression data confirmed the expression of these PABs in other reproductive tissues together with pollen .
Conversely, transcription factors showed more diverse spectra of expression profiles including early, constitutive and late. There was a considerable variation in the expression profiles of individual transcription factor families. The most over-represented was the C3H family, members of which are known to have roles in lignin and other phenylpropanoid pathways in plants . Although sporopollenin synthesis is believed to be under strict sporophytic control (see ), the diversity of gametophytic C3H transcription factors might suggest a function for these genes in regulating chemical interactions between phenylpropanoid precursors secreted by the tapetum. One candidate is the At1g74990 gene encoding a putative RING finger protein, which is abundantly and preferentially expressed at UNM and BCP stages.
The majority of core translation factors belonged to the early gene clusters. In contrast, a significant number of transcription-factor genes were strongly expressed during pollen maturation. These data alone did not obviously support the fact that pollen germination and early tube growth in many species are largely independent of transcription, but vitally dependent on translation . Similarly, we found that Arabidopsis pollen germination and tube growth were relatively independent of transcription, and that active pollen-tube growth, and to a lesser extent pollen germination, were dependent upon protein synthesis. It is known for some plant species that mRNAs and rRNAs accumulate during pollen maturation and are stored for use during pollen germination [62, 63]. Our results show that Arabidopsis pollen is charged with a diverse complement of stored mRNAs that could be used to support pollen germination and pollen tube growth. Moreover, the early synthesis of mRNAs encoding translation factors strongly suggests that these are preformed and stored in mature pollen grains to support rapid activation upon hydration and germination. We also suggest that some abundant late transcription factors could regulate maturation-associated genes or act as repressors of inappropriate transcription in growing pollen tubes.
The key impact of this work is that it provides a genome-wide view of the complexity of gene expression during single cell development in plants. Analysis of the male gametophytic transcriptome provides comprehensive and unequivocal evidence for the unique state of differentiation that distinguishes the developing male gametophyte from the sporophyte. Male gametogenesis is accompanied by large-scale repression of gene expression that is associated with the termination of cell proliferation and the selective activation of new groups of genes involved in maturation and post-pollination events. Development is associated with major early and late transcriptional programs and the expression of about 600 putative transcription factors that are potential regulators of these developmental programs. This wealth of information lays the foundation for new genomic-led studies of cellular functions and the identification of regulatory networks that operate to specify male gametophyte development and functions.
Materials and methods
Plant material and spore isolation
Arabidopsis thaliana ecotype Landsberg erecta plants were grown in controlled-environment cabinets at 21°C under illumination of 150 μmol/m2/sec with a 16-h photoperiod. Isolated spores from three stages of immature male gametophyte were obtained by modification of the protocol of Kyo and Harada [38, 39]. After removal of open flowers, inflorescences from 400 plants were collected and gently ground using a mortar and pestle in 0.3 M mannitol. The slurry was filtered through 100 μM and 53 μM nylon mesh. Mixed spores were concentrated by centrifugation (50 ml Falcon tubes, 450 g, 3 min, 4°C). Concentrated spores were loaded onto the top of 25%/45%/80% Percoll step gradient in a 10-ml centrifuge tube and centrifuged (450 g, 5 min, 4°C). Three fractions were obtained containing: (1) microspores mixed with tetrads; (2) microspores mixed with bicellular pollen; and (3) tricellular pollen (Figure 1). Fraction 2 was diluted with one volume of 0.3 M mannitol loaded onto the top of a 25%/30%/45% Percoll step gradient and centrifuged again under the same conditions. Three subfractions of immature pollen were obtained: (2.1) microspores; (2.2) microspores and bicellular pollen mixture; and (2.3) bicellular pollen. Spores in each fraction were concentrated by centrifugation (eppendorf tubes, 2,000 g, 1 min, 4°C) and stored at -80°C. The purity of isolated fractions was determined by light microscopy and 4',6-diaminophenylindole (DAPI) staining according to . Viability was assessed by fluorescein 3',6'-diacetate (FDA) treatment . Mature pollen was isolated as described previously . Pollen tubes were cultivated in vitro for 4 h according to . Pollen was scored as germinated when pollen tubes were longer than half a pollen grain diameter. Pollen-tube growth was scored by counting those with tubes longer than two pollen grain diameters.
RNA extraction, probe preparation and DNA chip hybridization
Total RNA was extracted from 50 mg of isolated spores at each developmental stage using the RNeasy Plant Kit (Qiagen) according to the manufacturer's instructions. The yield and RNA purity was determined spectrophotometrically and using an Agilent 2100 Bioanalyzer at the NASC.
Biotinylated target RNA was prepared from 20 μg of total RNA as described in the Affymetrix GeneChip expression analysis technical manual. Double-stranded cDNA was synthesized using SuperScript Choice System (Life Technologies) with oligo(dT)24 primer fused to T7 RNA polymerase promoter. Biotin-labeled target cRNA was prepared by cDNA in vitro transcription using the BioArray High-Yield RNA Transcript Labeling Kit (Enzo Biochem) in the presence of biotinylated UTP and CTP.
Arabidopsis ATH1 Genome Arrays were hybridized with 15 μg labeled target cRNA for 16 h at 45°C. Microarrays were stained with streptavidin-phycoerythrin solution and scanned with an Agilent 2500A GeneArray Scanner.
Sporophytic data from public baseline GeneChip experiments used for comparison with the pollen transcriptome were downloaded from the NASC website [41, 64]. The list of dataset codes was as follows: COT (three replicates), Cornah_A4-cornah-wsx_SLD_REP1-3; LEF (three replicates), A4-LLOYD-CON_REP1-3; PET (three replicates), Millenaar_A1-MILL-AIR-REP1-3; STM (two replicates), Turner_A-7-Turne-WT-Base1-2_SLD; ROT (two replicates), Sophie_A1-Fille-WT-nodex_SLD, Sophie_A5-Fille-WT-nodex_SLD; RHR (two replicates), Jones_A1-jones-WT1, SLD, Jones_A1-jones-WT2_SLD; SUS (three replicates), A1-WILLA-CON-REP1-3.
All gametophytic and sporophytic datasets were normalized using freely available dChip 1.3 software . The reliability and reproducibility of analyses was ensured by the use of duplicates or triplicates in each experiment, the normalization of all 26 arrays to the median probe intensity level and the use of normalized CEL intensities of all arrays for the calculation of model-based gene-expression values based on the Perfect Match-only model [66, 67]. A given gene was scored as 'expressed' when it gave a reliable expression signal in all replicates. Expression signal value '0' means that the detection call value was not 'present' in all replicates provided. All raw and dChip-normalized gametophytic datasets are available at the Institute of Experimental Botany AS CR website . Although a RT-PCR validation of microarray data was not performed specifically for the purpose of this publication, our confidence in the quality of the data presented is based on our previously published RT-PCR validation of the expression of 70 genes [21, 35, 41].
Microsoft Excel was used to manage and filter the microarray data. For annotation of genes present on the ATH1 Array, the Arabidopsis Genome Annotation Release 3.0 published by The Institute for Genomic Research  was used. Genes were sorted into functional categories created according to data mined from the Munich Information Center for Protein Sequences Arabidopsis thaliana Database , Kyoto Encyclopedia of Genes and Genomes  and TAIR . Hierarchical clustering of expressed genes was performed using expression-profile data clustering and analysis software EPCLUST , with correlation measure based distance and average linkage clustering methods.
Additional data files
The following additional data is available with the online version of this article: Additional data file 1 is an Excel file containing the following items. The table Data contains the complete transcriptomic datasets used. Data were normalized using dChip 1.3 as described in Materials and methods. Expression signal value '0' means that the detection call value for particular gene was not 'present' in all replicates provided. In the column 'Cluster', the appropriate cluster for each male gametophyte-expressed gene is shown. The table Clusters gives the number of genes comprising all 37 clusters of genes coexpressed during male gametophyte development. The table Cell-cycle data lists core cell-cycle genes showing their expression values in male gametophytic datasets. Genes were defined according to . The chart shows expression profiles of male gametophyte-expressed core cell-cycle genes. The table Transcription data lists transcription-factor genes, showing their expression values in male gametophytic datasets. Genes comprising Arabidopsis transcription factor families were derived by compilation of data available at the Ohio State University Arabidopsis Gene Regulatory Information Server , data published in  and database homology searches. MADS-box and bHLH gene families were defined according to  and , respectively. The table Translation data lists core translation-factor genes showing their expression values in male gametophytic datasets. Genes were defined according to the FIAT database . The chart shows expression profiles of male gametophyte-expressed core translation-factor genes. The Transcription table summarizes transcription factor gene families showing the number of genes expressed during male gametophyte development. Additional data file 2 lists a complete set of 39 clusters of genes coexpressed during male gametophyte development. Clusters were determined using EPCLUST software with a threshold value of 0.05. The list of genes comprising each cluster is given in Additional data file 1. Additional data file 3 gives the expression profiles of male gametophyte-expressed transcription factors sorted into individual gene families. Expression data are given in Additional data file 1.
We gratefully acknowledge support from the BBSRC and the GARNet transcriptomic centre at NASC for performing pollen microarray hybridizations. We thank Andy Johnson for help with microspore extraction, John Okyere for advice on microarray normalization protocols and all members of the Twell laboratory for helpful comments on the manuscript. D.H. was supported through a Royal Society/NATO Fellowship, a Ministry of Education of the Czech Republic Grant 1K03018 and a Grant Agency of the ASCR grant KJB6038409.
- Burton GR, Guan Y, Nagarajan R, McGehee RE: Microarray analysis of gene expression during early adipocyte differentiation. Gene. 2002, 293: 21-31. 10.1016/S0378-1119(02)00726-6.PubMedView ArticleGoogle Scholar
- Klebes A, Biehs B, Cifuentes F, Kornberg TB: Expression profiling of Drosophila imaginal discs. Genome Biol. 2002, 3: research0038.1-0038.16. 10.1186/gb-2002-3-8-research0038.View ArticleGoogle Scholar
- Stathopoulos A, Van Drenth M, Erives A, Markstein M, Levine M: Whole-genome analysis of dorsal-ventral patterning in the Drosophila embryo. Cell. 2002, 111: 687-701. 10.1016/S0092-8674(02)01087-5.PubMedView ArticleGoogle Scholar
- Chiang MK, Melton DA: Single-cell transcript analysis of pancreas development. Dev Cell. 2003, 4: 383-393. 10.1016/S1534-5807(03)00035-2.PubMedView ArticleGoogle Scholar
- Breyne P, Dreesen R, Vandepoele K, De Veylder L, Van Breusegem F, Callewaert L, Rombauts S, Raes J, Cannoot B, Engler G: Transcriptome analysis during cell division in plants. Proc Natl Acad Sci USA. 2002, 99: 14825-14830. 10.1073/pnas.222561199.PubMedPubMed CentralView ArticleGoogle Scholar
- Menges M, Hennig L, Gruissem W, Murray JA: Cell cycle-regulated gene expression in Arabidopsis. J Biol Chem. 2002, 277: 41987-42002. 10.1074/jbc.M207570200.PubMedView ArticleGoogle Scholar
- Brandt S, Kloska S, Altmann T, Kehr J: Using array hybridization to monitor gene expression at the single cell level. J Exp Bot. 2002, 53: 2315-2323. 10.1093/jxb/erf093.PubMedView ArticleGoogle Scholar
- Leonhardt N, Kwak JM, Robert N, Waner D, Leonhardt G, Schroeder JI: Microarray expression analyses of Arabidopsis guard cells and isolation of a recessive abscisic acid hypersensitive protein phosphatase 2C mutant. Plant Cell. 2004, 16: 596-615. 10.1105/tpc.019000.PubMedPubMed CentralView ArticleGoogle Scholar
- Demura T, Tashiro G, Horiguchi G, Kishimoto N, Kubo M, Matsuoka N, Minami A, Nagata-Hiwatashi M, Nakamura K, Okamura Y, et al: Visualization by comprehensive microarray analysis of gene expression programs during transdifferentiation of mesophyll cells into xylem cells. Proc Natl Acad Sci USA. 2002, 99: 15794-15799. 10.1073/pnas.232590499.PubMedPubMed CentralView ArticleGoogle Scholar
- Milioni D, Sado PE, Stacey NJ, Roberts K, McCann MC: Early gene expression associated with the commitment and differentiation of a plant tracheary element is revealed by cDNA-amplified fragment length polymorphism analysis. Plant Cell. 2002, 14: 2813-2824. 10.1105/tpc.005231.PubMedPubMed CentralView ArticleGoogle Scholar
- Phillips J, Eberwine JH: Antisense RNA amplification: A linear amplification method for analyzing the mRNA population from single living cells. Methods. 1996, 10: 283-288. 10.1006/meth.1996.0104.PubMedView ArticleGoogle Scholar
- Luo L, Salunga RC, Guo H, Bittner A, Joy KC, Galindo JE, Xiao H, Rogers KE, Wan JS, Jackson MR, et al: Gene expression profiles of laser-captured adjacent neuronal subtypes. Nat Med. 1999, 5: 117-122. 10.1038/4806.PubMedView ArticleGoogle Scholar
- Birnbaum K, Shasha DE, Wang JY, Jung JW, Lambert GM, Galbraith DW, Benfey PN: A gene expression map of the Arabidopsis root. Science. 2003, 302: 1956-1960. 10.1126/science.1090022.PubMedView ArticleGoogle Scholar
- McCormick S: Male gametophyte development. Plant Cell. 1993, 5: 1265-1275. 10.1105/tpc.5.10.1265.PubMedPubMed CentralView ArticleGoogle Scholar
- Twell D, Park SK, Lalanne E: Asymmetric division and cell-fate determination in developing pollen. Trends Plant Sci. 1998, 3: 305-310. 10.1016/S1360-1385(98)01277-1.View ArticleGoogle Scholar
- Twell D: Pollen developmental biology. In Plant Reproduction. Annual Plant Reviews. Edited by: O'Neil SD, Roberts JA. 2002, Sheffield: Sheffield Academic Press, 6: 86-153.Google Scholar
- McCormick S: Control of male gametophyte development. Plant Cell. 2004, 16 (Suppl): S142-S153. 10.1105/tpc.016659.PubMedPubMed CentralView ArticleGoogle Scholar
- Grini PE, Schnittger A, Schwarz H, Zimmermann I, Schwab B, Jurgens G, Hulskamp M: Isolation of ethyl methanesulfonate-induced gametophytic mutants in Arabidopsis thaliana by a segregation distortion assay using the multimarker chromosome 1. Genetics. 1999, 151: 849-863.PubMedPubMed CentralGoogle Scholar
- Howden R, Park SK, Moore JM, Orme J, Grossniklaus U, Twell D: Selection of T-DNA-tagged male and female gametophytic mutants by segregation distortion in Arabidopsis. Genetics. 1998, 149: 621-631.PubMedPubMed CentralGoogle Scholar
- Lalanne E, Twell D: Genetic control of male germ unit organization in Arabidopsis. Plant Physiol. 2002, 129: 865-875. 10.1104/pp.003301.PubMedPubMed CentralView ArticleGoogle Scholar
- Lalanne E, Honys D, Johnson A, Borner GHH, Lilley KS, Dupree P, Grossniklaus U, Twell D: SETH1 and SETH2, two components of the glycosylphosphatidylinositol anchor biosynthetic pathway, are required for pollen germination and tube growth in Arabidopsis. Plant Cell. 2004, 16: 229-240. 10.1105/tpc.014407.PubMedPubMed CentralView ArticleGoogle Scholar
- Lalanne E, Michaelidis C, Moore JM, Gagliano W, Johnson A, Patel R, Howden R, Vielle-Calzada J-P, Grossniklaus U, Twell D: Analysis of transposon insertion mutants highlights the diversity of mechanisms underlying male progamic development in Arabidopsis. Genetics. 2004, 167: 1975-1986. 10.1534/genetics.104.030270.PubMedPubMed CentralView ArticleGoogle Scholar
- Lobstein E, Guyon A, Férault M, Twell D, Pelletier G, Bonhomme S: The putative Arabidopsis homolog of yeast Vps52p is required for pollen tube elongation, localizes to Golgi and might be involved in vesicle trafficking. Plant Physiol. 2004, 135: 1480-1490. 10.1104/pp.103.037747.PubMedPubMed CentralView ArticleGoogle Scholar
- Oh S-A, Park S-K, Jang I, Howden R, Moore JM, Grossniklaus U, Twell D: Halfman, an Arabidopsis male gametophytic mutant associated with a 150 kb chromosomal deletion at the site of transposon insertion. Sex Plant Reprod. 2003, 16: 99-102. 10.1007/s00497-003-0181-9.View ArticleGoogle Scholar
- Park SK, Howden R, Twell D: The Arabidopsis thaliana gametophytic mutation gemini pollen1 disrupts microspore polarity, division asymmetry and pollen cell fate. Development. 1998, 125: 3789-3799.PubMedGoogle Scholar
- Park S-K, Rahman D, Oh S-A, Twell D: Gemini pollen 2, a male and female gametophytic cytokinesis defective mutation. Sex Plant Reprod. 2004, 17: 63-70. 10.1007/s00497-004-0216-x.PubMedPubMed CentralView ArticleGoogle Scholar
- Procissi A, de Laissardiere S, Ferault M, Vezon D, Pelletier G, Bonhomme S: Five gametophytic mutations affecting pollen development and pollen tube growth in Arabidopsis thaliana. Genetics. 2001, 158: 1773-83.PubMedPubMed CentralGoogle Scholar
- Procissi A, Guyon A, Pierson ES, Giritch A, Knuiman B, Grandjean O, Tonelli C, Derksen J, Pelletier G, Bonhomme S: KINKY POLLEN encodes a SABRE-like protein required for tip growth in Arabidopsis and conserved among eukaryotes. Plant J. 2003, 36: 894-904. 10.1046/j.1365-313X.2003.01933.x.PubMedView ArticleGoogle Scholar
- Twell D, Park SK, Hawkins TJ, Schubert D, Schmidt R, Smertenko A, Hussey PJ: MOR1/GEM1 plays an essential role in the plant-specific cytokinetic phragmoplast. Nat Cell Biol. 2002, 4: 711-714. 10.1038/ncb844.PubMedPubMed CentralView ArticleGoogle Scholar
- Engel ML, Chaboud A, Dumas C, McCormick S: Sperm cells of Zea mays have a complex complement of mRNAs. Plant J. 2003, 34: 697-707. 10.1046/j.1365-313X.2003.01761.x.PubMedView ArticleGoogle Scholar
- Xu H, Swoboda I, Bhalla P, Singh MB: Male gametic cell-specific gene expression in flowering plants. Proc Natl Acad Sci USA. 1999, 96: 2554-2558. 10.1073/pnas.96.5.2554.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu H, Weterings K, Vriezen W, Feron R, Xue Y, Derksen J, Mariani C: Isolation and characterization of male-germ-cell transcripts in Nicotiana tabacum. Sex Plant Reprod. 2002, 14: 339-346. 10.1007/s00497-002-0128-6.View ArticleGoogle Scholar
- Zhang Z, Xu H, Singh MB, Russell SD: Isolation and collection of two populations of viable sperm cells from the pollen of Plumbago zeylanica. Zygote. 1998, 6: 295-298. 10.1017/S0967199498000240.PubMedView ArticleGoogle Scholar
- Lee JY, Lee DH: Use of serial analysis of gene expression technology to reveal changes in gene expression in Arabidopsis pollen undergoing cold stress. Plant Physiol. 2003, 132: 517-529. 10.1104/pp.103.020511.PubMedPubMed CentralView ArticleGoogle Scholar
- Honys D, Twell D: Comparative analysis of the Arabidopsis pollen transcriptome. Plant Physiol. 2003, 132: 640-652. 10.1104/pp.103.020925.PubMedPubMed CentralView ArticleGoogle Scholar
- Becker JD, Boavida LC, Carneiro J, Haury M, Feijo JA: Transcriptional profiling of Arabidopsis tissues reveals the unique characteristics of the pollen transcriptome. Plant Physiol. 2003, 133: 713-725. 10.1104/pp.103.028241.PubMedPubMed CentralView ArticleGoogle Scholar
- Da Costa-Nunes JA, Grossniklaus U: Unveiling the gene-expression profile of pollen. Genome Biol. 2003, 5: 205-10.1186/gb-2003-5-1-205.PubMedPubMed CentralView ArticleGoogle Scholar
- Kyo M, Harada H: Studies on conditions for cell division and embryogenesis in isolated pollen culture of Nicotiana rustica. Plant Physiol. 1985, 79: 90-94.PubMedPubMed CentralView ArticleGoogle Scholar
- Kyo M, Harada H: Control of the developmental pathway of tobacco pollen in vitro. Planta. 1986, 168: 427-432.PubMedView ArticleGoogle Scholar
- NASCArrays. [http://affymetrix.arabidopsis.info/narrays/experimentbrowse.pl]
- NASC: The European Arabidopsis Stock Centre. [http://Arabidopsis.info]
- Boyes DC, Zayed AM, Ascenzi R, McCaskill AJ, Hoffman NE, Davis KR, Gorlach J: Growth stage-based phenotypic analysis of Arabidopsis : a model for high throughput functional genomics in plants. Plant Cell. 2001, 13: 1499-1510. 10.1105/tpc.13.7.1499.PubMedPubMed CentralView ArticleGoogle Scholar
- Sze H, Padmanaban S, Cellier F, Honys D, Cheng N-H, Bock KW, Conéjéro G, Li X, Twell D, Ward JM, Hirschi KD: Expression patterns of a novel AtCHX gene family highlight potential roles in osmotic adjustment and K+ homeostasis in pollen development. Plant Physiol. 2004, 136: 2532-2547. 10.1104/pp.104.046003.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamada K, Lim J, Dale JM, Chen H, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M, et al: Empirical analysis of transcriptional activity in the Arabidopsis genome. Science. 2003, 302: 842-846. 10.1126/science.1088305.PubMedView ArticleGoogle Scholar
- Vandepoele K, Raes J, De Veylder L, Rouze P, Rombauts S, Inze D: Genome-wide analysis of core cell cycle genes in Arabidopsis. Plant Cell. 2002, 14: 903-916. 10.1105/tpc.010445.PubMedPubMed CentralView ArticleGoogle Scholar
- TAIR homepage. [http://www.arabidopsis.org]
- Arabidopsis thaliana web page. [http://arabidopsis.med.ohio-state.edu]
- Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.PubMedView ArticleGoogle Scholar
- Parenicova L, de Folter S, Kieffer M, Horner DS, Favalli C, Busscher J, Cook HE, Ingram RM, Kater MM, Davies B, et al: Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis : new openings to the MADS world. Plant Cell. 2003, 15: 1538-1551. 10.1105/tpc.011544.PubMedPubMed CentralView ArticleGoogle Scholar
- Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 15: 1749-1770. 10.1105/tpc.013839.PubMedPubMed CentralView ArticleGoogle Scholar
- FIAT: factors in Arabidopsis translation. [http://www.cm.utexas.edu/browning/db]
- TIGR FTP site. [ftp://ftp.tigr.org/pub/data/a_thaliana/ath1]
- Carpenter JL, Ploense SE, Snustad DP, Silflow CDP: Preferential expression of an alpha-tubulin gene of Arabidopsis in pollen. Plant Cell. 1992, 4: 557-571. 10.1105/tpc.4.5.557.PubMedPubMed CentralGoogle Scholar
- Roberts MR, Foster GD, Blundell RP, Robinson SW, Kumar A, Draper J, Scott RP: Gametophytic and sporophytic expression of an anther-specific Arabidopsis thaliana gene. Plant J. 1993, 3: 111-20. 10.1046/j.1365-313X.1993.t01-5-00999.x.PubMedView ArticleGoogle Scholar
- Truernit E, Stadler R, Baier K, Sauer NP: A male gametophyte-specific monosaccharide transporter in Arabidopsis. Plant J. 1999, 17: 191-201. 10.1046/j.1365-313X.1999.00372.x.PubMedView ArticleGoogle Scholar
- Hennig L, Menges M, Murray JAH, Gruissem W: Arabidopsis transcript profiling on Affymetrix GeneChip arrays. Plant Mol Biol. 2003, 53: 457-465. 10.1023/B:PLAN.0000019069.23317.97.PubMedView ArticleGoogle Scholar
- Twell D: The diversity and regulation of gene expression in the pathway of male gametophyte development. Soc ExpBiol Semin. 1994, Ser 55: 83-135.Google Scholar
- Eady C, Lindsey K, Twell D: The significance of microspore division and division symmetry for vegetative cell-specific transcription and generative cell differentiation. Plant Cell. 1995, 7: 65-74. 10.1105/tpc.7.1.65.PubMedPubMed CentralView ArticleGoogle Scholar
- Belostotsky DA: Unexpected complexity of poly(A)-binding protein gene families in flowering plants: three conserved lineages that are at least 200 million years old and possible auto- and cross-regulation. Genetics. 2003, 163: 311-319.PubMedPubMed CentralGoogle Scholar
- Franke R, Humphreys JM, Hemm MR, Denault JW, Ruegger MO, Cusumano JC, Chapple C: The Arabidopsis REF8 gene encodes the 3-hydroxylase of phenylpropanoid metabolism. Plant J. 2002, 30: 33-45. 10.1046/j.1365-313X.2002.01266.x.PubMedView ArticleGoogle Scholar
- Hoekstra FA, Bruinsma J: Protein synthesis of binucleate and trinucleate pollen and its relationship to tube emergence and growth. Planta. 1979, 146: 559-566.PubMedView ArticleGoogle Scholar
- Schrauwen JAM, De Groot PFM, Van Herpen MMA, Van Den Lee T, Reynen WH, Weterings KAP, Wullems GJ: Stage-related expression of mRNAs during pollen development in lily and tobacco. Planta. 1990, 182: 298-304. 10.1007/BF00197125.PubMedView ArticleGoogle Scholar
- Honys D, Combe JP, Twell D, Čapková V: The translationally repressed pollen-specific ntp303 mRNA is stored in non-polysomal mRNPs during pollen maturation. Sex Plant Reprod. 2000, 13: 135-144. 10.1007/s004970000047.View ArticleGoogle Scholar
- Craigon DJ, James N, Okyere J, Higgins J, Jotham J, May S: NASCArrays: a repository for microarray data generated by NASC's transcriptomics service. Nucleic Acids Res. 2004, 32 Database issue: D575-D577. 10.1093/nar/gkh133.View ArticleGoogle Scholar
- DNA-Chip Analyzer (dChip). [http://www.dchip.org]
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98: 31-36. 10.1073/pnas.011404098.PubMedPubMed CentralView ArticleGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: Model validation, design issues and standard error application. Genome Biol. 2001, 2: research0032.1-0032.11. 10.1186/gb-2001-2-8-research0032.Google Scholar
- LBP IEB AS CR website. [http://www.ueb.cas.cz/laboratory_of_pollen_biology/frames_3.htm]
- MatDB: Arabidopsis thaliana. [http://mips.gsf.de/proj/thal/db/index.html]
- KEGG: Kyoto Encyclopedia of Genes and Genomes. [http://www.genome.ad.jp/kegg]
- EPCLUST - clustering, visualization, and analysis. [http://ep.ebi.ac.uk/EP/EPCLUST]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.