Genome-wide transcriptome analysis of the transition from primary to secondary stem development in Populus trichocarpa
© Dharmawardhana et al; licensee BioMed Central Ltd. 2010
Received: 24 November 2008
Accepted: 4 March 2010
Published: 4 March 2010
With its genome sequence and other experimental attributes, Populus trichocarpa has become the model species for genomic studies of wood development. Wood is derived from secondary growth of tree stems, and begins with the development of a ring of vascular cambium in the young developing stem. The terminal region of the developing shoot provides a steep developmental gradient from primary to secondary growth that facilitates identification of genes that play specialized functions during each of these phases of growth.
Using a genomic microarray representing the majority of the transcriptome, we profiled gene expression in stem segments that spanned primary to secondary growth. We found 3,016 genes that were differentially expressed during stem development (Q-value ≤ 0.05; >2-fold expression variation), and 15% of these genes encode proteins with no significant identities to known genes. We identified all gene family members putatively involved in secondary growth for carbohydrate active enzymes, tubulins, actins, actin depolymerizing factors, fasciclin-like AGPs, and vascular development-associated transcription factors. Almost 70% of expressed transcription factors were upregulated during the transition to secondary growth. The primary shoot elongation region of the stem contained specific carbohydrate active enzyme and expansin family members that are likely to function in primary cell wall synthesis and modification. Genes involved in plant defense and protective functions were also dominant in the primary growth region.
Our results describe the global patterns of gene expression that occur during the transition from primary to secondary stem growth. We were able to identify three major patterns of gene expression and over-represented gene ontology categories during stem development. The new regulatory factors and cell wall biogenesis genes that we identified provide candidate genes for further functional characterization, as well as new tools for molecular breeding and biotechnology aimed at improvement of tree growth rate, crown form, and wood quality.
With the completion of its genome sequence Populus trichocarpa, the genus Populus (poplar) has become the consensus model taxon for woody plant genomics and biotechnology. Many attributes of Populus makes it well suited for genomic research, including facile transformation of selected genotypes; abundant natural genetic diversity; ease of vegetative propagation; rapid growth; interspecific hybrids and mapping pedigrees; availability of genome scale microarrays; and abundant EST resources (reviewed in [1–3]). Although herbaceous short generation models such as Arabidopsis have powerful advantages for studying plant development, including fundamental aspects of lignocellulosic tissue development [4–6], poplar enables a detailed, integrated analysis of secondary growth and other developmental processes such as seasonal dormancy induction whose expression is not fully developed in annual herbaceous species.
The perennial stem growth habit of most tree species is characterized by secondary growth that results in a cumulative increase in girth during each growth cycle. This is achieved by cell division activity of the vascular cambium, with subsequent differentiation of secondary xylem towards the inside of the cambium and secondary phloem towards the outside. This coordinated and environmentally responsive program of meristematic differentiation and tissue patterning involves multiple developmental processes with interacting regulatory mechanisms (reviewed in [7, 8, 4]).
With the increasing availability of genomic tools, many studies have been conducted to help understand these developmental processes in more depth. Wood development has been characterized by expression profiling in aspen [9–11], Pinus[12, 13], black locust , Eucalyptus[15, 16] and spruce . These studies, and investigations on vascular cell development in herbaceous species Arabidopsis and Zinnia, have identified several classes of important transcription regulators and structural genes involved in secondary cell wall biogenesis (reviewed in [8, 18, 7, 19]). The developmental transition from primary to secondary growth that occurs along the stem has been used as a unique system to identify processes specific to secondary growth in sitka spruce  and aspen [20, 21]. The aspen studies, based on cDNA-AFLP, characterized 85 and 271 transcript derived fragments, respectively, as differentially expressed during secondary growth, and includes many regulatory and structural genes involved in secondary growth and secondary wall biogenesis. Our study extends this line of work to comprehensively cover the transcriptome of Populus during primary to secondary growth transition. We detected more than 36,000 Populus trichocarpa genes, and found that more than 3,000 were differentially expressed during stem development. We also provide an in depth look at the expression of several classes of transcription regulators and structural gene families whose functions are closely linked to stem and wood development.
Results and discussion
Poplar stem development
To determine transcriptional changes that occur during primary and secondary growth phases of stem development, we analyzed stem segments from successive internodes below the apex (plastochron indices 2, 3, 4, and 5: ), and from internode (IN) 9 from further down the stem. IN1 was not analyzed as it was too short to sample without including tissue from the nodes and axillary meristems.
Secondary Xylem Marker Genes are upregulated in IN9
Global expression profile during stem development
Total RNA derived from the 5 stem internode development series described above was used for expression profiling with a Nimblegen microarray representing the whole annotated transcriptome at the time that the array was created. Within the 5 stem internode samples we could detect the expression of 43,283 genes out of the 65,915 represented on the array (9,995 unigenes derived from aspen EST sequences: ). To identify statistically significant, differentially expressed genes within the developmental series we used the EDGE software package [26, 27]. A total of 4,279 Populus trichocarpa genes were identified as differentially expressed within the stem developmental series, when a Q value cutoff of 0.05 was employed; 87 genes were differentially expressed at a Q value cutoff of 0.01 (Additional file 2). To achieve a more stringent gene list we further filtered the 4,279 genes to remove weakly expressed genes (below 1.5 X the mean intensity of controls plus one standard error), and also removed genes with a fold-change below two. The resulting 3,016 genes were considered as differentially expressed in the manuscript, and used for further analysis and clustering. More than 15% of these genes were novel and had no significant BLAST returns ('No hit' cutoff e-10) in TAIR Arabidopsis or NCBI non-redundant sequence databases (June 2008). As the tissue sampling for this experiment was carried out at one time point during the active growth period, the set of differentially expressed genes might not reflect additional genes that have highly specific patterns of diurnal or seasonal regulation. To validate our differentially expressed gene list, we conducted qRT-PCR on three additional biological replicates using a subset of genes that were at the FDR cutoff of 0.05. Except for one out of the 9 genes tested, all showed a positive correlation between microarray and qRT-PCR expression patterns. The microarray and qRT-PCR derived fold-changes between elongation zone (IN2) and secondary growth dominant (IN9 or IN5) zones are illustrated in Additional File 3. The weakly regulated gene number one, identified as differentially expressed in the microarray dataset, did not show differential expression between IN2 and IN9 in the qRT-PCR analysis. As is commonly reported for microarray and qRT-PCR comparisons, the qRT-PCR expression fold-differences for all of the genes, except number 6 and 1, were considerably higher than that estimated from microarray data.
Based on gene expression, the IN5 sample largely clustered with IN9 (Figure 4A), rather than with the IN4 sample that is adjacent in the stem. This correlated well with the switch to stem radial growth following activation of secondary growth in IN5. Thus, most genes are likely to be involved in processes specific to vascular cambium activity and secondary growth. Highly expressed genes in these 2 clusters include lignin biosynthetic genes, cell wall carbohydrate metabolic genes, arabinogalactan proteins, cytoskeleton-associated genes, and a number of transcription factor genes (discussed below). Other very active genes during secondary growth included transporter genes (peptide and amino acid transporters, aquaporins), proteases, lipid transfer protein and metallothionin-like genes.
The upper quarter of cluster 3 shows a sub-cluster of genes that are upregulated mainly in IN5. The latter group of genes could be associated predominantly with early stages of secondary growth activity of the cambium. As evident in IN5 stem cross section (Figure 1), the secondary wall development and lignification of xylem is just beginning in this sample, whereas well developed and lignified walls appear mainly in IN9 xylem. A highly expressed xyloglucan endotrans glycosylase belonging to this group has been reported to be highly expressed in aspen tension wood and is suggested to be involved in maintaining crosslinks between G and S layers of highly cellulosic tension wood fibers . Several pectin methyl esterases and a UDP-glucose dehydrogenase that are likely involved in pectin and hemicellulose formation are also present within this group.
miRNAs have recently been recognized as critical post-transcriptional regulators of plant development and stress response. Our array contains only 47 miRNA precursors based on knowledge at the time of its creation; however, 234 miRNAs are now reported for poplar (miRBase miRNA database 2008 ). We were able to detect expression for 21 of these but only a few showed distinct regulation. miRNA164, 162a and 168 showed secondary development zone dominant expression. miRNA168 shows the strongest upregulation and is known to control Argonaute expression in Arabidopsis and poplar [30–32], which in turn is required for stem cell functioning, organ polarity, and root initiation in Arabidopsis [30, 33].
Highly expressed genes within cluster number 2 include genes associated with primary cell wall synthesis, cell wall loosening/extension, and plant defense. Putative cell elongation-related genes within this group include expansins, extensins, a cell elongation protein similar to DWARF1, and an apyrase. DWARF1  is a Ca/Calmodulin-regulated brassinosteroid biosynthetic gene required for cell elongation in Arabidopsis . Apyrases are documented to be involved in cell growth regulation and accumulate in rapidly growing tissues with high auxin concentrations . Thaumatin and osmotin-like proteins can have β-glucanase functions , and are also represented in the elongation zone dominant group. Biosynthetic genes that take part in plant secondary metabolite production are also well represented in this group and include several polyphenol oxidases, chlalcone synthases, and flavonol 3-hydroxylases. Phenolic and other defense-related genes have also been reported to be upregulated in young Sitka spruce stems . In addition to the IN9-upregulated lipid transfer protein, several others are upregulated in the elongation zone. Lipid transfer proteins are abundant lipid binding-capable proteins with diverse proposed (but mostly unconfirmed) functions . These include plant defense, signaling, cuticle biosynthesis, and cell wall loosening [40, 41]. Most of the lipid transfer proteins upregulated in developing Sitka spruce stem were also found in the primary growth region .
Gene function during early xylem maturation
To help assess the biological function of the regulated genes we identified during the early stages of cambial zone differentiation into xylem, we compared our results to previous studies of gene expression in aspen [42, 43, 9]. These studies used high resolution tangential cambial zone cryosectioning, but employed cDNA based microarrays that represented only ~30% of predicted gene models in poplar. Nonetheless, as discussed below, many genes were in common with the genes we identified.
We compared the genes showing the highest expression ratio in xylem maturation zone cryosections  to the differentially expressed gene set in our study. Of the 330 genes with the highest expression ratio (log2 ratio >2) in the initial maturation zone (MX1- ), we found 137 counterparts in our differentially expressed gene set (Figure 4A). Out of the 137 overlapping genes, 131 grouped in the IN5 and IN9 upregulated gene clusters (Figure 4A: Cluster 1 and 3), supporting their involvement in xylem maturation. The genes' closest hits and Arabidopsis homologs are listed in Additional file 4 and shows a dominant presence of genes associated with lignin biosynthesis, cell wall carbohydrate biosynthesis, and cell wall-associated proteins. Other major genes in common included transcription regulators (bZIP, MYB, ZF-C3HC4, LIM), hormone response genes (IAA, ethylene), and signaling genes (kinases, calmodulin binding), and cytoskeleton-associated genes (tubulins, microtubule-associated proteins).
Functional roles of differentially expressed genes during stem development
Cluster 3 genes also showed upregulation in both IN5 and IN9, but as opposed to cluster 1 peak expression was mainly in IN5. This contrast was visible in gene enrichment categories with primary cell wall and early stage secondary cell wall development-related categories showing dominance. The carbohydrate metabolic process and catalytic activity groups include cellulose biosynthesis, hemicelluloses/xylan biosynthetic genes as well as glycosyl group transferase activity and ubiquitin ligase activity-associated genes.
Cluster 2 genes, mainly upregulated in primary growth dominant IN2 and IN3 showed enrichment for many categories of biological processes, molecular functions and cellular components. Dominant biological process functional categories include active cell growth, development/morphogenesis/differentiation, and metabolic process-associated genes that reflect activities related to primary growth. Enrichment of nucleotide metabolism and DNA replication related categories could also reflect the high cell division activity during the initiation of cambium in the IN3 region. The secondary metabolism-enriched group includes flavonoid and anthocyanin biosynthetic genes.
Transcription factor and regulatory genes
There are 2,576 transcription factors and regulatory genes in two poplar databases (http://dptf.cbi.pku.edu.cn/ and http://plntfdb.bio.uni-potsdam.de/v2.0/), and 2,542 of these were represented on the array. We found that 1801 genes were expressed at detectable levels in at least one of the 5 internode samples. Figure 4B shows the clustering of 439 transcription factors that are differentially regulated during stem development (see additional file 5 for gene identities). A recent related study on transcriptome changes during shoot organogenesis in Populus identified a comparable number (588) of regulated transcription factors . We could recognize clusters upregulated in IN5-IN9, and in IN2-IN3. The majority of regulated transcription factors (69%) that were upregulated correlated with the onset of secondary growth.
With 297 identified members, MYB and MYB-like genes are the largest transcription factor family in poplar. Out of the 187 genes expressed in developing stems, more than 50% show secondary growth-associated expression patterns. MYBs are implicated as critical regulators of vascular differentiation and phenylpropanoid metabolism in plants, including the tree genera Picea, Pinus, Eucalyptus and Populus[45–48]. Overexpression of Pinus MYBs have resulted in ectopic lignification in tobacco . The secondary wall biosynthesis-regulating Arabidopsis MYB46 has also been demonstrated to be a direct target of the secondary wall-associated NAC domain transcription factor SND1 , further elaborating the network of interactions. The Populus trichocarpa homologue of Arabidopsis MYB46, gw1.IX.3536.1 is also uprgulated in IN5, IN9 samples (Figure 6, Additional file 7). The complete group of secondary growth-associated MYB family members is identified in Figure 6.
Similar to MYBs, nearly 50% of expressed NACs, one of the largest plant-specific transcription factor families, were upregulated during secondary growth (Figure 6). Arabidopsis NACs, VND(Vascular-related NAC-domain), and NST(NAC secondary wall thickening promoting factor) groups of genes have preferred expression in vascular tissue, and have been identified as regulators of xylem vessel element differentiation and secondary wall formation in 'secondary growth' tissue of Arabidopsis [51, 50, 52]. The ectopic expression of VND6, VND7, NST1, NST2 and NST3/SND1 have induced transdifferentiation of various cells to produce thick secondary walls [51, 53, 51]. The 7 VNDs and 3 NSTs fall into separate phylogenetic branches but are grouped in the same subfamily of NAC proteins . Although all putative Populus homologues of VND and NST are not regulated in our dataset, VND1 and VND2 homologues gw1.VII.3881.1 and gw184.108.40.206-JGI and NST3/SND1 homologue gw1.I.5485.1 are upregulated in IN5 and IN9 samples (Figure 6, Additional file 7). These genes were also specifically upregulated in xylem tissue in a Populus tissue-type dataset (results not shown), further supporting their role in vascular differentiation. Figure 6 illustrates more undescribed Populus NAC gene family members that show strong secondary growth-associated upregulation.
The homeodomain-containing superfamily of transcription factors participate in a wide variety of plant developmental processes. The leucine zipper-associated HD-Zip, wushel-related WOX, knotted-related KNOX, and zinc finger-associated ZF-HD homeodomain containing transcription factors have been associated with processes related to meristem function, organ polarity, and vascular development in several species . Ko et al.  have shown that a hybrid aspen class III HD-Zip protein, Pt-HB1 was closely associated with wood formation and that it was tightly regulated over developmental and seasonal gradients, with highest expression during the active growing season. A Populus trichocara homologue of aspen Pt-HB1, estExt_Genewise1_v1.C_660759, was within the IN5, IN9 upregulated cluster (Figure 6, Additional file 7), and we also found additional homeodomain superfamily members (HB and ZF-HD in Figure 6) that were strongly upregulated during secondary growth.
bZIP transcription factors play critical roles in regulating plant defense responses and development . In addition to being upregulated during Eucalyptus xylem development , the bZIP family has been identified as one of the transcription factor families most regulated during Arabidopsis stem development , with two specific candidates associated with xylem fiber differentiation. Poplar homologs of these two bZIP transcription factors are also represented within the many secondary growth-associated clusters of bZIP genes in Figure 6.
Although WRKY transcription factors are mainly implicated in regulating defense signaling , they have also been identified to be upregulated in Arabidopsis stem secondary growth and xylem tissue [59, 60] and in aspen tension wood . The WRKY family is highly diversified in plants. There were 21 poplar WRKY members that were upregulated during secondary growth, of a total of 104 total members in poplar. The functional study of this subset would help to expand knowledge of the diversity of WRKY developmental functions in plants.
There were 11 LIM genes associated with secondary development (Figure 6). LIM domain proteins in eukaryotes act as regulators of transcription, or in organization of the cytoskeleton. Tobacco NtWLIM1 binds F-actin and has been suggested to be involved in actin cytoskeleton stabilization . In sunflower HaPLIM1 is found in actin-enriched cells, and exhibits non-specific DNA and RNA binding activity . Tobacco NtLIM1 has been demonstrated to act as a transcription factor and is able to bind the conserved PAL-box found in promoters of many phenylpropanoid pathway genes . Down regulation of NtLIM1 in tobacco also results in reduced lignin and decreased expression of phenylpropanoid genes. In poplar some LIM genes have been reported to be upregulated in tension wood . Both the cytoskeletal component binding and transcription regulating activity of these genes can influence the lignin biosynthetic pathway, and thus influence secondary cell wall properties that affect wood quality and warrants further study into the secondary development associated LIM genes illustrated in Figure 6.
The few MADS box gene family members that were upregulated during secondary growth are illustrated in Figure 6. MADS box genes are well known for their involvement in reproductive development in many plants, and have also been shown to be important in vegetative development and dormancy in poplar [65, 66]. The seasonal expression pattern of aspen MADS-box5 (PTM5) that is also represented in the secondary growth-associated cluster of poplar genes (Figure 6) is specific to developing vascular tissue, including vascular cambium. It interacts with other MADS box proteins and an actin depolymerizing factor, to perhaps play a major in control of woody tissue development .
An ethylene response element binding factor-like gene that belongs to AP2/ERF domain transcription factor class is one of the most strongly regulated genes in secondary growth dominant cluster 1 (Figure 4B). It has also been identified in aspen  as differentially expressed and phloem localized in secondary tissue. AP2/ERF family members are known to be involved in integration of jasmonic acid and ethylene signals in plant defense , but also have members (BOLITA, BABYBOOM) that affect cell expansion, proliferation and differentiation pathways in Arabidopsis [68–71].
Carbohydrate active enzymes (CAZymes)
Carbohydrate catalyzing enzymes (CAZymes) expressed during stem development.
Major CAZyme Families
Total gene models (represented on poplar array)
Number of expressed transcripts
Glycoside hydrolases (GHs)
Glycosyl transferases (GTs)
Polysaccharide lyases (PLs)
Carbohydrate esterases (CEs)
Pectins include chemically diverse polymers that are major components of cell walls. They are important for wall plasticity and cell adhesion, and are major constituent of primary walls of poplar  that undergo substantial remodeling during development . In addition to some GHs, members of CE and PE gene families are involved in pectin de-esterification/acetylation and pectin degradation. Siedlecka et al.  have shown that aspen pectin methyl transferases are involved in intrusive and symplastic cell growth in developing wood cells. The small clusters of PL and CE gene family members that are preferentially expressed in the elongation zone (Additional file 8) are possible candidates for control of pectin restructuring during stem elongation. The secondary growth-associated small clusters of PL and CE gene families are likely to be involved in pectin remodeling during the early stages of secondary wall biosynthesis and radial and intrusive growth of cambial derivatives. This would be applicable to the CE family members that show highest expression in IN5.
Expansins and extensins
Extensins are cell wall structural proteins that take part in plant pathogen response and plant cell wall assembly, acting as scaffolds . During poplar stem development a large proportion of extensin family members are preferentially expressed in the elongation zone, with IN3 sample showing the highest expression (Figure 8B).
Phenylpropanoid and lignin-related genes
We have documented the transcriptome changes that occur during the transition from primary to secondary growth using sequential stem developmental stages. More than 3,000 genes were identified as differentially expressed, which is likely to be a conservative estimate due to the use of two biological replicates. Our results provide an extensive catalog of regulatory factors and cell wall biogenesis genes involved in stem elongation and secondary growth phase of tree development, whose potential applications include candidate genes for study of natural polymorphisms, targets for reverse genetic studies of tree development, and tools for marker-aided breeding or genetic engineering of tree growth and wood quality. The transcriptome patterns provided will also aid in making functional interpretations for the vast number of genes that have poorly known biochemical and developmental functions.
Two healthy, rapidly growing, three-year-old Populus trichocarpa trees (clone Nisqually-1) at a field site in Corvallis, Oregon (USA) were selected for the study. The sample collection was performed on a clear day between 10 am and 11 am in August 2005. RNA from each of the two trees was used for separate array hybridizations, serving as separate biological replicates. Technical replications of hybridizations were not performed due to high platform technical consistency as determined in a previous study. Following measurement of internode lengths, the internodes (excluding the nodes) from the apical bud to the base of the shoot were excised into liquid nitrogen in the field and stored at -80°C. A 2 mm long segment of each internode was fixed in FAA for microscopy. Following dehydration these tissue segments were embedded in wax for sectioning (WAX-IT Histology Services, Vancouver). De-waxed stem sections were stained with Toluidine Blue-O  and photographed.
RNA was isolated using Qiagen RNeasy kit using the manufacturers protocols, incorporating the Qiagen "on column" DNAse I treatment to remove contaminating genomic DNA. The quality of total RNA was assessed by absorption at 260 nm/280 nm, gel electrophoresis, and via the Agilent Bioanalyzer (Agilent Technologies, USA).
A poplar Nimblegen microarray http://www.nimblegen.com/ targeting 55,794 predicted genes from the Populus trichocarpa genome sequencing project, 126 mitochondrial and chloroplast gene models and 9995 unigenes derived from aspen EST sequences  were utilized for the experiment. Each gene model is represented by 3 replicated 60mer isothermal probes on the array. A total of 20 ug of quality tested total RNA from each sample was used for labeling. Array hybridization, quality control and data extraction were carried out by Nimblegen using their established microarray processing pipeline . The complete microarray dataset has been deposited at NCBI GEO database http://ncbi.nlm.nih.gov/geo/ with Accession # GES13043
Microarray data were normalized across all the arrays used in the experiment using the Bioconductor - Robust Multiarray Averaging (RMA) protocol incorporated in the NimbleScan (Nimblegen Inc.) software. Duplicate readings were averaged following normalization. Background hybridization intensity for determining expressed genes was estimated using signal intensity of negative control probes on the array (Escherichia coli DNA sequence derived Ambion 'ArrayControl' spike probes, #1-AE000151, #4-AE000185, #6-M34825, #8-AE000184). The mean of normalized fluorescence intensity value plus one standard error of the negative controls of all arrays in the experiment was defined as the expression threshold. Gene Spring 7.3 (Agilent Technologies, USA), EDGE  and Microarray Software Suite-MeV4.0 http://www.tm4.org/ software were utilized for further data analysis and graphical display of results. Normalized data were analyzed using EDGE  to generate a list of differentially expressed genes with a false discovery rate cut off of 0.05. EDGE is a cross-platform, compatible open-source software running on the R statistical software package, and is based on a novel statistical methods that includes an Optimal Discovery Procedure and a time course methodology [101, 26]. Clustering of gene expression was performed using hierarchical algorithm based on Pearson correlations and within-gene averaging. Populus trichocarpa genes homologous to Arabidopsis were identified using JGI annotation or Orthomap search tool http://orthomap.cgrb.oregonstate.edu. Functional enrichment/over-representation analysis was carried out using the network visualization program Cytoscape with GO plugins. For over representation determination, the Benjamini and Hochberg FDR-adjusted significance level cutoff was set at 0.05.
Array quality and data precision
The A260/A280 ratios of RNA samples used for hybridizations ranged from 1.7 to 2.0. The absence of contaminating genomic DNA and integrity of RNA samples were examined by Agilent 2100 Bioanalyzer. RNA Integrity Numbers (RIN) of samples used for hybridization ranged from 8.1 to 10.0. The correlation coefficients (r2) of pair wise comparisons of normalized intensity values between the biological replicates ranged from 0.90 to 0.97. Data precision was estimated by calculating the mean coefficient of variation for the two biological replicates for the first internode sample (IN2) for all expressed genes (36,284). The mean standard deviation was 0.71 (14.8%), and the mean standard error was 0.50 (10.5%). To evaluate the correspondence between the biological replicates, in additional file 10 we present scatter plots of intensity values for all differentially expressed genes (FDR 0.05, 2 fold) between sample 1 (tree 1) and sample 2 (tree 2) at each internode. The graph shows that 99% of the genes fall within the 2-fold demarcation lines. The biological accuracy and precision of data generated from our arrays are further supported by the quantitative correspondence of secondary xylem marker gene expression from previously published literature described above (Figure 3).
Real-time quantitative reverse transcription PCR (qRT-PCR) was employed to validate a select set of differentially expressed genes observed in our microarray study. Internode plant material was collected as described earlier from three individual trees different from what had been used for the microarray study. The sampling date was two days following the microarray study sampling date. Nine genes at the FDR cutoff (.049-.05) with low expression and fold-changes, and falling into the three main expression patterns (Figure 4A), were selected for verification. Elongation factor 1-beta (eugene3.00091463) was selected as the internal control gene  based on its consistency of expression level (mean SD = 11%) in the 5 stem tissue stages used in the study. RNA extraction, DNAse treatment, and quality control were performed as described above. cDNA synthesis from 1 ug of total RNA was carried out using a SuperScript III first-strand synthesis kit (Invitrogen) with oligo dT primers, as in the manufacturer's instructions. Primers were designed using Primer3 software (sequences in Additional file 11). Standard two-step RT-PCR and gel electrophoresis was used to evaluate primer pair efficacy and absence of genomic DNA contamination. Triplicate quantitative PCR assays per sample were performed using a Platinum SYBR Green Super mix (Invitrogen) with ROX reference dye according to the manufacturer's protocol. Dissociation curve analysis was performed at the end of each reaction to detect primer-dimers and multiple products. For relative quantification and comparison, we used the delta-delta-Ct method  with elongation factor 1-beta as the normalizing internal control gene. A 100% efficient amplification, resulting in a doubling of amplification product per cycle, was assumed (a requirement for the use of delta-delta-Ct method). This assumption was justified as only primer pairs showing an efficiency of 100% ± 5%, based on standard curves (pooled tissue type cDNA) were used for the experiments.
We thank the Office of Science (BER), U.S. Department of Energy Grant No. DE-FG02-04ER63788 for major support, and the industrial members of the TBGRC Research Cooperative at Oregon State University for partial support. We thank the Center for Genome Research and Biocomputing at Oregon State University for providing access to computing software.
- Brunner AM, Busov VB: Poplar genome sequence: functional genomics in an ecologically dominant plant species. Trends Plant Sci. 2004, 9: 49-56. 10.1016/j.tplants.2003.11.006.PubMedView Article
- Tuskan GA, DiFazio S: The genome of black cottonwood, Populus trichocarpa. Science. 2006, 313: 1596-1604. 10.1126/science.1128691.PubMedView Article
- Jansson S, Douglas CJ: Populus: a model system for plant biology. Ann Rev of Plant Biol. 2007, 58: 3485-3501.View Article
- Baucher M, El Jaziri M: From primary to secondary growth: origin and development of the vascular system. J Exp Bot. 2007, 58: 3485-3501. 10.1093/jxb/erm185.PubMedView Article
- Yokoyama R, Nishitani K: Identification and characterization of Arabidopsis thaliana genes involved in xylem secondary cell walls. J Plant Res. 2006, 119: 189-94. 10.1007/s10265-006-0261-7.PubMedView Article
- Yong W, Link B: Genomics of plant cell wall biogenesis. Planta. 2005, 221: 747-51. 10.1007/s00425-005-1563-z.PubMedView Article
- Groover A, Robischon M: Developmental mechanisms regulating secondary growth in woody plants. Curr, Opin in Plant Biol. 2006, 9: 55-58. 10.1016/j.pbi.2005.11.013.View Article
- Demura T, Fukuda H: Transcriptional regulation in wood formation. Trends Plant Sci. 2006, 12: 1360-1385.
- Schrader J, Nilsson J: A high-resolution transcript profile across the wood-forming meristem of poplar identifies potential regulators of cambial stem cell identity. Plant Cell. 2004, 16: 2278-92. 10.1105/tpc.104.024190.PubMed CentralPubMedView Article
- Israelsson MM, Eriksson ME: Changes in gene expression in the wood-forming tissue of transgenic hybrid aspen with increased secondary growth. Plant Mol Biol. 2003, 52: 893-903. 10.1023/A:1025097410445.PubMedView Article
- Sterky F, Regan S: Gene discovery in the wood-forming tissues of poplar: analysis of 5, 692 expressed sequence tags. Proc Natl Acad Sci USA. 95: 13330-5. 10.1073/pnas.95.22.13330.
- Paiva JA, Garces M: Molecular and phenotypic profiling from the base to the crown in maritime pine wood-forming tissue. New Phytol. 2008, 178: 283-301. 10.1111/j.1469-8137.2008.02379.x.PubMedView Article
- Nairn CJ, Lennon DM: Carbohydrate-related genes and cell wall biosynthesis in vascular tissues of loblolly pine (Pinus taeda). Tree Physiol. 2008, 28: 1099-110.PubMedView Article
- Yang J, Kamdem DP: Seasonal changes in gene expression at the sapwood-heartwood transition zone of black locust (Robinia pseudoacacia) revealed by cDNA microarray analysis. Tree Physiol. 2004, 24: 461-74.PubMedView Article
- Foucart C, Paux E: Transcript profiling of a xylem vs phloem cDNA subtractive library identifies new genes expressed during xylogenesis in Eucalyptus. New Phytol. 2006, 170: 739-52. 10.1111/j.1469-8137.2006.01705.x.PubMedView Article
- Paux E, Carocha V: Transcript profiling of Eucalyptus xylem genes during tension wood formation. New Phytol. 2005, 167: 89-100. 10.1111/j.1469-8137.2005.01396.x.PubMedView Article
- Pavy N, Boyle B: Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses. New Phytol. 2008, 180: 766-786. 10.1111/j.1469-8137.2008.02615.x.PubMedView Article
- Ye Z, York WS: Important new players in secondary wall synthesis. Trends Plant Sci. 2006, 11: 162-164. 10.1016/j.tplants.2006.02.001.PubMedView Article
- Sieburth LE, Deyholos MK: Vascular development: the long and winding road. Curr Opin Plant Biol. 2006, 9: 48-54. 10.1016/j.pbi.2005.11.008.PubMedView Article
- van Raemdonck D, Pesquet E: Molecular changes associated with the setting up of secondary growth in aspen. J Exp Bot. 2005, 56: 2211-2227. 10.1093/jxb/eri221.PubMedView Article
- Prassinos C, Ko JH: Transcriptome profiling of vertical stem segments provides insights into the genetic regulation of secondary growth in hybrid aspen trees. Plant Cell Physiol. 2005, 46: 1213-1225. 10.1093/pcp/pci130.PubMedView Article
- Larson PR, Isebrands JG: The plastochron index as applied to developmental studies of cottonwood. Can J Forest Res. 1971, 1: 1-11. 10.1139/x71-001.View Article
- Ko JH, Prassinos : Developmental and seasonal expression of PtaHB1, a Populus gene encoding a class III HD-Zip protein, is closely associated with secondary growth and inversely correlated with the level of microRNA (miR166). New Phytol. 2006, 169: 469-78. 10.1111/j.1469-8137.2005.01623.x.PubMedView Article
- Brunner AM, Beers EP: Transcription factors in poplar growth and development. Poplar. In the series: Genomics of Industrial Crop Plants, Series. Edited by: Joshi CP, DiFazio S, Kole C. 2010, Science Publishers, Enfield, New York,
- Sterky F, Bhalerao RR: A Populus EST resource for plant functional genomics. Proc Natl Acad Sci USA. 2004, 101: 13951-6. 10.1073/pnas.0401641101.PubMed CentralPubMedView Article
- Leek JT, Monsen E: EDGE: extraction and analysis of differential gene expression. Bioinformatics. 2006, 22: 507-8. 10.1093/bioinformatics/btk005.PubMedView Article
- Perelman E, Ploner A, Calza S, Pawitan Y: Detecting differential expression in microarray data: comparison of optimal procedures. BMC Bioinformatics. 2007, 8: 1-10. 10.1186/1471-2105-8-28.View Article
- Nishikubo N, Awano T: Xyloglucan Endo-transglycosylase (XET) functions in gelatinous layers of tension wood fibers in poplar-a glimpse into the mehanism of the balancing act of trees. Plant Cell Physiol. 2007, 48: 843-855. 10.1093/pcp/pcm055.PubMedView Article
- The miRNA database. [http://microrna.sanger.ac.uk/]
- Kidner CA, Martienssen RA: The developmental role of microRNA in plants. Curr Opin Plant Biol. 2005, 8: 38-44. 10.1016/j.pbi.2004.11.008.PubMedView Article
- Lu S, Sun YH: Stress-responsive microRNAs in Populus. Plant J. 2008, 55: 131-151. 10.1111/j.1365-313X.2008.03497.x.PubMedView Article
- Lu S, Sun YH: Novel and mechanical stress-responsive MicroRNAs in Populus trichocarpa that are absent from Arabidopsis. Plant Cell. 2005, 17: 2186-203. 10.1105/tpc.105.033456.PubMed CentralPubMedView Article
- Sorin C, Bussell JD: Auxin and light control of adventitious rooting in Arabidopsis require ARGONAUTE1. Plant Cell. 2005, 17: 1343-59. 10.1105/tpc.105.031625.PubMed CentralPubMedView Article
- Klahre U, Noguchi T: The Arabidopsis DIMINUTO/DWARF1 gene encodes a protein involved in steroid synthesis. Plant Cell. 10: 1677-90. 10.1105/tpc.10.10.1677.
- Du L, Poovaiah BW: Ca2+/calmodulin is critical for brassinosteroid biosynthesis and plant growth. Nature. 2005, 437: 741-5. 10.1038/nature03973.PubMedView Article
- Wu J, Steinebrunner I: Apyrases (nucleoside triphosphate-diphosphohydrolases) play a key role in growth control in Arabidopsis. Plant Physiol. 2007, 144: 961-75. 10.1104/pp.107.097568.PubMed CentralPubMedView Article
- Grenier J, Potvin C: Some thaumatin-like proteins hydrolyse polymeric beta-1,3-glucans. Plant J. 1999, 19: 473-80. 10.1046/j.1365-313X.1999.00551.x.PubMedView Article
- Friedmann M, Ralph SG: Microarray gene expression profiling of developmental transitions in Sitka spruce (Picea sitchensis) apical shoots. J Exp Bot. 2007, 58: 593-614. 10.1093/jxb/erl246.PubMedView Article
- Yeats TH, Rose JKC: The biochemistry and biology of extracellular plant lipid-transfer proteins (LTPs). Protein Sci. 2008, 17: 191-198. 10.1110/ps.073300108.PubMed CentralPubMedView Article
- Nieuwland J, Feron R: Lipid transfer proteins enhance cell wall extension in tobacco. Plant Cell. 2005, 17: 2009-2019. 10.1105/tpc.105.032094.PubMed CentralPubMedView Article
- Cameron KD, Teece MA: Increased accumulation of cuticular wax and expression of lipid transfer protein in response to periodic drying events in leaves of tree tobacco. Plant Physiol. 2006, 140: 176-83. 10.1104/pp.105.069724.PubMed CentralPubMedView Article
- Moreau CL, Pesquet E: A unique program for cell death in xylem fibers of Populus stem. Plant J. 2009, 58: 260-274. 10.1111/j.1365-313X.2008.03777.x.View Article
- Druart N, Johansson A: Environmental and hormonal regulation of the activity-dormancy cycle in the cambial meristem involves stage-specific modulation of transcriptional and metabolic networks. Plant J. 2007, 50: 557-573. 10.1111/j.1365-313X.2007.03077.x.PubMedView Article
- Bao Y, Dharmawardhana P, Mockler T, Strauss SH: Genome scale transcriptome analysis of shoot organogenesis in Populus. BMC Plant Biol. 2009, 9: 132-10.1186/1471-2229-9-132.PubMed CentralPubMedView Article
- Bedon F, Grima-Pettenati J: Conifer R2R3-MYB transcription factors: sequence analyses and gene expression in wood-forming tissues of white spruce (Picea glauca). BMC Plant Biol. 2007, 7: 17-10.1186/1471-2229-7-17.PubMed CentralPubMedView Article
- Zhong R, Richardson EA: The MYB46 transcription factor is a direct target of SND1 and regulates secondary wall biosynthesis in Arabidopsis. Plant Cell. 2007, 19: 2776-2792. 10.1105/tpc.107.053678.PubMed CentralPubMedView Article
- Karpinska B, Karlsson M: MYB transcription factors are differentially expressed and regulated during secondary vascular tissue development in hybrid aspen. Plant Mol Biol. 2004, 56: 255-270. 10.1007/s11103-004-3354-5.PubMedView Article
- Goicoechea M, Lacombe E: EgMYB2, a new transcriptional activator from Eucalyptus xylem, regulates secondary cell wall formation and lignin biosynthesis. Plant J. 2005, 43: 553-67. 10.1111/j.1365-313X.2005.02480.x.PubMedView Article
- Patzlaff A, McInnis S: Characterisation of a pine MYB that regulates lignification. Plant J. 2003, 36: 743-54. 10.1046/j.1365-313X.2003.01916.x.PubMedView Article
- Zhong R, Richardson EA: Two NAC domain transcription factors, SND1 and NST1, function redundantly in regulation of secondary wall synthesis in fibers of Arabidopsis. Planta. 2007, 225: 1603-1611. 10.1007/s00425-007-0498-y.PubMedView Article
- Mitsuda N, Iwase A: NAC transcription factors, NST1 and NST3, are key regulators of the formation of secondary walls in woody tissues of Arabidopsis. Plant Cell. 2007, 19: 270-280. 10.1105/tpc.106.047043.PubMed CentralPubMedView Article
- Yamaguchi M, Kubo M: Vascular-related NAC-DOMAIN7 is involved in the differentiation of all types of xylem vessels in Arabidopsis roots and shoots. Plant J. 2008, 55: 652-64. 10.1111/j.1365-313X.2008.03533.x.PubMedView Article
- Zhong R, Demura T: SND1, a NAC domain transcription factor, is a key regulator of secondary wall synthesis in fibers of Arabidopsis. Plant Cell. 2006, 18: 3158-70. 10.1105/tpc.106.047399.PubMed CentralPubMedView Article
- Ariel FD, Manavella PA: The true story of the HD-Zip family. Trends Plant Sci. 2007, 12: 419-426. 10.1016/j.tplants.2007.08.003.PubMedView Article
- Jakoby M, Weisshaar B: bZIP transcription factors in Arabidopsis. Trends in Plant Sci. 2002, 7: 106-11. 10.1016/S1360-1385(01)02223-3.View Article
- Paux E, Tamasloukht M: Identification of genes preferentially expressed during wood formation in Eucalyptus. Plant Mol Biol. 2004, 55: 263-80. 10.1007/s11103-004-0621-4.PubMedView Article
- Ehlting J, Mattheus N: Global transcript profiling of primary stems from Arabidopsis thaliana identifies candidate genes for missing links in lignin biosynthesis and transcriptional regulators of fiber differentiation. Plant J. 2005, 42: 618-640. 10.1111/j.1365-313X.2005.02403.x.PubMedView Article
- Eulgem T, Somssich IE: Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007, 10: 366-71. 10.1016/j.pbi.2007.04.020.PubMedView Article
- Ko JH, Han KH: Arabidopsis whole-transcriptome profiling defines the features of coordinated regulations that occur during secondary growth. Plant Mol Biol. 2004, 55: 433-53. 10.1007/s11103-004-1051-z.PubMedView Article
- Zhao C, Craig JC: The xylem and phloem transcriptomes from secondary tissues of the Arabidopsis root-hypocotyl. Plant Physiol. 2005, 138: 803-18. 10.1104/pp.105.060202.PubMed CentralPubMedView Article
- Andersson-Gunneras S, Mellerowicz EJ: Biosynthesis of cellulose-enriched tension wood in Populus: global analysis of transcripts and metabolites identifies biochemical and developmental regulators in secondary wall biosynthesis. Plant J. 2006, 45: 144-165. 10.1111/j.1365-313X.2005.02584.x.PubMedView Article
- Thomas C, Hoffmann C: Tobacco WLIM1 is a novel F-actin binding protein involved in actin cytoskeleton remodeling. Plant Cell. 2006, 18: 2194-206. 10.1105/tpc.106.040956.PubMed CentralPubMedView Article
- Mundel C, Baltz R: A LIM-domain protein from sunflower is localized to the cytoplasm and/or nucleus in a wide variety of tissues and is associated with the phragmoplast in dividing cells. Plant Mol Biol. 2000, 42: 291-302. 10.1023/A:1006333611189.PubMedView Article
- Kawaoka A, Kaothien P: Functional analysis of tobacco LIM protein Ntlim1 involved in lignin biosynthesis. Plant J. 2000, 22: 289-301. 10.1046/j.1365-313x.2000.00737.x.PubMedView Article
- Hoenicka H, Nowitzki O: Heterologous overexpression of the birch FRUITFULL-like MADS-box gene BpMADS4 prevents normal senescence and winter dormancy in Populus tremula L. Planta. 2008, 227: 1001-11. 10.1007/s00425-007-0674-0.PubMedView Article
- Cseke LJ, Ravinder N: Identification of PTM5 protein interaction partners, a MADS-box gene involved in aspen tree vegetative development. Gene. 2007, 391: 209-22. 10.1016/j.gene.2006.12.033.PubMedView Article
- Pre M, Atallah M: The AP2/ERF domain transcription factor ORA59 integrates jasmonic acid and ethylene signals in plant defense. Plant Physiol. 2008, 147: 1347-57. 10.1104/pp.108.117523.PubMed CentralPubMedView Article
- Marsch-Martinez N, Greco R: BOLITA, an Arabidopsis AP2/ERF-like transcription factor that affects cell expansion and proliferation/differentiation pathways. Plant Mol Biol. 2006, 62: 825-43. 10.1007/s11103-006-9059-1.PubMedView Article
- Passarinho P, Ketelaar T: BABY BOOM target genes provide diverse entry points into cell proliferation and cell growth pathways. Plant Mol Biol. 2008, 68: 225-37. 10.1007/s11103-008-9364-y.PubMedView Article
- Srinivasan C, Liu Z: Heterologous expression of the BABY BOOM AP2/ERF transcription factor enhances the regeneration capacity of tobacco (Nicotiana tabacum L.). Planta. 2007, 225: 341-51. 10.1007/s00425-006-0358-1.PubMedView Article
- Boutilier K, Offringa R: Ectopic expression of BABY BOOM triggers a conversion from vegetative to embryonic growth. Plant Cell. 2002, 14: 1737-1749. 10.1105/tpc.001941.PubMed CentralPubMedView Article
- Aspeborg H, Schrader J: Carbohydrate-active enzymes involved in the secondary cell wall biogenesis in hybrid aspen. Plant Physiol. 2005, 137: 983-997. 10.1104/pp.104.055087.PubMed CentralPubMedView Article
- Geisler-Lee J, Geisler M: Poplar carbohydrate-active enzymes. Gene identification and expression analyses. Plant Physiol. 2006, 140: 946-962. 10.1104/pp.105.072652.PubMed CentralPubMedView Article
- Zhong R, Pena MJ: Arabidopsis fragile fiber8, which encodes a putative glucuronyltransferase, is essential for normal secondary wall synthesis. Plant Cell. 2005, 17 (12): 3390-408. 10.1105/tpc.105.035501.PubMed CentralPubMedView Article
- Pena MJ, Zhong R: Arabidopsis irregular xylem8 and irregular xylem9: implications for the complexity of glucuronoxylan biosynthesis. Plant Cell. 2007, 19: 549-63. 10.1105/tpc.106.049320.PubMed CentralPubMedView Article
- Zhou GK, Zhong R: Molecular characterization of PoGT8D and PoGT43B, two secondary wall-associated glycosyltransferases in poplar. Plant Cell Physiol. 2007, 48: 689-699. 10.1093/pcp/pcm037.PubMedView Article
- Mellerowicz EJ, Baucher M: Unravelling cell wall formation in the woody dicot stem. Plant Mol Biol. 2001, 47: 239-74. 10.1023/A:1010699919325.PubMedView Article
- Micheli F: Pectin methylesterases: cell wall enzymes with important roles in plant physiology. Trends Plant Sci. 2001, 6: 414-9. 10.1016/S1360-1385(01)02045-3.PubMedView Article
- Siedlecka A, Wiklund S: Pectin methyl esterase inhibits intrusive and symplastic cell growth in developing wood cells of Populus. Plant Physiol. 2008, 146: 554-65. 10.1104/pp.107.111963.PubMed CentralPubMedView Article
- Suzuki S, Li L: The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in Populus trichocarpa. Plant Physiol. 2006, 142: 1233-1245. 10.1104/pp.106.086678.PubMed CentralPubMedView Article
- Djerbi S, Lindskog M: The genome sequence of black cottonwood (Populus trichocarpa) reveals 18 conserved cellulose synthase (CesA) genes. Planta. 2005, 221: 739-746. 10.1007/s00425-005-1498-4.PubMedView Article
- Kalluri UC, Joshi CP: Differential expression patterns of two cellulose synthase genes are associated with primary and secondary cell wall development in aspen trees. Planta. 2004, 220: 47-55. 10.1007/s00425-004-1329-z.PubMedView Article
- Cosgrove DJ: Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005, 6: 850-61. 10.1038/nrm1746.PubMedView Article
- Sampedro J, Cosgrove DJ: The expansin superfamily. Genome Biol. 2005, 6: 242-10.1186/gb-2005-6-12-242.PubMed CentralPubMedView Article
- Gray-Mitsumune M, Mellerowicz EJ: Expansins abundant in secondary xylem belong to subgroup A of the alpha-expansin gene family. Plant Physiol. 2004, 135: 1552-64. 10.1104/pp.104.039321.PubMed CentralPubMedView Article
- Gray-Mitsumune M, Blomquist K: Ectopic expression of a wood-abundant expansin PttEXPA1 promotes cell expansion in primary and secondary tissues in aspen. Plant Biotechnol J. 2008, 6: 62-72.PubMed
- Cannon MC, Terneus K: Self-assembly of the plant cell wall requires an extensin scaffold. Proc Natl Acad Sci USA. 2008, 105: 2226-31. 10.1073/pnas.0711980105.PubMed CentralPubMedView Article
- Wightman R, Turner SR: The roles of the cytoskeleton during cellulose deposition at the secondary cell wall. Plant J. 2008, 54: 794-805. 10.1111/j.1365-313X.2008.03444.x.PubMedView Article
- Spokevicius AV, Southerton SG: B-tubulin affects cellulose microfibril orientation in plant secondary fibre cell walls. Plant J. 2007, 51: 717-726. 10.1111/j.1365-313X.2007.03176.x.PubMedView Article
- Oakley RV, Wang YS: Differential expansion and expression of alpha- and beta-tubulin gene families in Populus. Plant Physiol. 2007, 145: 961-73. 10.1104/pp.107.107086.PubMed CentralPubMedView Article
- Bhandari S, Fujino T: Xylem-specific and tension stress-responsive coexpression of KORRIGAN endoglucanase and three secondary wall-associated cellulose synthase genes in aspen trees. Planta. 2006, 224: 828-837. 10.1007/s00425-006-0269-1.PubMedView Article
- Gaspar Y, Johnson KL: The complex structures of arabinogalactan-proteins and the journey towards understanding function. Plant Mol Biol. 2001, 47: 161-76. 10.1023/A:1010683432529.PubMedView Article
- Zhang Y, Brown G: An arabinogalactan protein associated with secondary cell wall formation in differentiating xylem of loblolly pine. Plant Mol Biol. 2003, 52: 91-102. 10.1023/A:1023978210001.PubMedView Article
- Loopstra CA, Puryear JD: Purification and cloning of an arabinogalactan-protein from xylem of loblolly pine. Planta. 2000, 210 (4): 686-9. 10.1007/s004250050061.PubMedView Article
- Lafarguette F, Leple JC: Poplar genes encoding fasciclin-like arabinogalactan proteins are highly expressed in tension wood. New Phytol. 2004, 164: 107-121. 10.1111/j.1469-8137.2004.01175.x.View Article
- Tsai C, Harding SA: Genome-wide analysis of the structural genes regulating defense phenylpropanoid metabolism in Populus. New Phytol. 2006, 172: 47-62. 10.1111/j.1469-8137.2006.01798.x.PubMedView Article
- Peters DJ, Constabel CP: Molecular analysis of herbivore-induced condensed tannin synthesis: cloning and expression of dihydroflavonol reductase from trembling aspen (Populus tremuloides). Plant J. 2002, 32: 701-12. 10.1046/j.1365-313X.2002.01458.x.PubMedView Article
- Olsen KM, Lea US: Differential expression of four Arabidopsis PAL genes; PAL1 and PAL2 have functional specialization in abiotic environmental-triggered flavonoid synthesis. J Plant Physiol. 2008, 165: 1491-1499. 10.1016/j.jplph.2007.11.005.PubMedView Article
- Jensen WA: Botanical Histochemistry: Principles and Practice. 1962, San Francisco and London, Freeman and Company
- Wang X, He H: NMPP: a user-customized NimbleGen microarray data processing pipeline. Bioinformatics. 2006, 22 (23): 2955-7. 10.1093/bioinformatics/btl525.PubMedView Article
- Storey JD, Xiao W: Significance analysis of time course microarray experiments. Proc Natl Acad Sci USA. 2006, 102: 12837-42. 10.1073/pnas.0504609102.View Article
- Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC Plant Biology. 2004, 4: 14-10.1186/1471-2229-4-14.PubMed CentralPubMedView Article
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView Article
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.