NorWood: a gene expression resource for evo‐devo studies of conifer wood development

Summary The secondary xylem of conifers is composed mainly of tracheids that differ anatomically and chemically from angiosperm xylem cells. There is currently no high‐spatial‐resolution data available profiling gene expression during wood formation for any coniferous species, which limits insight into tracheid development. RNA‐sequencing data from replicated, high‐spatial‐resolution section series throughout the cambial and woody tissues of Picea abies were used to generate the NorWood.conGenIE.org web resource, which facilitates exploration of the associated gene expression profiles and co‐expression networks. Integration within PlantGenIE.org enabled a comparative regulomics analysis, revealing divergent co‐expression networks between P. abies and the two angiosperm species Arabidopsis thaliana and Populus tremula for the secondary cell wall (SCW) master regulator NAC Class IIB transcription factors. The SCW cellulose synthase genes (CesAs) were located in the neighbourhoods of the NAC factors in A. thaliana and P. tremula, but not in P. abies. The NorWood co‐expression network enabled identification of potential SCW CesA regulators in P. abies. The NorWood web resource represents a powerful community tool for generating evo‐devo insights into the divergence of wood formation between angiosperms and gymnosperms and for advancing understanding of the regulation of wood development in P. abies.

[1]  R. Zhong,et al.  Global analysis of direct targets of secondary wall NAC master switches in Arabidopsis. , 2010, Molecular plant.

[2]  Dirk Inzé,et al.  CORNET 2.0: integrating plant coexpression, protein-protein interactions, regulatory interactions, gene associations and functional annotations. , 2012, The New phytologist.

[3]  A. Barabasi,et al.  Network biology: understanding the cell's functional organization , 2004, Nature Reviews Genetics.

[4]  Partners in Time: EARLY BIRD Associates with ZEITLUPE and Regulates the Speed of the Arabidopsis Clock1[W][OA] , 2011, Plant Physiology.

[5]  Isabelle Giguére,et al.  Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses. , 2008, The New phytologist.

[6]  Yves Van de Peer,et al.  PLAZA 3.0: an access point for plant comparative genomics , 2014, Nucleic Acids Res..

[7]  A. Nekrutenko,et al.  Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences , 2010, Genome Biology.

[8]  G. Taylor,et al.  A transcriptomic approach to identify genes associated with wood density in Picea sitchensis , 2011 .

[9]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[10]  Weiping Song,et al.  Bioenergy and biofuels: History, status, and perspective , 2015 .

[11]  T. Lumley,et al.  gplots: Various R Programming Tools for Plotting Data , 2015 .

[12]  Douglas G. Scofield,et al.  The Norway spruce genome sequence and conifer genome evolution , 2013, Nature.

[13]  Monica C Munoz-Torres,et al.  Web Apollo: a web-based genomic annotation editing platform , 2013, Genome Biology.

[14]  Gary D. Bader,et al.  Cytoscape Web: an interactive web-based network browser , 2010, Bioinform..

[15]  A. Myburg,et al.  SND2, a NAC transcription factor gene, regulates genes involved in secondary cell wall development in Arabidopsis fibres and increases fibre cell area in Eucalyptus , 2011, BMC Plant Biology.

[16]  Tim Kovacs,et al.  TreeVector: Scalable, Interactive, Phylogenetic Trees for the Web , 2010, PloS one.

[17]  Nicholas J. Provart,et al.  An “Electronic Fluorescent Pictograph” Browser for Exploring and Analyzing Large-Scale Biological Data Sets , 2007, PloS one.

[18]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[19]  R. Funada,et al.  Differences in patterns of cell death between ray parenchyma cells and ray tracheids in the conifers Pinus densiflora and Pinus rigida , 2008, Trees.

[20]  H. Tuominen,et al.  A unique program for cell death in xylem fibers of Populus stem. , 2009, The Plant journal : for cell and molecular biology.

[21]  Ari Pekka Mähönen,et al.  Cytokinin and Auxin Display Distinct but Interconnected Distribution and Signaling Profiles to Stimulate Cambial Activity , 2016, Current Biology.

[22]  T. Demura,et al.  VASCULAR-RELATED NAC-DOMAIN7 directly regulates the expression of a broad range of genes for xylem vessel formation. , 2011, The Plant journal : for cell and molecular biology.

[23]  H. Fukuda,et al.  TDIF Peptide Signaling Regulates Vascular Stem Cell Proliferation via the WOX4 Homeobox Gene in Arabidopsis[W] , 2010, Plant Cell.

[24]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[25]  B. Sundberg,et al.  Auxin as a positional signal in pattern formation in plants. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[26]  David M. Goodstein,et al.  Phytozome: a comparative platform for green plant genomics , 2011, Nucleic Acids Res..

[27]  A. Séguin,et al.  Large-scale screening of transcription factor–promoter interactions in spruce reveals a transcriptional network involved in vascular development , 2014, Journal of experimental botany.

[28]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[29]  Victor X. Jin,et al.  Computational analysis reveals a correlation of exon-skipping events with splicing, transcription and epigenetic factors , 2013, Nucleic acids research.

[30]  Xinguo Li,et al.  Transcriptome profiling of Pinus radiata juvenile wood with contrasting stiffness identifies putative candidate genes involved in microfibril orientation and cell wall mechanics , 2011, BMC Genomics.

[31]  Y. van de Peer,et al.  The Plant Genome Integrative Explorer Resource: PlantGenIE.org. , 2015, The New phytologist.

[32]  B. Sundberg,et al.  A transcriptional roadmap to wood formation , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[33]  K. Shinozaki,et al.  The NAC Transcription Factors NST1 and NST2 of Arabidopsis Regulate Secondary Wall Thickenings and Are Required for Anther Dehiscencew⃞ , 2005, The Plant Cell Online.

[34]  Stefan R. Henz,et al.  A gene expression map of Arabidopsis thaliana development , 2005, Nature Genetics.

[35]  Kengo Kinoshita,et al.  ATTED-II in 2016: A Plant Coexpression Database Towards Lineage-Specific Coexpression , 2015, Plant & cell physiology.

[36]  Steven J. M. Jones,et al.  Insights into Conifer Giga-Genomes1 , 2014, Plant Physiology.

[37]  J. Bohlmann,et al.  Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression , 2012, BMC Genomics.

[38]  P. Rigault,et al.  A White Spruce Gene Catalog for Conifer Genome Analyses1[W][OA] , 2011, Plant Physiology.

[39]  Torgeir R. Hvidsten,et al.  Guidelines for RNA-Seq data analysis , 2014 .

[40]  T. Demura,et al.  Vascular-related NAC-DOMAIN7 is involved in the differentiation of all types of xylem vessels in Arabidopsis roots and shoots. , 2008, The Plant journal : for cell and molecular biology.

[41]  J. Grima-Pettenati,et al.  Conifer R2R3-MYB transcription factors: sequence analyses and gene expression in wood-forming tissues of white spruce (Picea glauca) , 2007, BMC Plant Biology.

[42]  I. Longden,et al.  EMBOSS: the European Molecular Biology Open Software Suite. , 2000, Trends in genetics : TIG.

[43]  K. Shinozaki,et al.  NAC Transcription Factors, NST1 and NST3, Are Key Regulators of the Formation of Secondary Walls in Woody Tissues of Arabidopsis[W][OA] , 2007, The Plant Cell Online.

[44]  David Levine,et al.  A high-performance computing toolset for relatedness and principal component analysis of SNP data , 2012, Bioinform..

[45]  B. Usadel,et al.  PlaNet: Combined Sequence and Expression Comparisons across Plant Networks Derived from Seven Species[W][OA] , 2011, Plant Cell.

[46]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[47]  L. Rieseberg,et al.  Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii). , 2014, The New phytologist.

[48]  Jason A. Corwin,et al.  An Arabidopsis Gene Regulatory Network for Secondary Cell Wall Synthesis , 2014, Nature.

[49]  M. Sipi,et al.  Distributions of tracheid cross-sectional dimensions in different parts of Norway spruce stems , 2008 .

[50]  R. Zhong,et al.  A Battery of Transcription Factors Involved in the Regulation of Secondary Cell Wall Biosynthesis in Arabidopsis , 2008, The Plant Cell Online.

[51]  Nobutaka Mitsuda,et al.  Reconstitution of a Secondary Cell Wall in a Secondary Cell Wall-Deficient Arabidopsis Mutant , 2014, Plant & cell physiology.

[52]  L. MacNeil,et al.  Gene regulatory networks and the role of robustness and stochasticity in the control of gene expression. , 2011, Genome research.

[53]  J. Grima-Pettenati,et al.  Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis , 2008, Journal of experimental botany.

[54]  Jin-Hua Ran,et al.  Evolution and biogeography of gymnosperms. , 2014, Molecular phylogenetics and evolution.

[55]  Isabelle Giguére,et al.  Modular organization of the white spruce (Picea glauca) transcriptome reveals functional organization and evolutionary signatures , 2015, The New phytologist.

[56]  J. Paiva,et al.  Plasticity of maritime pine (Pinus pinaster) wood-forming tissues during a growing season. , 2008, The New phytologist.

[57]  Francisco R. Cantón,et al.  Reprogramming of gene expression during compression wood formation in pine: Coordinated modulation of S-adenosylmethionine, lignin and lignan related genes , 2012, BMC Plant Biology.

[58]  C. R. Becer,et al.  Lignocellulosic biomass: a sustainable platform for the production of bio-based chemicals and polymers , 2015, 1602.01684.

[59]  P. Nilsson,et al.  A High-Resolution Transcript Profile across the Wood-Forming Meristem of Poplar Identifies Potential Regulators of Cambial Stem Cell Identityw⃞ , 2004, The Plant Cell Online.

[60]  J. Collins,et al.  Large-Scale Mapping and Validation of Escherichia coli Transcriptional Regulation from a Compendium of Expression Profiles , 2007, PLoS biology.

[61]  Tetsuro Mimura,et al.  Transcription switches for protoxylem and metaxylem vessel formation. , 2005, Genes & development.

[62]  P. Nilsson,et al.  A genomic approach to investigate developmental cell death in woody tissues of Populus trees , 2005, Genome Biology.

[63]  R. Zhong,et al.  Two NAC domain transcription factors, SND1 and NST1, function redundantly in regulation of secondary wall synthesis in fibers of Arabidopsis , 2007, Planta.

[64]  G. Berlyn,et al.  Botanical Microtechnique and Cytochemistry , 1991 .

[65]  A. Myburg,et al.  Navigating the transcriptional roadmap regulating plant secondary cell wall deposition , 2013, Front. Plant Sci..

[66]  Rodrigo Lopez,et al.  Analysis Tool Web Services from the EMBL-EBI , 2013, Nucleic Acids Res..

[67]  Xinguo Li,et al.  Seasonal reorganization of the xylem transcriptome at different tree ages reveals novel insights into wood formation in Pinus radiata. , 2010, The New phytologist.

[68]  C. Nairn,et al.  Three loblolly pine CesA genes expressed in developing xylem are orthologous to secondary cell wall CesA genes of angiosperms. , 2005, The New phytologist.

[69]  Carsten O. Daub,et al.  Estimating mutual information using B-spline functions – an improved similarity measure for analysing gene expression data , 2004, BMC Bioinformatics.

[70]  Hajime Ohyanagi,et al.  Plant Omics Data Center: An Integrated Web Repository for Interspecies Gene Expression Networks with NLP-Based Curation , 2014, Plant & cell physiology.

[71]  O. Gascuel,et al.  New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. , 2010, Systematic biology.

[72]  Staffan Persson,et al.  The cell biology of cellulose synthesis. , 2014, Annual review of plant biology.

[73]  Ge Gao,et al.  PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants , 2016, Nucleic Acids Res..

[74]  T. Demura,et al.  Multiple classes of transcription factors regulate the expression of VASCULAR-RELATED NAC-DOMAIN7, a master switch of xylem vessel differentiation. , 2015, Plant & cell physiology.

[75]  T. Demura,et al.  Contribution of NAC Transcription Factors to Plant Adaptation to Land , 2014, Science.

[76]  David Sundell,et al.  ComPlEx: conservation and divergence of co-expression networks in A. thaliana, Populus and O. sativa , 2014, BMC Genomics.

[77]  M. Mutwil,et al.  Tools of the trade: studying molecular networks in plants. , 2016, Current opinion in plant biology.

[78]  H. Endo,et al.  NAC-MYB-based transcriptional regulation of secondary cell wall biosynthesis in land plants , 2015, Front. Plant Sci..

[79]  C. Ritland,et al.  Global Transcriptome Analysis of Constitutive Resistance to the White Pine Weevil in Spruce , 2011, Genome biology and evolution.