Transcriptional programming during cell wall maturation in the expanding Arabidopsis stem
© Hall and Ellis; licensee BioMed Central Ltd. 2013
Received: 20 November 2012
Accepted: 21 January 2013
Published: 25 January 2013
Plant cell walls are complex dynamic structures that play a vital role in coordinating the directional growth of plant tissues. The rapid elongation of the inflorescence stem in the model plant Arabidopsis thaliana is accompanied by radical changes in cell wall structure and chemistry, but analysis of the underlying mechanisms and identification of the genes that are involved has been hampered by difficulties in accurately sampling discrete developmental states along the developing stem.
By creating stem growth kinematic profiles for individual expanding Arabidopsis stems we have been able to harvest and pool developmentally-matched tissue samples, and to use these for comparative analysis of global transcript profiles at four distinct phases of stem growth: the period of elongation rate increase, the point of maximum growth rate, the point of stem growth cessation and the fully matured stem. The resulting profiles identify numerous genes whose expression is affected as the stem tissues pass through these defined growth transitions, including both novel loci and genes identified in earlier studies. Of particular note is the preponderance of highly active genes associated with secondary cell wall deposition in the region of stem growth cessation, and of genes associated with defence and stress responses in the fully mature stem.
The use of growth kinematic profiling to create tissue samples that are accurately positioned along the expansion growth continuum of Arabidopsis inflorescence stems establishes a new standard for transcript profiling analyses of such tissues. The resulting expression profiles identify a substantial number of genes whose expression is correlated for the first time with rapid cell wall extension and subsequent fortification, and thus provide an important new resource for plant biologists interested in gene discovery related to plant biomass accumulation.
KeywordsCell wall Anisotropy Growth kinematic profiling Transcriptome Microarray Arabidopsis Inflorescence stem
Directional (anisotropic) cell wall expansion is an integral part of most plant developmental processes, where it facilitates the structural changes necessary for proper cell and organ morphogenesis. The initial expansive growth phase, which requires both addition of new extracellular polymers and remodeling of existing components in the primary cell walls, is often succeeded by cell wall thickening and rigidification processes to create secondary cell walls that enhance the structural integrity of the organ, but also curtail further wall extension. These sequential processes require a high degree of dynamic, context-specific coordination of cell wall building, reconstruction and fortification in order to harness the underlying driving force of turgor pressure in a spatially-defined manner (reviewed in [1, 2]). Consistent with such developmental complexity, at least one thousand genes in Arabidopsis have been shown to have some association with cell wall synthesis and remodeling .
The gene expression patterns associated with cell wall expansion and/or secondary cell wall formation in Arabidopsis have been analyzed in several studies in efforts to identify participating genes and understand the biological roles of their products [4–8]. To specifically address cell expansion processes, for example, transcript profiles have been compared with protein accumulation profiles in Arabidopsis seedling hypocotyls that were undergoing rapid cell elongation without significant cell division . Transcript profiling has also been conducted with in vitro cultured Zinnia mesophyll cells [10, 11] and Arabidopsis subcultured cells  that were induced to trans-differentiate into tracheary element-like cells, a process that is accompanied by deposition of distinctive secondary cell wall thickenings on top of the original primary cell walls. In another approach, large-scale correlation analysis of public microarray data enabled the in silico identification of genes whose expression is strongly aligned with expression of specific members of the Arabidopsis cellulose synthase (CesA) gene family that are believed to be predominantly involved in either primary or secondary cell wall biogenesis .
The Arabidopsis inflorescence stem represents an attractive experimental system for such gene discovery studies since it provides more substantial amounts of tissue for analyses, and its tissue architecture is largely established prior to bolting, which means that stem expansion is primarily the product of cell elongation, rather than division. However, the growing stem also presents a continuum of developmental states along the organ as its component cells transition from early anisotropic expansion growth to growth cessation, and finally to cell wall fortification. Integrated with these changes in cell expansion activity are additional changes associated with the differentiation and maturation of the discrete tissues that comprise the stem.
In order to accurately monitor the gene expression changes, or transcriptional programming, that accompany these various stem growth transitions, it is essential to sample and pool stem tissues that are verifiably associated with specific stages of development. A number of inflorescence stem profiling studies have attempted to compare the global transcriptional changes occurring between specific developmental stages [4, 5, 7, 14], but the experimental strategies employed have typically compared tissues from visually selected regions of multiple plants, and have operated under two untested assumptions: 1) that the pooled plants whose stems are being sampled all have similar developmental proportioning, and 2) that the sampling guidelines for the harvested plants, derived from destructive analysis of a different set of plants, accurately associate features such as appearance of lignification in the interfascicular fibres of the stem  with specific developmental growth stages. Contrary to these assumptions, growth kinematic profiling of expanding inflorescence stems of individual Arabidopsis plants has recently demonstrated that stem growth profiles actually vary widely from plant to plant, even within genetically homogeneous populations . As a consequence, data obtained from indirectly selected and harvested stem regions are likely to be relatively poorly correlated with onset of processes involved in cell wall extension or modification events associated with specific growth stages.
In this study, we have applied growth kinematic profiling (GKP) to a series of individual Arabidopsis inflorescence stems as described in Hall and Ellis (2012) , and used the resulting growth rate profiles to generate pooled samples of stem tissue that accurately represent four discrete growth stages along the cell wall expansion developmental continuum. Our use of GKP-based sampling was expected to reduce the biological noise associated with indirect sampling strategies used in previous studies, and thus increase the sensitivity (power to detect actual differential expression) and accuracy of transcript profiling. Microarray-based assessment of global transcript abundance in these GKP-matched stem samples then enabled us to generate transcriptomic datasets that can be positioned with confidence within a validated developmental context of cell wall expansion performance (relative elemental growth rates). The resulting gene expression profiles demonstrate the participation of many genes that had earlier been linked to primary or secondary cell wall synthesis, but they also highlight expression changes in a range of unique genes whose role(s) in cell wall maturation or stem expansion have yet to be assessed.
Stage-specific transcriptome analysis
To maximize both statistical power and flexibility in analysis, our experiment directly co-hybridized all pair-wise combinations of developmental stage samples (‘complete factorial’ experimental design) and utilized a ‘mixed effects model’ analysis  to compare the four developmental stages on the basis of six biological replicates, each pooled from a common set of thirty-four randomly-assigned plants.
Examination of differential expression between stages
The goal of this study was to identify genes whose expression in the inflorescence stem differs most strongly between different growth stages, since these are expected to represent the loci most actively involved in the accompanying transcriptional reprogramming. The mixed effects model-based analysis of this experiment generates six possible pair-wise comparisons between the four stages, for which the complete statistical analysis is presented in Additional files 2, 3 and 4: Table S1.
Twenty most differentially expressed genes with higher expression in YNG stage relative to MGR stage
elongation factor 1-alpha / EF-1-alpha
spliceosome associated protein-related
elongation factor 1-alpha / EF-1-alpha
ARABIDOPSIS THALIANA HOMEOBOX PROTEIN 33 (AtHB33)
MATERNAL EFFECT EMBRYO ARREST 48 (MEE48)
aspartyl protease family protein located in membrane, plant-type cell wall
chaperonin, putative with domain Cpn60/TCP-1 (InterPro:IPR002423)
chalcone and stilbene synthase family protein involved in phenylpropanoid biosynthetic process
RECEPTOR-LIKE KINASE 902 (RLK902)
myrosinase binding protein, putative
GDSL-MOTIF LIPASE 7 (GLIP7)
EMBRYO DEFECTIVE 2296 (EMB2296)
protease inhibitor/seed storage/lipid transfer protein (LTP) family protein
CLAVATA3/ESR-RELATED 16 (CLE16)
LIPOXYGENASE 2 (LOX2)
leucine-rich repeat transmembrane protein kinase
ARABIDOPSIS NAC DOMAIN CONTAINING PROTEIN 63 (anac063)
Twenty most differentially expressed genes with higher expression in MGR stage relative to YNG stage
ARABINOGALACTAN PROTEIN 12 (AGP12)
putative polygalacturonase (pectinase)
MYB DOMAIN PROTEIN 61 (MYB61)
RARE-COLD-INDUCIBLE 2A (RCI2A)
similar to 1-aminocyclopropane-1-carboxylate oxidase
Auxin responsive SAUR protein
SOUL heme-binding family protein
ARABINOGALACTAN PROTEIN 13 (AGP13)
BRANCHED-CHAIN AMINOTRANSFERASE4 (BCAT4)
basic helix-loop-helix (bHLH) family protein
Overexpression leads to PEL (Pseudo-Etiolation in Light) phenotype
putative gibberellin-responsive protein (GASA6)
MINI ZINC FINGER 1 (MIF1)
protease inhibitor/seed storage/lipid transfer protein (LTP) family protein
peroxidase 64 (PER64) located in plant-type cell wall
Twenty most differentially expressed genes with higher expression in MGR stage relative to CSS stage
MLP-LIKE PROTEIN 423 (MLP423)
PHLOEM PROTEIN 2-B7 (AtPP2-B7)
LIPID TRANSFER PROTEIN 1 (LP1)
VEGETATIVE STORAGE PROTEIN 1 (VSP1)
SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 3 (SPL3)
LI-TOLERANT LIPASE 1 (LTL1)
chloroplast 60 kDa chaperonin beta subunit
GERMIN 3 (GER3)
radical SAM domain-containing protein
bundle-sheath defective protein 2 family / bsd2 family
GAST1 PROTEIN HOMOLOG 4 (GASA4)
GLUTAMINE-DEPENDENT ASPARAGINE SYNTHETASE (ASN1)
CHAPERONIN 20 (CPN20)
protease inhibitor/lipid transfer protein (LTP) family protein
40S ribosomal protein S19 (RPS19C)
elongation factor P (EF-P) family protein
F-box family protein (FBW1)
LOW-MW CYSTEINE-RICH 68 (LCR68)(PDF2.3)
Twenty most differentially expressed genes with higher expression in CSS stage relative to MGR stage
CBL-INTERACTING PROTEIN KINASE 25 (CIPK25)(SnRK3.25)
XYLOGLUCAN ENDOTRANSGLUCOSYLASE (XTH19)
UDP-glucuronic acid decarboxylase 3 (UXS3)
LACCASE 4 (IRX12)
IRREGULAR XYLEM 9 (IRX9)
REDUCED WALL ACETYLATION 1 (RWA1)
S/T protein kinase
TRICHOME BIREFRINGENCE-LIKE 3 (TBL3)
Nucleotide/sugar transporter family protein
plastocyanin-like domain-containing protein
glucuronic acid substitution of xylan1 (GUX1)
CHITINASE-LIKE PROTEIN 2 (CTL2)
FRAGILE FIBER 8 (FRA8)
aspartyl protease family protein
Twenty most differentially expressed genes with higher expression in CSS stage relative to OLD stage
GERMIN 3 (GER3)
pectate lyase family protein
GERMIN-LIKE PROTEIN 1 (GER1)
ARABIDOPSIS THALIANA METHIONINE GAMMA-LYASE (ATMGL)
glycoside hydrolase family 28 protein / polygalacturonase (pectinase) family protein
hydrolase, alpha/beta-fold family protein
ribulose bisphosphate carboxylase small chain 1B / RuBisCO small subunit 1B (RBCS-1B) (ATS1B)
PLASMA MEMBRANE INTRINSIC PROTEIN 2E (PIP2E)
LIPID TRANSFER PROTEIN 1 (LTP1)
DELTA TONOPLAST INTEGRAL PROTEIN (DELTA-TIP)
coproporphyrinogen oxidase activity in porphyrin biosynthetic process within chloroplast
copper-binding family protein in metal ion transport
family II extracellular lipase 3 (EXL3), carboxylesterase activity, acyltransferase activity
glycosyl hydrolase family 17 protein
ribulose bisphosphate carboxylase small chain 2B / RuBisCO small subunit 2B (RBCS-2B) (ATS2B)
ALPHA-XYLOSIDASE 1 (XYL1)
DNA-DAMAGE REPAIR/TOLERATION 100 (DRT100)
Twenty most differentially expressed genes with higher expression in OLD stage relative to CSS stage
EXTENSIN 3 (ATEXT3)
proline-rich extensin-like family protein, structural constituent of cell wall
MYB DOMAIN PROTEIN 49 (MYB49)
Bet v I allergen family
GLUTATHIONE S-TRANSFERASE F3 (ATGSTF3)
acyl-transferase family protein
GLUTATHIONE S-TRANSFERASE 6 (GSTF6)
DEFECTIVELY ORGANIZED TRIBUTARIES 1 (DOT1)
LOW-MOLECULAR-WEIGHT CYSTEINE-RICH 67 (LCR67)
PLANT DEFENSIN 1.2B (PDF1.2b)
Wax ester synthase homologue
PLANT DEFENSIN 1.2 (PDF1.2)
PLANT DEFENSIN 1.3 (PDF1.3)
PLANT DEFENSIN 1.2C (PDF1.2c)
FATTY ACID REDUCTASE 6 (FAR6)
To establish the status of the stem transcriptome prior to the period at which the maximal elongation rate has been achieved, gene expression profiling was conducted on the top 1 cm segment of the stem (with flowers removed), and this YNG transcript profile was then compared to the profile generated from segments representing the GKP-identified maximum growth rate (MGR) phase. The forty genes exhibiting either high transcript copy number in the YNG stage sample relative to the MGR stage (positive fold-change values), or in the MGR stage sample relative to the YNG stage (negative fold-change values) are listed in Tables 1 and 2, respectively Among the most up-regulated genes in the YNG-dominant set are two members of the four-member ELONGATION FACTOR 1-ALPHA gene family  (EF-1-α A2 (At1g07930) and EF-1-α A3 (At1g07940)). This association with active protein synthesis is also reflected in the GO term enrichment analysis of this set (Additional file 5: Table S2), which indicates >80-fold enrichment in GO terms containing ‘translation’. Also found within this ‘YNG-up-regulated’ list are genes related to signaling (RLK902; CLE16(CLAVATA3 homologue); LOX2; At1g62950, a LRR protein kinase), as well as transcription factors (ZF-HD class AtHB33; NAC063). Cell-cell communication mediated by peptides derived from CLE gene products, acting together with cognate receptor kinases, represents part of the elaborate signaling network that helps guide plant development . While known cell wall-associated genes are not notably over-represented within the ‘YNG-up-regulated’ list, one gene encoding a putative glucan endo-β(1→3)-glucosidase (At4g14080)  is up-regulated 15-fold over the MGR stage.
Although cell wall expansion is expected to be taking place in both the YNG and the MGR stage tissues, genes whose products are uniquely required for rapid expansion should be relatively more highly expressed in the latter. The most strongly differentially up-regulated (~40-fold) gene in the MGR tissues relative to YNG is a peroxidase (PER64) that has been previously reported to be up-regulated in stems in response to mechanical load . The peroxidase gene family in Arabidopsis is large, and its members play a number of roles in cellular metabolism, including modulation of reactive oxygen species accumulation  and the oxidative coupling of aromatic metabolites such as the monolignols that serve as precursors for the lignin polymer [26, 27]. The expression of PER64 in Arabidopsis has been shown to be concentrated in the protoxylem , where lignification of patterned secondary cell wall thickenings contributes to cell wall stabilization during vascular elongation, a spatial specificity that is consistent with the strong PER64 expression in MGR tissues.
The MGR up-regulated list also contains several genes more directly related to primary cell wall formation and re-modeling, including a xyloglucan endotransglycosylase/hydrolase MERISTEM-5 (MERI5B/XTH24), a putative pectinase (At1g80170), two arabinogalactan proteins (AGP12, AGP13) and a MYB transcription factor (MYB61) that has recently been shown to contribute to both cell wall synthesis and regulation of plant carbon allocation [29–32]. In addition, several genes encoding proteins associated with phytohormone signalling are more highly expressed in the MGR tissues, including a putative ACC oxidase, the GA-responsive MINI ZINC FINGER 1 (MIF1)  and another gene GASA6 (At1g74670) reported to be GA-responsive .
Comparison of the MGR stage gene expression patterns to those observed at the more mature CSS stage provides another view of those genes that are most relevant to active stem expansion, by contrasting their performance in the rapidly expanding MGR tissues with that seen in the CSS tissues where cell wall expansion has ceased. Interestingly, the list of twenty genes whose expression is ‘Higher in MGR relative to CSS’ (Table 3) is led, not by genes known to be associated with cell wall synthesis or modification, but by MAJOR LATEX PROTEIN 423 (MLP423), a member of the BET V1 class of allergens that exhibits sequence homology to ABA- and stress-responsive proteins from various plant species (EMBL-EBI database information). MLP423 is accompanied by two members of the large (108-member) GDSL-type lipase homologue gene family, and by other genes associated with lipid metabolism/transport, but few, if any, genes known to be directly involved in cell wall synthesis are included. This profile implies that the genes populating the ‘Higher in MGR relative to CSS’ list are primarily those whose expression is relatively strongly reduced as the cells make their transition from rapid anisotropic expansion to maturation.
The ‘Higher in CSS relative to MGR’ gene list (Table 4), on the other hand, would be expected to capture those genes that make a major contribution to the re-programming associated with transition to a phase of cell wall stabilization and rigidification. Consistent with this prediction, this list is dominated by genes associated with formation of non-expanding walls, including all three of the cellulose synthase genes believed to be involved in cellulose microfibril deposition during secondary cell wall biosynthesis (CESA4/IRX5, CESA7/IRX3, and CESA8/IRX1) [14, 35–37], and CHITINASE-LIKE PROTEIN 2 (CTL2)  whose loss-of-function mutant displays cellulose biosynthesis defects . Also strongly represented are genes required for xylan biosynthesis/modification, including a UDP-GLUCURONIC ACID DECARBOXYLASE 3/UXS3 that provides UDP-xylose for xylan backbone synthesis, IRREGULAR XYLEM 9 (IRX9)  and FRAGILE FIBER 8 (FRA8)  whose encoded proteins build and extend the glucuronosylxylan polymer, and two xylan modification genes: a xyloglucan-specific endotransglycosylase/hydrolase 19 (XTH19) , and REDUCED WALL ACETYLATION 1 (RWA1). Other cell wall modification genes are present, including two pectinesterases (At2g43050, At2g45220), one of which is the most strongly differentially-expressed gene in the list. The prominence of these pectin de-methylating enzymes in the MGR→CSS transition list is consistent with a current model for plant cell wall rigidification in which a reduction in the levels of pectin methylesterification leads to enhanced calcium ion cross-linking and wall stiffening [45–48].
In addition to genes whose encoded products affect cell wall polysaccharide biosynthesis, the list includes IRREGULAR XYLEM 12 (IRX12/LAC4). Laccases are thought to contribute to polymerization of lignin in secondary walls, and LAC4 expression has previously been shown to be specific to xylary and interfascicular fibres in the Arabidopsis stem. Lignin deposition is largely unaffected in the lac4 loss-of-function mutant, but is strongly reduced in the lac4/lac17 double loss-of-function mutant . It is noteworthy that we observed no significant difference in expression of LAC17 between the CSS and MGR stages (1.3-fold differential, CSS/MGR). Overall, nine of the twenty genes featured in this list also occur among a set of ‘xylem-specific’ Arabidopsis genes identified through analysis of public datasets , consistent with a metabolic commitment in CSS tissues to cell wall rigidification in xylem fibres and tracheary elements once stem expansion ceases.
While growth kinematic data cannot precisely position the base of the stem along the developmental continuum (growth kinematic profiling can only distinguish stem regions on the basis of their rates of expansion), it is clear from previous microscopic analysis [4, 51] that the OLD stage tissue displays an advanced phase of organ growth and cell wall maturation in the 10-15 cm tall Columbia plants examined in this study. Based on our present understanding of the stem maturation process, the CSS and OLD samples are expected to contain tissues actively engaged in earlier and later stages of secondary cell wall formation and reinforcement, respectively. Tables 5 and 6 present the twenty genes whose expression is ‘Higher in CSS relative to OLD’ and the twenty genes whose expression is ‘Higher in OLD relative to CSS’, respectively.
Displaying high expression in the CSS relative to OLD samples are GERMIN-LIKE PROTEIN 3 (GER3/GLP3)(At5g20630) and GERMIN-LIKE PROTEIN 1 (GER1/GLP1) (At1g72610). GER3 also appeared in the list of genes more highly expressed in MGR tissues than in CSS (Table 3), indicating that expression of this member of the GER gene family follows a steeply declining trajectory during the stem maturation process. While specific developmental roles for GLP1 and 3 have yet to be identified, GER proteins are apoplastic glycoproteins that have been widely associated with plant disease resistance and ROS modulation, particularly in the cereals . Interestingly, another Arabidopsis GER homologue (GLP10, At3G62020) whose expression was previously found to be highly correlated with secondary cell wall-associated CESAs (CesA4, 7 and 8) in regression analysis of public microarray datasets , also displayed elevated expression at both the CSS and OLD stages in our study (Additional files 2, 3 and 4: Table S1).
Also more highly expressed at this earlier stage of cell wall maturation are two pectate lyases (polygalacturonases), At3g07010 and At3g15720, previously associated with cell separation , and ALPHA-XYLOSIDASE 1/AXY3 (At1g68560), an exoglycosylase that acts specifically on non-fucosylated xyloglucans  and is essential for apoplastic xyloglucan modification . Several other up-regulated genes are less clearly linked to cell wall processes, but the functions of their encoded proteins may be related to the over-representation of ‘turgor pressure’ in the GO term enrichment analysis for this gene set (Additional file 7: Table S4).
The list of genes most highly expressed in OLD tissues relative to CSS tissues (Table 6) is particularly striking: six of the eight most highly up-regulated genes encode PLANT DEFENSIN (PDF) proteins, small cysteine-rich peptides homologous to anti-microbial peptides that are widely distributed within the eukaryotes . Since both CSS and OLD tissues were harvested only seconds apart, an artifactual pattern of wounding-induced gene induction is not likely. Instead, it appears that accumulation of the products of such classical “defense” genes may form an integral part of the normal maturation of the inflorescence stem, perhaps reflecting a commitment to protection of these tissues until fertilization and seed dispersal are successfully completed.
Relatively few cell wall-specific genes appear in the ‘higher in OLD than in CSS’ short list, with the exception of EXTENSIN 3/RSH and another proline-rich extensin-like family protein. EXT3/RSH plays an essential role in cell wall deposition through formation of EXTENSIN protein scaffolds that cross-link other cell wall constituents, thereby contributing to cell wall rigidification [57, 58]. The most up-regulated of all the genes at the OLD stage relative to the CSS stage is the chloroplast-localized FATTY ACID REDUCTASE 6 (FAR6). A similar pattern of elevated FAR6 expression was earlier observed in microarray analysis of epidermal peels from the stem base  as well as in stem sections harvested from the base of mature Arabidopsis Col-0 plants . Accompanying FAR6 in this list of most highly expressed genes is a wax synthase homologue (At5g22490), a co-occurrence pattern consistent with epidermal cells in fully mature stems actively synthesizing both their cuticle polyester network and the associated wax matrix. The modest representation of explicitly cell wall-associated genes in this CSS-to-OLD transition list implies that the CSS and OLD stage tissues share quite similar transcriptional profiles in terms of secondary cell wall formation processes, and that the metabolic commitment to cell wall fortification in stem tissues does not change dramatically after cessation of active elongation.
Stage-specific, whole-genome co-expression analysis
While differential gene expression datasets contrasting discrete growth stages provide initial insights into the biology underlying specific developmental transitions, potential functional relationships between gene products can also be revealed by considering transcript abundances across all the sampled developmental stages. The underlying rationale is that genes co-expressed at one stage and exhibiting similar association patterns across a broader developmental range may represent a subset of genes involved in specific biological processes.
The ‘mixed effects model’ approach used in this study allowed us to generate developmental stage 'estimates' from two-channel arrays, which can be expressed as mean fold-change values (biological replicates = 6) of transcript abundance at one stage relative to a hypothetical mean value of zero across the entire experiment. It should be noted that these 'estimates' can be computed with the same statistical power as applies to the log2 differential expression ratios reported in Tables 1, 2, 3, 4, 5, 6. This treatment provides a more intuitive means of visualizing gene expression trajectories, and provides the basis for formal co-expression analysis. The genes associated with each co-expression set (cluster) are identified in a filterable column within the full-genome dataset (Additional files 2, 3 and 4: Table S1), and their AGI codes are also listed separately in Additional file 8: Table S5 for easier access.
To test for functional relatedness of the genes populating these clusters, we examined gene ontology enrichment within several sub-clusters that exhibited marked patterns of coordinated up- or down-regulation specific to single developmental stages (Figure 2; sub-clusters 1.1, 2.1, 5.1, 5.2, and Cluster 8). Although box-plotting of the mean 'estimates' for these sub-clusters clearly demonstrates the extent to which these co-expression sets are synchronously up- or down-regulated with respect to the other stages (Additional file 9: Figure S2), gene ontology analysis revealed only modest GO term enrichment within these subsets of co-expressed genes (Additional file 10: Figure S3, Additional file 11: Table S6), indicating that, despite their shared expression pattern, the genes in each sub-cluster do not display obvious functional relatedness.
While gene expression profiling has been applied previously to the expanding inflorescence stem in Arabidopsis, those studies have all suffered from various limitations that impede our ability to accurately align the resulting expression profiles with the developmental state of the tissue being sampled. To address this issue, we have made use of growth kinematic profiling (GKP) to establish, for each plant being sampled, the precise state of growth extension that each stem section represents. Collectively, these sections span the growth and maturation states of the stem, and they also represent a cell wall development continuum. Thus, GKP-guided pooling of sections associated with discrete zones along that continuum (e.g. the point at which extension growth ceases) makes it possible to generate tissue samples whose gene expression profiles can be confidently aligned with specific developmental states.
While the growing stem is a complex organ consisting of multiple tissues, the common denominator across all of these cell types is the coordinated and initially rapid anisotropic expansion of their cell walls along the axis of growth. This expansion ultimately comes to a halt as the walls of some tissues, most notably the vascular tissues and supporting fibres, are reinforced with non-extensible secondary wall layers. Thus, while the biological processes being probed in these samples are not restricted to cell wall expansion/modification, the latter processes can be expected to dominate the broader landscape of transcriptional changes that accompany the maturation of the stem and its component tissues.
‘Young’ stage tissue displays a complex transcriptional profile
The YNG stage sampled in this study captures the top 1 cm of the Arabidopsis stem and so encompasses a developmentally complex region containing the shoot apical meristem, and up to twenty short internodes which bore flowers and/or siliques prior to harvest. Since multiple tissue development trajectories are being initiated within YNG samples, it is not surprising that cell wall-forming/modifying processes do not dominate the gene expression profile of YNG stage tissue in comparison with the MGR stage. It is interesting that the up-regulated MEE48 endo-β(1→3)-glucosidase (Figure 2) was previously characterized as an anther-specific gene whose proposed function involved callose degradation during pollen exine formation . Since all floral tissue had been deliberately removed from the YNG tissue at the time of sampling, MEE48 must play additional roles in development. The prominence of MEE48 expression in the YNG transcriptome may be related to the importance of callose hydrolysis for the development of new cell plate structures during cytokinesis , a process that is actively underway in the apical meristem.
Rapid tissue extension is associated with a unique transcriptional signature
Both the YNG/MGR comparison and MGR/CSS comparison gene lists provide a perspective on the genes up-regulated in cells undergoing extension growth at their maximum rate (MGR). Genes associated with gibberellic acid (GA)-mediated elongation are prominently expressed at the MGR stage, consistent with the known role of gibberellic acid as an effector of directional cell growth . The redox-associated cysteine-rich signal peptide, GA-STIMULATED ARABIDOPSIS 4 (GASA4) , is 3.4-fold up-regulated in MGR relative to CSS (Table 3), while the GA-responsive transcription factor, MINI ZINC FINGER 1 (MIF1) , is up-regulated 12-fold in the MGR stage relative to the YNG stage. Loss-of-function at the MIF1 locus results in unresponsiveness to GA and inflorescence stem dwarfism . MERISTEM-5 (MERI5B, XTH24; At4g30270) encodes a Group 2 xyloglucan endotransglycosylase/hydrolase , a group of proteins that facilitate the remodelling of hemicellulose to allow cellulose microfibril separation and ‘creep’ during anisotropic cell wall expansion . MERI5B expression is also elevated in MGR relative to YNG tissues, and clusters with MIF1 across all the developmental stages studied (Figure 2). In addition to displaying co-expression, MIF1 and MERI5B co-cluster with the arabinogalactan protein, AGP13, and the MYB61 transcription factor in a set that exhibits peak expression somewhat later, at the onset of cessation (Figure 2, cluster 7).
In contrast to XTH24, expression of another XTH, ENDOXYLOGLUCAN TRANSFERASE A1 (EXGT-A1), required for normal cell wall expansion  is restricted to the MGR stage (Figure 2, cluster 5), where it clusters with another member of the AGP family, AGP12, whose expression is >6-fold higher in the MGR stage relative to the YNG stage. An additional five AGPs (AGP14, 21, 22, 24 and FLA13) were found to be significantly up-regulated (q-value<0.05) in the MGR stage relative to YNG stage (Additional files 2, 3 and 4: Table S1), suggesting that a suite of AGPs may be contributing to the unique structural dynamics of rapidly expanding cells at the MGR stage.
‘Cessation’ stage gene expression is dominated by secondary cell wall processes
The composition of the shortlist of twenty genes most up-regulated at the CSS stage relative to the MGR stage (Table 4) is particularly striking, since thirteen appear to be functionally related to secondary cell wall biosynthesis. It is also notable that the population of this list is completely distinct from those genes whose expression is dominant in CSS tissues relative to OLD tissues, suggesting that these are genes whose expression becomes elevated as cell expansion slows and then remains elevated through ensuing stem maturation.
Central among these cell wall-associated genes are the three cellulose synthases that are essential for secondary cell wall formation [13, 35]. CESA8 is thought to belong to the same multi-protein biosynthetic complex as CESA4 and CESA7 , which have similar contributions to secondary cell wall synthesis , although their relative proportions in the CESA complex remain unknown. In our co-expression analysis, CESA8 clusters differently from CESA4 and CESA7, primarily due to increasingly elevated expression of CESA8 in the OLD tissue sample (Figure 2, cluster 8). In contrast, the expression of CESA4 and CESA7 (Figure 2, cluster 6) does not change significantly from the CSS to OLD stages. Since the relative stability and turnover rates for the three CESA proteins are unknown, these differences in gene expression do not necessarily conflict with the predicted abundance of their cognate proteins in the plasma membrane. It is also possible that the relative proportions of secondary cell wall-associated CESAs within cellulose synthase complexes do not remain fixed throughout the period of secondary cell wall formation.
By definition, only primary cell walls are capable of expanding , and the great majority of this expansion would be occurring above the point of cessation being sampled in this study. We therefore anticipated that cellulose synthase genes associated with primary cell wall formation (notably, CESA1, 3 and 6) would be up-regulated in the YNG and MGR stages relative to their expression in the CSS and OLD stages. Instead, CESA1 and 6 are not significantly differentially expressed between pre- and post-cessation stages, and CESA3 is actually more highly expressed in the CSS and OLD stages than in the YNG and MGR stages (Figure 2, Additional files 2, 3 and 4: Table S1). This is not an isolated example of the behavior of the CESA3 gene; Ko and Han (2004)  had earlier observed elevated CESA3 expression in the base of fully mature Arabidopsis Col-0 stems (>25 cm in height) relative to the bases of less mature stems (5 and 10 cm in height), while CESA1 and CESA6 expression declined with stem maturation. Ehlting et al. (2005)  also detected higher levels of CESA3 expression at the mid-point of 10 cm stems of Landsberg erecta (Ler) plants than in the top 3 cm stem region of those plants. Another study found that, while AtCESA4, 7 and 8 were up-regulated in mature Col-0 stems, none of the canonical ‘primary cell wall’ CESA genes (CESA1, 3, 6) were present in the list of CESAs significantly up-regulated in actively elongating tissues , likely due to the persistent expression of CESA1, 3 and 6 at later, non-elongating stages as well. In the present study, CESA6, which is generally considered to be important for cellulose deposition in primary cell walls , is most highly up-regulated in the MGR and CSS stage tissues (Figure 2, Cluster 5), where it is co-expressed with CESA2 (previously associated with radial cell wall reinforcement ) and with CESA10.
Collectively, these data suggest that deployment of particular cellulose synthases in plant cells does not follow a pattern of exclusive association with either actively elongating tissues (i.e. with primary cell wall synthesis) or post-elongation tissues (i.e. secondary cell wall synthesis). Instead, a more diverse co-occurring set of cell wall-forming/modifying processes may recruit distinct sub-sets of CESA and CESA-LIKE family members for specific developmental programming (e.g. intrusive growth of interfascicular fibres).
Another strong indication that the CSS tissue sample accurately captures the transition from primary to secondary cell wall formation is the presence in the ‘higher-in-CSS’ gene lists of a suite of genes specifically associated with accumulation of glucuronylarabinoxylans (‘xylans’), including XYLOSE SYNTHASE 3, FRA8, IRX9, RWA1, XTH19 and GUX1 (Figure 2, cluster 7). Other xylan-related genes (XYLOSE SYNTHASE 6; XTR4; BXL1; BXL2; EXGT-A1; XTH18) are also up-regulated in the CSS tissue (Additional files 2, 3 and 4: Table S1), but failed to qualify for the “top twenty” short-list of genes more highly expressed in CSS than in OLD.
A possible positive regulator of secondary cell wall development, the MYB61 transcription factor gene, is also strongly up-regulated in the CSS stage (Figures 2, S7). MYB61 has been proposed to promote cell wall lignification , and more specifically to regulate three cell wall–associated genes encoding the KNAT7 transcription factor, the lignin biosynthesis enzyme CAFFEOYL-COA O-METHYLTRANSFERASE 7, and a pectin methylesterase (At2g45220) . While our expression data confirm an association of MYB61 with secondary cell wall formation, such involvement is likely conditioned by other factors since MYB61 activity has also been associated with a wide range of biological processes in plants, including seed coat mucilage production , stomatal closure , and pleiotropic control of photosynthate partitioning .
Arabinogalactan proteins (AGP) form a largely uncharacterized class of proteoglycans that likely play structural and/or signaling roles in cell wall development. Indeed, a number of AGP family members exhibit significant modulation of expression at the onset of secondary cell wall formation (Additional files 2, 3 and 4: Table S1). For example, expression of AGP18, a member of a lysine-rich, GPI-anchored sub-family that includes AGP17 and AGP19, is down-regulated at the OLD stage relative to the MGR stage, and the loss-of-function agp18 mutant possesses shortened inflorescence stems , indicating a possible role for AGP18 in promoting cell wall expansion. AGP12 also shows a significant drop in expression coincident with the onset of cessation (Additional files 2, 3 and 4: Tables S1), consistent with a similar functional association.
FASCICLIN-LIKE 8 (FLA8), on other hand, is up-regulated at the OLD stage relative to the MGR stage (Additional files 2, 3 and 4: Table S1). FLA8/AGP8 belongs to sub-family of AGPs that contain a fasciclin domain and often possess a glycosyl phosphatidyl inositol (GPI) anchor . A poplar homologue of AtFLA8 was observed to be up-regulated significantly in tension wood, but not in opposite wood, in poplar stems, when compared to its expression in differentiating xylem . AGP21 also appears significantly up-regulated in the OLD stage relative to the MGR stage (Additional files 2, 3 and 4: Table S1). Interestingly, AGP21, similar in sequence to AGP12 and AGP14, is down-regulated ~4-fold upon silencing of the transcription factor PRODUCTION OF ANTHOCYANIN PIGMENT 1 (PAP1/MYB75), coincident with increased cell wall thickness in xylary and interfascicular fibres . The members of the large AGP family thus appear to have functionally diverged, as revealed in part through differences in spatiotemporal regulation of their expression .
Protection and fortification are hallmarks of OLD stem tissue
The base of the stem of 10-15 cm Columbia plants (OLD tissue samples) contains highly contrasting tissues, including live, photosynthetically active cells located adjacent to thick-walled, highly lignified fibres of the interfascicular region that are presumably in the advanced stages of programmed cell death.
In general, however, genes related to secondary cell wall synthesis are most active in this region of the lower stem (Figure 2; clusters 2, 6, 7 and 8). Interestingly, CESA8, which appears in cluster 8, exhibits its highest level of expression at this stage, as does XYLOGLUCAN ENDOTRANSGLYCOSYLASE/HYDROLASE 18 (XTH18). This expression data is consistent with other results linking xyloglucan deposition with late stages of secondary cell wall synthesis. For instance, incorporation of xyloglucan has been observed to continue in cotton fibres after cessation of wall extension , and PttXET16 activity was associated with secondary vasculature of poplar .
In the final stages of fibre secondary cell wall maturation, the polysaccharide matrix is typically impregnated with the phenylpropanoid polymer, lignin (reviewed in ). Several genes whose products are associated with the shikimic acid and phenylpropanoid pathways, and lignification exhibit corresponding expression patterns within this dataset, although differences in their clustering suggest subtle distinctions in the timing of their expression (Figure 2). For instance, PAL1, a member of the PHENYLALANINE AMMONIA- LYASE (PAL) gene family whose activity is required for phenylalanine allocation to phenylpropanoid metabolism, appears highly expressed in both the CSS and OLD stages, coincident with up-regulation of 3-DEOXY-D-ARABINO-HEPTULOSONATE 7-PHOSPHATE (DAHP) SYNTHASE 3 (DHS3), which regulates the intake of carbon into the shikimate pathway. Maximum expression of PAL4, on the other hand, occurs in the OLD stage, suggesting that different PAL family members may be playing distinct roles.
The large plant-to-plant variation in stem growth kinematic profiles that we identified earlier  makes it clear that earlier studies in which stems from multiple plants have been pooled to create biological replicate samples are inevitably compromised in their ability to accurately place cellular changes in an elongative development context. By contrast, the concordance of our GKP-guided gene expression data sets with current knowledge of cell wall biology provides strong evidence of the ability of this approach to capture development stage-specific information. At the same time, those known players in our data sets are accompanied by numerous genes of currently unknown biological function, which makes them high priority candidates for further research into the processes underpinning both plant cell expansion and deposition of the cellulose-rich cell walls that comprise plant biomass.
Plant growth, growth kinematic profiling and sampling
Plant growth and imaging was conducted as described in Hall & Ellis (2012) , using applied paper tags as synthetic optical markers for growth kinematic profiling. Tagged and imaged plants were harvested sequentially between 1 and 3 pm (mid-day where daylight cycle occurs between 6 am and 10 pm on a 16hL:8hD regime) at 20-minute intervals. Stem segments (~1 cm) were immediately snap-frozen in liquid nitrogen and deposited into 0.2 mL PCR tubes for −80°C storage. Segments were subsequently pooled on the basis of growth kinematic profiling data (shown in Additional file 1: Figure S1) and experimental design objectives, as outlined in Figure 1.
Whole stem segments (pooled according to growth kinematic profile equivalence) were homogenized in liquid nitrogen with a pre-cooled mortar and pestle. The frozen powder was then transferred to 1.5 mL microcentrifuge tubes, weighed and combined with TRIzol™ reagent (cat#15596-026, Invitrogen)(1 mL TRIzol per 100 mg tissue), vortexed and incubated at room temperature for 5 minutes. Chloroform (0.2 mL for each 1 mL of TRIzol) was added, vortexed for 15 seconds, incubated for 1 minute at room temperature, and centrifuged at 15000 g for 10 minutes at 4°C. The aqueous phase was transferred to fresh RNAse-free tubes and then combined with an equal volume of isopropanol and incubated 20 minutes on ice. RNAtotal was pelleted by centrifugation at 15000 g for 10 minutes, and pellets were washed with 1 mL 75% ethanol in RNAse-free water. Following a 5-minute pellet drying phase, pellets were resuspended in 25 μL RNAse-free water and incubated on ice for 1 hour. Each resuspension was treated with 1/10th volume (~2 μL) 10X DNAse I buffer and 1 μL 10X DNAse I (from RNAqueous® Micro kit; cat#AM1931, Ambion) for 20 minutes at 37°C followed by addition of 2 μL DNAse inactivation reagent (also from RNAqueous® Micro kit) and incubated at room temperature for 2 minutes. Samples were centrifuged at 13000 g for 1.5 minutes and the supernatant transferred to RNAse-free tubes and stored at −80°C.
Reverse transcription and labelling
For each biological replicate, approximately 20 μg RNAtotal was incubated with 2 ug oligo(dT) primer (cat#18418-012, 0.5 ug/ul, Invitrogen) in a 22.5 μL volume of RNA-primer mix and denatured at 70°C for 10 minutes. Reaction mix was prepared such that each sample contained 9 μL 5X First Strand buffer (supplied with Superscript II, Invitrogen, cat#18064-014), 0.23 μL each of 0.1 mM dATP (cat#10216-108, Invitrogen), dCTP (cat#10217-016, Invitrogen), and dGTP (cat#10218-014,Invitrogen), as well as 0.045 μL dTTP (cat#10219-012, Invitrogen) for a total reaction mix volume of 18.5 μL. This reaction mix was combined with 18.5 μL RNA-primer mix along with 1.5 μL (1.5 moles/μL) of the appropriate Cyanine dye; Cy5-dUTP (cat#45-000-740, Fisher) or Cy3-dUTP (cat#45-000-738, Fisher). After incubation at 42°C for 2 minutes, 1 uL 40U/uL RNAase Inhibitor (cat#10777-019, 40U/ul, Invitrogen) and 1.2 ul 200 U/μL Superscript II (cat#18064-014,200 U/ul, Invitrogen) were added for a final volume of 45 μL which was incubated at 42°C for 2.5 hours, then deactivated with 0.5 M NaH/50 mM EDTA at 65°C for 15 minutes. The reaction was neutralized with 7.1 uL 1 M Tris–HCl (pH7.5). Samples were cleaned of unlabeled probe via centrifugal filtration using Amicon 0.5-Ultra 30 kDa filters (cat#UFC503096, Millipore) prior to array hybridization.
For transcript profiling, we employed custom two-channel microarrays spotted with 26 929 70-mer oligonucleotides originally synthesized on the basis of ‘The Arabidopsis Information Resource’ (TAIR) ‘5’ release of the Arabidopsis genome (http://www.Arabidopsis.org), with gene annotations subsequently updated to the current genome release (TAIR10) . The microarray slides were first pre-conditioned by incubating them in Coplin jars with 50°C 2X SSC for 20 minutes, followed by room temperature washes with 0.2X SSC and ddH2O for 5 and 3 minutes, respectively using an Advawash AV400 machine (Advalytix/Beckman-Coulter). Pre-hybridization solution of 1X formamide-based hybridization buffer (pre-warmed to 80°C) from Vial 7 of the 3DNA Array 350 kit (cat#W300130, Genisphere) was then added to the gap between each slide and a pre-placed m-Series lifterslip (cat#48382-251, VWR) within the Slidebooster (Advalytix/Beckman-Coulter) hybridization chamber and subsequently incubated for 1–1.5 hours at 50°C with sonication (power=15, pulse=1 second 'on', 9 seconds 'off'). Slides were then washed in 2X SSC (0.2% SDS) for 15 minutes at 65°C followed by room temperature washes in 2X SSC and 0.2X SSC for 10 minutes each, and centrifuged at 700 rpm until dry in Advatubes (cat# OAX05216, Advalytix/Beckman-Coulter). Equal volumes of labeled Cy3 and Cy5 mixes (12.5 μL each) were combined with 25 μL 2X formamide buffer (Vial 7, 3DNA Array350 kit) and the 50 μL hybridization mix added to the gap between the 42°C pre-warmed slides and the pre-placed m-Series lifter-slips. Slides were then incubated at 42°C for 16–18 hours with sonication (power=15, pulse=1 second 'on', 9 seconds 'off'). Post-hybridization washing was carried out in reduced lighting with 42°C 2X SSC (0.2% SDS) for 15 minutes followed by room temperature washes with 2X SSC and 0.2X SSC for 15 minutes each, then centrifuged until dry at 700 rpm. Slides were stored in a light-proof desiccating chamber until fluorescence scanning.
Microarray scanning and spot quantification
Hybridized arrays were scanned with a ScanArray Express HT (Perkin-Elmer) scanner and associated software, using 543 nm laser irradiation for Cy3, and 633 m laser for Cy5 fluorescence. Laser power was adjusted for each slide individually within the range of 95-100% such that ~1-2% of spotted probes (presumed positive controls) yielded saturated signals. PMT gain ranged from 60-95%, set for each slide such that fluorescence intensity of sub-grid regions surrounding spots did not exceed 400 (16-bit grayscale). TIFF images of array scans were imported into Imagene (Biodiscovery Software) and grid templates were roughly placed before applying the 'auto-adjust' function to best fit the subgrids on a per-spot basis, allowing spot size variation from 15-20 μm. Median pixel intensities computed from spot regions were used to represent spot intensity in subsequent analyses.
Microarray data analysis
Data analysis was carried out in the statistical programming environment R (cran.r-project.org/) using custom scripts and contributed packages. To remove local background noise, the mean signal intensity of the dimmest five percent of spots within each of 48 subgrids was subtracted from each array element using a custom script, then variance stabilization normalization (VSN) was applied to each channel to normalize for non-linearity in variance across spot intensities  using the function ‘vsn’ (‘vsn’ package, Bioconductor). Normalized intensities were then fit to the mixed effects model  using the ‘lme’ function (‘nlme’ package), and all pairwise differential expressions for array elements were computed as the log2 intensity difference values between treatment class intensities. Associated measures of significance (p-values relative to null hypothesis, log2 difference equals zero) were corrected for false-discovery rate using a custom script based upon standard q-value calculation , and ‘estimates’ were computed as the log2 intensity difference of each treatment class to the mean of all treatment classes (normalized to zero). Associated measures of significance (p-values relative to null hypothesis of log2 difference = 0) were also corrected for false-discovery rate as described above. Raw and output data were exported along with TAIR10 annotations in the supplemental data (Additional files 2, 3 and 4: Table S1). For hierarchical clustering, dissimilarity matrices were computed from filtered datasets using the ‘diana’ function (‘cluster’ package) and rendered as dendrograms using the ‘dendro’ function (‘cluster’ package). Heatmaps were generated using the 'heatmap.2' function ('gplots' package).
Availability of supporting data
Gene annotations, raw expression data, statistical analysis, mean differentials, mean estimates, and gene categorization for the full genome are provided in Additional files 2, 3 and 4: Table S1, and have been deposited with ArrayExpress following the MIAME conventions , as accession E-MEXP-3525.
We would like to acknowledge invaluable discussions with Dr. Jürgen Ehlting (University of Victoria) regarding microarray design and application, and inflorescence stem sampling challenges. We are greatly indebted to Anne Haegert and Dr. Stéphane Le Bihan at The Jack Bell Prostate Centre, Vancouver, for microarray printing and quality control. This project would not have been possible without the resources of the Michael Smith Laboratories, UBC, and financial support from Genome Canada, and the Natural Sciences and Engineering Research Council of Canada. Finally, we are grateful to Rick White (SCARL, University of British Columbia) for statistical consultation.
- Cosgrove DJ: Assembly and enlargement of the primary cell wall in plants. Annu Rev Cell Dev Biol. 1997, 13: 171-201. 10.1146/annurev.cellbio.13.1.171.PubMedView Article
- Szymanski DB, Cosgrove DJ: Dynamic coordination of cytoskeletal and cell wall systems during plant cell morphogenesis. Curr Biol. 2009, 19: R800-11. 10.1016/j.cub.2009.07.056.PubMedView Article
- Somerville CR, Bauer S, Brininstool G, Facette M, Hamann T, Milne J, Osborne E, Paredez A, Persson S, Raab T, Vorwerk S, Youngs H: Toward a systems approach to understanding plant cell walls. Science. 2004, 306 (5705): 2206-2211. 10.1126/science.1102765.PubMedView Article
- Ehlting J, Mattheus N, Aeschliman DS, Li E, Hamberger B, Cullis IF, Zhuang J, Kaneda M, Mansfield SD, Samuels L, Ritland K, Ellis BE, Bohlmann J, Douglas CJ: 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
- Imoto K, Yokoyama R, Nishitani K: Comprehensive approach to genes involved in cell wall modifications in Arabidopsis thaliana. Plant Mol Biol. 2005, 58 (2): 177-192. 10.1007/s11103-005-5344-7.PubMedView Article
- Ko J, Han KH: Arabidopsis whole-transcriptome profiling defines the features of coordinated regulations that occur during secondary growth. Plant Mol Biol. 2004, 55 (3): 433-453. 10.1007/s11103-004-1051-z.PubMedView Article
- Minic Z, Jamet E, San-Clemente H, Pelletier S, Renou JP, Rihouey C, Okinyo DP, Proux C, Lerouge P, Jouanin L: Transcriptomic analysis of Arabidopsis developing stems: a close-up on cell wall genes. BMC Plant Biol. 2009, 9 (1): 6-10.1186/1471-2229-9-6.PubMedPubMed CentralView Article
- Oh S, Park S, Han KH: Transcriptional regulation of secondary growth in Arabidopsis thaliana. J Exp Bot. 2003, 54 (393): 2709-2722. 10.1093/jxb/erg304.PubMedView Article
- Jamet E, Roujol D, San-Clemente H, Irshad M, Soubigou-Taconnat L, Renou J, Pont-Lezica R: Cell wall biogenesis of Arabidopsis thaliana elongating cells: transcriptomics complements proteomics. BMC Genomics. 2009, 10 (1): 505-10.1186/1471-2164-10-505.PubMedPubMed CentralView Article
- Milioni D, Sado P, Stacey N, Domingo C, Roberts K, McCann M: Differential expression of cell-wall-related genes during the formation of tracheary elements in the Zinnia mesophyll cell system. Plant Mol Biol. 2001, 47 (1–2): 221-238.PubMedView Article
- Demura T, Tashiro G, Horiguchi G, Kishimoto N, Kubo M, Matsuoka N, Minami A, Nagata-Hiwatashi M, Nakamura K, Okamura Y, Sassa N, Suzuki S, Yazaki J, Kikuchi S, Fukuda H: Visualization by comprehensive microarray analysis of gene expression programs during transdifferentiation of mesophyll cells into xylem cells. Proc Natl Acad Sci USA. 2002, 99 (24): 15794-15799. 10.1073/pnas.232590499.PubMedPubMed CentralView Article
- Ohashi-Ito K, Oda Y, Fukuda H: Arabidopsis VASCULAR-RELATED NAC-DOMAIN6 directly regulates the genes that govern programmed cell death and secondary wall formation during xylem differentiation. Plant Cell. 2010, 22 (10): 3461-3473. 10.1105/tpc.110.075036.PubMedPubMed CentralView Article
- Persson S, Wei H, Milne J, Page GP, Somerville CR: Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets. Proc Natl Acad Sci USA. 2005, 102: 8633-8638. 10.1073/pnas.0503392102.PubMedPubMed CentralView Article
- Brown DM, Zeef LAH, Ellis J, Goodacre R, Turner SR: Identification of novel genes in Arabidopsis involved in secondary cell wall formation using expression profiling and reverse genetics. Plant Cell. 2005, 17 (8): 2281-2295. 10.1105/tpc.105.031542.PubMedPubMed CentralView Article
- Hall HC, Ellis BE: Developmentally equivalent tissue sampling based on growth kinematic profiling of Arabidopsis inflorescence stems. New Phytol. 2012, 194 (1): 287-296. 10.1111/j.1469-8137.2012.04060.x.PubMedView Article
- Tempelman RJ: Assessing statistical precision, power, and robustness of alternative experimental designs for two color microarray platforms based on mixed effects models. Vet Immunol Immunopathol. 2005, 105 (3–4): 175-186.PubMedView Article
- Vandepoele K, Quimbaya M, Casneuf T, De Veylder L, Van de Peer Y: Unraveling transcriptional control in Arabidopsis using cis-regulatory elements and coexpression networks. Plant Physiol. 2009, 150 (2): 535-546. 10.1104/pp.109.136028.PubMedPubMed CentralView Article
- Axelos M, Bardet C, Liboz T, Le Van Thai A, Curie C, Lescure B: The gene family encoding the Arabidopsis thaliana translation elongation factor EF-1 alpha: molecular cloning, characterization and expression. Mol Gen Genet. 1989, 219: 106-112.PubMedView Article
- Tarutani Y, Morimoto T, Sasaki A, Yasuda M, Nakashita H, Yoshida S, Yamaguchi I, Suzuki Y: Molecular characterization of two highly homologous receptor-like kinase genes, RLK902 and RKL1, in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2004, 68: 1935-1941. 10.1271/bbb.68.1935.PubMedView Article
- Cock JM, McCormick S: A large family of genes that share homology with CLAVATA3. Plant Physiol. 2001, 126: 939-942. 10.1104/pp.126.3.939.PubMedPubMed CentralView Article
- Bell E, Creelman RA, Mullet JE: A chloroplast lipoxygenase is required for wound-induced jasmonic acid accumulation in Arabidopsis. Proc Natl Acad Sci. 1995, 92 (19): 8675-8679. 10.1073/pnas.92.19.8675.PubMedPubMed CentralView Article
- Jun J, Fiume E, Roeder AHK, Meng L, Sharma VK, Osmont KS, Baker C, Ha CM, Meyerowitz EM, Feldman LJ, Fletcher JC: Comprehensive analysis of CLE polypeptide signaling gene expression and overexpression activity in Arabidopsis. Plant Physiol. 2010, 154 (4): 1721-1736. 10.1104/pp.110.163683.PubMedPubMed CentralView Article
- Dobritsa AA, Geanconteri A, Shrestha J, Carlson A, Kooyers N, Coerper D, Urbanczyk-Wochniak E, Bench BJ, Sumner LW, Swanson R, Preuss D: A large-scale genetic screen in Arabidopsis to identify genes involved in pollen exine production. Plant Physiol. 2011, 157 (2): 947-970. 10.1104/pp.111.179523.PubMedPubMed CentralView Article
- Koizumi K, Yokoyama R, Nishitani K: Mechanical load induces upregulation of transcripts for a set of genes implicated in secondary wall formation in the supporting tissue of Arabidopsis thaliana. J Plant Res. 2009, 122: 651-659. 10.1007/s10265-009-0251-7.PubMedView Article
- Apel K, Hirt H: Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol. 2004, 55: 373-399. 10.1146/annurev.arplant.55.031903.141701.PubMedView Article
- Boerjan W, Ralph J, Baucher M: Lignin biosynthesis. Annu Rev Plant Biol. 2003, 54: 519-546. 10.1146/annurev.arplant.54.031902.134938.PubMedView Article
- Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W: Lignin biosynthesis and structure. Plant Physiol. 2010, 153 (3): 895-905. 10.1104/pp.110.155119.PubMedPubMed CentralView Article
- Wenzel CL, Hester Q, Mattsson J: Identification of genes expressed in vascular tissues using NPA-induced vascular overgrowth in Arabidopsis. Plant Cell Physiol. 2008, 49 (3): 457-468. 10.1093/pcp/pcn023.PubMedView Article
- Romano JM, Dubos C, Prouse MB, Wilkins O, Hong H, Poole M, Kang K, Li E, Douglas CJ, Western TL, Mansfield SD, Campbell MM: AtMYB61, an R2R3-MYB transcription factor, functions as a pleiotropic regulator via a small gene network. New Phytol. 2012, 195 (4): 774-786. 10.1111/j.1469-8137.2012.04201.x.PubMedView Article
- Liang Y, Dubos C, Dodd IC, Holroyd GH, Hetherington AM, Campbell MM: AtMYB61, an R2R3-MYB transcription factor controlling stomatal aperture in Arabidopsis thaliana. Curr Biol. 2005, 15 (13): 1201-1206. 10.1016/j.cub.2005.06.041.PubMedView Article
- Newman LJ, Perazza DE, Juda L, Campbell MM: Involvement of the R2R3-MYB, AtMYB61, in the ectopic lignification and dark-photomorphogenic components of the det3 mutant phenotype. Plant J. 2004, 37 (2): 239-250. 10.1046/j.1365-313X.2003.01953.x.PubMedView Article
- Penfield S, Meissner RC, Shoue DA, Carpita NC, Bevan MW: MYB61 is required for mucilage deposition and extrusion in the Arabidopsis seed coat. Plant Cell. 2001, 13 (12): 2777-2791.PubMedPubMed CentralView Article
- Hu W, Ma H: Characterization of a novel putative zinc finger gene MIF1: involvement in multiple hormonal regulation of Arabidopsis development. Plant J. 2006, 45 (3): 399-422. 10.1111/j.1365-313X.2005.02626.x.PubMedView Article
- Lin P, Pomeranz MC, Jikumaru Y, Kang SG, Hah C, Fujioka S, Kamiya Y, Jang J: The Arabidopsis tandem zinc finger protein AtTZF1 affects ABA- and GA-mediated growth, stress and gene expression responses. Plant J. 2011, 65 (2): 253-268. 10.1111/j.1365-313X.2010.04419.x.PubMedView Article
- Gardiner JC, Taylor NG, Turner SR: Control of cellulose synthase complex localization in developing xylem. Plant Cell. 2003, 15 (8): 1740-1748. 10.1105/tpc.012815.PubMedPubMed CentralView Article
- Hamann T, Osborne E, Young HL, Misson J, Nussaum L, Somerville CR: Global expression analysis of CESA and CSL genes in Arabidopsis. Cellulose. 2004, 11: 279-286.View Article
- Holland N, Holland D, Helentjaris T, Dhugga KS, Xoconostle-Cazares B, Delmer DP: A comparative analysis of the plant cellulose synthase (CesA) gene family. Plant Physiol. 2000, 123 (4): 1313-1324. 10.1104/pp.123.4.1313.PubMedPubMed CentralView Article
- Hossain MA, Noh H, Kim K, Koh E, Wi S, Bae H, Lee H, Hong S: Mutation of the chitinase-like protein-encoding AtCTL2 gene enhances lignin accumulation in dark-grown Arabidopsis seedlings. J Plant Physiol. 2010, 167 (8): 650-658. 10.1016/j.jplph.2009.12.001.PubMedView Article
- Sánchez-Rodríguez C, Bauer S, Hématy K, Saxe F, Ibanez AB, Vodermaier V, Konlechner C, Sampathkumar A, Rüggeberg M, Aichinger E, Neumetzler L, Burgert I, Somerville C, Hauser M, Persson S: CHITINASE-LIKE1/POM-POM1 and its homolog CTL2 Are glucan-interacting proteins important for cellulose biosynthesis in Arabidopsis. Plant Cell. 2012, 24 (2): 589-607. 10.1105/tpc.111.094672.PubMedPubMed CentralView Article
- Pattathil S, Harper AD, Bar-Peled M: Biosynthesis of UDP-xylose: characterization of membrane-bound AtUxs2. Planta. 2005, 221 (4): 538-548. 10.1007/s00425-004-1471-7.PubMedView Article
- Wu A, Hörnblad E, Voxeur A, Gerber L, Rihouey C, Lerouge P, Marchant A: Analysis of the Arabidopsis IRX9/IRX9-L and IRX14/IRX14-L pairs of glycosyltransferase genes reveals critical contributions to biosynthesis of the hemicellulose glucuronoxylan. Plant Physiol. 2010, 153 (2): 542-554. 10.1104/pp.110.154971.PubMedPubMed CentralView Article
- Lee C, Teng Q, Huang W, Zhong R, Ye Z: The F8H glycosyltransferase is a functional paralog of FRA8 involved in glucuronoxylan biosynthesis in Arabidopsis. Plant Cell Physiol. 2009, 50 (4): 812-827. 10.1093/pcp/pcp025.PubMedView Article
- Maris A, Kaewthai N, Eklöf JM, Miller JG, Brumer H, Fry SC, Verbelen J, Vissenberg K: Differences in enzymic properties of five recombinant xyloglucan endotransglucosylase/hydrolase (XTH) proteins of Arabidopsis thaliana. J Exp Bot. 2011, 62 (1): 261-271. 10.1093/jxb/erq263.PubMedView Article
- Lee C, Teng Q, Zhong R, Ye Z: The four Arabidopsis REDUCED WALL ACETYLATION genes are expressed in secondary wall-containing cells and required for the acetylation of xylan. Plant Cell Physiol. 2011, 52: 1289-1301. 10.1093/pcp/pcr075.PubMedView Article
- Derbyshire P, McCann MC, Roberts K: Restricted cell elongation in Arabidopsis hypocotyls is associated with a reduced average pectin esterification level. BMC Plant Biol. 2007, 7: 31-10.1186/1471-2229-7-31.PubMedPubMed CentralView Article
- Hongo S, Sato K, Yokoyama R, Nishitani K: Demethylesterification of the primary wall by PECTIN METHYLESTERASE35 provides mechanical support to the Arabidopsis stem. Plant Cell. 2012, 24 (6): 2624-34. 10.1105/tpc.112.099325.PubMedPubMed CentralView Article
- Kim JB, Carpita NC: Changes in esterification of the uronic acid groups of cell wall polysaccharides during elongation of maize coleoptiles. Plant Physiol. 1992, 98: 646-653. 10.1104/pp.98.2.646.PubMedPubMed CentralView Article
- Siedlecka A, Wiklund S, Peronne M, Micheli F, Lesniewska J, Sethson I, Edlund U, Richard L, Sundberg B, Mellerowicz EJ: Pectin methyl esterase inhibits intrusive and symplastic cell growth in developing wood cells of Populus. Plant Physiol. 2008, 146: 554-565.PubMedPubMed CentralView Article
- Berthet S, Demont-Caulet N, Pollet B, Bidzinski P, Cézard L, Le Bris P, Borrega N, Hervé J, Blondet E, Balzergue S, Lapierre C, Jouanin L: Disruption of LACCASE4 and 17 results in tissue-specific alterations to lignification of Arabidopsis thaliana stems. Plant Cell. 2011, 23 (3): 1124-1137. 10.1105/tpc.110.082792.PubMedPubMed CentralView Article
- Ko J, Beers E, Han K: Global comparative transcriptome analysis identifies gene network regulating secondary xylem development in Arabidopsis thaliana. Mol Genet Genomics. 2006, 276 (6): 517-531. 10.1007/s00438-006-0157-1.PubMedView Article
- Ko J, Han K, 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 (2): 1069-1083. 10.1104/pp.104.038844.PubMedPubMed CentralView Article
- Manosalva PM, Bruce M, Leach JE: Rice 14-3-3 protein (GF14e) negatively affects cell death and disease resistance. Plant J. 2011, 68 (5): 777-787. 10.1111/j.1365-313X.2011.04728.x.PubMedView Article
- Sun L, van Nocker S: Analysis of promoter activity of members of the PECTATE LYASE-LIKE (PLL) gene family in cell separation in Arabidopsis. BMC Plant Biol. 2010, 10 (1): 152-10.1186/1471-2229-10-152.PubMedPubMed CentralView Article
- Sampedro J, Pardo B, Gianzo C, Guitián E, Revilla G, Zarra I: Lack of α-xylosidase activity in Arabidopsis alters xyloglucan composition and results in growth defects. Plant Physiol. 2010, 154 (3): 1105-1115. 10.1104/pp.110.163212.PubMedPubMed CentralView Article
- Günl M, Pauly M: AXY3 encodes a α-xylosidase that impacts the structure and accessibility of the hemicellulose xyloglucan in Arabidopsis plant cell walls. Planta. 2011, 233 (4): 707-719. 10.1007/s00425-010-1330-7.PubMedPubMed CentralView Article
- Thomma B, Cammue B, Thevissen K: Plant defensins. Planta. 2002, 216: 193-202. 10.1007/s00425-002-0902-6.PubMedView Article
- Rapaport H: Ordered peptide assemblies at interfaces. Supramol Chem. 2006, 18 (5): 445-454. 10.1080/10610270600665905.View Article
- Cannon MC, Terneus K, Hall Q, Tan L, Wang Y, Wegenhart BL, Chen L, Lamport DTA, Chen Y, Kieliszewski MJ: Self-assembly of the plant cell wall requires an extensin scaffold. Proc Natl Acad Sci. 2008, 105 (6): 2226-2231. 10.1073/pnas.0711980105.PubMedPubMed CentralView Article
- Suh MC, Samuels AL, Jetter R, Kunst L, Pollard M, Ohlrogge J, Beisson F: Cuticular lipid composition, surface structure, and gene expression in Arabidopsis stem epidermis. Plant Physiol. 2005, 139 (4): 1649-1665. 10.1104/pp.105.070805.PubMedPubMed CentralView Article
- Chen X, Kim J: Callose synthesis in higher plants. Plant Signal Behav. 2009, 4 (6): 489-492. 10.4161/psb.4.6.8359.PubMedPubMed CentralView Article
- Schwechheimer C: Understanding gibberellic acid signaling - are we there yet?. Curr Opin Plant Biol. 2008, 11 (1): 9-15. 10.1016/j.pbi.2007.10.011.PubMedView Article
- Rubinovich L, Weiss D: The Arabidopsis cysteine-rich protein GASA4 promotes GA responses and exhibits redox activity in bacteria and in planta. Plant J. 2010, 64 (6): 1018-1027. 10.1111/j.1365-313X.2010.04390.x.PubMedView Article
- Rose JKC, Braam J, Fry SC, Nishitani K: The XTH family of enzymes involved in xyloglucan endotransglucosylation and endohydrolysis: Current perspectives and a new unifying nomenclature. Plant Cell Physiol. 2002, 43 (12): 1421-1435. 10.1093/pcp/pcf171.PubMedView Article
- Akamatsu T, Hanzawa Y, Ohtake Y, Takahashi T, Nishitani K, Komeda Y: Expression of endoxyloglucan transferase genes inacaulis mutants of Arabidopsis. Plant Physiol. 1999, 121 (3): 715-722. 10.1104/pp.121.3.715.PubMedPubMed CentralView Article
- Endler A, Persson S: Cellulose Synthases and synthesis in Arabidopsis. Mol Plant. 2011, 4 (2): 199-211. 10.1093/mp/ssq079.PubMedView Article
- McNeil M, Darvill AG, Fry SC, Albersheim P: Structure and function of the primary cell walls of plants. Annu Rev Biochem. 1984, 53 (1): 625-663. 10.1146/annurev.bi.53.070184.003205.PubMedView Article
- Doering A, Lathe R, Persson S: An update on xylan synthesis. Mol Plant. 2012, 5 (4): 769-771. 10.1093/mp/sss049.PubMedView Article
- Yang J, Sardar HS, McGovern KR, Zhang Y, Showalter AM: A lysine-rich arabinogalactan protein in Arabidopsis is essential for plant growth and development, including cell division and expansion. Plant J. 2007, 49 (4): 629-640. 10.1111/j.1365-313X.2006.02985.x.PubMedView Article
- MacMillan CP, Mansfield SD, Stachurski ZH, Evans R, Southerton SG: Fasciclin-like arabinogalactan proteins: specialization for stem biomechanics and cell wall architecture in Arabidopsis and eucalyptus. Plant J. 2010, 62: 689-703. 10.1111/j.1365-313X.2010.04181.x.PubMedView Article
- Lafarguette F, Leplé J, Déjardin A, Laurans F, Costa G, Lesage-Descauses M, Pilate G: Poplar genes encoding fasciclin-like arabinogalactan proteins are highly expressed in tension wood. New Phytol. 2004, 164 (1): 107-121. 10.1111/j.1469-8137.2004.01175.x.View Article
- Ito S, Suzuki Y, Miyamoto K, Ueda J, Yamaguchi I: AtFLA11, a fasciclin-like arabinogalactan-protein, specifically localized in sclerenchyma cells. Biosci Biotechnol Biochem. 2005, 69 (10): 1963-1969. 10.1271/bbb.69.1963.PubMedView Article
- Bhargava A, Mansfield SD, Hall HC, Douglas CJ, Ellis B: MYB75 functions in regulation of secondary cell wall formation in the Arabidopsis inflorescence stem. Plant Physiol. 2010, 154 (3): 1428-1438. 10.1104/pp.110.162735.PubMedPubMed CentralView Article
- Tokumoto H, Wakabayashi K, Kamisaka S, Hoson T: Changes in the sugar composition and molecular mass distribution of matrix polysaccharides during cotton fiber development. Plant Cell Physiol. 2002, 43 (4): 411-418. 10.1093/pcp/pcf048.PubMedView Article
- Bourquin V, Nishikubo N, Abe H, Brumer H, Denman S, Eklund M, Christiernin M, Teeri TT, Sundberg B, Mellerowicz EJ: Xyloglucan endotransglycosylases have a function during the formation of secondary cell walls of vascular tissues. Plant Cell. 2002, 14 (12): 3073-3088. 10.1105/tpc.007773.PubMedPubMed CentralView Article
- Humphreys JM, Chapple C: Rewriting the lignin roadmap. Curr Opin Plant Biol. 2002, 5 (3): 224-229. 10.1016/S1369-5266(02)00257-1.PubMedView Article
- Douglas CJ, Ehlting J: Arabidopsis thaliana full genome longmer microarrays: a powerful gene discovery tool for agriculture and forestry. Transgenic Res. 2005, 14: 551-561. 10.1007/s11248-005-8926-x.PubMedView Article
- Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18 (suppl 1): S96-S104. 10.1093/bioinformatics/18.suppl_1.S96.PubMedView Article
- Storey J: A direct approach to false discovery rates. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2002, 64: 479-10.1111/1467-9868.00346.View Article
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.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.