<title>Background</title> <p>Organ formation is not limited to embryonic development, but can occur throughout the lifetime of a plant. The potential to develop new organs post-embryonically is attributed to meristems located at the growing tips of the plants, with the root apical meristem generating the underground part of the plant and the shoot apical meristem (SAM) giving rise to the entire shoot system after seed germination.</p> <p>Like its root counterpart, the SAM contains a pool of pluripotent stem cells that can self-maintain as well as produce the cells that can differentiate into multiple cell and tissue types [reviewed in [##REF##16724254##1##]]. While lateral organs such as the leaves are initiated from the peripheral regions of the SAM, the basal regions of the SAM contribute to the formation of the stem. The stem cells of the SAM must thus replenish areas where cells have been recruited and at the same time maintain the population of stem cells. This is generally attributed to an active process of communication among neighbouring SAM cells in the microenvironment of the stem cells [##REF##16315035##2##, ####REF##16616580##3##, ##REF##17499544##4####17499544##4##].</p> <p>Elegant genetic work carried out in the model plant, <italic>Arabidopsis thaliana</italic>, has enhanced our understanding of this vital developmental process [reviewed by [##REF##18031217##5##]]. This is exemplified by the identification in <italic>Arabidopsis </italic>of WUSCHEL (WUS), a homeodomain transcription factor essential for maintaining the pools of stem cells in an undifferentiated state [##REF##8565856##6##], and the <italic>CLAVATA </italic>group of genes that act together to restrict the proliferation of stem cells [##REF##12927973##7##]. While the <italic>Arabidopsis </italic>genus provides invaluable model plants for enhancing our understanding of plant biology, it does not represent all the diverse developmental, environmental and physiological processes operating in the plant kingdom. There thus remains a need to extend the knowledge gained to other plant species especially crop plants.</p> <p>Applying modern genomics research techniques to improving crops requires new knowledge and the development of new genomics resources. Legume species belonging to the family Fabaceae are cultivated for seeds rich in proteins, and represent important components of the diet in many parts of the world, especially pea, lentil and soybean. Legumes have attracted the attention of biologists because of their unique characteristics such as nitrogen fixation, secondary metabolism and pod development, and these are among the various processes that cannot be studied in <italic>Arabidopsis </italic>species. On the other hand, garden pea (<italic>Pisum sativum</italic>) has been a classic model species used in genetics [##UREF##0##8##,##UREF##1##9##] and plant-development studies. Its extensive use in studies of flowering initiation and development has provided important insights into this transition process [##REF##12114569##10##, ####REF##14563931##11##, ##REF##16096961##12##, ##UREF##2##13####2##13##]. In addition, the availability of various developmental and flowering pea mutants [##REF##14563931##11##,##UREF##3##14##, ####UREF##4##15##, ##REF##10948249##16##, ##REF##15749758##17##, ##UREF##5##18####5##18##] makes this tractable for genomics studies.</p> <p>In this study, we applied a transcriptomics approach to investigate the gene expression profiles associated with the SAM of the garden pea, an agriculturally and commercially important model legume. We also investigated the use of micro-dissected SAMs in unravelling the transcriptome profile of the SAM. By identifying genes that exhibit differential expression between SAMs and non-meristematic tissues (NM), we aimed to elucidate the transcriptional signature of the SAM and thereby identify genes that might play important roles in regulating SAM activity.</p> <p>To this end, three directional cDNA libraries were constructed using SAMs that were carefully micro-dissected from garden peas. These libraries comprised the standard cDNA library plus a normalized library and a subtracted library, in order to increase the likelihood of recovering rare cDNAs, allowing the sampling of the wide diversity of genes expressed in the pea SAM. The EST sequences derived from these libraries were used in the subsequent design of a CombiMatrix CustomArray™ 4 × 2 K oligonucleotide array that was representative of the gene content of the SAM.</p> <p>In this paper, we present the EST and transcriptional profiling data from this genomics project. The transcriptional profiling experiment represents the first analysis of genes that exhibit differential expression between the pea SAM and NM. The data reveal that transcripts putatively annotated as being associated with cell division and proliferation, epigenetic regulation, auxin-mediated responses and microRNA (miRNA) regulation are more abundant in the SAM than in the NM. In contrast, sequences related to photosynthesis, abiotic or biotic stress responses, reactive oxygen species (ROS) homeostasis and general cell-wall maintenance are down-regulated in the SAM.</p>
<title>Methods</title> <title>Plant materials and cDNA libraries synthesis</title> <p>Garden pea (<italic>Pisum sativum</italic>) cultivar Torsdag was grown in a greenhouse located at the University of Melbourne, Australia. SAMs were micro-dissected from 10-day-old peas under the dissecting microscope at 40× magnification. Any leaf primordia were excluded in order to create a meristem-enriched tissue collection and the location of tissue sample is indicated in Figure ##FIG##4##5##. Dissected samples were quickly frozen in liquid nitrogen and stored at -80°C until used for RNA extraction and cDNA library synthesis [##UREF##9##49##]. For the subtracted library, the driver sequences were derived from an equal mix of RNA extracted from non-meristematic (NM) tissue consisted of primary stem (without axillary meristems), mature leaf lamina, primary roots (without root apical meristem and root hairs), whereas the tester sequences consisted of RNA harvested from dissected SAM. All cDNAs were cloned into pBlueScriptIISK+ plasmid vector.</p> <title>EST sequencing</title> <p>A total of about 500, 1000 and 1500 randomly picked clones from each of the standard (C), subtracted (S) and normalized (N) libraries were sequenced at Australian Genome Research Facility (AGRF), Australia and subsequently at Macrogen Korea using T7 primer. These sequences have been deposited in GenBank under the accession numbers <ext-link ext-link-type="gen" xlink:href="EX568682">EX568682</ext-link> to <ext-link ext-link-type="gen" xlink:href="EX571416">EX571416</ext-link>.</p> <title>Sequence analysis and annotation</title> <p>Sequence data were trimmed off vector, adaptor and low quality sequences using SEQTools [##UREF##10##50##]. Trimmed sequences that were shorter than 100 basepairs were excluded from further analysis. Blast score-based clustering method with a score cut-off of 0.6 from SEQTools was then used to assemble the sequences. All clusters and singletons resulting from this automated clustering were considered to be the best estimation of a minimal gene set for our EST library and we have called this set as "unigenes".</p> <p>All sequences were then imported into Blast2GO, a web-based Gene Ontology (GO) annotation and analysis tool [##REF##16081474##21##] for subsequent analysis. This involved automated retrieval of GO terms associated to the hits obtained after a BLASTX search of the corresponding unigene sequence against NCBI (nr) protein database. The e-value cut off was set at 1e-5.</p> <title>RNA Extraction for microarray experiments</title> <p>Total RNA was extracted from dissected SAM (approximately 80 SAMs per extraction) or other plant parts (primary stem, primary roots and mature leaves) using Qiagen RNeasy Mini Kit. Four independent tissue collections and RNA extractions (designated A, B, C and D) were performed for each of the microarray hybridization experiment.</p> <title>Design of pea SAM Combimatrix CustomArray™ 4X2K</title> <p>The Combimatrix arrays are semiconductor-based oligonucleotide microarrays and are generated based on CombiMatrix technology (hyperlink) of <italic>in situ </italic>synthesis [##UREF##11##51##]. The CustomArray™ 4x2k is a microarray that is divided into 4 sectors, each of which can contain up to 2,240 different oligonucleotide probes (spots) and can be hybridized individually with different targets using a provided sectored hybridization cap. A total of 1686 pea sequences (290 sequences from C library, 300 sequences from S library and 1086 ESTs from the N library) together with 500 pea sequences randomly selected from GenBank pea protein database were submitted for probe design using the open source CombiMatrix probe design system. The length of probes ranged from 35–40 bases in length. A variety of control elements were also arrayed on the slide and these include blank spot, housekeeping genes (actin) as well as non-plant transgenes.</p> <title>Target preparation and hybridization to microarray</title> <p>Target preparation and hybridization were performed in Australian Genome Research Facility Ltd (AGRF) according to the standard CombiMatrix protocol described in detail at <ext-link ext-link-type="uri" xlink:href="http://www.combimatrix.com/docs/PTL005_00_4x2K_%20Hyb_%20Imaging.pdf"/>. One microgram of total RNA from SAM or NM was labelled using the Kreatech's ULS™ RNA ampULSe kit to generate Cy5 or Cy3 labelled targets. The Cy5- or Cy3-labelled cDNA was then hybridized to different sector of the chip according to a balanced block design dual label experiment scheme [##REF##12424114##23##]:</p> <p>Sector 1: Cy3-SAM A vs Cy5-NM A</p> <p>Sector 2: Cy5-SAM B vs Cy3-NM B</p> <p>Sector 3: Cy3-SAM C vs Cy5-NM C</p> <p>Sector 4: Cy5-SAM D vs Cy3-NM D</p> <title>Image acquisition and data analysis</title> <p>The Cy5- and Cy3-hybridized chip was then scanned using Genepix 4000B microarray scanner (Axon Instruments, CA, USA) according to manufacturer's instructions. The generated tiff image files were then imported into Combimatrix Microarray Imager to produce intensity data. The LimmaGUI software, which is an implementation of the Empirical Bayes linear modelling approach, was used for subsequent statistical analysis of the resulting data [##REF##15297296##24##]. A robust spline method was chosen for within array normalization and a least-square linear model fit was computed with the p-value adjusted using the Benjamini-Hochberg procedure. The entire differentially-expressed transcripts (p value &lt; 0.05) in the SAM are listed in Additional File ##SUPPL##0##1## and Additional File ##SUPPL##1##2##. Microarray data have been deposited in the Gene Expression Omnibus database [##UREF##6##20##] under accession number <ext-link ext-link-type="gen" xlink:href="GSE9278">GSE9278</ext-link>.</p> <title>RT-PCR analysis</title> <p>The one-tube, two enzyme Access RT-PCR system (Promega, Annandale, New South Wales, Australia) was used according to manufacturer's instructions in all RT-PCR analysis. Ten ng of RNA isolated from the SAM, mature stem, mature leaf and primary root of 10-day-old pea seedlings were used as a template in a 10 μl reaction volume. The pea actin gene was used as an internal control. The number of cycles used for the transcripts investigated was routinely between 25–28 and 80 % of the PCR reaction was separated on 1% agarose gel containing 0.1 μg/μl ethidium bromide and visualized under UV light.</p>
<title>Results and Discussion</title> <title>Features of generated ESTs</title> <p>A total of 3000 clones from three cDNA libraries were single-pass sequenced from their 5' ends. Sequence cleaning processes as outlined in the Materials and Methods yielded 2735 ESTs. These sequences had an average trimmed length of 519 base pairs and were assembled into 348 clusters and 1332 singlets, resulting in the final annotation of 1686 unigenes. Clusters ranged in membership from 67 ESTs (one) to two ESTs (253). The redundancy levels were 15.0%, 20.3% and 62.8% within the normalized, standard and subtracted libraries, respectively. A high redundancy level of cDNA libraries constructed using a similar subtraction method from <italic>Thellungiella </italic>plants under abiotic stress conditions has also been observed previously [##REF##16021339##19##].</p> <p>The translated products of the 1686 unigenes were searched against the non-redundant protein database provided by GenBank [##UREF##6##20##] to putatively assign their functions. At the time of writing, 918 (54.4%) unigenes showed significant similarity (<italic>E </italic>value cut-off of 10<sup>-5</sup>) to genes of known or putative function, whereas 549 (32.6%) ESTs were assigned to transcripts with unknown function, which includes hypothetical genes predicted in genomes of model organisms (Table ##TAB##0##1##). The remaining 219 ESTs (13.0%) had no homologues in the public protein sequence database. A further BLASTN search against GenBank EST collections revealed that 62 of the 219 ESTs were likely to be novel sequences.</p> <title>Overview of putative genes represented in SAM-derived cDNA libraries</title> <p>More than 50% of the SAM unigenes could be assigned as genes with a known or putative function based on sequence similarity. However, the lack of information on the encoded products meant that many of these transcripts could not be functionally categorized according to the Gene Ontology Consortium (GO). Using BLAST2GO [##REF##16081474##21##], we successfully classified 593 and 557 unigenes in terms of GO molecular functions (Figure ##FIG##0##1A##) and biological processes (Figure ##FIG##0##1B##), respectively. A single gene product might be assigned to more than one GO term, and hence the total number of GO mappings in each of the ontologies exceeded the number of ESTs.</p> <p>The successfully classified unigenes cover a broad range of GO functional categories (Figure ##FIG##0##1##). Under a molecular-function classification, most of the genes (53%) were assigned to the "binding" class (Figure ##FIG##0##1A##). This class includes sequences with putative involvement in mainly nucleic acid binding, a substantial number of which are predicted to encode transcription factors that are known to be essential to the regulation of plant development. Unigenes predicted to encode histone subunits and histone-modification proteins, chromatin remodelling factors and DNA methyltransferases represent another group of sequences linked to nucleic acid binding.</p> <p>To investigate the different types of transcription-factor families represented by our EST libraries, a search using the best-matching <italic>Arabidopsis </italic>locus for SAM unigenes (based on a BLASTX search against the TAIR <italic>Arabidopsis </italic>protein database) was performed at the <italic>Arabidopsis </italic>Gene Regulatory Information Server [##UREF##7##22##]. There are 50 families of transcription factors currently listed in the database, 19 of these are represented by the SAM unigenes (Table ##TAB##1##2##). Six members of the family of homeobox transcription factors that play key roles in the regulation of development may represent interesting candidate genes for further studies.</p> <p>In the biological-processes category, 40% and 35% of the unigenes were involved with cellular physiological processes and metabolism, respectively (Figure ##FIG##0##1B##). The former contains gene products that play an important role in cell organization and biogenesis, while the latter has sequences related mainly to protein metabolism. The sequences relevant to protein metabolism ranged from those associated with protein biosynthesis, such as different ribosomal subunits, to those that modify and degrade proteins, including various sequences involved in the ubiquitin-proteasome pathway.</p> <title>Detection of differentially expressed transcripts in the SAM using the pea 2 K array</title> <p>An oligonucleotide microarray has been developed using our EST resource and about 500 pea sequences randomly retrieved from the GenBank pea protein database, which contains approximately 2000 entries (see Materials and Methods). We utilized this array to compare the transcriptional profiles of the SAM and NM. Four independent replications of balanced-block-design dual-label experiments were performed [Materials and methods, [##REF##12424114##23##]] and the resulting data were depicted in Figure ##FIG##1##2##.</p> <p>Differentially expressed genes that were detected using LimmaGUI [##REF##15297296##24##] (at an adjusted probability value of &lt; 0.05) were subjected to further selection based on the relative change, with up- and down-regulated transcripts defined using cut-offs of greater than 1.3-fold or less than 0.7-fold, respectively. Based on these criteria, we identified 184 and 175 transcripts that were significantly up- or down-regulated in the pea SAM relative to the NM. These transcripts were annotated based on the best BLASTX match, and exhibited changes relative to the NM of 0.1- to 10.7-fold (see Additional File ##SUPPL##0##1## &amp;##SUPPL##1##2##). The change of 0.1-fold was for a gene encoding a type 1 metallothionein (<ext-link ext-link-type="gen" xlink:href="AB176564">AB176564</ext-link>), while that of 10.7-fold was for a sequence annotated as <italic>vegetative lectin </italic>(<ext-link ext-link-type="gen" xlink:href="AAA33691">AAA33691</ext-link>). A study of similar genes encoding type 1 metallothionein in cotton revealed their abundant (although not exclusive) expression in the root [##REF##8790303##25##]. Meanwhile, the high expression of a sequence encoding vegetative lectin in the pea apex has been observed [##REF##16666629##26##]. These and other studies listed in Table ##TAB##2##3## provide independent verification of our microarray data. We also performed RT-PCR analysis on five selected transcripts and as shown in Figure ##FIG##2##3##, the outcome is generally in good agreement with the microarray data.</p> <p>Further comparison of the functional categories identified to be differentially regulated in maize data [##REF##8790303##25##] and this study (Figure ##FIG##3##4##) revealed that transcripts associated with the categories of transcription, chromatin and cell division are highly represented in the list of up-regulated genes, while sequences categorized as being related to metabolism, stress and photosynthesis are down-regulated in the SAM. However, our data also highlight differences between the transcriptional repertoires of the monocot SAM [maize, [##REF##17764504##27##]] and the dicot SAM (pea, present study). In particular, the high retrotransposon-related transcriptional activity reported in maize [##REF##17764504##27##] does not appear to be a conserved feature of the SAM, since it was not present in the pea SAM ESTs.</p> <p>Since the SAM we investigated consists of distinct functional zones, the averaging of signals across the whole SAM probably attenuated signals associated with any given stem-cell region. However, our experiment showed that a sequence annotated as encoding PROTODERMAL FACTOR 1 was up-regulated 1.8-fold in the SAM relative to the NM (Table ##TAB##2##3##). An <italic>Arabidopsis </italic>counterpart gene has been found to be exclusively expressed in the outermost (L1) layer of the SAM [##REF##10483121##28##]. This indicates the ability to identify transcripts that are specifically expressed in certain domains of the SAM, though the magnitude of the relative change might have been reduced and low abundance transcript could have been missed detection.</p> <p>A closer inspection of the differentially expressed transcripts revealed that about 60% of the identified down-regulated sequences were derived from GenBank (Table ##TAB##3##4##). In contrast, only 12% of the genes with higher expression in SAMs were retrieved from GenBank, while the other up-regulated ESTs were from clones derived from our library collection. This is not surprising since the sequences from GenBank were generated primarily from tissues other than the SAM, whereas our EST collection was derived from dissected SAMs. This indicates the significance of constructing a library from specialized tissues and further suggests the utility of our EST collection as a valuable resource in studying the molecular processes underlying the functions of the plant meristem. Meanwhile, 30% of the differentially expressed genes detected were not identified, consisting of genes annotated as hypothetical or expressed protein as well as sequences that have no BLASTX matches in the public database. This list of unknowns could represent intriguing candidate genes for functional analysis.</p> <title>Cell division and proliferation in the SAM</title> <p>Putative functions could be assigned to 116 of the up-regulated transcripts based on protein sequence similarity. Manual inspection of the corresponding transcripts revealed a high representation of ESTs predicted to encode proteins associated with cell division and proliferation (Table ##TAB##3##4##; Figure ##FIG##3##4##).</p> <p>Table ##TAB##3##4## indicates that ESTs encoding all five subtypes of histones, minichromosome-maintenance proteins and cell-proliferating nuclear antigens were among the transcripts whose expression was higher in the SAM than in the NM. The expression of these genes is known to be associated with DNA synthesis and cell proliferation, and they are thus abundant in the meristematic tissue [##REF##11089867##29##,##REF##1709277##30##] since this region consists of actively dividing cells.</p> <p>In the same category, there is a transcript predicted to encode HIGH MOBILITY GROUP (HMG) protein. HMG are proteins that are known to play an architecture role in modifying DNA conformation to facilitate the assembly of multiprotein-DNA complexes. They may serve only to maintain physical orders but the involvement of these proteins in the network regulating SAM activity seems plausible. This is in view of recent studies that demonstrate the binding of HMG proteins to functionally important regions of plant gene promoter and stimulate transcriptions [##REF##11389716##31##].</p> <p>Other up-regulated sequences included transcripts predicted to encode cyclin D (<ext-link ext-link-type="gen" xlink:href="AB008188">AB008188</ext-link>) and cyclin B (<ext-link ext-link-type="gen" xlink:href="EX570164">EX570164</ext-link>), which are involved in the progression of the cell cycle. D-cyclins are one of the main rate-limiting factors for cell proliferation, and several of them may play key roles in the association between the cell cycle and meristem function, in particular primordia formation [reviewed in [##REF##16724263##32##]].</p> <title>Epigenetic regulation of the SAM activity</title> <p>Increases in transcript abundance were also found for genes predicted to encode histone-modification protein (<ext-link ext-link-type="gen" xlink:href="EX569897">EX569897</ext-link>), chromatin remodelling factors (<ext-link ext-link-type="gen" xlink:href="EX570366">EX570366</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570306">EX570306</ext-link> and <ext-link ext-link-type="gen" xlink:href="DQ026703">DQ026703</ext-link>) and DNA methyltransferases (<ext-link ext-link-type="gen" xlink:href="EX568764">EX568764</ext-link> and <ext-link ext-link-type="gen" xlink:href="AF034419">AF034419</ext-link>; Table ##TAB##3##4##). These proteins are known to play a role in the epigenetic regulation of gene expression by participating in mechanisms that alter chromatin structure so as to activate or repress particular sets of genes [reviewed in [##REF##15862099##33##]].</p> <p>There is an emerging recognition of the significance of the chromatin remodelling process in regulating the activity of plant stem cells [reviewed in [##REF##15953752##34##]]. For instance, mutation of the <italic>Arabidopsis FASCIATA1 </italic>(<italic>FAS1</italic>) and <italic>FASCIATA2 </italic>(<italic>FAS2</italic>) genes that encode subunits of the chromatin assembly factor leads to dysfunction of the SAM [##REF##11163246##35##]. This was associated with the down-regulation of <italic>WUS </italic>gene expression in both mutants, suggesting that regulation of meristem function and organogenesis by chromatin remodelling factors is primarily achieved through regulation of the expression of the homeobox transcription factor [##REF##11163246##35##]. However, direct links between chromatin remodelling factors and the regulation of the expression of key meristem genes remain to be established. Nonetheless, the up-regulation of the expression of these sequences in the SAM suggests that the mechanism of epigenetic regulation is important to maintaining the identity of stem cells in plants, as has been reported in animals [##REF##15281084##36##].</p> <title>Transcription factors and hormonal regulation of SAM activity</title> <p>We also identified a few putative transcription factors among the genes whose expression was higher in the SAM than in the NM. These included an EST (<ext-link ext-link-type="gen" xlink:href="EX570025">EX570025</ext-link>) annotated as a putative plant-specific transcription factor from a YABBY-family protein – members of this family are reportedly involved in the abaxial cell-fate specification in lateral organs of <italic>Arabidopsis </italic>[##REF##10457020##37##]. Also on the list was a sequence predicted to encode AUXIN RESPONSE FACTOR 9 (ARF9; <ext-link ext-link-type="gen" xlink:href="EX570238">EX570238</ext-link>). Similar members of the auxin response factor group of proteins are known to regulate auxin-mediated transcript activation or repression. For example, the expression of several genes, such as those encoding members of LATERAL ORGAN BOUNDARIES domain proteins and AUXIN-REGULATED GENE INVOLVED IN ORGAN SIZE, are disrupted in the double mutant of ARF7 and ARF19, implicating their roles in auxin-mediated plant development [##REF##15659631##38##].</p> <p>Other gene products that might be related to auxin include a putative auxin-efflux carrier (PIN1), auxin-conjugating protein (DWARF IN LIGHT1) and a ribosomal protein L24B (Additional File ##SUPPL##0##1##, Table ##TAB##3##4##). Similar genes in <italic>Arabidopsis </italic>have been implicated in auxin-mediated developmental regulation [##REF##11169197##39##, ####REF##16227452##40##, ##REF##12495745##41####12495745##41##]. For example, DWARF IN LIGHT1 is involved in auxin signal transduction, and inhibits shoot and hypocotyl-cell elongation [##REF##11169197##39##]. The presence of several auxin-related transcripts on our list reflects the well-established roles of auxin in organ initiation and positioning at the meristem [##REF##10760240##42##,##REF##14628043##43##].</p> <p>Intriguingly, miRNAs that are endogenous 21-nucleotide riboregulators have been shown to target several mRNAs implicated in auxin responses, including DWARF IN LIGHT1 [##REF##15829600##44##]. There is increasing amount of evidence indicating that miRNA-mediated repression plays an important role in the spatial expression of plant cell-fate regulatory genes [e.g. [##UREF##8##45##]]. A protein called ARGONAUTE is known to function as a catalytic component of the RNA-induced silencing complex, which targets mRNA for degradation using miRNA as a guide [##REF##18073770##46##]. Nevertheless, the precise role for the up-regulation of a transcript annotated as encoding a similar ARGONAUTE in our dataset (Table ##TAB##3##4##) awaits further study.</p> <title>Stress responses in the SAM</title> <p>Many of the down-regulated transcripts were potentially associated with biotic and abiotic stress responses (Table ##TAB##4##5##), including transcripts predicted to encode dehydrin-related protein (<ext-link ext-link-type="gen" xlink:href="AY065655">AY065655</ext-link>), pathogenesis-related protein (<ext-link ext-link-type="gen" xlink:href="AJ586324">AJ586324</ext-link>), disease-resistance-response protein (<ext-link ext-link-type="gen" xlink:href="AF139018">AF139018</ext-link>), antimicrobial defensin (<ext-link ext-link-type="gen" xlink:href="AF525685">AF525685</ext-link>) and chitinase (<ext-link ext-link-type="gen" xlink:href="AB037832">AB037832</ext-link>). This might be attributable to the stems, leaves and roots generally being exposed to greater biotic or abiotic stress than the well-shielded SAM, with the former therefore requiring the constitutive presence of these gene products at a higher level than in the SAM in order to maintain successful defence responses. However, we found at least one other sequence (<ext-link ext-link-type="gen" xlink:href="EX571142">EX571142</ext-link>) potentially related to stress responses whose expression was higher in the SAM (Table ##TAB##3##4##). Although the molecular basis for this is unknown, it is possible that the corresponding encoded product plays dual roles in both stress responses and development. This is supported by a recent study finding a network of rice genes associated with stress responses and seed development [##REF##12684538##47##].</p> <title>ROS homeostasis in the SAM</title> <p>Surprisingly, the down-regulated transcripts included various sequences encoding proteins that scavenge or generate ROS, such as thioredoxin (<ext-link ext-link-type="gen" xlink:href="AJ319808">AJ319808</ext-link>), catalase (<ext-link ext-link-type="gen" xlink:href="EX568770">EX568770</ext-link>), Cu-Zn superoxide dismutase (<ext-link ext-link-type="gen" xlink:href="AB189165">AB189165</ext-link>), copper amine oxidase (<ext-link ext-link-type="gen" xlink:href="AB026253">AB026253</ext-link>) and peroxidase (<ext-link ext-link-type="gen" xlink:href="AB087838">AB087838</ext-link>). This implies that the level of ROS is lower in the SAM than in the NM, and hence the expression of genes encoding ROS scavengers is lower in the former. The absence of photosynthesis (a source of ROS) in the SAM might explain the lower ROS level therein. However, there is increasing evidence that plants use ROS as signalling molecules for regulating development and various physiological responses, and for mediating abscisic-acid-induced stomatal closure, as well as in auxin signalling and gravitropism in roots [##REF##15377225##48##]. Whether this indicates that the regulatory role of ROS is less prominent in the SAM than in the NM awaits further investigation.</p> <title>Photosynthesis- and cell-wall-related transcripts</title> <p>Several of the genes more abundant in the NM than in the SAM were predicted to be related to photosynthesis, including subunits of photosystem I and II and chlorophyll-a/b-binding protein. It is well known that meristematic cells do not contain differentiated plastids, which may explain the lower expression level of genes associated with photosynthesis in the SAM.</p> <p>Some of the sequences predicted to encode products that play roles in cellulose synthesis, such as cellulose synthase (<ext-link ext-link-type="gen" xlink:href="EX569643">EX569643</ext-link>), pectin methylesterase (<ext-link ext-link-type="gen" xlink:href="AF056493">AF056493</ext-link>), xyloglucan endotransglycosylase (<ext-link ext-link-type="gen" xlink:href="AB042531">AB042531</ext-link>), were found to be down-regulated in the SAM. This probably reflects the cell-wall structure in the stem and leaf being more complex that the thin primary cell wall in meristematic cells [##REF##17499544##4##].</p>
<title>Results and Discussion</title> <title>Features of generated ESTs</title> <p>A total of 3000 clones from three cDNA libraries were single-pass sequenced from their 5' ends. Sequence cleaning processes as outlined in the Materials and Methods yielded 2735 ESTs. These sequences had an average trimmed length of 519 base pairs and were assembled into 348 clusters and 1332 singlets, resulting in the final annotation of 1686 unigenes. Clusters ranged in membership from 67 ESTs (one) to two ESTs (253). The redundancy levels were 15.0%, 20.3% and 62.8% within the normalized, standard and subtracted libraries, respectively. A high redundancy level of cDNA libraries constructed using a similar subtraction method from <italic>Thellungiella </italic>plants under abiotic stress conditions has also been observed previously [##REF##16021339##19##].</p> <p>The translated products of the 1686 unigenes were searched against the non-redundant protein database provided by GenBank [##UREF##6##20##] to putatively assign their functions. At the time of writing, 918 (54.4%) unigenes showed significant similarity (<italic>E </italic>value cut-off of 10<sup>-5</sup>) to genes of known or putative function, whereas 549 (32.6%) ESTs were assigned to transcripts with unknown function, which includes hypothetical genes predicted in genomes of model organisms (Table ##TAB##0##1##). The remaining 219 ESTs (13.0%) had no homologues in the public protein sequence database. A further BLASTN search against GenBank EST collections revealed that 62 of the 219 ESTs were likely to be novel sequences.</p> <title>Overview of putative genes represented in SAM-derived cDNA libraries</title> <p>More than 50% of the SAM unigenes could be assigned as genes with a known or putative function based on sequence similarity. However, the lack of information on the encoded products meant that many of these transcripts could not be functionally categorized according to the Gene Ontology Consortium (GO). Using BLAST2GO [##REF##16081474##21##], we successfully classified 593 and 557 unigenes in terms of GO molecular functions (Figure ##FIG##0##1A##) and biological processes (Figure ##FIG##0##1B##), respectively. A single gene product might be assigned to more than one GO term, and hence the total number of GO mappings in each of the ontologies exceeded the number of ESTs.</p> <p>The successfully classified unigenes cover a broad range of GO functional categories (Figure ##FIG##0##1##). Under a molecular-function classification, most of the genes (53%) were assigned to the "binding" class (Figure ##FIG##0##1A##). This class includes sequences with putative involvement in mainly nucleic acid binding, a substantial number of which are predicted to encode transcription factors that are known to be essential to the regulation of plant development. Unigenes predicted to encode histone subunits and histone-modification proteins, chromatin remodelling factors and DNA methyltransferases represent another group of sequences linked to nucleic acid binding.</p> <p>To investigate the different types of transcription-factor families represented by our EST libraries, a search using the best-matching <italic>Arabidopsis </italic>locus for SAM unigenes (based on a BLASTX search against the TAIR <italic>Arabidopsis </italic>protein database) was performed at the <italic>Arabidopsis </italic>Gene Regulatory Information Server [##UREF##7##22##]. There are 50 families of transcription factors currently listed in the database, 19 of these are represented by the SAM unigenes (Table ##TAB##1##2##). Six members of the family of homeobox transcription factors that play key roles in the regulation of development may represent interesting candidate genes for further studies.</p> <p>In the biological-processes category, 40% and 35% of the unigenes were involved with cellular physiological processes and metabolism, respectively (Figure ##FIG##0##1B##). The former contains gene products that play an important role in cell organization and biogenesis, while the latter has sequences related mainly to protein metabolism. The sequences relevant to protein metabolism ranged from those associated with protein biosynthesis, such as different ribosomal subunits, to those that modify and degrade proteins, including various sequences involved in the ubiquitin-proteasome pathway.</p> <title>Detection of differentially expressed transcripts in the SAM using the pea 2 K array</title> <p>An oligonucleotide microarray has been developed using our EST resource and about 500 pea sequences randomly retrieved from the GenBank pea protein database, which contains approximately 2000 entries (see Materials and Methods). We utilized this array to compare the transcriptional profiles of the SAM and NM. Four independent replications of balanced-block-design dual-label experiments were performed [Materials and methods, [##REF##12424114##23##]] and the resulting data were depicted in Figure ##FIG##1##2##.</p> <p>Differentially expressed genes that were detected using LimmaGUI [##REF##15297296##24##] (at an adjusted probability value of &lt; 0.05) were subjected to further selection based on the relative change, with up- and down-regulated transcripts defined using cut-offs of greater than 1.3-fold or less than 0.7-fold, respectively. Based on these criteria, we identified 184 and 175 transcripts that were significantly up- or down-regulated in the pea SAM relative to the NM. These transcripts were annotated based on the best BLASTX match, and exhibited changes relative to the NM of 0.1- to 10.7-fold (see Additional File ##SUPPL##0##1## &amp;##SUPPL##1##2##). The change of 0.1-fold was for a gene encoding a type 1 metallothionein (<ext-link ext-link-type="gen" xlink:href="AB176564">AB176564</ext-link>), while that of 10.7-fold was for a sequence annotated as <italic>vegetative lectin </italic>(<ext-link ext-link-type="gen" xlink:href="AAA33691">AAA33691</ext-link>). A study of similar genes encoding type 1 metallothionein in cotton revealed their abundant (although not exclusive) expression in the root [##REF##8790303##25##]. Meanwhile, the high expression of a sequence encoding vegetative lectin in the pea apex has been observed [##REF##16666629##26##]. These and other studies listed in Table ##TAB##2##3## provide independent verification of our microarray data. We also performed RT-PCR analysis on five selected transcripts and as shown in Figure ##FIG##2##3##, the outcome is generally in good agreement with the microarray data.</p> <p>Further comparison of the functional categories identified to be differentially regulated in maize data [##REF##8790303##25##] and this study (Figure ##FIG##3##4##) revealed that transcripts associated with the categories of transcription, chromatin and cell division are highly represented in the list of up-regulated genes, while sequences categorized as being related to metabolism, stress and photosynthesis are down-regulated in the SAM. However, our data also highlight differences between the transcriptional repertoires of the monocot SAM [maize, [##REF##17764504##27##]] and the dicot SAM (pea, present study). In particular, the high retrotransposon-related transcriptional activity reported in maize [##REF##17764504##27##] does not appear to be a conserved feature of the SAM, since it was not present in the pea SAM ESTs.</p> <p>Since the SAM we investigated consists of distinct functional zones, the averaging of signals across the whole SAM probably attenuated signals associated with any given stem-cell region. However, our experiment showed that a sequence annotated as encoding PROTODERMAL FACTOR 1 was up-regulated 1.8-fold in the SAM relative to the NM (Table ##TAB##2##3##). An <italic>Arabidopsis </italic>counterpart gene has been found to be exclusively expressed in the outermost (L1) layer of the SAM [##REF##10483121##28##]. This indicates the ability to identify transcripts that are specifically expressed in certain domains of the SAM, though the magnitude of the relative change might have been reduced and low abundance transcript could have been missed detection.</p> <p>A closer inspection of the differentially expressed transcripts revealed that about 60% of the identified down-regulated sequences were derived from GenBank (Table ##TAB##3##4##). In contrast, only 12% of the genes with higher expression in SAMs were retrieved from GenBank, while the other up-regulated ESTs were from clones derived from our library collection. This is not surprising since the sequences from GenBank were generated primarily from tissues other than the SAM, whereas our EST collection was derived from dissected SAMs. This indicates the significance of constructing a library from specialized tissues and further suggests the utility of our EST collection as a valuable resource in studying the molecular processes underlying the functions of the plant meristem. Meanwhile, 30% of the differentially expressed genes detected were not identified, consisting of genes annotated as hypothetical or expressed protein as well as sequences that have no BLASTX matches in the public database. This list of unknowns could represent intriguing candidate genes for functional analysis.</p> <title>Cell division and proliferation in the SAM</title> <p>Putative functions could be assigned to 116 of the up-regulated transcripts based on protein sequence similarity. Manual inspection of the corresponding transcripts revealed a high representation of ESTs predicted to encode proteins associated with cell division and proliferation (Table ##TAB##3##4##; Figure ##FIG##3##4##).</p> <p>Table ##TAB##3##4## indicates that ESTs encoding all five subtypes of histones, minichromosome-maintenance proteins and cell-proliferating nuclear antigens were among the transcripts whose expression was higher in the SAM than in the NM. The expression of these genes is known to be associated with DNA synthesis and cell proliferation, and they are thus abundant in the meristematic tissue [##REF##11089867##29##,##REF##1709277##30##] since this region consists of actively dividing cells.</p> <p>In the same category, there is a transcript predicted to encode HIGH MOBILITY GROUP (HMG) protein. HMG are proteins that are known to play an architecture role in modifying DNA conformation to facilitate the assembly of multiprotein-DNA complexes. They may serve only to maintain physical orders but the involvement of these proteins in the network regulating SAM activity seems plausible. This is in view of recent studies that demonstrate the binding of HMG proteins to functionally important regions of plant gene promoter and stimulate transcriptions [##REF##11389716##31##].</p> <p>Other up-regulated sequences included transcripts predicted to encode cyclin D (<ext-link ext-link-type="gen" xlink:href="AB008188">AB008188</ext-link>) and cyclin B (<ext-link ext-link-type="gen" xlink:href="EX570164">EX570164</ext-link>), which are involved in the progression of the cell cycle. D-cyclins are one of the main rate-limiting factors for cell proliferation, and several of them may play key roles in the association between the cell cycle and meristem function, in particular primordia formation [reviewed in [##REF##16724263##32##]].</p> <title>Epigenetic regulation of the SAM activity</title> <p>Increases in transcript abundance were also found for genes predicted to encode histone-modification protein (<ext-link ext-link-type="gen" xlink:href="EX569897">EX569897</ext-link>), chromatin remodelling factors (<ext-link ext-link-type="gen" xlink:href="EX570366">EX570366</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570306">EX570306</ext-link> and <ext-link ext-link-type="gen" xlink:href="DQ026703">DQ026703</ext-link>) and DNA methyltransferases (<ext-link ext-link-type="gen" xlink:href="EX568764">EX568764</ext-link> and <ext-link ext-link-type="gen" xlink:href="AF034419">AF034419</ext-link>; Table ##TAB##3##4##). These proteins are known to play a role in the epigenetic regulation of gene expression by participating in mechanisms that alter chromatin structure so as to activate or repress particular sets of genes [reviewed in [##REF##15862099##33##]].</p> <p>There is an emerging recognition of the significance of the chromatin remodelling process in regulating the activity of plant stem cells [reviewed in [##REF##15953752##34##]]. For instance, mutation of the <italic>Arabidopsis FASCIATA1 </italic>(<italic>FAS1</italic>) and <italic>FASCIATA2 </italic>(<italic>FAS2</italic>) genes that encode subunits of the chromatin assembly factor leads to dysfunction of the SAM [##REF##11163246##35##]. This was associated with the down-regulation of <italic>WUS </italic>gene expression in both mutants, suggesting that regulation of meristem function and organogenesis by chromatin remodelling factors is primarily achieved through regulation of the expression of the homeobox transcription factor [##REF##11163246##35##]. However, direct links between chromatin remodelling factors and the regulation of the expression of key meristem genes remain to be established. Nonetheless, the up-regulation of the expression of these sequences in the SAM suggests that the mechanism of epigenetic regulation is important to maintaining the identity of stem cells in plants, as has been reported in animals [##REF##15281084##36##].</p> <title>Transcription factors and hormonal regulation of SAM activity</title> <p>We also identified a few putative transcription factors among the genes whose expression was higher in the SAM than in the NM. These included an EST (<ext-link ext-link-type="gen" xlink:href="EX570025">EX570025</ext-link>) annotated as a putative plant-specific transcription factor from a YABBY-family protein – members of this family are reportedly involved in the abaxial cell-fate specification in lateral organs of <italic>Arabidopsis </italic>[##REF##10457020##37##]. Also on the list was a sequence predicted to encode AUXIN RESPONSE FACTOR 9 (ARF9; <ext-link ext-link-type="gen" xlink:href="EX570238">EX570238</ext-link>). Similar members of the auxin response factor group of proteins are known to regulate auxin-mediated transcript activation or repression. For example, the expression of several genes, such as those encoding members of LATERAL ORGAN BOUNDARIES domain proteins and AUXIN-REGULATED GENE INVOLVED IN ORGAN SIZE, are disrupted in the double mutant of ARF7 and ARF19, implicating their roles in auxin-mediated plant development [##REF##15659631##38##].</p> <p>Other gene products that might be related to auxin include a putative auxin-efflux carrier (PIN1), auxin-conjugating protein (DWARF IN LIGHT1) and a ribosomal protein L24B (Additional File ##SUPPL##0##1##, Table ##TAB##3##4##). Similar genes in <italic>Arabidopsis </italic>have been implicated in auxin-mediated developmental regulation [##REF##11169197##39##, ####REF##16227452##40##, ##REF##12495745##41####12495745##41##]. For example, DWARF IN LIGHT1 is involved in auxin signal transduction, and inhibits shoot and hypocotyl-cell elongation [##REF##11169197##39##]. The presence of several auxin-related transcripts on our list reflects the well-established roles of auxin in organ initiation and positioning at the meristem [##REF##10760240##42##,##REF##14628043##43##].</p> <p>Intriguingly, miRNAs that are endogenous 21-nucleotide riboregulators have been shown to target several mRNAs implicated in auxin responses, including DWARF IN LIGHT1 [##REF##15829600##44##]. There is increasing amount of evidence indicating that miRNA-mediated repression plays an important role in the spatial expression of plant cell-fate regulatory genes [e.g. [##UREF##8##45##]]. A protein called ARGONAUTE is known to function as a catalytic component of the RNA-induced silencing complex, which targets mRNA for degradation using miRNA as a guide [##REF##18073770##46##]. Nevertheless, the precise role for the up-regulation of a transcript annotated as encoding a similar ARGONAUTE in our dataset (Table ##TAB##3##4##) awaits further study.</p> <title>Stress responses in the SAM</title> <p>Many of the down-regulated transcripts were potentially associated with biotic and abiotic stress responses (Table ##TAB##4##5##), including transcripts predicted to encode dehydrin-related protein (<ext-link ext-link-type="gen" xlink:href="AY065655">AY065655</ext-link>), pathogenesis-related protein (<ext-link ext-link-type="gen" xlink:href="AJ586324">AJ586324</ext-link>), disease-resistance-response protein (<ext-link ext-link-type="gen" xlink:href="AF139018">AF139018</ext-link>), antimicrobial defensin (<ext-link ext-link-type="gen" xlink:href="AF525685">AF525685</ext-link>) and chitinase (<ext-link ext-link-type="gen" xlink:href="AB037832">AB037832</ext-link>). This might be attributable to the stems, leaves and roots generally being exposed to greater biotic or abiotic stress than the well-shielded SAM, with the former therefore requiring the constitutive presence of these gene products at a higher level than in the SAM in order to maintain successful defence responses. However, we found at least one other sequence (<ext-link ext-link-type="gen" xlink:href="EX571142">EX571142</ext-link>) potentially related to stress responses whose expression was higher in the SAM (Table ##TAB##3##4##). Although the molecular basis for this is unknown, it is possible that the corresponding encoded product plays dual roles in both stress responses and development. This is supported by a recent study finding a network of rice genes associated with stress responses and seed development [##REF##12684538##47##].</p> <title>ROS homeostasis in the SAM</title> <p>Surprisingly, the down-regulated transcripts included various sequences encoding proteins that scavenge or generate ROS, such as thioredoxin (<ext-link ext-link-type="gen" xlink:href="AJ319808">AJ319808</ext-link>), catalase (<ext-link ext-link-type="gen" xlink:href="EX568770">EX568770</ext-link>), Cu-Zn superoxide dismutase (<ext-link ext-link-type="gen" xlink:href="AB189165">AB189165</ext-link>), copper amine oxidase (<ext-link ext-link-type="gen" xlink:href="AB026253">AB026253</ext-link>) and peroxidase (<ext-link ext-link-type="gen" xlink:href="AB087838">AB087838</ext-link>). This implies that the level of ROS is lower in the SAM than in the NM, and hence the expression of genes encoding ROS scavengers is lower in the former. The absence of photosynthesis (a source of ROS) in the SAM might explain the lower ROS level therein. However, there is increasing evidence that plants use ROS as signalling molecules for regulating development and various physiological responses, and for mediating abscisic-acid-induced stomatal closure, as well as in auxin signalling and gravitropism in roots [##REF##15377225##48##]. Whether this indicates that the regulatory role of ROS is less prominent in the SAM than in the NM awaits further investigation.</p> <title>Photosynthesis- and cell-wall-related transcripts</title> <p>Several of the genes more abundant in the NM than in the SAM were predicted to be related to photosynthesis, including subunits of photosystem I and II and chlorophyll-a/b-binding protein. It is well known that meristematic cells do not contain differentiated plastids, which may explain the lower expression level of genes associated with photosynthesis in the SAM.</p> <p>Some of the sequences predicted to encode products that play roles in cellulose synthesis, such as cellulose synthase (<ext-link ext-link-type="gen" xlink:href="EX569643">EX569643</ext-link>), pectin methylesterase (<ext-link ext-link-type="gen" xlink:href="AF056493">AF056493</ext-link>), xyloglucan endotransglycosylase (<ext-link ext-link-type="gen" xlink:href="AB042531">AB042531</ext-link>), were found to be down-regulated in the SAM. This probably reflects the cell-wall structure in the stem and leaf being more complex that the thin primary cell wall in meristematic cells [##REF##17499544##4##].</p>
<title>Conclusion</title> <p>The development of our EST collection from the pea SAM represents an important advance towards understanding SAM function and maintenance, especially due to the under-representation of SAM-related transcripts in the public database as demonstrated in this study. Subsequent transcriptional profiling experiments using the microarray constructed from ESTs yielded the transcriptional signatures of the pea SAM, and we have reported a repertoire of transcripts with putative or unknown functions that are differentially regulated in the SAM. <italic>In silico </italic>analysis of the predicted gene products has implicated several processes in the complex molecular network that regulates this developmental process. Future studies of these genes should attempt to reveal how they interact in the complex molecular network that maintains and regulates the dynamics of the SAM.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Despite the importance of the shoot apical meristem (SAM) in plant development and organ formation, our understanding of the molecular mechanisms controlling its function is limited. Genomic tools have the potential to unravel the molecular mysteries of the SAM, and legume systems are increasingly being used in plant-development studies owing to their unique characteristics such as nitrogen fixation, secondary metabolism, and pod development. Garden pea (<italic>Pisum sativum</italic>) is a well-established classic model species for genetics studies that has been used since the Mendel era. In addition, the availability of a plethora of developmental mutants makes pea an ideal crop legume for genomics studies. This study aims to utilise genomics tools in isolating genes that play potential roles in the regulation of SAM activity.</p> <title>Results</title> <p>In order to identify genes that are differentially expressed in the SAM, we generated 2735 ESTs from three cDNA libraries derived from freshly micro-dissected SAMs from 10-day-old garden peas (<italic>Pisum sativum </italic>cv Torsdag). Custom-designed oligonucleotide arrays were used to compare the transcriptional profiles of pea SAMs and non-meristematic tissues. A total of 184 and 175 transcripts were significantly up- or down-regulated in the pea SAM, respectively. As expected, close to 61% of the transcripts down-regulated in the SAM were found in the public database, whereas sequences from the same source only comprised 12% of the genes that were expressed at higher levels in the SAM. This highlights the under-representation of transcripts from the meristematic tissues in the current public pea protein database, and demonstrates the utility of our SAM EST collection as an essential genetic resource for revealing further information on the regulation of this developmental process. In addition to unknowns, many of the up-regulated transcripts are known to encode products associated with cell division and proliferation, epigenetic regulation, auxin-mediated responses and microRNA regulation.</p> <title>Conclusion</title> <p>The presented data provide a picture of the transcriptional profile of the pea SAM, and reveal possible roles of differentially expressed transcripts in meristem function and maintenance.</p>
<title>Abbreviations</title> <p>SAM: shoot apical meristem; NM: non-meristem; EST: expressed sequence tag.</p> <title>Authors' contributions</title> <p>CEW carried out the EST analysis, participated in the microarray experiment and RT-PCR analysis, and drafted the manuscript. PLB and MBS were responsible for design of the project, standardisation &amp; organization of meristem micro-dissection &amp; making EST libraries as well as overall coordination of the experiments and manuscript editing. HO contributed to EST library characterization and EST sequence analysis. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Andrea Merrell and Cathy Jensen for their assistance in the pea meristem dissections. This work was supported by the Australian Research Council Centre of Excellence (grant no. CEO348212) to the University of Melbourne Node of the Centre of Excellence for Integrative Legume Research.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Categorization of SAM unigenes according to the Gene Ontology (GO).</bold> Unigenes with a BLASTX score of &lt; 10<sup>-5 </sup>were classified using the BLAST2GO automated system [##REF##16081474##21##]. Note that a single gene can be assigned to more than one category in the GO classification system.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Results of transcript profiling experiments using the custom designed pea SAM Combimatrix 2 K chip.</bold> Expression profiles of SAM were compared with those of NM tissues. Average plot were generated from the experimental data of four independent biological replicates after normalization (see Materials and Methods for details). Black dotsrepresents data generated from the sequences derived from this study while red applies to those from the GenBank. Green corresponds to empty spot on the array. log(ratio) = log<sub>10</sub>(I<sub>SAM</sub>/I<sub>NM</sub>); log(intensity) = 0.5 log<sub>10</sub>(I<sub>SAM</sub>.I<sub>NM</sub>) where I<sub>SAM </sub>and I<sub>NM </sub>are signal intensities for a transcript in the SAM or the NM tissues, respectively.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Verification of microarray data using RT-PCR analysis.</bold> RT-PCR analysis was carried out under linear amplification conditions for five randomly selected transcripts as indicated. The <italic>actin </italic>gene was used as an internal control. 1, SAM; 2, Leaf; 3, Stem; 4, Root.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Functional categories of transcripts differentially expressed in the pea SAM.</bold> The number of genes differentially expressed between the pea SAM and NM (p &lt; 0.05) with changes in expression higher than indicated cut offs of the fold change are plotted with different functional categories highlighted by colour coding. Transcripts corresponding to cell division and proliferation, epigenetic regulation, and gene silencing are among the up-regulated categories while stress response and the metabolism class of genes are in down-regulated categories.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Sampling strategy for the pea SAM.</bold> Micrograph of a representative pea SAM, with the box showing the location of the tissue sample. Scale bar = 100 μm.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>A summary of the results obtained from BLASTX search against the NCBI (nr) protein database with expect value cut-off at 1e-5.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Category</bold></td><td align="left"><bold>Number of unique ESTs</bold></td></tr></thead><tbody><tr><td align="left">Match to genes with known or putative function</td><td align="left">918 (54.4%)</td></tr><tr><td align="left">Match to unknown or hypothetical genes</td><td align="left">549 (32.6%)</td></tr><tr><td align="left">No hits found</td><td align="left">219 (13.0%)</td></tr><tr><td colspan="2"><hr/></td></tr><tr><td align="left">Total</td><td align="left">1686</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>ESTs related to the ontology of nucleic acid binding.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>SAM Clone</bold></td><td align="left"><bold>Family</bold></td><td align="left"><bold>Matching Arabidopsis Locus<sup>a</sup></bold></td></tr></thead><tbody><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568781">EX568781</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569978">EX569978</ext-link></td><td align="left">CCAAT-HAP5</td><td align="left">At1g07980, At1g08970</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569532">EX569532</ext-link></td><td align="left">MYB-Related</td><td align="left">At1g09770</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569985">EX569985</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569919">EX569919</ext-link></td><td align="left">TUB</td><td align="left">At1g16070, At2g18280</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569797">EX569797</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569977">EX569977</ext-link></td><td align="left">GRAS</td><td align="left">At1g21450, At5g66770</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569951">EX569951</ext-link></td><td align="left">AP2-EREBP</td><td align="left">At1g28360</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571097">EX571097</ext-link></td><td align="left">NAC</td><td align="left">At1g28470</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570087">EX570087</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570238">EX570238</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569491">EX569491</ext-link></td><td align="left">ARF</td><td align="left">At1g30330, At4g23980, At5g20730</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571252">EX571252</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569634">EX569634</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570163">EX570163</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569945">EX569945</ext-link></td><td align="left">WRKY</td><td align="left">At1g30650, At2g03340, At2g47260, At2g24570</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570518">EX570518</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569818">EX569818</ext-link></td><td align="left">bZIP</td><td align="left">At1g42990, At4g34590</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570905">EX570905</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569205">EX569205</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570291">EX570291</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX568744">EX568744</ext-link></td><td align="left">C2H2</td><td align="left">At1g55110, At1g75710, At2g23740, At2g27100</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568999">EX568999</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX571110">EX571110</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570832">EX570832</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX571147">EX571147</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569211">EX569211</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570070">EX570070</ext-link></td><td align="left">Homeobox</td><td align="left">At1g62360, At2g27990, At4g36870, At4g40060, At4g32980, At5g03790</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570236">EX570236</ext-link></td><td align="left">G2-like</td><td align="left">At1g69580</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569752">EX569752</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569559">EX569559</ext-link></td><td align="left">ARID</td><td align="left">At1g76110, At1g76510</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571125">EX571125</ext-link></td><td align="left">BHLH</td><td align="left">At2g27230</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568890">EX568890</ext-link></td><td align="left">CCAAT-HAP3</td><td align="left">At2g37060</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569425">EX569425</ext-link></td><td align="left">MADS</td><td align="left">At2g45660</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569354">EX569354</ext-link></td><td align="left">C2C2-Gata</td><td align="left">At3g06740</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569807">EX569807</ext-link></td><td align="left">SBP</td><td align="left">At3g60030</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Differentially-regulated transcripts with corresponding orthologues known to be highly expressed in shoot apical meristems (SAMs) or non-meristem (NM) tissues.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Probe ID<sup>a</sup></bold></td><td align="left"><bold>Annotation</bold></td><td align="left"><bold>Fold change</bold></td><td align="left"><bold>References</bold></td></tr></thead><tbody><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB176564">AB176564</ext-link></td><td align="left">Metallothionein</td><td align="left">0.1</td><td align="left">Hudspeth et al., 1996</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF029242">AF029242</ext-link></td><td align="left">Dormancy associated gene 1 (DRM1)</td><td align="left">0.4</td><td align="left">Stafstrom et al., 1998</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570325">EX570325</ext-link></td><td align="left">MERISTEMATIC RECEPTOR-LIKE KINASE</td><td align="left">1.4</td><td align="left">Fujita et al., 2003</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568912">EX568912</ext-link></td><td align="left">F-box protein (STAMINA PISTILLOIDA)</td><td align="left">1.8</td><td align="left">Taylor et al., 2001</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570203">EX570203</ext-link></td><td align="left">Protodermal factor 1</td><td align="left">1.8</td><td align="left">Abe et al., 2001</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568701">EX568701</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX571325">EX571325</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX569084">EX569084</ext-link>, <ext-link ext-link-type="gen" xlink:href="EX570634">EX570634</ext-link></td><td align="left">Histone subunits</td><td align="left">2.0</td><td align="left">References in Meshi et al., 2000</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570428">EX570428</ext-link></td><td align="left">Mini-chromosome maintenance proteins</td><td align="left">2.2</td><td align="left">Stevens et al., 2002</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570270, AB008186">EX570270, AB008186</ext-link></td><td align="left">Proliferating cell nuclear antigen</td><td align="left">4.8, 3.8</td><td align="left">Kosugi et al., 1991</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB031227">AB031227</ext-link></td><td align="left">PsAD1</td><td align="left">5.3</td><td align="left">Madoka &amp;Mori, 2000</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570531">EX570531</ext-link></td><td align="left">Vegetative lectin</td><td align="left">10.7</td><td align="left">Dobres &amp; Thompson, 1988</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Representative transcripts that are detected to be significantly up-regulated in the pea shoot apical meristem (SAM) in comparison to the non-meristem (NM) tissues.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Probe<sup>a</sup></bold></td><td align="left"><bold>Fold Change</bold></td><td align="left"><bold>Annotation<sup>b</sup></bold></td></tr></thead><tbody><tr><td align="left" colspan="3"><bold>Cell division and proliferation</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570197">EX570197</ext-link></td><td align="left">2.7</td><td align="left">High mobility group protein (HMGI/Y)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571067">EX571067</ext-link></td><td align="left">1.6</td><td align="left">SAR-DNA binding protein 2</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571065">EX571065</ext-link></td><td align="left">1.7</td><td align="left">SAR-DNA binding protein 1</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570428">EX570428</ext-link></td><td align="left">1.8</td><td align="left">Mini-chromosome maintenance protein 6</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569064">EX569064</ext-link></td><td align="left">2.0</td><td align="left">Small nuclear ribonucleoprotein associated protein B</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568701">EX568701</ext-link></td><td align="left">2.0</td><td align="left">Histone H1</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568991">EX568991</ext-link></td><td align="left">2.1</td><td align="left">Delta DNA polymerase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570486">EX570486</ext-link></td><td align="left">2.1</td><td align="left">Helicase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570226">EX570226</ext-link></td><td align="left">2.2</td><td align="left">Phosphoesterase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571345">EX571345</ext-link></td><td align="left">3.3</td><td align="left">Germinal histone H4</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569084">EX569084</ext-link></td><td align="left">3.7</td><td align="left">Histone H2a.1</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570634">EX570634</ext-link></td><td align="left">4.2</td><td align="left">Histone H3</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570270">EX570270</ext-link></td><td align="left">4.8</td><td align="left">Proliferating cell nuclear antigen 2</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568755">EX568755</ext-link></td><td align="left">5.9</td><td align="left">Histone H4</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB008188">AB008188</ext-link></td><td align="left">1.9</td><td align="left">Cyclin D3.1 protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570164">EX570164</ext-link></td><td align="left">1.8</td><td align="left">Mitotic cyclin B1-1</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570388">EX570388</ext-link></td><td align="left">1.5</td><td align="left">Cell cycle protein kinase</td></tr><tr><td align="left" colspan="3"><bold>Protein synthesis and turnover</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569344">EX569344</ext-link></td><td align="left">1.4</td><td align="left">Transducin family of protein (SLOW WALKER 1)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570594">EX570594</ext-link></td><td align="left">1.9</td><td align="left">60S ribosomal protein L18a</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569298">EX569298</ext-link></td><td align="left">1.9</td><td align="left">Ribosomal protein L30</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569956">EX569956</ext-link></td><td align="left">2.0</td><td align="left">40S ribosomal protein S18</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568807">EX568807</ext-link></td><td align="left">2.1</td><td align="left">Fibrillarin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570044">EX570044</ext-link></td><td align="left">2.1</td><td align="left">60S ribosomal protein L44</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570565">EX570565</ext-link></td><td align="left">2.1</td><td align="left">Ribosomal protein L23</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568742">EX568742</ext-link></td><td align="left">2.2</td><td align="left">Ribosomal protein L24</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570468">EX570468</ext-link></td><td align="left">2.2</td><td align="left">60S ribosomal protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568861">EX568861</ext-link></td><td align="left">2.4</td><td align="left">40S ribosomal protein S17</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569048">EX569048</ext-link></td><td align="left">2.6</td><td align="left">Ribosomal protein S2</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568908">EX568908</ext-link></td><td align="left">2.9</td><td align="left">Ribosomal protein S15-like</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568954">EX568954</ext-link></td><td align="left">3.0</td><td align="left">Ribosomal protein S3</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB021873">AB021873</ext-link></td><td align="left">3.2</td><td align="left">Ribosome sedimenting protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570400">EX570400</ext-link></td><td align="left">5.0</td><td align="left">Ribosomal protein S4</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569708">EX569708</ext-link></td><td align="left">2.3</td><td align="left">Ubiquitin extension protein 2</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568912">EX568912</ext-link></td><td align="left">1.8</td><td align="left">F-box protein (STAMINA PISTILLOIDA)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570030">EX570030</ext-link></td><td align="left">1.9</td><td align="left">Chaperonin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568849">EX568849</ext-link></td><td align="left">1.6</td><td align="left">Elongation factor 1-beta/EF-1-beta</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570437">EX570437</ext-link></td><td align="left">1.4</td><td align="left">Cyclophilin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568710">EX568710</ext-link></td><td align="left">1.4</td><td align="left">Eukaryotic initiation factor</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570401">EX570401</ext-link></td><td align="left">1.4</td><td align="left">Translation initiation protein</td></tr><tr><td align="left" colspan="3"><bold>Transcription factors or hormonal regulation</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569978">EX569978</ext-link></td><td align="left">1.5</td><td align="left">Putative Hap5 transcription factor</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570025">EX570025</ext-link></td><td align="left">1.5</td><td align="left">YABBY family transcription factor</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570238">EX570238</ext-link></td><td align="left">1.5</td><td align="left">Auxin response factor 9</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570380">EX570380</ext-link></td><td align="left">3.4</td><td align="left">Putative transcriptional co-activator (KIWI)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571172">EX571172</ext-link></td><td align="left">1.6</td><td align="left">Transcription factor</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570262">EX570262</ext-link></td><td align="left">2.1</td><td align="left">Argonaute protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570187">EX570187</ext-link></td><td align="left">2.9</td><td align="left">Indole-3-acetic acid amido synthetase (DWARF IN LIGHT 1)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF325121">AF325121</ext-link></td><td align="left">1.6</td><td align="left">Brassinosteroid biosynthetic protein</td></tr><tr><td align="left" colspan="3"><bold>Epigenetic regulation-related</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568764">EX568764</ext-link></td><td align="left">1.4</td><td align="left">Chromomethylase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF034419">AF034419</ext-link></td><td align="left">1.4</td><td align="left">DNA methyltransferase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569897">EX569897</ext-link></td><td align="left">1.6</td><td align="left">H3 lysine-9 specific SUVH4</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570366">EX570366</ext-link></td><td align="left">1.9</td><td align="left">SWI/SNF-like ATPase subunits, DDM1</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570306">EX570306</ext-link></td><td align="left">2.0</td><td align="left">WD-40 repeat protein (MSI3)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="DQ026703">DQ026703</ext-link></td><td align="left">1.5</td><td align="left">WD-40 repeat protein (MSI1)</td></tr><tr><td align="left" colspan="3"><bold>Other developmental regulation-related</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569130">EX569130</ext-link></td><td align="left">2.0</td><td align="left">Mandelonitrile lyase family of FAD containing oxidoreductases (HOTHEAD)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY343326">AY343326</ext-link></td><td align="left">1.5</td><td align="left">Late-flowering gene</td></tr><tr><td align="left" colspan="3"><bold>Unclassified/Unknown/No BLASTX match</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571192">EX571192</ext-link></td><td align="left">9.0</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569831">EX569831</ext-link></td><td align="left">3.3</td><td align="left">Proline-rich protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570071">EX570071</ext-link></td><td align="left">3.7</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571142">EX571142</ext-link></td><td align="left">3.6</td><td align="left">Pathogenesis-related group 5 protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX571013">EX571013</ext-link></td><td align="left">2.8</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570256">EX570256</ext-link></td><td align="left">2.4</td><td align="left">Unnamed protein product</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569233">EX569233</ext-link></td><td align="left">2.1</td><td align="left">Serine/threonine dehydratase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569274">EX569274</ext-link></td><td align="left">4.5</td><td align="left">Lipid transfer protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570925">EX570925</ext-link></td><td align="left">3.1</td><td align="left">Kunitz inhibitor ST1-like</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB032830">AB032830</ext-link></td><td align="left">3.3</td><td align="left">Endo-1,4-beta-glucanase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570531">EX570531</ext-link></td><td align="left">10.7</td><td align="left">Vegetative lectin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569068">EX569068</ext-link></td><td align="left">2.4</td><td align="left">Expressed protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569663">EX569663</ext-link></td><td align="left">2.1</td><td align="left">Expressed protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570209">EX570209</ext-link></td><td align="left">2.2</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570712">EX570712</ext-link></td><td align="left">2.0</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570576">EX570576</ext-link></td><td align="left">1.9</td><td align="left">Cytochrome P450 monoxygenase(CYP78A8)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569596">EX569596</ext-link></td><td align="left">1.8</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570325">EX570325</ext-link></td><td align="left">1.5</td><td align="left">MERISTEMATIC RECEPTOR-LIKE KINASE 2</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570536">EX570536</ext-link></td><td align="left">2.4</td><td align="left">Unknown protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570361">EX570361</ext-link></td><td align="left">1.5</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570202">EX570202</ext-link></td><td align="left">1.4</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570695">EX570695</ext-link></td><td align="left">1.4</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569023">EX569023</ext-link></td><td align="left">1.6</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570276">EX570276</ext-link></td><td align="left">1.5</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569420">EX569420</ext-link></td><td align="left">1.4</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568900">EX568900</ext-link></td><td align="left">1.6</td><td align="left">No BLASTX match</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>Representative transcripts that are detected to be significantly down-regulated in the pea shoot apical meristem (SAM) in comparison to the non-meristem (NM) tissues.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Probe<sup>a</sup></bold></td><td align="left"><bold>Fold Change</bold></td><td align="left"><bold>Annotation<sup>b</sup></bold></td></tr></thead><tbody><tr><td align="left" colspan="3"><bold>Stress responses</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY065655">AY065655</ext-link></td><td align="left">0.1</td><td align="left">Ultraviolet B-repressible dehydrin-related protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY065659">AY065659</ext-link></td><td align="left">0.1</td><td align="left">Ultraviolet B- inducible protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ586324">AJ586324</ext-link></td><td align="left">0.1</td><td align="left">Putative basic PR1 protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568748">EX568748</ext-link></td><td align="left">0.3</td><td align="left">Dehydrin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF139018">AF139018</ext-link></td><td align="left">0.3</td><td align="left">Disease resistance protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF525685">AF525685</ext-link></td><td align="left">0.4</td><td align="left">Antimicrobial defensin</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ278699">AJ278699</ext-link></td><td align="left">0.5</td><td align="left">Protease</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB087832">AB087832</ext-link></td><td align="left">0.4</td><td align="left">Class 1 chitinase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF175278">AF175278</ext-link></td><td align="left">0.7</td><td align="left">Wound-inducible P450 hydrolase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF137351">AF137351</ext-link></td><td align="left">0.2</td><td align="left">Pathogenesis-related protein 4</td></tr><tr><td align="left" colspan="3"><bold>Reactive oxygen species homeostasis</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568770">EX568770</ext-link></td><td align="left">0.3</td><td align="left">Catalase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB026253">AB026253</ext-link></td><td align="left">0.5</td><td align="left">Copper amine oxidase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ50832">AJ50832</ext-link></td><td align="left">0.4</td><td align="left">Germine-like protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB189165">AB189165</ext-link></td><td align="left">0.4</td><td align="left">Copper zinc superoxide dismutase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB087837">AB087837</ext-link></td><td align="left">0.7</td><td align="left">Glutathione-S-transferase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ319808">AJ319808</ext-link></td><td align="left">0.3</td><td align="left">Thioredoxin H</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB087838">AB087838</ext-link></td><td align="left">0.5</td><td align="left">Peroxidase</td></tr><tr><td align="left" colspan="3"><bold>Photosynthesis</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569880">EX569880</ext-link></td><td align="left">0.1</td><td align="left">Light harvesting chlorophyll a/b binding protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY845255">AY845255</ext-link></td><td align="left">0.2</td><td align="left">Light harvesting chlorophyll a/b binding protein 3</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569551">EX569551</ext-link></td><td align="left">0.2</td><td align="left">Light harvesting chlorophyll a/b binding protein type 1 (CAB)</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY292531">AY292531</ext-link></td><td align="left">0.1</td><td align="left">Oxygen-evolving enhancer protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569675">EX569675</ext-link></td><td align="left">0.3</td><td align="left">Type II chlorophyll a/b binding protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569989">EX569989</ext-link></td><td align="left">0.4</td><td align="left">PSI-K subunit of Photosystem I</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY007467">AY007467</ext-link></td><td align="left">0.4</td><td align="left">Photosystem II CP47 protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570832">EX570832</ext-link></td><td align="left">0.1</td><td align="left">Ribulose 1,5-biphosphate carboxylase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY065656">AY065656</ext-link></td><td align="left">0.2</td><td align="left">RUBISCO activase</td></tr><tr><td align="left" colspan="3"><bold>Metabolism</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY112702">AY112702</ext-link></td><td align="left">0.6</td><td align="left">Vacuolar acid invertase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ012080">AJ012080</ext-link></td><td align="left">0.3</td><td align="left">Sucrose synthase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="Y08728">Y08728</ext-link></td><td align="left">0.6</td><td align="left">ADP-glucose phosphorylase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569667">EX569667</ext-link></td><td align="left">0.4</td><td align="left">Ribulose-5-phosphate-3-epimerase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570956">EX570956</ext-link></td><td align="left">0.7</td><td align="left">Phosphofructokinase</td></tr><tr><td align="left" colspan="3"><bold>Cell wall modification</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF056493">AF056493</ext-link></td><td align="left">0.7</td><td align="left">Pectin methylesterase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB042531">AB042531</ext-link></td><td align="left">0.2</td><td align="left">Xyloglucan endotransglycosylase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569643">EX569643</ext-link></td><td align="left">0.6</td><td align="left">Cellulose synthase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ621355">AJ621355</ext-link></td><td align="left">0.4</td><td align="left">KORRIGAN</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB015428">AB015428</ext-link></td><td align="left">0.6</td><td align="left">Endoxyloglucan transferase 1</td></tr><tr><td align="left" colspan="3"><bold>Membrane- or transport-related</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ243307">AJ243307</ext-link></td><td align="left">0.1</td><td align="left">Putative plasma membrane intrinsic protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AJ243309">AJ243309</ext-link></td><td align="left">0.3</td><td align="left">Putative tonoplast intrinsic protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AB027616">AB027616</ext-link></td><td align="left">0.5</td><td align="left">Apyrase</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF109922">AF109922</ext-link></td><td align="left">0.7</td><td align="left">Sucrose transport protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570713">EX570713</ext-link></td><td align="left">0.6</td><td align="left">Aquaporin-like transmembrane protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570839">EX570839</ext-link></td><td align="left">0.7</td><td align="left">Outward-rectifying potassium channel</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569706">EX569706</ext-link></td><td align="left">0.7</td><td align="left">Sulfate transporter</td></tr><tr><td align="left" colspan="3"><bold>Unclassified/Unknown/No BLASTX match</bold></td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AY371200">AY371200</ext-link></td><td align="left">0.2</td><td align="left">Ripening-related protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF369889">AF369889</ext-link></td><td align="left">0.5</td><td align="left">Embryo-abundant protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX568773">EX568773</ext-link></td><td align="left">0.5</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569332">EX569332</ext-link></td><td align="left">0.6</td><td align="left">Hypothetical protein</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX569635">EX569635</ext-link></td><td align="left">0.5</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="EX570920">EX570920</ext-link></td><td align="left">0.7</td><td align="left">No BLASTX match</td></tr><tr><td align="left"><ext-link ext-link-type="gen" xlink:href="AF515795">AF515795</ext-link></td><td align="left">0.1</td><td align="left">Dormancy-associated protein 3</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Transcripts identified to be up-regulated in the pea shoot apical meristem (SAM) in comparison to the non-meristem (NM) tissues.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Transcripts identified to be down-regulated in the pea shoot apical meristem (SAM) in comparison to the non-meristem (NM) tissues.</p></caption></supplementary-material>
<table-wrap-foot><p>Arabidopsis transcription factor family represented by the shoot apical meristem (SAM) unigenes. The matching Arabidopsis locus for each clone is indicated.</p><p>a. Annotation is based on the best BLASTX match against TAIR Arabidopsis protein database (Evalue &lt; 1e-5)</p></table-wrap-foot> <table-wrap-foot><p>The data reported is in good agreement with the transcript profiling data.</p><p>a. GenBank accession number that begins with EX is derived from this study.</p></table-wrap-foot> <table-wrap-foot><p>a, Probe with GenBank accession number that begins with EX is derived from this study.</p><p>b. Annotation is based on the expect value of BLASTX hit that is set at 1e-5</p></table-wrap-foot> <table-wrap-foot><p>a, Probe with GenBank accession number that begins with EX is derived from this study.</p><p>b. Annotation is based on the expect value of BLASTX hit that is set at 1e-5</p></table-wrap-foot>
<graphic xlink:href="1471-2229-8-73-1"/> <graphic xlink:href="1471-2229-8-73-2"/> <graphic xlink:href="1471-2229-8-73-3"/> <graphic xlink:href="1471-2229-8-73-4"/> <graphic xlink:href="1471-2229-8-73-5"/>
<media xlink:href="1471-2229-8-73-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2229-8-73-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Blixt", "RC K"], "given-names": ["S"], "article-title": ["The pea"], "source": ["The Handbook of Genetics"], "year": ["1974"], "volume": ["Vol 2"], "publisher-name": [" New York: Plenum Press"], "fpage": ["181"], "lpage": ["221"]}, {"surname": ["Weeden", "Muffet"], "given-names": ["NF", "M"], "article-title": ["Identification of genes affecting root mass and root/shoot ratio in a JI 1794 x \u2018Slow\u2019 RIL population"], "source": ["Pisum Genetics"], "year": ["2002"], "volume": ["34"], "fpage": ["28"], "lpage": ["31"]}, {"surname": ["Murfet", "Reid", "Hebblethwaite PD, Heath MC, Dawkins TCK"], "given-names": ["IC", "JB"], "article-title": ["The control of flowering and internode length in Pisum."], "source": ["The pea crop"], "year": ["1985"], "publisher-name": [" London:  Butterworth"], "fpage": ["67"], "lpage": ["80"]}, {"surname": ["Singer", "Maki", "Sollinger", "Plotz", "Fitzgerald", "Fishbach", "Mullen"], "given-names": ["S", "S", "J", "J", "K", "J", "H"], "article-title": ["Flower development in pea: Role of Proliferating Inflorescence Meristem (PIM), an AP1 homolog"], "source": ["Developmental Biology"], "year": ["2002"], "volume": ["247"], "fpage": ["518"], "lpage": ["518"]}, {"surname": ["Sollinger", "Singer"], "given-names": ["JD", "SR"], "article-title": ["Unraveling the flower with pea developmental mutants-homologies and hidden potentials"], "source": ["Developmental Biology"], "year": ["2002"], "volume": ["247"], "fpage": ["444"], "lpage": ["445"]}, {"surname": ["Yaxley", "Jablonski", "Reid"], "given-names": ["JL", "W", "JB"], "article-title": ["Leaf and flower development in pea ("], "italic": ["Pisum sativum", "cochleata", "unifoliata"], "source": ["Annals of Botany"], "year": ["2001"], "volume": ["88"], "fpage": ["225"], "lpage": ["234"], "pub-id": ["10.1006/anbo.2001.1448"]}, {"article-title": ["National Centre for Biotechnology Information"]}, {"article-title": ["Arabidopsis Gene Regulatory Information Server"]}, {"surname": ["Nagasaki", "Itoh", "Hayashi", "Hibara", "Satoh-Nagasawa", "Nosaka", "Mukouhata", "Ashikari", "Kitano", "Matsuoka", "Nagato", "Sato"], "given-names": ["H", "J", "K", "K", "N", "M", "M", "M", "H", "M", "Y", "Y"], "article-title": ["The small interfering RNA production pathway is required for shoot meristem initiation in rice"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["2007"], "volume": ["104"], "fpage": ["14867"], "lpage": ["14871"], "pub-id": ["10.1073/pnas.0704339104"]}, {"article-title": ["VERTIS"]}, {"surname": ["Rasmussen"], "given-names": ["SW"], "article-title": [" SEQtools, a software package for analysis of nucleotide and protein sequences"], "year": ["2006"]}, {"article-title": ["Combimatrix Corporation"]}]
BMC Plant Biol. 2008 Jun 30; 8:73