A genomic approach to investigate developmental cell death in woody tissues of Populustrees
© Moreau et al.; licensee BioMed Central Ltd. 2005
Received: 17 November 2004
Accepted: 21 February 2005
Published: 22 March 2005
Poplar (Populus sp.) has emerged as the main model system for molecular and genetic studies of forest trees. A Populus expressed sequence tag (EST) database (POPULUSDB) was previously created from 19 cDNA libraries each originating from different Populus tree tissues, and opened to the public in September 2004. We used this dataset for in silico transcript profiling of a particular process in the woody tissues of the Populus stem: the programmed death of xylem fibers.
One EST library in POPULUSDB originates from woody tissues of the Populus stem where xylem fibers undergo cell death. Analysis of EST abundances and library distribution within the POPULUSDB revealed a large number of previously uncharacterized transcripts that were unique in this library and possibly related to the death of xylem fibers. The in silico analysis was complemented by a microarray analysis utilizing a novel Populus cDNA array with a unigene set of 25,000 sequences.
In silico analysis, combined with the microarray analysis, revealed the usefulness of non-normalized EST libraries in elucidating transcriptional regulation of previously uncharacterized physiological processes. The data suggested the involvement of two novel extracellular serine proteases, nodulin-like proteins and an Arabidopsis thaliana OPEN STOMATA 1 (AtOST1) homolog in signaling fiber-cell death, as well as mechanisms responsible for hormonal control, nutrient remobilization, regulation of vacuolar integrity and autolysis of the dying fibers.
The woody tissues of angiosperm trees, the xylem fibers and vessels, are formed from the lateral meristem of the stem, the vascular cambium. In contrast to vessel elements, which differentiate very rapidly close to the vascular cambium, fiber differentiation is a relatively slow process involving initial expansion of the cells in both the radial and longitudinal dimensions, followed by extensive synthesis of the secondary cell walls. The final phase in maturation of both vessel elements and fibers is cell death and autolysis of the cell contents.
Xylem-cell death involves a range of morphological and nuclear changes in a strictly spatially and temporally coordinated and programmed fashion [1, 2]. The programmed cell death (PCD) of xylem has been analyzed in detail in an in vitro system of Zinnia elegans, in which mesophyll cells of Zinnia transdifferentiate into xylem vessels commonly called as tracheary elements in a semi-synchronized manner . In Zinnia cells, irreversible differentiation into tracheary elements is marked by the accumulation of hydrolytic enzymes in the vacuole and deposition of the secondary cell walls, followed by tonoplast disruption, release of the vacuolar proteases and nucleases into the cytoplasm, and finally the autolytic loss of cell contents [2, 4]. Several different types of proteases have been detected in Zinnia [5, 6], and an S1-type nuclease, capable of hydrolyzing both DNA and RNA, seems to control nuclear DNA degradation in the Zinnia tracheary elements . Even though the chain of events during tracheary element PCD is well characterized in the Zinnia system, very little is known about the regulation of this process in intact plants. In addition, the Zinnia system has not allowed analysis of the different cell types of the xylem, such as the fibers.
Programmed cell death also occurs in plants in response to external factors, such as avirulent pathogens, giving rise to the so-called hypersensitive response (HR) and in response to shortening daylength - manifested in the senescence of leaves. HR cell death is usually fast and it shares certain features with the apoptotic death of animal cells, such as nuclear shrinkage and fragmentation of DNA into oligonucleosomal multiples of 180-bp fragments . Senescence-induced cell death is a much slower process, involving nuclear degradation, DNA fragmentation and thorough proteolytic degradation of the cellular contents and controlled remobilization of the nutrients . The death of the xylem elements is different from HR and senescence-related PCD in that the organellar structure remains intact until vacuolar collapse and the oligonucleosomal DNA fragmentation does not precede cell death . Whether these processes are related at the molecular level is unknown, but the differences in temporal and spatial regulation, and in cellular morphology, suggest that there are significant differences not only in the early regulation, but also in the execution of the various plant PCD processes.
The genus Populus has emerged as the main model system for trees, because of its amenability for genomic and molecular analyses . Populus is also suitable for analysis of xylem development [11, 12]. A Populus expressed sequence tag (EST) database (POPULUSDB) was created from 19 different cDNA libraries . The database consists of 102,019 ESTs, assembled into a unigene set of 11,885 clusters and 12,759 non-clustered singletons corresponding altogether to 24,644 unique sequences or transcripts . The great diversity of the tissue types giving rise to the different cDNA libraries enables digital analysis of gene expression by comparison of the EST frequencies in the different libraries. One of the libraries was produced from Populus woody tissues composed of xylem fibers undergoing cell death. In this work, we studied gene expression in the process of fiber death by in silico analysis of this 'fiber death library' and by a microarray analysis with a novel Populus 25K cDNA microarray. In addition to its economic importance as one of the processes that regulate wood quality, fiber-cell death is an interesting biological process that as yet is poorly understood. Our analysis identified several novel candidate regulatory genes for xylem PCD.
Results and discussion
The unique characteristics of fiber-cell death in Populuswood
The most abundant transcripts during fiber-cell death
The most abundant transcript in the fiber death library was glycine hydroxymethyltransferase (GHMT; POPLAR.161). GHMT was also highly abundant in several other libraries derived from xylem-containing tissues, such as the tension wood and the roots (Figure 3). GHMT has also been identified as one of the most abundant proteins in Populus xylem , Arabidopsis roots  and loblolly pine xylem . GHMT is a key enzyme in one-carbon metabolism, catalyzing reversible conversion between serine and glycine to produce 5,10-methylenetetrahydrofolate, which can be used to recover methionine from 5-methyl-tetrahydrofolate and homocysteine . One-carbon metabolism is known to be active in photorespiration, but its preferential expression in the late maturing fibers suggests that the Populus GHMT is also involved in some other process(es). Also, three other enzymes that participate in one-carbon metabolism were all highly abundant in the fiber death library. 5-Methyltetrahydropteroyltriglutamate-homocysteine S-methyltransferase (POPLAR.649) catalyzes biosynthesis of methionine, while S-adenosylmethionine synthetase (POPLAR.155) and adenosylhomocysteinase (POPLAR.147) are involved in methionine catabolism (Figure 3). It is possible that one-carbon metabolism is required during xylem maturation for the production of glycine, which is abundant in the cell wall proteins. S-adenosylmethionine, which is synthesized from methionine, can also be used for methylation reactions which occur during secondary wall formation. It is, however, unlikely that these enzymes are only needed for fiber-cell death, because they are also highly abundant in libraries derived from the cambial zone and tension wood (Figure 3).
Cell-wall proteins were highly abundant in the fiber death library (Figure 3). Arabinogalactan proteins (AGP) belong to a superfamily of proteoglycans encompassing several subclasses, such as classical AGPs and fasciclin-like AGPs. Both classical AGPs (POPLAR.862 and POPLAR.999) and a fasciclin-like AGP (POPLAR.3203) were highly abundant in the fiber death library (Figure 3). AGPs are a major proteinaceous constituent of the cell walls, but their function has remained unclear. A recent report on a hybrid protein containing AGP domains supports the hypothesis that AGPs have a critical function in mediating cell-cell interactions during vascular differentiation . The fasciclin-like AGPs, on the other hand, seem to be involved in the control of cell adhesion [19, 20]. AGPs have also been implicated in PCD [21, 22]. AGP genes in clusters POPLAR.999 and POPLAR.3203 were also highly expressed in other libraries, especially the library from the tension wood-forming xylem tissues, suggesting that they have a more general function during cell-wall formation (Figure 3). However, the POPLAR.862 AGP was enriched in the fiber death library. POPLAR.862 has also been shown to be suppressed in a microarray analysis of tension wood where fiber-cell death is inhibited (S. Andersson-Gunnerås, E. Mellerowicz and B. Sundberg, personal communication), strongly indicating a role for this AGP in the stimulation of fiber-cell death. Two other cell-wall proteins, a glycine-rich protein (POPLAR.1776) and an extensin-like protein (POPLAR.9554), were also highly abundant and overrepresented in the fiber death library (Figure 3), suggesting that these proteins also have functions during the late maturation of xylem fibers.
The fiber death library-specific transcripts are putatively novel regulators of cell death
The unique transcripts in the fiber death library within the POPULUSDB
Number of ESTs
Most similar Arabidopsis gene
Expression ratio B/A
Glycosyl hydrolase family 1 protein
Proline-rich protein family
Oligopeptide transport family protein
Transferase family protein
Basic helix-loop-helix family protein
PHD-type zinc finger protein
None of the transcripts that were shown to be unique to the fiber death library has previously been implicated in the regulation of cell death. The oligopeptide transporter (POPLAR.11639) is involved in amino-acid metabolism related to remobilization of nutrients from dying tissues, but the reason for the specific expression pattern of the other transcripts is not clear. Several of them seem to be membrane proteins or targeted to the endomembrane system, as predicted on the basis of their closest Arabidopsis homologs (Table 1). A glycosyl hydrolase of family 1 (POPLAR.11628) is the most abundant unique transcript in the fiber death library, but the substrates of this glucosidase and its exact role in fiber-cell death remain to be elucidated. Interestingly, the closest homologs of this protein are cyanogenic in nature, and it is tempting to speculate that cyanide production by this enzyme could participate in regulation of cell death in xylem fibers. Reverse transcription PCR (RT-PCR) analysis of five unique transcripts confirmed that POPLAR.11628 (nine ESTs), POPLAR.11639 (five ESTs) and POPLAR.11646 (four ESTs) were indeed very specifically expressed in xylem fibers undergoing cell death (Figure 4). POPLAR.11658 (six ESTs) and POPLAR.11624 (two ESTs) showed highest expression in the dying xylem fibers, but were also expressed in other types of Populus tree tissues (Figure 4).
Combination of the in silicoanalysis with a microarray analysis refines selection of candidate regulatory genes
Among the most upregulated transcripts that were unique or highly enriched in the fiber death library, there are transcripts that participate in amino-acid metabolism and transport (X024D04, POPLAR.872), proteolysis (POPLAR.11632, POPLAR.4995, POPLAR.10724, X077D02), and also transcripts, such as kinases (X053F08, POPLAR.9347), nodulin-like proteins (X002G05, POPLAR.4667) and a CACTIN-like protein (POPLAR.11667), that seem to have signaling functions rather than mere cellular disintegration of the dying fibers (Figure 5). High expression of nodulin-like proteins in dying xylem fibers suggests that nodulins, which regulate nodule formation in response to Rhizobium infection, also have an important function in the maturation of xylem fibers. Interestingly, certain nodulins have been shown to regulate accumulation of reactive oxygen species (ROS) [24, 25], and it is possible that nodulin-like proteins regulate ROS accumulation that occurs during the late maturation or cell death of xylem fibers. ROS accumulation is also implied by the high expression levels of peroxidases (POPLAR.11659 and 11669) and a protein kinase (singleton X053F08) that seems to encode the Populus ortholog of Arabidopsis thaliana OPEN STOMATA 1 (OST1; Figure 5). OST1 regulates abscisic acid (ABA)-mediated accumulation of ROS related to stomatal closure in Arabidopsis, and it has also been shown to be expressed in vascular tissues of Arabidopsis leaves and roots . Other important proteins in the early signaling of fiber-cell death could include the basic helix-loop-helix transcription factor V031H02 and the C2H2 class zinc-finger protein POPLAR.10810. Interestingly, POPLAR.11144 is most similar to the Arabidopsis gene VACUOLELESS1, which is required for proper vacuole formation and autophagy . In the Z. elegans cell culture system, the permeability and integrity of the vacuolar membrane regulates initiation of xylem-cell death , and our results suggest that VACUOLELESS1 could be involved in this regulation during the cell death of fibers.
Combining the in silico expression analysis in 19 different tissue types with a focused microarray analysis is expected to facilitate the identification of candidate genes better than either method alone. While the microarray analysis allowed selection of genes with high expression levels in xylem fibers undergoing cell death compared to early developing fibers, the in silico analysis of POPULUSDB facilitated exclusion of those genes that were highly abundant in other types of Populus tree tissues. The candidate genes selected here are potential novel regulators of cell death in xylem fibers.
Comparison of proteases in three different cell death processes
The expression of serine, cysteine and aspartic proteases was analyzed in detail in the fiber death library, because proteases have well established roles in the control of cell death in both animal cells and plants [28, 29]. We also compared their expression in the fiber death library and three other cDNA libraries: the leaf senescence library (library I) and the virus/fungal infected leaf library (library Y), which are expected to be enriched in transcripts related to cell death, and the tension wood library (library G), which should theoretically be devoid of fiber-death-related transcripts.
Expression of cysteine, serine and aspartic proteases in POPULUSDB and in a microarray analysis
Most similar Arabidopsis gene
Number of ESTs in each library
Expression ratio B/A
Cysteine endopeptidase CysEP
Papain-like cysteine proteinase RD19A
Papain-like cysteine proteinase
Vacuolar processing enzyme
Papain-like cysteine endopeptidase XCP2
Papain-like cysteine proteinase SAG12
Papain-like cysteine proteinase AALP
Papain-like cysteine proteinase RD21A
Subtilisin-like serine protease
Serine carboxypeptidase III-like protein
Subtilisin-like serine protease
Aspartic-type endopeptidase (phytepsin)
Aspartic-type endopeptidase (family A1)
Aspartic-type endopeptidase (family A1)
Aspartic-type endopeptidase (phytepsin)
In contrast to the cysteine proteases, the serine proteases displayed higher specificity to given cell death processes (Table 2). This conforms well to the reported role of serine proteases in the early signaling of apoptotic cell death in virus-infected animal cells . In plant cells, serine proteases have been implicated in pathogen-induced cell death [35, 36], and serine protease activities have also been demonstrated during xylem-cell death [5, 6]. Groover and Jones  showed that serine protease activity stimulated, and was necessary for, the death of Z. elegans tracheary elements grown in vitro. In accordance with these findings, a 40 kDa serine protease was shown to accumulate in the culture medium of in vitro tracheary elements . Our data revealed the identities of two serine proteases possibly related to the regulation of xylem-cell death. The serine proteases POPLAR.10138 and POPLAR.4995 were enriched in the fiber death library, and also highly upregulated in the microarray analysis of the xylem fibers undergoing cell death (Table 2). A PSORT analysis  of protein sequences derived from ESTs and manually assembled genomic sequences supported extracellular locations for both POPLAR.10138 and POPLAR.4995 (data not shown).
The targets of these serine proteases are not known. Groover and Jones  suggested that the extracellular 40 kDa serine protease was responsible for activation of Ca2+ channels, which is a prerequisite of xylem-cell death. Other possible targets are membrane-bound leucine-rich-repeat-containing proteins, which have been shown to interact with serine proteases during hypersensitive cell death . One such target could be a plasma-membrane localized leucine-rich-repeat-containing receptor kinase that was unique to the fiber death library (singleton X021F12) and significantly upregulated during fiber death (PU27994; Additional data file 2). However, regardless of their targets, it seems evident from our data that modification of the extracellular matrix occurs during fiber-cell death by two extracellular serine proteases, probably as part of the early signal transduction process.
An aspartic protease has been localized specifically in the tracheary elements undergoing cell death . In this analysis, we found no evidence for any fiber-death specific aspartic proteases. All four aspartic proteases identified in POPULUSDB were either present in several cell death libraries or also in the tension wood library, where cell-death-related transcripts should not accumulate (Table 2).
Hormonal control of fiber maturation
Arabidopsis has been used to analyze transcriptomes during xylem development [40, 41]. In Arabidopsis, formation of the vascular cambium, giving rise to the so-called secondary growth, takes place in the hypocotyls after two to three months of growth requiring continuous removal of the inflorescence stems . The early development and maturation of xylem vessels and fibers during the secondary growth of the hypocotyl is similar to Populus except for that the fibers in Arabidopsis do not die in the highly coordinated manner as in Populus, even after an extended growth period of three months (our unpublished work). It is therefore possible to analyze similarities between Populus and Arabidopsis in the process of fiber maturation but not cell death.
To identify common patterns in the transcriptomes of Populus and Arabidopsis during fiber maturation, we identified Arabidopsis homologs to the Populus genes that were upregulated at least two times (P < 0.001 and B > 0) in the fiber cell death sample (late maturing fibers; sample B) compared to the early fiber development (sample A). This dataset, denoted as 'Populus B/A', was compared to two previously published Arabidopsis datasets . The first Arabidopsis datasets, denoted as 'treatment', consists of genes that were upregulated at least two times during secondary growth (9 weeks of growth) compared to the primary growth of the hypocotyls (3 weeks of growth), and is therefore expected to enrich transcripts related to secondary growth including fiber maturation. The second Arabidopsis dataset, denoted as 'xylem', consists of genes that were upregulated more than two times in secondary xylem tissues compared to bark tissues of hypocotyls, and is expected to enrich transcripts related to all aspects of xylem development during secondary growth including death of the xylem vessels and fiber maturation. The common features of these datasets are shown in Additional data file 3. Cell-death-related transcripts were rarely shared by the different datasets due to the sampling method for the comparisons 'treatment' and 'xylem'. However, a large number of transcription factors and plant hormone-related transcripts were often shared by the Arabidopsis and the Populus datasets. We will discuss here the latter group of transcripts as very little previous knowledge exists on the role of plant hormones in the late maturation events of xylem fibers.
Indole-3-acetic acid (IAA) is an auxin-type plant hormone which regulates several different developmental processes in plants and also the activity of the vascular cambium and xylem formation . In Populus stems, IAA shows a steep radial concentration gradient across the stem with the highest concentration close to the vascular cambium and a very low concentration in the late-maturing xylem . Several transcriptional regulators of IAA-related gene expression, such as the auxin response factor (ARF) family proteins as well as the auxin/indoleacetic acid (IAA) family proteins, were observed in all the datasets (Additional data file 3), suggesting involvement of these proteins not only in the early development of secondary xylem, as suggested by the high concentration of IAA in these tissues, but also during late maturation of the xylem fibers. The role of IAA during late maturation of fibers is not clear, but it is possible that functional IAA signaling is required for the survival of the cells. At2g21620 encodes a universal stress protein (USP) family member, which is regulated by auxin . USPs have been implicated in the regulation of stress-related metabolic shifts, and activation of the auxin-regulated USP during secondary growth in Arabidopsis and during late fiber maturation in Populus (Additional data file 3), suggests a role for IAA in late maturing xylem fibers that experience some kind of stress. This could be osmotic stress due to condensation of the cytoplasmic contents or nutrient depletion due to the increasing distance from the nutrient-transporting phloem cells. The role of IAA during cellular stress is supported also by the fact that IAA induces expression of enzymes involved in the biosynthesis of ethylene , which is a plant hormone that is produced in response to several different stress conditions.
Ethylene does not seem to be needed for normal xylem development on the basis of the fact that ethylene-insensitive genotypes in Arabidopsis and in other species grow normally. However, ethylene biosynthesis is activated in gravitationally stimulated Populus stems, that is, when the stem is displaced from its vertical position, resulting in the production of tension wood . By analogy with several other cell-death processes in plants , it is expected that ethylene is also involved in regulation of xylem cell death. This is supported by the activation of several ethylene-related transcripts in both Arabidopsis and the late-maturing xylem fibers in Populus (Additional data file 3). ETHYLENE-INSENSITIVE 3 (EIN3), EIN3-BINDING F-BOX 1 (EBF1) and ETHYLENE RESPONSE SENSOR 1 (ERS1) mediate known parts of ethylene signal transduction . EBF1 and ERS1 are also induced by ethylene. Because both EBF1 and ERS1 function to suppress ethylene signaling, it seems that the ethylene signal required for xylem maturation needs to be suppressed or only transiently activated.
Comparison of the datasets suggests involvement of two additional plant hormones. PATHOGENESIS-RELATED PROTEIN and PHENYLALANINE AMMONIA-LYASE 2 are related to salicylic acid (SA) signaling and biosynthesis, respectively, and CORONATINE INSENSITIVE 1 to jasmonic acid (JA) signaling. Both JA and SA control stress-related processes, such as cell death and other defense responses to pathogens . The transcriptional activation of these genes in the Arabidopsis xylem tissues and in the late maturing Populus fibers suggests that both JA and SA are involved in the regulation of cell death also in xylem fibers (Additional data file 3).
Even though Arabidopsis is a suitable model system for studying early vascular development and primary growth, it cannot be easily used for studying secondary growth of the stem. Xylem fibers are formed in Arabidopsis only after two to three months of growth, which is equivalent to the time required to grow Populus trees to a size that allows collection of large amounts of woody tissues from the stem. In addition, because of the larger diameter growth of the tree stem, tissues can be easily collected from the various developmental phases without mixing the different tissue types . Therefore, Populus has several advantages over Arabidopsis for wood analysis. Development of appropriate genomic and bioinformatic tools has further strengthened Populus as the main model system for wood formation.
This analysis fully exploited the advantages of the Populus system. Combining in silico analysis of the POPULUSDB with a microarray analysis using the novel 25K Populus microarrays revealed a number of candidate genes that were unique and highly abundant in the late-developing woody tissues, and possibly related to the regulation of cell death in xylem fibers. We found fiber death-specific transcription factors, as well as nodulins and subtilisin-like serine proteases, all of which are expected to play a role in the early signaling of fiber-cell death. Sequences corresponding to the Populus homolog of Arabidopsis OPEN STOMATA 1 and peroxidases suggest involvement of ROS and ABA in fiber-cell death signaling as well. Specific expression patterns of downstream components, such as the Populus homolog of Arabidopsis XYLEM CYSTEINE PEPTIDASE 2 and the oligopeptide transporters, support the importance of these proteins in protein degradation and remobilization of nutrients in the dying fibers. Interestingly, two proteins that are known to regulate vacuolar assembly and vacuole integrity in other cellular processes, the Populus homolog of Arabidopsis VACUOLELESS1 and a VACUOLAR PROCESSING ENZYME, were specifically expressed during fiber-cell death. Permeability of the vacuolar membrane and vacuolar integrity are known to regulate death of the xylem cells, and we can now for the first time suggest candidate regulatory proteins for this process. Taken together, our analyses have identified a number of candidate genes that may be important in the regulation of fiber-cell death. Overexpression and transcriptional suppression experiments should be undertaken to elucidate the function of these genes in the process of fiber-cell death and their impact on economically important traits of woody tissues.
Materials and methods
Computational biology and database analyses
Construction of the fiber death library and the POPULUSDB  is described in . Data were analyzed using mysql, PHP, C++ and Filemaker software. The gene ontology (GO) comparison was made by listing all GO terms  in the hierarchy for each clone and selecting clones from appropriate GO hierarchy levels. Around 10% of the clones were assigned to more than one class. The EST clone distribution within POPULUSDB was performed using mysql and Filemaker, and visualized with the R software package  according to the program code described in . Annotations were derived from BLASTX analysis against the Arabidopsis proteome or the Swiss-Prot database in cases where no Arabidopsis proteins with sufficient homology were identified according to an annotation pipeline described in .
Preparation of the 25K PopuluscDNA microarray
The microarrays used here constitute the second generation of the global Populus cDNA microarrays and contain in total 24,735 different cDNA fragments. This array is based on the first generation 13K Populus array  with clones from seven cDNA libraries, representing: the cambial zone (AB), young leaves (C), floral buds (F), tension wood (G), senescing leaves (I), dormant cambium (UA), and active cambium (UB). The 25K array contains clones from 12 additional cDNA libraries, representing the apical shoot (K), cold-stressed leaves (L), roots (R), bark (N), shoot meristem (T), male catkins (V), dormant buds (Q), female catkins (M), petioles (P), fiber death (X), imbibed seeds (S) and virus/fungal infected leaves (Y). For a detailed description of the construction and sequencing of the cDNA libraries, see . All clones in the unigene set were resequenced from both the 5'-end, as were the original EST sequences, and the 3'-end. Each unigene clone is defined by a PU number on the microarray. Sequence information of the clones can be found in the POPULUSDB .
Plasmid preparations were made in a 96-plate format from bacterial cell suspensions with a Montage 96 Plasmid Prep kit (Millipore) using a Biorobot 8000 molecular biology workstation (Qiagen). PCR amplifications were done in 100 μl reaction volumes, and purified in Montage PCR384 Filter Plates (Millipore) with a Biorobot 8000 workstation (Qiagen). The purified PCR products were dissolved in 40 μl of 30% DMSO and split between a storage plate and a printing plate.
The microarrays were printed with a QArray (Genetix) instrument with 24 SMP2.5 pins (Telechem) on Ultra GAPS slides (Corning). The 24,735 cDNA fragments, together with eight copies each of the 23 different Lucidea Universal Scorecard controls (Amersham Biosciences), were spotted with a feature center-to-center distance of 180 μm. The quality of the spotted slides was assessed by staining with Syto61 (Molecular Probes) and by hybridization with random nonamers. The slides were UV cross-linked at 250 mJ/cm2 followed by baking at 75°C for 2 h.
Samples for RNA isolation were collected from the base of the stem of a 6-month-old hybrid aspen tree (Populus tremula x tremuloides) grown in a greenhouse. The bark was peeled off, and the developmental phases of the xylem were identified by light microscopy and by the texture and color of the different tissue types (Figure 1). The A sample, consisting of the remains of the cambial zone, expanding xylem and secondary cell wall depositing xylem, was collected by scraping the part of the xylem that was yellowish in color and still relatively soft with a knife. The B sample, consisting of the thick-cell-walled and late maturing xylem fibers approaching cell death, was light yellow in color and relatively difficult to scrape with the knife. The B sample was scraped to the border of the dead wood, which was discernible by its completely white color, dryness and hard texture. Total RNA was prepared from the two samples according to , and mRNA was prepared from 1 μg total RNA using paramagnetic oligo(dT) beads (Dynabeads, Dynal Biotech) in a 10 μl elution volume.
For cDNA synthesis, the mRNA sample was sheared by repeated suction into a 10 μl syringe with a beveled, non-coring needle point (Hamilton Bonaduz AG). First-strand cDNA was prepared with 1 μg oligo-dT15 primer and 200 U SuperScript II RNase H- reverse transcriptase (Invitrogen). Second-strand cDNA was prepared in a final volume of 150 μl at 16°C for 2 h with 10 U Escherichia coli DNA ligase (Invitrogen), 40 U E. coli DNA polymerase I (MBI Fermentas), 2 U RNase H (USB, GE Healthcare Bio-Sciences) in a buffer containing 0.2 mM of each dNTP, 19 mM Tris-HCl (pH 6.9), 91 mM KCl, 4.5 mM MgCl2 and 10 mM (NH4)2SO4. A further incubation was performed at 16°C for 20 min with 2 U T4 DNA polymerase (Invitrogen). The reaction was stopped with 10 μl 0.5 M EDTA, and the cDNA was recovered by phenol and chloroform:isoamylalcohol extraction followed by ethanol precipitation. The precipitates were dissolved in 2 mM Tris-HCl (pH 7.5) and purified using Autoseq G-50 columns (Amersham Biosciences). A blunt-end adapter was created by annealing two single-stranded oligonucleotides, (0.5 mmol MarAdLong: 5'-CTAATACGACTCACTATAGGGCTCGAGCGGCCGCCCGGGCAGGT-3' and 0.5 mmol MarAdShort 5'-PO4- ACCTGCCC- NH4-3') in a ligase buffer containing 66 mM Tris-HCl (pH 7.6), 5 mM MgCl2, 5 mM DTT and 7.5 μg BSA in a total volume of 15 μl with a 0.7°C/min temperature gradient from 50°C to 10°C. The adapter was ligated to the 5' template strand with 10 U T4 DNA ligase (Invitrogen) in a ligase buffer containing 33.8 mM Tris-HCl (pH 7.6), 2.56 mM MgCl2, 2.56 mM DTT and 24 μg BSA in a total volume of 48.8 μl at room temperature. The cDNA was purified with a QIAquick PCR purification kit (Qiagen).
The cDNA was amplified by PCR in a 100-μl volume containing 0.2 mM of each dNTP, 0.75 μM MaraAP1 (5'-CCATCCTAATACGACTCACTATAGGGC-3'), 0.75 μM oligo-dT15, 67 mM Tris-HCl (pH 8.8), 4 mM MgCl2, 16 mM (NH4)2SO4, 3 μg BSA and 0.6 μl AmpliTaq DNA polymerase (Applied Biosystems). The PCR procedure was 95°C for 1 min, 72°C for 5 min, addition of the AmpliTaq, and 17 to 29 cycles of 95°C for 1 min, 50°C for 1 min and 72°C for 2 min. The appropriate cycle number was defined as two cycles before saturation of the PCR product as detected by gel electrophoresis. The PCR product was purified using a QIAquick PCR purification kit and its concentration was measured by spectrophotometry (NanoDrop ND-1000, NanoDrop Technologies).
Labeling of the amplified cDNA samples was performed by direct incorporation of 3 μl Cy3-dUTP or Cy5-dUTP (Amersham Biosciences) in an asymmetric PCR reaction with 100 ng cDNA, 1 μM MaraAP1 primer, 0.6 μl AmpliTaq, 67 mM Tris-HCl (pH 8.8), 4 mM MgCl2, 16 mM (NH4)2SO4, 80 μM of each dATP, dCTP, dGTP and 20 μM dTTP in a total reaction volume of 50 μl. The PCR conditions were 95°C for 1 min followed by nine cycles of 95°C for 30 sec, 50°C for 30 sec and 72°C for 10 min. The PCR product was purified with a QIAquick PCR purification kit and eluted twice in 30 μl 4 mM KPO4 buffer (pH 8.5). The final volume was decreased to 41 μl.
Microarray hybridizations, as well as scanning of the slides and image analysis were done according to . The two samples A and B were hybridized against each other five times (including dye-swaps). The microarray raw data, including tiff and gpr files from scanning and image analysis, is deposited to the UPSC-BASE microarray database . For statistics, B-values based on Bayesian statistics  and parametric t-tests were obtained with the program R, version 1.8.1 . The microarray data shown in Additional data files 1 and 2 includes the log2 differential expression ratio (M) between the two samples B and A, P-value from the t-test and the B-value.
Samples were collected from ten different tissue types of a 6-month-old hybrid aspen (Populus tremula × tremuloides) tree grown in the greenhouse. Stem tissues; cortex, phloem, expanding xylem, secondary cell wall-depositing xylem, and fiber-cell death tissues, were collected by scraping with a knife from the base of the stem (see Figure 1). The developmental phase of the different xylem tissues was verified in transverse sections of the stem. The xylem fiber-cell death sample corresponded to the tissues collected for the B sample in the microarray analysis (see above) and the tissues prepared for construction of the fiber death cDNA library analyzed in this study (for a description of the library see ). Other tissues collected for the RT-PCR analysis were the apical shoot, root tips, young leaves, old leaves and male catkins.
RNA was prepared according to . The samples obtained were treated with RQ1 DNase (Promega), and cDNA was produced using SuperScript II RNase H- (Invitrogen) and random decamers (Ambion), followed by RNase H treatment. Two microliters of cDNA was PCR-amplified using gene-specific primers and an appropriate mixture of the universal 18S gene-specific primers and the 18S competimers (Ambion). The gene-specific primers were ACAGATGAAGCGTTCAGCAA and CAGTGCAGCACACCTAGGAA for cluster POPLAR.11628, CAGGAAAGCTCTCCGTTTCTT and TCAAAGCTCTTCCCTTCTGC for POPLAR.11658, AATCCCATGAATATTACCCCTAGA and TCTCTTGCATGGGTAGACATTTT for POPLAR.11639, ACCTCCATAGCCACCCAAG and CTGCAAGCTGATGCAGAAGT for POPLAR.11646, TGAACAGCAGGAGGTGTGAG and AACAGGTGTCCCCATCTGAG for POPLAR.11624, and TGGCAACTCCAATGAAGAAC and CACCAACAGTTTATTTATTATTCAGATG for POPLAR.9335.
Analysis of proteases in the POPULUSDB
A list of Populus proteases was compiled, based on known protease sequences of A. thaliana collected from the following databases: the Arabidopsis Information Resource , TrEMBL , Merops , InterPro , TIGR  and Munich Information Center for Protein Sequences . Sequences of 602 Arabidopsis proteases were found. The contig sequences of the POPULUSDB were then queried with BLASTX against the Arabidopsis proteases. Only the Populus sequences that obtained a BLASTX score higher than 100 and were annotated as proteases or as coding for proteins with a biological function related to proteolysis or peptidolysis in the TIGR database were accepted. The corresponding clusters, having three or more ESTs in any of the libraries of interest (tension wood, senescing leaves, fiber death and virus/fungal infected leaves) in POPULUSDB were chosen for the analysis.
Additional data files
The following additional data is available with the online version of this article. Additional data file 1 is a table containing the complete list of fiber death library-specific transcripts and the corresponding gene expression data from the microarray analysis. Additional data file 2 is a table containing the complete gene expression data in xylem fibers undergoing cell death using the 25K Populus cDNA array and the POPULUSDB. Additional data file 3 is a table containing a comparison of the microarray data from xylem tissues of Populus and Arabidopsis. Additional data file 4 is a table showing GenBank accession(s) for the EST sequences(s) in each POPULUSDB unigene cluster and singleton.
We are grateful to Andreas Sjödin for statistical analysis of the microarray results and for visualization of the data and to Edouard Pesquet for data analysis. The work was supported by the Wallenberg foundation, the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning, the Swedish Graduate School in Genomics and Bioinformatics and the Kempe foundation.
- Fukuda H: Programmed cell death of tracheary elements as a paradigm in plants. Plant Mol Biol. 2000, 44: 245-253. 10.1023/A:1026532223173.PubMedView ArticleGoogle Scholar
- Fukuda H: Xylogenesis: initiation, progression, and cell death. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 299-325. 10.1146/annurev.arplant.47.1.299.PubMedView ArticleGoogle Scholar
- Fukuda H, Komamine A: Establishment of an experimental system for the tracheary element differentiation from single cells isolated from the mesophyll of Zinnia elegans. Plant Physiol. 1980, 128: 84-94.Google Scholar
- Groover A, Jones AM: Tracheary element differentiation uses a novel mechanism coordinating programmed cell death and secondary cell wall synthesis. Plant Physiol. 1999, 119: 375-384. 10.1104/pp.119.2.375.PubMedPubMed CentralView ArticleGoogle Scholar
- Beers EP, Freeman TB: Proteinase activity during tracheary element differentiation in Zinnia mesophyll cultures. Plant Physiol. 1997, 113: 873-880.PubMedPubMed CentralGoogle Scholar
- Ye ZH, Varner JE: Induction of cysteine and serine proteases during xylogenesis in Zinnia elegans. Plant Mol Biol. 1996, 30: 1233-1246. 10.1007/BF00019555.PubMedView ArticleGoogle Scholar
- Ito J, Fukuda H: ZEN1 is a key enzyme in the degradation of nuclear DNA during programmed cell death of tracheary elements. Plant Cell. 2002, 14: 3201-3211. 10.1105/tpc.006411.PubMedPubMed CentralView ArticleGoogle Scholar
- Heath MC: Hypersensitive response-related death. Plant Mol Biol. 2000, 44: 321-334. 10.1023/A:1026592509060.PubMedView ArticleGoogle Scholar
- Thomas H, Ougham HJ, Wagstaff C, Stead AD: Defining senescence and death. J Exp Bot. 2003, 54: 1127-1132. 10.1093/jxb/erg133.PubMedView ArticleGoogle Scholar
- Bradshaw HD, Ceulemans R, Davis J, Stettler R: Emerging model systems in plant biology: Poplar (Populus) as a model forest tree. J Plant Growth Regul. 2000, 19: 306-313. 10.1007/s003440000030.View ArticleGoogle Scholar
- Tuominen H, Puech L, Fink S, Sundberg B: A radial concentration gradient of indole-3-acetic acid is related to secondary xylem development in hybrid aspen. Plant Physiol. 1997, 115: 577-585.PubMedPubMed CentralGoogle Scholar
- Hertzberg M, Aspeborg H, Schrader J, Andersson A, Erlandsson R, Blomqvist K, Bhalerao R, Uhlen M, Teeri TT, Lundeberg J, et al: A transcriptional roadmap to wood formation. Proc Natl Acad Sci USA. 2001, 98: 14732-14737. 10.1073/pnas.261293398.PubMedPubMed CentralView ArticleGoogle Scholar
- Sterky F, Bhalerao RR, Unneberg P, Segerman B, Nilsson P, Brunner AM, Charbonnel-Campaa L, Lindvall JJ, Tandre K, Strauss SH, et al: A Populus EST resource for plant functional genomics. Proc Natl Acad Sci USA. 2004, 101: 13951-13956. 10.1073/pnas.0401641101.PubMedPubMed CentralView ArticleGoogle Scholar
- POPULUSDB. [http://www.populus.db.umu.se]
- Vander Mijnsbrugge K, Meyermans H, Van Montagu M, Bauw G, Boerjan W: Wood formation in poplar: identification, characterization, and seasonal variation of xylem proteins. Planta. 2000, 210: 589-598. 10.1007/s004250050048.PubMedView ArticleGoogle Scholar
- McClung CR, Hsu M, Painter JE, Gagne JM, Karlsberg SD, Salome PA: Integrated temporal regulation of the photorespiratory pathway. Circadian regulation of two Arabidopsis genes encoding serine hydroxymethyltransferase. Plant Physiol. 2000, 123: 381-392. 10.1104/pp.123.1.381.PubMedPubMed CentralView ArticleGoogle Scholar
- Allona I, Quinn M, Shoop E, Swope K, St Cyr S, Carlis J, Riedl J, Retzel E, Campbell MM, Sederoff R, et al: Analysis of xylem formation in pine by cDNA sequencing. Proc Natl Acad Sci USA. 1998, 95: 9693-9698. 10.1073/pnas.95.16.9693.PubMedPubMed CentralView ArticleGoogle Scholar
- Motose H, Sugiyama M, Fukuda H: A proteoglycan mediates inductive interaction during plant vascular development. Nature. 2004, 429: 873-878. 10.1038/nature02613.PubMedView ArticleGoogle Scholar
- Majewska-Sawka A, Nothnagel EA: The multiple roles of arabinogalactan proteins in plant development. Plant Physiol. 2000, 122: 3-10. 10.1104/pp.122.1.3.PubMedPubMed CentralView ArticleGoogle Scholar
- Shi H, Kim Y, Guo Y, Stevenson B, Zhu JK: The Arabidopsis SOS5 locus encodes a putative cell surface adhesion protein and is required for normal cell expansion. Plant Cell. 2003, 15: 19-32. 10.1105/tpc.007872.PubMedPubMed CentralView ArticleGoogle Scholar
- Schindler T, Bergfeld R, Schopfer P: Arabinogalactan proteins in maize coleoptiles: developmental relationship to cell death during xylem differentiation but not to extension growth. Plant J. 1995, 7: 25-36. 10.1046/j.1365-313X.1995.07010025.x.PubMedView ArticleGoogle Scholar
- Gao M, Showalter AM: Yariv reagent treatment induces programmed cell death in Arabidopsis cell cultures and implicates arabinogalactan protein involvement. Plant J. 1999, 19: 321-331. 10.1046/j.1365-313X.1999.00544.x.PubMedView ArticleGoogle Scholar
- Brodersen P, Petersen M, Pike HM, Olszak B, Skov S, Odum N, Jorgensen LB, Brown RE, Mundy J: Knockout of Arabidopsis accelerated-cell-death11 encoding a sphingosine transfer protein causes activation of programmed cell death and defense. Genes Dev. 2002, 16: 490-502. 10.1101/gad.218202.PubMedPubMed CentralView ArticleGoogle Scholar
- Shaw SL, Long SR: Nod factor inhibition of reactive oxygen efflux in a host legume. Plant Physiol. 2003, 132: 2196-2204. 10.1104/pp.103.021113.PubMedPubMed CentralView ArticleGoogle Scholar
- Mohammad A, Miranda-Rios J, Estrada Navarrete G, Quinto C, Olivares JE, Garcia-Ponce B, Sanchez F: Nodulin 22 from Phaseolus vulgaris protects Escherichia coli cells from oxidative stress. Planta. 2004, 219: 993-1002. 10.1007/s00425-004-1303-9.PubMedView ArticleGoogle Scholar
- Mustilli AC, Merlot S, Vavasseur A, Fenzi F, Giraudat J: Arabidopsis OST1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell. 2002, 14: 3089-3099. 10.1105/tpc.007906.PubMedPubMed CentralView ArticleGoogle Scholar
- Rojo E, Gillmor CS, Kovaleva V, Somerville CR, Raikhel NV: VACUOLELESS1 is an essential gene required for vacuole formation and morphogenesis in Arabidopsis. Dev Cell. 2001, 1: 303-310. 10.1016/S1534-5807(01)00024-7.PubMedView ArticleGoogle Scholar
- Beers EP, Woffenden BJ, Zhao C: Plant proteolytic enzymes: possible roles during programmed cell death. Plant Mol Biol. 2000, 44: 399-415. 10.1023/A:1026556928624.PubMedView ArticleGoogle Scholar
- Danial NN, Korsmeyer SJ: Cell death: critical control points. Cell. 2004, 116: 205-19. 10.1016/S0092-8674(04)00046-7.PubMedView ArticleGoogle Scholar
- Minami A, Fukuda H: Transient and specific expression of a cysteine endopeptidase associated with autolysis during differentiation of Zinnia mesophyll cells into tracheary elements. Plant Cell Physiol. 1995, 36: 1599-1606.PubMedGoogle Scholar
- Woffenden BJ, Freeman TB, Beers EP: Proteasome inhibitors prevent tracheary element differentiation in zinnia mesophyll cell cultures. Plant Physiol. 1998, 118: 419-430. 10.1104/pp.118.2.419.PubMedPubMed CentralView ArticleGoogle Scholar
- Funk V, Kositsup B, Zhao C, Beers EP: The Arabidopsis xylem peptidase XCP1 is a tracheary element vacuolar protein that may be a papain ortholog. Plant Physiol. 2002, 128: 84-94. 10.1104/pp.128.1.84.PubMedPubMed CentralView ArticleGoogle Scholar
- Kinoshita T, Yamada K, Hiraiwa N, Kondo M, Nishimura M, Hara-Nishimura I: Vacuolar processing enzyme is up-regulated in the lytic vacuoles of vegetative tissues during senescence and under various stressed conditions. Plant J. 1999, 19: 43-53. 10.1046/j.1365-313X.1999.00497.x.PubMedView ArticleGoogle Scholar
- Hatsugai N, Kuroyanagi M, Yamada K, Meshi T, Tsuda S, Kondo M, Nishimura M, Hara-Nishimura I: A plant vacuolar protease, VPE, mediates virus-induced hypersensitive cell death. Science. 2004, 305: 855-8. 10.1126/science.1099859.PubMedView ArticleGoogle Scholar
- Jorda L, Coego A, Conejero V, Vera P: A genomic cluster containing four differentially regulated subtilisin-like processing protease genes is in tomato plants. J Biol Chem. 1999, 274: 2360-2365. 10.1074/jbc.274.4.2360.PubMedView ArticleGoogle Scholar
- Coffeen WC, Wolpert TJ: Purification and characterization of serine proteases that exhibit caspase-like activity and are associated with programmed cell death in Avena sativa. Plant Cell. 2004, 16: 857-873. 10.1105/tpc.017947.PubMedPubMed CentralView ArticleGoogle Scholar
- PSORT. [http://psort.nibb.ac.jp]
- Tornero P, Mayda E, Gomez MD, Canas L, Conejero V, Vera P: Characterization of LRP, a leucine-rich repeat (LRR) protein from tomato plants that is processed during pathogenesis. Plant J. 1996, 10: 315-330. 10.1046/j.1365-313X.1996.10020315.x.PubMedView ArticleGoogle Scholar
- Runeberg-Roos P, Saarma M: Phytepsin, a barley vacuolar aspartic proteinase, is highly expressed during autolysis of developing tracheary elements and sieve cells. Plant J. 1998, 15: 139-145. 10.1046/j.1365-313X.1998.00187.x.PubMedView ArticleGoogle Scholar
- Ko JH, Han KH, Park S, Yang J: Plant body weight-induced secondary growth in Arabidopsis and its transcription phenotype revealed by whole-transcriptome profiling. Plant Physiol. 2004, 135: 1069-1083. 10.1104/pp.104.038844.PubMedPubMed CentralView ArticleGoogle Scholar
- Oh S, Park S, Han KH: Transcriptional regulation of secondary growth in Arabidopsis thaliana. J Exp Bot. 2003, 54: 2709-2722. 10.1093/jxb/erg304.PubMedView ArticleGoogle Scholar
- Mellerowicz EJ, Baucher M, Sundberg B, Boerjan W: Unravelling cell wall formation in the woody dicot stem. Plant Mol Biol. 2001, 47: 239-274. 10.1023/A:1010699919325.PubMedView ArticleGoogle Scholar
- Uggla C, Moritz T, Sandberg G, Sundberg B: Auxin as a positional signal in pattern formation in plants. Proc Natl Acad Sci USA. 1996, 93: 9282-9286. 10.1073/pnas.93.17.9282.PubMedPubMed CentralView ArticleGoogle Scholar
- Kerk D, Bulgrien J, Smith DW, Gribskov M: Arabidopsis proteins containing similarity to the universal stress protein domain of bacteria. Plant Physiol. 2003, 131: 1209-1219. 10.1104/pp.102.016006.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsuchisaka A, Theologis A: Unique and overlapping expression patterns among the Arabidopsis 1-amino-cyclopropane-1-carboxylate synthase gene family members. Plant Physiol. 2004, 136: 2982-3000. 10.1104/pp.104.049999.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersson-Gunneras S, Hellgren JM, Bjorklund S, Regan S, Moritz T, Sundberg B: Asymmetric expression of a poplar ACC oxidase controls ethylene production during gravitational induction of tension wood. Plant J. 2003, 34: 339-349. 10.1046/j.1365-313X.2003.01727.x.PubMedView ArticleGoogle Scholar
- Overmyer K, Brosche M, Kangasjarvi J: Reactive oxygen species and hormonal control of cell death. Trends Plant Sci. 2003, 8: 335-342. 10.1016/S1360-1385(03)00135-3.PubMedView ArticleGoogle Scholar
- Guo H, Ecker JR: The ethylene signaling pathway: new insights. Curr Opin Plant Biol. 2004, 7: 40-49. 10.1016/j.pbi.2003.11.011.PubMedView ArticleGoogle Scholar
- Ihaka R, Gentleman R: A language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 299-314.Google Scholar
- Bhalerao R, Keskitalo J, Sterky F, Erlandsson R, Bjorkbacka H, Birve SJ, Karlsson J, Gardestrom P, Gustafsson P, Lundeberg J, et al: Gene expression in autumn leaves. Plant Physiol. 2003, 131: 430-442. 10.1104/pp.012732.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersson A, Keskitalo J, Sjodin A, Bhalerao R, Sterky F, Wissel K, Tandre K, Aspeborg H, Moyle R, Ohmiya Y, et al: A transcriptional timetable of autumn senescence. Genome Biol. 2004, 5: R24-10.1186/gb-2004-5-4-r24.PubMedPubMed CentralView ArticleGoogle Scholar
- Chang S, Puryear J, Cairney C: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993, 11: 113-116.View ArticleGoogle Scholar
- Smith CM, Rodriguez-Buye M, Karlsson J, Campbell MM: The response of the poplar transcriptome to wounding and subsequent infection by a viral pathogen. New Phytol. 2004, 164: 123-136. 10.1111/j.1469-8137.2004.01151.x.View ArticleGoogle Scholar
- UPSC-BASE 1.2.10+. [http://www.upscbase.db.umu.se]
- Lonnstedt I, Speed T: Replicated microarray data. Stat Sinica. 2002, 12: 31-46.Google Scholar
- TAIR. [http://www.arabidopsis.org]
- ExPASy - Swiss-Prot and TrEMBL. [http://us.expasy.org/sprot]
- MEROPS - The peptidase database. [http://merops.sanger.ac.uk]
- InterPro. [http://www.ebi.ac.uk/interpro]
- The Institute for Genomic Research. [http://www.tigr.org]
- Munich Information Center for Protein Sequences. [http://mips.gsf.de]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.