Transcriptome Analysis of Reaction Wood in Gymnosperms by Next-Generation Sequencing


Special xylem tissue called “compression wood” is formed on the lower side of inclined stems when gymnosperms grow on a slope. We investigated the molecular mechanism of compression wood formation. Transcriptome analysis by next-generation sequencing (NGS) was applied to the xylem of Chamaecyparis obtusa to develop a catalog of general gene expression in differentiating xylem during compression and normal wood formation. The sequencing output generated 234,924,605 reads and 40,602 contigs (mean size = 529 bp). Based on a sequence similarity search with known proteins, 54.2% (22,005) of the contigs showed homology with sequences in the databases. Of these annotated contigs, 19,293 contigs were assigned to Gene Ontology categories. Differential gene expression between the compression and normal wood libraries was analyzed by mapping the reads from each library to the assembled contigs. In total, 2875 contigs were identified as differentially expressed, including 1207 that were up-regulated and 1668 that were down-regulated in compression wood. We selected 30 genes and compared the transcript abundance between compression and normal wood by quantitative polymerase chain reaction analysis to validate the NGS results. We found that 27 of the 30 genes showed the same expression patterns as the original NGS results.

Share and Cite:

Sato, S. , Yoshida, M. , Hiraide, H. , Ihara, K. and Yamamoto, H. (2014) Transcriptome Analysis of Reaction Wood in Gymnosperms by Next-Generation Sequencing. American Journal of Plant Sciences, 5, 2785-2798. doi: 10.4236/ajps.2014.518295.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Timell, T.E. (1986) Compression Wood in Gymnosperms. 1, Springer-Verlag, Berlin.
[2] Yamamoto, H., Okuyama, T., Yoshida, M. and Sugiyama, K. (1991) Generation Process of Growth Stresses in Cell Walls III. Growth Stress in Compression Wood. MokuzaiGakkaishi, 37, 94-100.
[3] Plomion, C., Pionneau, C., Brach, J., Costa, P. and Baillères, H. (2000) Compression Wood-Responsive Proteins in Developing Xylem of Maritime Pine (Pinuspinaster Ait.). Plant Physiology, 123, 959-969.
[4] Gion, J.M., Lalanne, C., Le Provost, G., Ferry-Dumazet, H., Paiva, J., Chaumeil, P., Frigerio, J.M., Brach, J., Barré, A., de Daruvar, A., Claverol, S., Bonneu, M., Sommerer, N., Negroni, L. and Plomion, C. (2005) The Proteome of Maritime Pine Wood Forming Tissue. Proteomics, 5, 3731-3751.
[5] Yamashita, S., Yoshida, M., Yamamoto, H. and Okuyama, T. (2008) Screening Genes That Change Expression during Compression Wood Formation in Chamaecyparis obtusa. Tree Physiology, 28, 1331-1340.
[6] Ramos, P., Le Provost, G., Gantz, C., Plomion, C. and Herrera, R. (2012) Transcriptional Analysis of Differentially Expressed Genes in Response to Stem Inclination in Young Seedlings of Pine. Plant Biology, 14, 923-933.
[7] Villalobos, D.P., Díaz-Moreno, S.M., Said, E.S.S., Canas, R.A., Osuna, D., Van Kerckhoven, S.H.E., Bautista, R., Claros, M.G., Canovas, F.M. and Canton, F.R. (2012) Reprogramming of Gene Expression during Compression Wood Formation in Pine: Coordinated Modulation of S-Adenosylmethionine, Lignin and Lignin Related Genes. BMC Plant Biology, 12, 100.
[8] Marioni, J.C., Mason, C.E., Mane, S.M., Stephens, M. and Gilad, Y. (2008) RNA-seq: An Assessment of Technical Reproducibility and Comparison with Gene Expression Arrays. Genome Research, 18, 1509-1517.
[9] Agarwal, A., Koppstein, D., Rozowsky, J., Sboner, A., Habegger, L., Hillier, L.W., Sasidharan, R., Reinke, V., Waterston, R.H. and Gerstein, M. (2010) Comparison and Calibration of Transcriptome Data from RNA-Seq and Tiling Arrays. BMC Gemomics, 11, 383.
[10] Mizrachi, E., Hefer, C.A., Ranik, M., Joubert, F. and Myburg, A.A. (2010) De novo Assembled Expressed Gene Catalog of a Fast-Growing Eucalyptus Tree Produced by Illumina mRNA-Seq. BMC Genomics, 11, 681.
[11] Hao, D.C., Ge, G.B., Xiao, P.G., Zhang, Y.Y. and Yang, L. (2011) The First Insight into the Tissue Specific Taxus Transcriptome via Illumina Second Generation Sequencing. PLoS ONE, 6, Article ID: e21220.
[12] Li, D.J., Deng, Z., Qin, B., Liu, X.H. and Men, Z.H. (2012) De novo Assembly and Characterization of Bark Transcriptome Using Illumina Sequencing and Development of EST-SSR Markers in Rubber Tree (Hevea brasiliensis Muell. Arg.). BMC Genomics, 13, 192.
[13] Sun, X.D., Zhou, S.M., Meng, F.L. and Liu, S.Q. (2012) De novo Assembly and Characterization of the Garlic (Allium sativum) Bud Transcriptome by Illumina Sequencing. Plant Cell Reports, 31, 1823-1828.
[14] Dong, Q., Schlueter, S.D. and Brendel, V. (2004) PlantGDB, Plant Genome Database and Analysis Tools. Nucleic Acids Research, 32, D354-D359.
[15] Conesa, A., Gotz, S., Garcia-Gomez, J.M., Terol, J., Talon, M. and Robles, M. (2005) Blast2GO: A Universal Tool for Annotation, Visualization and Analysis in Functional Genomics Research. Bioinformatics, 21, 3674-3676.
[16] Iseli, C., Jongeneel, C.V. and Bucher, P. (1999) ESTScan: A Program for Detecting, Evaluating and Reconstructing Potential Coding Regions in EST Sequences. Proceedings of the 7th International Conference on Intelligent Systems for Molecular Biology, Heidelberg, 6-10 August 1999, 138-148.
[17] Rice, P., Longden, I. and Bleasby, A. (2000) EMBOSS: The European Molecular Biology Open Software Suite. Trends in Genetics, 16, 276-277.
[18] Tatusov, R.L., Fedorova, N.D., Jackson, J.D., Jacobs, A.R., Kiryutin, B., Koonin, E.V., Krylov, D.M., Mazumder, R., Mekhedov, S.L., Nikolskaya, A.N., Rao, B.S., Smirnov, S., Sverdlov, A.V., Vasudevan, S., Wolf, Y.I., Yin, J.J. and Natale, D.A. (2003) The COG Database: An Updated Version Includes Eukaryotes. BMC Bioinformatics, 4, 41.
[19] Kanehisa, M. and Goto, S. (2000) KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research, 28, 27-30.
[20] Wu, S.T., Zhu, Z.W., Fu, L.M., Niu, B.F. and Li, W.Z. (2011) WebMGA: A Customizable Web Server for Fast Metagenomic Sequence Analysis. BMC Genomics, 12, 444.
[21] Sun, J., Nishiyama, T., Shimizu, K. and Kadota, K. (2013) TCC: An R Package for Comparing Tag Count Data with Robust Normalization Strategies. BMC Bioinformatics, 14, 219.
[22] Untergasser, A., Nijveen, H., Rao, X., Bisseling, T., Geurts, R. and Leunissen, J.A. (2007) Primer3Plus, an Enhanced Web Interface to Primer3. Nucleic Acids Research, 35, W71-W74.
[23] Livak, K.J. and Schmittgen, T.D. (2001) Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2(-Delta delta C(T)) Method. Methods, 25, 402-408.
[24] Ujino-Ihara, T., Kanamori, H., Yamane, H., Taguchi, Y., Namiki, N., Mukai, Y., Yoshimura, K. and Tsumura, Y. (2005) Comparative Analysis of Expressed Sequence Tags of Conifers and Angiosperms Reveals Sequence Specifically Conserved in Conifers. Plant Molecular Biology, 59, 895-907.
[25] Xia, Z., Xu, H., Zhai, J., Li, D., Luo, H., He, C. and Huang, X. (2011) RNA-Seq Analysis and de novo Transcriptome Assembly of Hevea brasiliensis. Plant Molecular Biology, 77, 299-308.
[26] Gordo, S.M.C., Pinheiro, D.G., Moreira, E.C.O., Rodrigues, S.M., Poltronien, M.C., De Lemos, O.F., Tojal da Silva, I., Ramos, R.T.J., Silva, A., Schneider, H., Silva Jr., W.A., Sampaio, I. and Darnet, S. (2012) High-Throughput Sequencing of Black Pepper Root Transcriptome. BMC Plant Biology, 12, 168.
[27] Wang, S., Wang, X., He, Q., Liu, X., Xu, W., Li, L., Gao, J. and Wang, F. (2012) Transcriptome Analysis of the Roots at Early and Late Seedling Stages Using Illumina Paired-End Sequencing and Development of EST-SSR Markers in Radish. Plant Cell Reports, 31, 1437-1447.
[28] Bräutigam, A., Mullick, T., Schliesky, S. and Weber, A.P.M. (2011) Critical Assessment of Assembly Strategies for Non-Model Species mRNA-Seq Data and Application of Next-Generation Sequencing to the Comparison of C3 and C4 Species. Journal of Experimental Botany, 62, 3093-3102.
[29] Parchman, T.L., Geist, K.S., Grahnen, J.A., Benkman, C.W. and Buerkle, C.A. (2010) Transcriptome Sequencing in an Ecologically Important Tree Species: Assembly, Annotation, and Marker Discovery. BMC Genomics, 11, 180.
[30] Sui, C., Zhang, J., Wei, J.H., Chen, S.L., Li, Y., Xu, J.S., Jin, Y., Xie, C.X., Gao, Z.H., Chen, H.J., Yang, C.M., Zhang, Z. and Xu, Y.H. (2011) Transcriptome Analysis of Bupleurum chinese Focusing on Genes Involved in the Biosynthesis of Saikosaponins. BMC Genomics, 12, 539.
[31] Garg, R., Patel, R.K., Jhanwar, S., Priya, P., Bhattacharjee, A., Yadav, G., Bhatia, S., Chattopadhyay, D., Tyagi, A.K. and Jain, M. (2011) Gene Discovery and Tissue-Specific Transcriptome Analysis in Chickpea with Massively Parallel Pyrosequencing and Web Resource Development. Plant Physiology, 156, 1661-1678.
[32] Zhou, Y.J., Gao, F., Liu, R., Feng, J.C. and Li, H.J. (2012) De novo Sequencing and Analysis of Root Transcriptome Using 454 Pyrosequencing to Discover Putative Genes Associated with Drought Tolerance in Ammopiptanthus mongolicus. BMC Genomics, 13, 266.
[33] Marquez, Y., Brown, J.W.S., Simpson, C., Barta, A. and Kalyna, M. (2012) Transcriptome Survey Reveals Increased Complexity of the Alternative Splicing Landscape in Arabidopsis. Genome Research, 22, 1184-1195.
[34] Wang, X.W., Luan, J.B., Li, J.M., Bao, Y.Y., Zhang, C.X. and Liu, S.S. (2010) De novo Characterization of a Whitefly Transcriptome and Analysis of Its Gene Expression during Development. BMC Genomics, 11, 400.
[35] Dixon, R.A., Chen, F., Guo, D. and Parvathi, K. (2001) The Biosynthesis of Monolignols: A “Metabolic Grid”, or Independent Pathways to Guaiacyl and Syringyl Units? Phytochemistry, 57, 1069-1084.
[36] Bailly, A., Yang, H., Martinoia, E., Geisler, M. and Murphy, A.S. (2012) Plant Lessons: Exploring ABCB Functionality through Structural Modeling. Frontiers in Plant Science, 2, 108.
[37] Liu, C.J. (2012) Deciphering the Enigma of Lignification: Precursor Transport, Oxidation, and the Topochemistry of Lignin Assembly. Molecular Plant, 5, 304-317.
[38] Sanders, D., Pelloux, J., Brownlee, C. and Harper, J.F. (2002) Calcium at the Crossroads of Signaling. The Plant Cell, 14, S401-S417.
[39] Michalak, M., Milner, R.E., Burns, K. and Opas, M. (1992) Calreticulin. Biochemical Journal, 285, 681-692.
[40] Allona, I., Quinn, M., Shoop, E., Swope, K., St. Cyr, S., Carlis, J., Riedl, J., Retzel, E., Campbell, M.M., Sederoff, R. and Whetten, R.W. (1998) Analysis of Xylem Formation in Pine by cDNA Sequencing. Proceedings of the National Academy of Sciences of the United States of America, 95, 9693-9698.
[41] Du, S. and Yamamoto, F. (2003) A Study on the Role of Calcium in Xylem Development and Compression Wood Formation in Taxodium distichum Seedling. IAWA Journal, 24, 75-85.
[42] Kurusu, T., Kuchitsu, K., Nakano, M., Nakayama, Y. and Iida, H. (2013) Plant Mechanosensing and Ca2+ Transport. Trends in Plant Science, 18, 227-233.
[43] Kao, Y.L., Deavours, B.E., Phelps, K.K., Walker, R.A. and Reddy, A.S.N. (2000) Bundling of Microtubules by Motor and Tail Domains of a Kinesin-Like Calmodulin-Binding Protein from Arabidopsis: Regulation by Ca2+/Calmodulin. Biochemical and Biophysical Research Communications, 267, 201-207.
[44] Preuss, M.L., Delmer, D.P. and Liu, B. (2003) The Cotton Kinesin-Like Calmodulin-Binding Protein Associates with Cortical Microtubules in Cotton Fibers. Plant Physiology, 132, 154-160.
[45] Barnett, J.R. and Bonham, V.A. (2004) Cellulose Microfibril Angle in the Cell Wall of Wood Fibres. Biological Reviews, 79, 461-472.
[46] Bedon, F., Grima-Pettenati, J. and Mackay, J. (2007) Conifer R2R3-MYB Transcription Factors: Sequence Analyses and Gene Expression in Wood-Forming Tissues of White Spruce (Picea glauca). BMC Plant Biology, 7, 17.
[47] Li, X.G., Yang, X.H. and Wu, H.X. (2013) Transcriptome Profiling of Radiate Pine Branches Reveals New Insights into Reaction Wood Formation with Implications in Plant Gravitropism. BMC Genomics, 14, 768.
[48] Zhong, R.Q., Lee, C. and Ye, Z.H. (2010) Functional Characterization of Poplar Wood-Associated NAC Domain Transcription Factors. Plant Physiology, 152, 1044-1055.
[49] Choi, D., Kim, J.H. and Lee, Y. (2008) Expansins in Plant Development. Advances in Botanical Research, 47, 47-97.
[50] Sato, S., Hiraide, H., Yoshida, M. and Yamamoto, H. (2013) Changes in Xylem Tissue and Laccase Transcript Abundance Associated with Posture Recovery in Chamaecyparis obtusa Saplings Growing on an Incline. Functional Plant Biology, 40, 637-643.

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.