Transcript profiling of Populus tomentosa genes in normal, tension, and opposite wood by RNA-seq
© Chen et al.; licensee BioMed Central. 2015
Received: 26 October 2014
Accepted: 24 February 2015
Published: 10 March 2015
Wood formation affects the chemical and physical properties of wood, and thus affects its utility as a building material or a feedstock for biofuels, pulp and paper. To obtain genome-wide insights on the transcriptome changes and regulatory networks in wood formation, we used high-throughput RNA sequencing to characterize cDNA libraries of mature xylem from tension wood (TW), opposite wood (OW), and normal wood (NW), in the industrial tree species Populus tomentosa.
Our sequencing generated 140,978,316 (TW), 128,972,228 (OW), and 117,672,362 (NW) reads, corresponding to 10,127 (TW), 10,129 (OW), and 10,129 (NW) unique genes. Of these, 361 genes were differentially transcribed between TW and OW (log2FC ≥ 1 or ≤ -1, FDR < 0.05), 2,658 differed between OW and NW, and 2,417 differed between TW and NW. This indicates that NW differs significantly from the wood in branches; GO term analysis also indicated that OW experienced more transcriptome remodeling. The differentially expressed genes included 97 encoding transcription factors (TFs), 40 involved in hormone signal transduction, 33 in lignin biosynthesis, 21 in flavonoid biosynthesis, and 43 in cell wall metabolism, including cellulose synthase, sucrose synthase, and COBRA. More than half of the differentially expressed TF showed more than 4-fold lower transcript levels in NW compared with TW or OW, indicating that TF abundances differed dramatically in different wood types and may have important roles in the formation of reaction wood. In addition, transcripts of most of the genes involved in lignin biosynthesis were more abundant in OW compared with TW, consistent with the higher lignin content of OW. We constructed two transcriptomic networks for the regulation of lignin and cellulose biosynthesis, including TFs, based on the co-expression patterns of different genes. Lastly, we used reverse transcription quantitative PCR to validate the differentially expressed genes identified.
Here, we identified the global patterns and differences in gene expression among TW, OW, and NW, and constructed two transcriptomic regulatory networks involved in TW formation in P. tomentosa. We also identified candidate genes for molecular breeding of wood quality, and provided a starting point to decipher the molecular mechanisms of wood formation in Populus.
Trees have important functions in ecosystems, and in providing feedstocks for global industry. For example, wood has uses in construction, pulp and paper, and a major potential role in biofuels as a renewable, cost-effective alternative to fossil fuels. The complex chemical makeup of wood also makes it an ideal raw material for a potential future “ligno-chemical” industry that could replace the petrochemical industry, providing raw materials for the manufacture of plastic, chemicals, food, and textiles . However, our understanding of how wood develops remains far from complete, and very little is known about the cellular, molecular, and developmental processes that underlie wood formation .
Wood formation requires the coordinated expression of numerous structural genes involved in cell division, cell differentiation, programmed cell death, and heartwood formation and wood develops under the control of regulatory genes that remain virtually unknown . Reaction wood, which develops under stress caused by branch bending or gravity, has been used as a model system for the functional genomics of wood formation, particularly in exploration of carbon partitioning between lignin and cellulose in woody plants. Tension wood (TW), which develops under tension stress, has more cellulose and less lignin and hemicelluloses than normal wood (NW); by contrast, opposite wood (OW, also termed compression wood, CW, in gymnosperm species), which develops under compression stress, has more lignin . Thus, many studies have used microarray analysis of TW and OW in artificially or naturally bent trunks of Populus to identify key genes that contribute to the formation of reaction wood, including genes encoding plant hormones, lignin and cellulose biosynthetic enzymes, and transcription factors (TFs) and other potential regulators [3-6]. For example, a cDNA microarray analysis of TW induced by bending force in young Populus tremula (L.) demonstrated that TW formation involved reprogramming of carbohydrate metabolism, including increases of PttCesA8-2 and PttCesA3-2 transcripts and a decrease of transcripts involved in lignin biosynthesis . For lignin biosynthesis, 4-Coumarate-CoA Ligase (4CL) expression increases during CW formation in pine species  and transgenic suppression of 4CL in Pinus radiata decreased lignin contents . Plant hormones such as auxin, ethylene, and gibberellin also have important roles in the formation of TW . For example, transgenic approaches showed that the endogenous ethylene produced in leaning trees acts as a key regulator of the asymmetrical cambial growth in Populus TW . Recent work also reported that the formation of TW and stem gravitropism in Acacia mangium seedlings requires gibberellins .
The complex process of wood formation requires various genes and pathways; therefore, genome-wide transcriptome analysis, especially by high-throughput RNA sequencing (RNA-seq), provides a useful approach [3-6] to explore the mechanisms underlying wood formation. RNA-seq can detect rare transcripts, splice variants, and novel transcripts . Moreover, RNA-seq data provide absolute transcript levels, rather than relative measurements, thus overcoming many limitations of microarray analysis . To date, most studies have focused on the difference between TW and NW in artificially bent Populus trunks, and have used cDNA microarrays. However, little is known about transcription and regulation in Populus branches (TW and OW) under gravity stress, especially combined with analysis of NW using RNA-Seq.
To provide accurate and comprehensive genome-wide insights into the molecular mechanisms involved in the formation of TW, we used RNA-seq to reveal transcriptome changes in TW, OW, and NW in Populus tomentosa (Carr.), an important industrial species for pulp and paper in China. Our results improve our understanding of the formation of reaction wood in response to gravity, including identifying co-expression networks and TFs likely involved in the regulatory network controlling cellulose and lignin biosynthesis. To the best of our knowledge, this study is the first to characterize the xylem transcriptome of P. tomentosa using RNA-seq, and may serve as a foundation for further studies of wood formation, particularly the formation of special wood in Populus.
RNA isolation, library construction, and Solexa sequencing
As biological replicates, this study used three individuals, 30-year-old clones from one genotype of P. tomentosa. The three trees were planted in the national nursery of Guan Xian County (Shandong Province, China; N36°28′28″, E115°26′17″) and were not treated to induce any stress conditions. Mature xylem tissue samples from NW, TW, and OW (normal, tension, and opposite wood) were isolated from the same individual to enable comparisons in the same genetic background. From each plant, we chose three branches that are in the same side of the tree and similar diameter (about 8 cm), and were more than 10 years old. TW (upper side of the branch) and OW (lower side of the branch) were collected from the same branch, using a sharp chisel and after removing the bark of the sampling area, following the method described in Li et al. . NW, which represents the stem xylem, was isolated from the same side of the tree at breast height, approximately one meter above the ground. All the samples were about 2 cm × 1 cm and 4–5 mm deep. Samples were collected in the morning, immediately frozen in liquid nitrogen, and stored at −80°C for isolation of RNA. Total RNA was isolated by a modified CTAB method  with isopropanol instead of lithium chloride for RNA precipitation. DNase was applied to eliminate any genomic DNA in the total RNA. RNA quality was monitored using a NanoDrop ND-1000 and Agilent Bioanalyzer 2100. Extracted RNA was used for RNA library construction. First, the high-quality total RNA was purified with the RNeasy micro kit (Cat#74004, Qiagen). Purified total RNA was processed with the TruSeq RNA Sample Preparation Guide to build the cDNA library. The sequencing of cDNA was carried by Shanghai Biotechnology Corporation (Shanghai, China) using Illumina HiSeq 2000 following the cBot User Guide and HiSeq 2000 User Guide with the paired-end program.
Transcriptome mapping and differentially expressed transcripts
The fold change (FC) for NW vs TW, for example, equals the FPKM of NW divided by FPKM of TW, and so on. The differentially expressed genes were selected using log2FC ≥ 1 or FC ≤ -1 and FDR < 0.05 (false discovery rate control, q-value).
Gene annotation and construction of the co-expression network
Gene annotations were carried out using PopGenie (http://www.popgenie.org/)  and gene ontology terms were analyzed using agriGO (http://bioinfo.cau.edu.cn/agriGO/index.php) . The enriched GO categories were checked using an FDR-adjusted value of ≤0.05 as the cutoff for significant GO categories. The co-expression network analysis was performed in R using the weighted gene co-expression network analysis (WGCNA) package, as previously described . Briefly, only differentially expressed genes involved in cellulose and lignin biosynthesis, and TFs were used to build an unsupervised co-expression based similarity matrix using Pearson’s correlation coefficient. Then the R package WGCNA version 1.35 was used to create the networks , which were modeled with Cytoscape 3.2 .
Quantitative real time PCR (qRT-PCR)
qRT-PCR was performed as described , using the TaKaRa ExTaq R PCR Kit, SYBR green dye (TaKaRa, Dalian, China) and a DNA Engine Opticon 2 machine (MJ Research, Waltham, MA). Fifteen genes including cellulose and lignin biosynthesis genes (Pt-CESA2.1, Pt-ATH.2, Pt-GLAC90.1, Pt-PRX1.8, Pt-PAL1.2 and Pt-PAL1.3) were validated, and the primers are shown in Additional file 1. The efficiency of the primers was calculated by performing real-time PCR on several dilutions of first-strand cDNAs. Efficiencies of the different primer sets were similar. The specificity of each primer set was checked by sequencing PCR products. The reactions were carried out in a 20 μl volume containing 2 μl of diluted cDNA, 200 nM of each primer, and PCR Master Mix with the following conditions: 95°C for 30 s, and 45 cycles of 95°C for 5 s, 58°C for 15 s, and 72°C for 20 s. Then, a thermal denaturing cycle of 95°C for 15 s and 60°C for 1 min was applied to determine the dissociation curves, which were used to verify the specificity of PCR amplifications. All reactions were run in triplicate for each sample. Relative expression levels of candidate genes were calculated by the 2−ΔCt method. The results obtained for the different tissues were analyzed and standardized to the mRNA levels of poplar ACTINII-like (Accession number: EF145577), which shows stable expression.
Global transcriptome analysis of the RNA-seq data
Summary of RNA-seq for NW, TW, and OW
Number of raw reads
Number after intensity trimming
Number after quality trimming
Number after adaptor trimming
Mapping result statistics in NW, TW, and OW
Mapped pair reads
Mapped broken-pair reads
Mapped unique reads
Mapped multi reads
Functional annotation of differentially expressed genes
Notably, unlike the differentially expressed genes in NW vs OW, and NW vs TW, less information about GO terms was found for TW vs OW (Additional file 5). No significant terms were found in the up-regulated genes of TW vs OW. Only eight terms for “molecular function”, such as “cofactor binding” and “heme binding and tetrapyrrole binding”, were found in the down-regulated genes in TW vs OW (Additional file 5).
Differential expression of genes involved in stress responses and transcriptional regulation
In the “response to chemical stimulus”, “response to desiccation”, “response to abiotic stimulus”, and “response to oxidative stress” terms, 34 genes were up-regulated in NW compared with TW or OW (Additional file 6). Of these, only 15 genes were annotated, including 6 genes that encode peroxidases (POs). Importantly, our study identified a large proportion of differentially expressed genes (97) related to transcriptional regulation, based on the GO terms “transcription regulator activity”, “transcription factor activity”, and “transcription repressor activity” (Additional file 6). In particular, all the putative transcriptional regulators were down-regulated in NW compared with TW and OW and included genes encoding TFs of the families: MYB, NAC, heat shock TF (HSF), AP2 domain-containing TF family protein, bZIP TF, basic helix-loop-helix (bHLH) family protein, and WRKY TF (Additional file 6). These genes also showed strong changes in expression level; of these genes, 39 and 50 genes were down-regulated more than 4-fold in NW compared with TW or OW, respectively, and 6 genes were down-regulated more than 3-fold in TW compared with OW (Additional file 6). Of these 97 TFs, six and ten TF genes involved in lignin and cellulose biosynthesis were differentially expressed in TW, OW, and NW, including LIM, MYB, VNI, ASL, and KNAT family members. Although both Pt-LIM1 (NW > OW > TW) and Pt-LIM2 (NW > TW > OW) were up-regulated in NW, Pt-LIM1 was up-regulated in OW, and Pt-LIM2 was down-regulated. The four MYB TF genes (Pt-MYB090, Pt-MYB156, Pt-MYB167 and Pt-MYB221) were all up-regulated in OW compared with TW (OW > TW). All the genes encoding TFs that belong to VNI, ASL, and KNAT family members showed the same expression pattern (NW > TW > OW).
Differential expression of genes involved in hormone responses
In total, 40 genes related to plant hormone signal transduction showed different transcript levels in TW, OW, and NW (Additional file 6), including genes related to auxin, abscisic acid (ABA), ethylene, cytokinin, gibberellin, brassinosteroids, jasmonic acid (JA), and salicylic acid. Of these, eight auxin-related genes were up-regulated in NW compared with TW or OW, and two were down-regulated, Pt-SAUR1 encoding a SAUR family protein and POPTR_0005s25800 encoding an auxin response factor (Additional file 6). The ABA-related genes encoding serine/threonine-protein kinase (SRK2), protein phosphatase 2C (PP2C), ABA responsive element binding factor, and ABA receptor PYR/PYL family member were all down-regulated in NW compared with TW or OW (Additional file 6). Similar to ABA-related genes, ethylene- and JA-related genes encoding ethylene-insensitive protein, ethylene receptor, MYC2, JA-amino synthetase, and a jasmonate ZIM domain-containing protein were all down-regulated in NW compared with TW or OW (Additional file 6). Most plant hormone signal transduction genes were differentially expressed among TW, NW, and OW (Additional file 6), implying that genes related to hormone signal transduction are involved in the development of reaction wood. Also, three genes involved in hormone biosynthesis, such as Pt-ACO1.3 related to ethylene, cytokinin oxidase (POPTR_0001s05830) related to cytokinin, and gibberellin 2-oxidase (Pt-GA2.7) related to gibberellin, were all down-regulated in TW and NW compared with OW (Additional file 7). Similarly, Pt-PAP2.3 encoding an auxin-responsive protein and POPTR_0018s06080 encoding an ethylene-responsive transcriptional coactivator were repressed in TW compared with OW (Additional file 7).
Differential expression of genes involved in lignin, flavonoid, and cellulose biosynthesis
The differentially-expressed transcripts involved in lignin biosynthesis
log 2 FC(NW/OW)
log 2 FC(NW/TW)
Among the differentially expressed genes, 43 genes related to cell wall metabolism were identified, including cellulose synthase (CesA) (Additional file 6). The majority of these genes were up-regulated in NW compared with TW or OW (Additional file 6). Of these, the cell wall structural protein gene POPTR_0008s01310 encoding a fasciclin-like domain-containing protein, was down-regulated over 8-fold in TW and OW. Among the genes related to the cell wall, we also found two differentially expressed genes encoding sucrose synthase (SUS), Pt-SUS1 (NW > OW > TW) and Pt-SUS2 (OW > TW > NW), and their transcript levels indicated that different SUS genes have distinct expression patterns in NW and TW. In total, 11 and 2 genes encoding CesA and COBRA were found to be differently transcribed, such as Pt-CESA2.1, which was up-regulated over 26- and 68-fold in NW compared with TW and OW, respectively (Additional file 6). Also, transcripts of Pt-CESA2.6, encoding a protein similar to CesA7A-like, increased over 4-fold in NW compared with TW and OW. Although the transcripts of most of these CesA genes were more abundant in NW than TW or OW, the Pt-MANS.1 and Pt-MANS.2 transcripts were most abundant in TW (TW > NW > OW). POPTR_0014s01060 encoding a protein similar to an Arabidopsis thaliana rhamnogalacturonate lyase family protein was up-regulated 1,325 fold in NW compared with OW and did not express in TW, or had so little message that it was undetectable (Additional file 6). One of the two genes encoding COBRA proteins, which are involved in cell wall expansion and/or cellulose deposition, was up-regulated in TW, and the other was down-regulated compared with OW. In addition, 14 genes encoding proteins related to transporters were identified, including a gene encoding a sucrose transporter up-regulated in TW, and two aquaporin genes, Pt-TIP2.4 and Pt-PIP2.3. In particular, nine genes encoding proteins related to cell expansion (expansin-like and expansin-related protein) and cell wall modification (xyloglucan endotransglycosylases, XTHs), S-adenosylmethionine (SAM) were also differentially expressed, including Pt-EXLB1.1, XTH, and Pt-SAM1.2 (Additional file 7).
Co-expression networks underlying the regulation of lignin and cellulose biosynthesis
The cellulose biosynthesis network has 13 cellulose synthesis genes, including CesA and COBRA related genes, and ten TF genes orthologous to Arabidopsis VNI2 (TW > OW), MYB103 (TW > OW), ASL11 (TW > OW), and KNAT3 (OW > TW) (Figure 7). Since most of the cellulose synthesis genes were up-regulated in TW compared with OW (except POPTR_0002s25970, POPTR_0001s12940, POPTR_0004s11620, and Pt-CESA4.1), KNAT3 may have a negative role in regulating cellulose biosynthesis (Figures 5 and 7). VNI2, MYB103, and ASL11 were likely to be co-expressed with the up-regulated cellulose synthesis genes, indicating they might have a positive role in cellulose biosynthesis in TW (Figures 5 and 7).
qRT-PCR validation of RNA-seq data
Extensive transcriptome remodeling underlies the drastic differences between NW and branch wood
Earlier studies showed that reaction wood, abnormal wood induced by mechanical bending stress or gravity and other natural environmental factors, has many differences with NW [2,13]. Many studies have also uncovered the differences in expression of genes between CW or TW and NW [13,28,29]. These studies revealed important discoveries but also had limitations, in particular because they only considered the differences between NW and TW or between TW and OW [13,28,29]. In this paper, we identified genes differentially transcribed among NW, TW and OW and annotated their various functions in cell division, cell expansion, primary wall synthesis, secondary wall deposition, hormone signaling, transcription, and environmental stresses. Previous studies found drastic differences between NW and TW or TW and OW of Populus, indicating that gravity affects cell division, secondary wall deposition, cellulose microfibril orientation, and overall wood properties [2,30,31]. The differences between the TW and OW transcriptomes observed in our study were mostly in agreement with previous data derived from bent trunks of Populus . Importantly, we identified 2,658 and 2,417 differentially expressed genes in OW and TW, respectively, in comparison to NW. That shows that more genes were up- or down-regulated in OW than in TW. Thus, these findings support the idea that OW experienced more transcriptome remodeling than TW, resulting in greater phenotypic variation between OW and wood formed in normal conditions (NW), compared to the variation between TW and NW. Our study also detected many metabolic and biosynthetic processes among the significantly enriched GO terms for up-regulated genes in NW compared with TW or OW, suggesting that these processes were mainly repressed in TW or OW. Moreover, all of these results may indicate that the normal biological processes of TW and OW were dramatically affected, compared with NW. Interestingly, the GO terms analysis of differentially expressed genes indicates that the down-regulated genes in NW vs OW and NW vs TW function in similar processes (Additional file 5). In other words, the genes up-regulated in TW or OW compared with NW are more likely to have the same functions.
Another interesting discovery in our study was that the “molecular function” terms “transcription regulator activity” and “transcription factor activity” were identified among the down-regulated genes in NW compared with TW or OW. Also, more than half of the differentially expressed TF genes decreased over 4-fold in NW, indicating that TF abundances differed dramatically in these tissues; these TFs may have important roles in reaction wood (Additional file 6). Furthermore, the term “transcription repressor activity” was detected among the genes down-regulated in NW compared with TW, suggesting that repressors were up-regulated in TW and may repress the expression of their target genes . Previous work found many differentially regulated TFs in developing TW in P. tremula (L.) × tremuloides (Michx.) TW-forming tissues induced by bending field-grown trees compared to normal wood . Since TFs can recruit coactivator or corepressor proteins to a transcription factor DNA complex , the TFs up-regulated in reaction wood may play positive or negative roles in reaction wood formation. In addition, our results indicated that several classes of genes encoding HSPs were mostly up-regulated in TW. HSPs maintain protein homeostasis inside cells and promote the proper folding, stability, and degradation of polypeptides . Moreover, previous studies have suggested a role for HSPs in early xylem development in wood in P. tomentosa, thereby emphasizing the importance of HSP genes as candidate genes for the genetic improvement of wood quality . Extended cell wall thickening and a delay of programmed cell death were observed with the proteins, including HSPs . Here, the genes encode HSPs were up-regulated in TW, which may reflect stronger cell wall thickening and delays of programmed cell death in TW, compared with OW.
Ehylene and gibberellin have potential regulatory roles in reaction wood formation
The plant hormones have important roles in wood development and have been studied for decades [10,35]. A previous study used microarrays to examine the differentially transcribed genes in CW and OW in radiate pine branches and showed that many genes related to hormones and calcium signaling, as well as various environmental stresses, were exclusively up-regulated in CW . For Populus, the up-regulation of 1-Aminocyclopropane-1-carboxylate oxidase (ACO) in TW compared with NW was detected in hybrid aspen . In our study, Pt-ACO7, encoding a putative ACO, was down-regulated in NW compared with either TW or OW and was down-regulated in TW compared with OW (Additional file 7). It could be proposed that the abundance of ACO is highest in OW, then TW and NW (OW > TW > NW), which may indicate that ethylene participates in TW formation. Ethylene interacts synergistically with gibberellins and cytokinin in TW formation in Leucaena leucocephala  and the formation of TW in A. mangium seedlings requires gibberellin, especially gibberellin 2-oxidase . Consistent with this, Pt-GA2.7 encoding a protein similar to gibberellin 2-oxidase was up-regulated to its highest levels in OW (Additional file 7). In summary, we have identified many genes related to hormone signaling and biosynthesis, as well as environmental stresses, as differentially transcribed in NW, TW, and OW. This provides valuable clues to improve our understanding of reaction wood formation in response to gravity stimulus. Overall, our results that show ethylene and gibberellin have potential regulatory roles in reaction wood formation.
Transcriptomic networks underlying the regulation of lignin and cellulose biosynthesis in NW, TW, and OW
Previous studies revealed that TW has more cellulose and less lignin than OW or NW [30,31]; in accordance, genes in the phenylpropanoid pathway for biosynthesis of lignin monomers and flavonoids were greatly down-regulated [5,29]. In our study, transcript levels of 33 genes encoding 11 enzymes of the lignin biosynthetic pathway, including CCoAOMT, PAL, C4H, CAM, CCR, C3H, 4CL, PO, CAld5H, CAD, and LAC, were lower in TW compared with OW (Table 3). This is consistent with previous work showing that poplar TW has lower expression of PAL, 4CL, and CAD . The genes encoding LACs, which are involved in the polymerization of monolignols to produce lignin macromolecules, were also significantly down-regulated in TW . Similar results were found in our study, as five genes encoding proteins like LAC were all down-regulated in TW compared with NW; this also agrees with a report showing that expression of these genes decreased during TW formation in Populus . Our results confirm that genes in the phenylpropanoid pathway affecting biosynthesis of both lignin monomers and flavonoid were greatly down-regulated in TW compared with OW. This could explain the lower lignin content in TW compared with OW.
In recent years, many studies have examined the genes involved in lignin biosynthesis , and the regulatory network underlying lignin biosynthesis is gradually being revealed to include TFs, specific expression of gene family members, and control by pathway intermediates . Previous studies reported that TFs including LIM [26,27] and MYB  family members could regulate the genes involved in lignin biosynthesis and affect the lignin content or lignin monomer ratio. Further study revealed that plant LIM domain proteins may act as transcriptional activators of lignin biosynthesis, activating the expression of PAL, 4CL, and CAD and/or function as actin binding and bundling proteins . Also, the expression pattern of the twelve duplicated poplar Ptr-LIM (P. tremula × P. alba) genes investigated by qRT-PCR showed that poplar Ptr-XLIM1a, Ptr-XLIM1b, and Ptr-WLIM1b were preferentially expressed in the secondary xylem, including TW and OW, suggesting a specific function in wood formation and mechanical stress response . According to the differential expression of two LIM TFs observed in our study, Pt-LIM1 and Pt-LIM2, it is tempting to speculate that they may have roles in wood formation and may also function in the response to mechanical stress (Figure 6). For the MYB family, a recent study reported that when Pt-MYB216 (P. tomentosa) was overexpressed, PAL4, 4CL5, C3H3, and CCR2 were up-regulated and CAld5H was down-regulated, but no obvious changes occurred in transcript levels of CCoAOMT1, COMT2, and CAD1, demonstrating that Pt-MYB216 is a transcriptional activator of lignin biosynthesis during secondary wall formation . In our study, four MYB TFs (Pt-MYB090, Pt-MYB156, Pt-MYB167, and Pt-MYB221), which probably regulate lignin biosynthesis, were all up-regulated in OW compared with TW. However, few prior studies have constructed a possible regulatory network to examine lignin biosynthesis. Here, our study provided a transcriptomic regulatory network including TFs and enzymes of the lignin biosynthetic pathway (Figure 6) via a co-expression network among TW, OW, and NW; this network may enable the analysis of the regulation of lignin biosynthesis. The co-expression analysis found that Pt-MYB090, Pt-MYB156, Pt-MYB167, and Pt-MYB221 may co-express with Pt-PAL1, Pt-PAL2, Pt-PAL3, Pt-PAL4, Pt-4CL3, Pt-4CL5, Pt-CCR2, and Pt-C3H3, and have a negative expression relationship with Pt-CAld5H1 and Pt-CAld5H2. Also, we found that Pt-LIM1 co-expressed with Pt-PAL1, Pt-PAL2, Pt-PAL3, Pt-PAL4, Pt-4CL3, Pt-4CL5, and Pt-HCT1 and Pt-LIM2 co-expressed with Pt-HCT6 and Pt-CAD1. Together with previous studies [25-27], these findings indicate that the transcriptomic regulatory networks (Figure 6) provided in our study reflect the lignin biosynthetic pathway regulated by TFs.
The gradual exploration of cellulose biosynthesis has identified more and more genes affected cellulose biosynthesis, including CesAs , KORRIGAN , SUS , COBRA , and various TFs [45,46]. Previous studies showed the factors (hormones, light, mechanical stimuli, nutrition, and interactions with the cytoskeleton) that influence cellulose deposition by affecting the levels of substrate, and the abundance of cellulose synthase complexes (CSCs). Different CesA have different functions and different expression patterns [41,47], especially in the primary and secondary cell wall. In the tension stress-responsive genes examined in aspen, CesA genes, namely PtrCesA1, PtrCesA2, and PtrCesA3, which are closely associated with secondary cell wall development in the xylem cells, also exhibited similar tension stress-responsive behavior . In another study of TW genes in Populus by microarray and metabolite analysis, the PttCesA8-2 and PttCesA3-2 transcripts were found to be up-regulated in TW compared with NW . Interestingly, the “cellulose metabolic process” and “cellulose biosynthetic process” terms were also found in to be enriched in up-regulated genes in NW compared with TW or OW in our study. However, for the “molecular function” terms, “cellulose synthase (UDP-forming) activity” and “cellulose synthase activity” were only found to be enriched among the up-regulated genes in NW vs OW, showing that the genes related to CesA were down-regulated in OW. Importantly, nine CesA related genes were up-regulated in NW, and two genes (Pt-MANS.1 and Pt-MANS.2) were up-regulated in TW in our study. Together with previous studies [5,42], this supports the idea that different CesAs might be present in the biosynthetic complexes of TW and NW. Similar to CesAs, SUSs have an important role in cellulose synthesis, and overexpression of a cotton SUS gene in poplar resulted in increased cellulose synthesis, supporting the idea of a direct connection between sucrose supply, sucrose breakdown, and cellulose production through SUS . Interestingly, different SUSs have different effects, since overexpression of a bean SUS sequence in poplar did not increase the amount of cellulose , and another study in pea indicated that different SUSs are associated with different metabolic fates of sucrose . Consequently, the two differentially expressed SUS genes, Pt-SUS1 (NW > OW > TW) and Pt-SUS2 (OW > TW > NW), identified in our study have opposite expression patterns, which may indicate that different SUSs have different functions in cellulose biosynthesis. COBRAs are involved in cell wall expansion and cellulose deposition by affecting cellulose microfibril orientation during plant morphogenesis . AtCOBL4 is required for cellulose biosynthesis in the secondary wall and was identified by the similarity of its expression to Arabidopsis AtCesAs implicated in secondary cell wall synthesis [52,53]. However, little is known about COBRAs, and the differentially expressed COBRA related genes POPTR_0004s11620 and POPTR_0015s07110 may be useful in determining the function of COBRAs. In summary, our study identified a series of differentially expressed genes involved in cellulose biosynthesis and suggested more complicated functions of different members.
Taking advantage of genome-wide changes, transcriptome profiling of wood formation has identified a number of TFs that are preferentially expressed in developing wood and may regulate gene expression [25,54,55]. Among these, a set of MYB TFs that are functional orthologs of Arabidopsis TFs has been shown to be involved in the regulation of secondary wall biosynthesis during wood formation in Populus [25,54,55]. In detail, PtrMYB3 and PtrMYB20, functional orthologs of Arabidopsis MYB46 and MYB83, can activate the biosynthetic pathways of cellulose, xylan, and lignin, suggesting that they regulate the biosynthesis of all three major wood components in poplar . Our study found a positive relationship between Pt-MYB103 expression and expression of most of the cellulose synthesis genes, suggesting a positive role for Pt-MYB103. For NAC domain TFs, the positive master regulator of secondary cell wall biosynthesis in Arabidopsis including the functionally redundant VASCULAR-RELATED NAC-DOMAIN7 (VND7) and VND6 TFs function in secondary cell wall formation in protoxylem and metaxylem [56,57]. Previous work revealed that these NAC domain TFs can directly activate the expression of secondary wall-specific biosynthetic genes and activate the expression of several downstream transcription factor genes that also directly regulate genes involved in the biosynthesis of secondary wall components [55,58]. Our study found that the three TF orthologs of Arabidopsis VNI2 co-express with most of the cellulose synthesis machinery genes. However, few studies have explored the role of NAC domain and KNAT TFs in Populus wood formation . Considered negative regulators, the Class II KNOX gene KNAT7 negatively regulates secondary wall formation in Arabidopsis and is functionally conserved in Populus [45,60]. The KNAT3 related genes (Figure 7) were co-expressed with the down-regulated cellulose synthesis machinery genes, which supports the idea that they have a negative role in TW cellulose biosynthesis. For trees, few advances have been achieved in understanding the regulation of cellulose biosynthesis; this requires more studies to uncover the molecular mechanism . The possible regulatory network of cellulose biosynthesis genes and TFs described here may help us to understand the complicated regulation of cellulose biosynthesis (Figure 7). Understanding the complex process of wood formation will be important for optimizing the use of wood as a renewable energy source. Although significant progress has been made in understanding the cellulose and lignin synthesis machinery , many fundamental aspects remain to be addressed. Wood formation or secondary cell wall formation is an exciting area of study with many challenges and opportunities. Uncovering the transcriptional regulators controlling wood biosynthesis will provide novel tools to alter the biosynthetic pathways of wood components based on our needs.
In addition, apart from the genes directly involved in lignin and cellulose biosynthesis, many genes related to cell wall also have important roles in reaction wood formation, particularly several expansin, XTH, and SAM genes, which were differentially transcribed in our study. The quantity of wood formation is largely related to cell division and expansion during primary cell wall development. Expansins are involved in cell expansion in all tissues of the plant and have been isolated from the secondary xylem . Differential expression of expansins in TW has been reported . In our study, the genes encoding expansin-like and expansin-related protein were differentially expressed between TW and OW, for example, Pt-EXLB1.1 was up-regulated in OW (Additional file 7). An earlier study confirmed that XTHs function during the formation of secondary cell walls of vascular tissues and are believed to be important regulators of primary wall expansion . Our study showed the up-regulation of XTH in TW compared with OW (Additional file 7). Differential transcription of these genes between TW and OW could provide a molecular explanation for the similar tracheid diameters (either radial or tangential directions) in these two types of wood. A previous study found the transcripts of SAM in the developing xylem and showed that SAM is important in methylation reactions during biosynthesis of coniferyl and sinapyl alcohols . Differences in the availability of SAM may affect wood quality by altering lignin content and composition . Recent studies have revealed that SAMs were up-regulated in CW [13,29] in conifer. However, we have little information on SAMs in Populus. In our study, Pt-SAM1.2 mRNA levels increased in the OW, which may point to a role for SAMs in reaction wood formation. All together, genes involved in primary wall modification, wood growth, and tracheid dimensions were differentially expressed in branch wood of Populus, providing more clues to the molecular mechanisms affecting wood formation.
This genome-wide transcriptome profiling of branch (TW and OW) and normal stem wood (NW) from P. tomentosa provides more accurate insights into the molecular basis of TW formation in response to gravity stimulus by identifying differentially expressed genes using RNA-Seq. In total, 3,058 significantly differentially expressed genes were identified, including genes involved in secondary cell wall structure and wall modification and composition, such as cellulose and lignin biosynthetic genes. Two transcriptomic networks that underlie lignin and cellulose biosynthesis regulation were constructed, including TFs, based on the co-expression network, which may provide clues for understanding the regulation of cellulose and lignin biosynthesis. However, many key differentially expressed genes identified in our study were poorly annotated, which impedes our understanding of the molecular mechanism of wood formation. Importantly, the high-resolution expression patterns presented here improve our understanding of the molecular mechanisms and expand our knowledge of reaction wood formation in Populus. These findings may have potential applications for the improvement of wood properties in plants via genetic engineering.
Availability of supporting data
The Illumina reads have been deposited in the Sequence Read Archive at NCBI (http://www.ncbi.nlm.nih.gov/sra) and they are available under study accession number SRP040531. And the sequences of 44 genes involve in cellulose and lignin biosynthesis regulation pathways were uploaded in NCBI GenBank (http://www.ncbi.nlm.nih.gov/nuccore/) with accession number KP769977-KP770020.
This work was supported by the National ‘863’ Plan Project (No. 2013AA102702), the State Key Basic Research Program of China (No. 2012CB114506), the 111 Project (No. B13007), and the Fundamental Research Funds for the Central Universities (Project No. BLYJ201409).
- Plomion C, Leprovost G, Stokes A. Wood formation in trees. Plant Physiol. 2001;127:1513–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Kwon M. Review: tension wood as a model system to explore the carbon partitioning between lignin and cellulose biosynthesis in woody plants. J App Bio Chem. 2007;50:83–7.Google Scholar
- Jin H, Do J, Moon D, Noh EW, Kim W, Kwon M. EST analysis of functional genes associated with cell wall biosynthesis and modification in the secondary xylem of the yellow poplar (Liriodendron tulipifera) stem during early stage of tension wood formation. Planta. 2011;234:959–77.View ArticlePubMedGoogle Scholar
- Wang M, Qi X, Zhao S, Zhang S, Lu M-Z. Dynamic changes in transcripts during regeneration of the secondary vascular system in Populus tomentosa Carr. revealed by cDNA microarrays. BMC Genomics. 2009;10:215.View ArticlePubMed CentralPubMedGoogle Scholar
- Andersson-Gunnerås S, Mellerowicz EJ, Love J, Segerman B, Ohmiya Y, Coutinho PM, et al. 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–65.View ArticlePubMedGoogle Scholar
- Bao H, Li E, Mansfield SD, Cronk QC, El-Kassaby YA, Douglas CJ. The developing xylem transcriptome and genome-wide analysis of alternative splicing in Populus trichocarpa (black cottonwood) populations. BMC Genomics. 2013;14:359.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang X-H, Chiang VL. Molecular cloning of 4-coumarate: coenzyme A ligase in loblolly pine and the roles of this enzyme in the biosynthesis of lignin in compression wood. Plant Physiol. 1997;113:65–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Wagner A, Donaldson L, Kim H, Phillips L, Flint H, Steward D, et al. Suppression of 4-coumarate-CoA ligase in the coniferous gymnosperm Pinus radiata. Plant Physiol. 2009;149:370–83.View ArticlePubMed CentralPubMedGoogle Scholar
- Love J, Björklund S, Vahala J, Hertzberg M, Kangasjärvi J, Sundberg B. Ethylene is an endogenous stimulator of cell division in the cambial meristem of Populus. Proc Natl Acad Sci U S A. 2009;106:5984–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Nugroho WD, Yamagishi Y, Nakaba S, Fukuhara S, Begum S, Marsoem SN, et al. Gibberellin is required for the formation of tension wood and stem gravitropism in Acacia mangium seedlings. Ann Bot. 2012;110:887–95.View ArticlePubMed CentralPubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticlePubMedGoogle Scholar
- Zhao S, Fung-Leung W-P, Bittner A, Ngo K, Liu X. Comparison of RNA-Seq and microarray in transcriptome profiling of activated T cells. PLoS One. 2014;9:e78644.View ArticlePubMed CentralPubMedGoogle Scholar
- Li X, Yang X, Wu HX. Transcriptome profiling of radiata pine branches reveals new insights into reaction wood formation with implications in plant gravitropism. BMC Genomics. 2013;14:768.View ArticlePubMed CentralPubMedGoogle Scholar
- Chang S, Puryear J, Cairney J. A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993;11:113–6.View ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Sjödin A, Street NR, Sandberg G, Gustafsson P, Jansson S. The Populus Genome Integrative Explorer (PopGenIE): a new resource for exploring the Populus genome. New Phytol. 2009;182:1013–25.View ArticlePubMedGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.PubMedGoogle Scholar
- Saito R, Smoot ME, Ono K, Ruscheinski J, Wang PL, Lotia S, et al. A travel guide to Cytoscape plugins. Nat Methods. 2012;9:1069–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang D, Du Q, Xu B, Zhang Z, Li B. The actin multigene family in Populus: organization, expression and phylogenetic analysis. Mol Genet Genomics. 2010;284:105–19.View ArticlePubMedGoogle Scholar
- Geraldes A, Pang J, Thiessen N, Cezard T, Moore R, Zhao Y, et al. SNP discovery in black cottonwood (Populus trichocarpa) by population transcriptome resequencing. Mol Ecol Resour. 2011;11:81–92.View ArticlePubMedGoogle Scholar
- Petre B, Morin E, Tisserant E, Hacquard S, Da Silva C, Poulain J, et al. RNA-seq of early-infected poplar leaves by the rust pathogen Melampsora larici-populina uncovers PtSultr3;5, a fungal-induced host sulfate transporter. PLoS One. 2012;7:e44408.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu L, Filkov V, Groover A. Modeling transcriptional networks regulating secondary growth and wood formation in forest trees. Physiol Plant. 2014;151:156–63.View ArticlePubMedGoogle Scholar
- Wang S, Li E, Porth I, Chen JG, Mansfield SD, Douglas CJ. Regulation of secondary cell wall biosynthesis by poplar R2R3 MYB transcription factor PtrMYB152 in Arabidopsis. Sci Rep. 2014;4:5054.PubMed CentralPubMedGoogle Scholar
- Kawaoka A, Kaothien P, Yoshida K, Endo S, Yamada K, Ebinuma H. Functional analysis of tobacco LIM protein Ntlim1 involved in lignin biosynthesis. Plant J. 2000;22:289–301.View ArticlePubMedGoogle Scholar
- Rogers LA, Campbell MM. The genetic control of lignin deposition during plant growth and development. New Phytol. 2004;164:17–30.View ArticleGoogle Scholar
- Villalobos DP, Díaz-Moreno SM, Said E-SS, Cañas RA, Osuna D, Van Kerckhoven SH, et al. Reprogramming of gene expression during compression wood formation in pine: Coordinated modulation of S-adenosylmethionine, lignin and lignan related genes. BMC Plant Biol. 2012;12:100.View ArticlePubMed CentralPubMedGoogle Scholar
- Jin H, Kwon M. Mechanical bending-induced tension wood formation with reduced lignin biosynthesis in Liriodendron tulipifera. J Wood Sci. 2009;55:401–8.View ArticleGoogle Scholar
- Timell T. The chemical composition of tension wood. Svensk Papp Tidn. 1969;72:173–81.Google Scholar
- Pilate G, Déjardin A, Laurans F, Leplé JC. Tension wood as a model for functional genomics of wood formation. New Phytol. 2004;164:63–72.View ArticleGoogle Scholar
- Yamasaki K, Kigawa T, Seki M, Shinozaki K, Yokoyama S. DNA-binding domains of plant-specific transcription factors: structure, function, and evolution. Trends Plant Sci. 2013;18:267–76.View ArticlePubMedGoogle Scholar
- Xu L, Glass CK, Rosenfeld MG. Coactivator and corepressor complexes in nuclear receptor function. Curr Opin Genet Dev. 1999;9:140–7.View ArticlePubMedGoogle Scholar
- Tiroli-Cepeda AO, Ramos CH. An overview of the role of molecular chaperones in protein homeostasis. Protein Pept Lett. 2011;18:101–9.View ArticlePubMedGoogle Scholar
- Pramod S, Patel PB, Rao KS. Ethylene and its interaction with other hormones in tension wood formation in Leucaena leucocephala (Lam.) de Wit. J Ind Acad Wood Sci. 2012;9:130–9.View ArticleGoogle Scholar
- Andersson Gunnerås S. Wood formation and transcript analysis with focus on tension wood and ethylene biology. In: PhD thesis, Swedish university of agricultural sciences, department of forest genetics and plant physiology. 2005.Google Scholar
- Chen HC, Song J, Wang JP, Lin YC, Ducoste J, Shuford CM, et al. Systems biology of lignin biosynthesis in Populus trichocarpa: Heteromeric 4-Coumaric Acid:Coenzyme A Ligase protein complex formation, regulation, and numerical modeling. Plant Cell. 2014;26:876–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhao Q, Dixon RA. Transcriptional networks for lignin biosynthesis: more complex than we thought? Trends Plant Sci. 2011;16:227–33.View ArticlePubMedGoogle Scholar
- Arnaud D, Dejardin A, Leple JC, Lesage-Descauses MC, Boizot N, Villar M, et al. Expression analysis of LIM gene family in poplar, toward an updated phylogenetic classification. BMC Res Notes. 2012;5:102.View ArticlePubMed CentralPubMedGoogle Scholar
- Tian Q, Wang X, Li C, Lu W, Yang L, Jiang Y, et al. Functional characterization of the poplar R2R3-MYB transcription factor PtoMYB216 involved in the regulation of lignin biosynthesis during wood formation. PLoS One. 2013;8:e76369.View ArticlePubMed CentralPubMedGoogle Scholar
- Suzuki S, Li L, Sun YH, Chiang VL. The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in Populus trichocarpa. Plant Physiol. 2006;142:1233–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Bhandari S, Fujino T, Thammanagowda S, Zhang D, Xu F, Joshi CP. 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–37.View ArticlePubMedGoogle Scholar
- An X, Chen Z, Wang J, Ye M, Ji L, Liao W, et al. Identification and characterization of the Populus sucrose synthase gene family. Gene. 2014;539:58–67.View ArticlePubMedGoogle Scholar
- Liu L, Shang-Guan K, Zhang B, Liu X, Yan M, Zhang L, et al. Brittle Culm1, a COBRA-like protein, functions in cellulose assembly through binding cellulose microfibrils. PLoS Genet. 2013;9:22.Google Scholar
- Handakumbura PP, Hazen SP. Transcriptional regulation of grass secondary cell wall biosynthesis: playing catch-up with Arabidopsis thaliana. Front Plant Sci. 2012;3:74.View ArticlePubMed CentralPubMedGoogle Scholar
- Li E, Bhargava A, Qiang W, Friedmann MC, Forneris N, Savidge RA, et al. The Class II KNOX gene KNAT7 negatively regulates secondary wall formation in Arabidopsis and is functionally conserved in Populus. New Phytol. 2012;194:102–15.View ArticlePubMedGoogle Scholar
- Lei L, Li S, Gu Y. Cellulose synthase complexes: composition and regulation. Front Plant Sci. 2012;3:75.View ArticlePubMed CentralPubMedGoogle Scholar
- Coleman HD, Yan J, Mansfield SD. Sucrose synthase affects carbon partitioning to increase cellulose production and altered cell wall ultrastructure. Proc Natl Acad Sci U S A. 2009;106:13118–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Konishi T, Ohmiya Y, Hayashi T. Evidence that sucrose loaded into the phloem of a poplar leaf is used directly by sucrose synthase associated with various beta-glucan synthases in the stem. Plant Physiol. 2004;134:1146–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Barratt DH, Barber L, Kruger NJ, Smith AM, Wang TL, Martin C. Multiple, distinct isoforms of sucrose synthase in pea. Plant Physiol. 2001;127:655–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Roudier F, Fernandez AG, Fujita M, Himmelspach R, Borner GH, Schindelman G, et al. COBRA, an Arabidopsis extracellular glycosyl-phosphatidyl inositol-anchored protein, specifically controls highly anisotropic expansion through its involvement in cellulose microfibril orientation. Plant Cell. 2005;17:1749–63.View ArticlePubMed CentralPubMedGoogle Scholar
- 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:2281–95.View ArticlePubMed CentralPubMedGoogle Scholar
- 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 U S A. 2005;102:8633–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhong R, Lee C, Ye ZH. Evolutionary conservation of the transcriptional network regulating secondary cell wall biosynthesis. Trends Plant Sci. 2010;15:625–32.View ArticlePubMedGoogle Scholar
- McCarthy RL, Zhong R, Fowler S, Lyskowski D, Piyasena H, Carleton K, et al. The poplar MYB transcription factors, PtrMYB3 and PtrMYB20, are involved in the regulation of secondary wall biosynthesis. Plant Cell Physiol. 2010;51:1084–1090.56.View ArticlePubMedGoogle Scholar
- Yamaguchi M, Mitsuda N, Ohtani M, Ohme-Takagi M, Kato K, Demura T. VASCULAR-RELATED NAC-DOMAIN 7 directly regulates the expression of a broad range of genes for xylem vessel formation. Plant J. 2011;66:579–90.View ArticlePubMedGoogle Scholar
- Yamaguchi M, Kubo M, Fukuda H, Demura T. 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.View ArticlePubMedGoogle Scholar
- Zhong R, Lee C, Zhou J, McCarthy RL, Ye ZH. A battery of transcription factors involved in the regulation of secondary cell wall biosynthesis in Arabidopsis. Plant Cell. 2008;20:2763–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu R, Qi G, Kong Y, Kong D, Gao Q, Zhou G. Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa. BMC Plant Biol. 2010;10:145.View ArticlePubMed CentralPubMedGoogle Scholar
- Li E, Wang S, Liu Y, Chen JG, Douglas CJ. OVATE FAMILY PROTEIN4 (OFP4) interaction with KNAT7 regulates secondary cell wall formation in Arabidopsis thaliana. Plant J. 2011;67:328–41.View ArticlePubMedGoogle Scholar
- Gray-Mitsumune M, Mellerowicz EJ, Abe H, Schrader J, Winzéll A, Sterky F, et al. Expansins abundant in secondary xylem belong to subgroup A of the α-expansin gene family. Plant Physiol. 2004;135:1552–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Whetten R, Sun Y-H, Zhang Y, Sederoff R. Functional genomics and cell wall biosynthesis in loblolly pine. Plant Mol Biol. 2001;47:275–91.View ArticlePubMedGoogle Scholar
- Shen B, Li C, Tarczynski MC. High free-methionine and decreased lignin content result from a mutation in the Arabidopsis S-adenosyl-L-methionine synthetase 3 gene. Plant J. 2002;29:371–80.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.