Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Abstract

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant–pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.

References

[1] Schrader, J. et al. Cambial meristem dormancy in trees involves extensive remodelling of the transcriptome. Plant J. 40, 173-87 (2004).

[2] Qiu, Z. et al. Genome-wide analysis reveals dynamic changes in expression of microRNAs during vascular cambium development in Chinese fir, Cunninghamia lanceolata. J. Exp. Bot. 66, 3041-54 (2015).

[3] Bhalerao, R., Nilsson, O. & Sandberg, G. Out of the woods: forest biotechnology enters the genomic era. Curr. Opin. Biotechnol. 14, 206-213 (2003).

[4] Baucher, M., El Jaziri, M. & Vandeputte, O. From primary to secondary growth: origin and development of the vascular system. J. Exp. Bot. 58, 3485-501 (2007).

[5] Lopes, S. T. et al. Phellem versus xylem: genome-wide transcriptomic analysis reveals novel regulators of cork formation in cork oak. Tree Physiol. 40, 129-141 (2020).

[6] Demura, T. & Fukuda, H. Transcriptional regulation in wood formation. Trends Plant Sci. 12, 64-70 (2007).

[7] Oliveira, J. M., Santarosa, E., Pillar, V. D. & Roig, F. A. Seasonal cambium activity in the subtropical rain forest tree Araucaria angustifolia. Trees-Struct. Funct. 23, 107-115 (2009).

[8] Ding, Q., Zeng, J. & He, X. Q. Deep sequencing on a genome-wide scale reveals diverse stage-specific microRNAs in cambium during dormancy-release induced by chilling in poplar. BMC Plant Biol. 14, 267 (2014).

[9] Gricar, J. et al. Plastic and locally adapted phenology in cambial seasonality and production of xylem and phloem cells in Picea abies from temperate environments. Tree Physiol. 34, 869-81 (2014).

[10] Law, J. A. & Jacobsen, S. E. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat. Rev. Genet. 11, 204-20 (2010).

[11] Rodrigues, J. A. & Zilberman, D. Evolution and function of genomic imprinting in plants. Genes Dev. 29, 2517-31 (2015).

[12] Kumar, S., Kumari, R., Sharma, V. & Sharma, V. Roles, and establishment, maintenance and erasing of the epigenetic cytosine methylation marks in plants. J. Genet. 92, 629-66 (2013).

[13] Du, J. M., Johnson, L. M., Jacobsen, S. E. & Patel, D. J. DNA methylation pathways and their crosstalk with histone methylation. Nat. Rev. Mol. Cell Biol. 16, 519-532 (2015).

[14] Gouil, Q. & Baulcombe, D. C. DNA methylation signatures of the plant chromomethyltransferases. PLoS Genet. 12, e1006526 (2016).

[15] Choi, Y. et al. DEMETER, a DNA glycosylase domain protein, is required for endosperm gene imprinting and seed viability in Arabidopsis. Cell 110, 33-42 (2002).

[16] Xiao, W. Y. et al. DNA methylation is critical for Arabidopsis embryogenesis and seed viability. Plant Cell. 18, 805-814 (2006).

[17] Narsai, R. et al. Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination. Genome Biol. 18, 172 (2017).

[18] Chen, W. et al. Comparative WGBS identifies genes that influence non-ripe phenotype in tomato epimutant Colourless non-ripening. Sci. China Life Sci. 61, 244-252 (2018).

[19] Jullien, P. E., Susaki, D., Yelagandula, R., Higashiyama, T. & Berger, F. DNA methylation dynamics during sexual reproduction in Arabidopsis thaliana. Curr. Biol. 22, 1825-1830 (2012).

[20] Ibarra, C. A. et al. Active DNA demethylation in plant companion cells reinforces transposon methylation in gametes. Science 337, 1360-1364 (2012).

[21] Park, K. et al. DNA demethylation is initiated in the central cells of Arabidopsis and rice. Proc. Natl Acad. Sci. USA 113, 15138-15143 (2016).

[22] Deleris, A., Halter, T. & Navarro, L. DNA methylation and demethylation in plant immunity. Annu. Rev. Phytopathol. 54, 579-603 (2016).

[23] Hossain, M. S. et al. Divergent cytosine DNA methylation patterns in single-cell, soybean root hairs. New Phytol. 214, 808-819 (2017).

[24] Secco, D. et al. Spatio-temporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery. Plant Cell. 25, 4285-304 (2013).

[25] Bouyer, D. et al. DNA methylation dynamics during early plant life. Genome Biol. 18, 179 (2017).

[26] Zemach, A. et al. The Arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin. Cell 153, 193-205 (2013).

[27] Mirouze, M. et al. Selective epigenetic control of retrotransposition in Arabidopsis. Nature 461, 427-30 (2009).

[28] Qiu, Z. et al. The regulation of cambial activity in Chinese fir (Cunninghamia lanceolata) involves extensive transcriptome remodeling. New Phytol. 199, 708-19 (2013).

[29] Xu, H. et al. Gene expression and proteomic analysis of shoot apical meristem transition from dormancy to activation in Cunninghamia lanceolata (Lamb.) Hook. Sci. Rep. 6, 19938 (2016).

[30] Chano, V., Collada, C. & Soto, A. Transcriptomic analysis of wound xylem formation in Pinus canariensis. BMC Plant Biol. 17, 234 (2017).

[31] Hao, Z., Liu, S., Hu, L., Shi, J. & Chen, J. Transcriptome analysis and metabolic profiling reveal the key role of carotenoids in the petal coloration of Liriodendron tulipifera. Hortic. Res. 7, 70 (2020).

[32] Druart, N. et al. Environmental and hormonal regulation of the activity-dormancy cycle in the cambial meristem involves stage-specific modulation of transcriptional and metabolic networks. Plant J. 50, 557-73 (2007).

[33] Ruttink, T. et al. A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell 19, 2370-90 (2007).

[34] Kirst, M. et al. Apparent homology of expressed genes from wood-forming tissues of loblolly pine (Pinus taeda L.) with Arabidopsis thaliana. Proc. Natl Acad. Sci. USA 100, 7383-8 (2003).

[35] Foucart, C. et al. Transcript profiling of a xylem vs phloem cDNA subtractive library identifies new genes expressed during xylogenesis in Eucalyptus. New Phytol. 170, 739-52 (2006).

[36] Larisch, C. et al. Poplar wood rays are involved in seasonal remodeling of tree physiology. Plant Physiol. 160, 1515-29 (2012).

[37] Reyes, J. L. & Chua, N. H. ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 49, 592-606 (2007).

[38] Li, D. et al. Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening. Hortic. Res. 6, 26 (2019).

[39] Wang, Y., Ding, Y., Yu, D., Xue, W. & Liu, J. High-throughput sequencing-based genome-wide identification of microRNAs expressed in developing cotton seeds. Sci. China Life Sci. 58, 778-786 (2015).

[40] Cai, H. et al. MiRNA-target pairs regulate adventitious rooting in Populus: a functional role for miR167a and its target Auxin response factor 8. Tree Physiol. 39, 1922-1936 (2019).

[41] Fan, D. et al. MicroRNA6443-mediated regulation of FERULATE 5-HYDROXYLASE gene alters lignin composition and enhances saccharification in Populus tomentosa. New Phytol. 226, 410-425 (2020).

[42] Song, Y., Tian, M., Ci, D. & Zhang, D. Methylation of microRNA genes regulates gene expression in bisexual flower development in andromonoecious poplar. J. Exp. Bot. 66, 1891-905 (2015).

[43] Li, B., Qin, Y., Duan, H., Yin, W. & Xia, X. Genome-wide characterization of new and drought stress responsive microRNAs in Populus euphratica. J. Exp. Bot. 62, 3765-79 (2011).

[44] Li, B. et al. Global identification of miRNAs and targets in Populus euphratica under salt stress. Plant Mol. Biol. 81, 525-39 (2013).

[45] Chen, L. et al. Genome-wide profiling of novel and conserved Populus microRNAs involved in pathogen stress response by deep sequencing. Planta 235, 873-83 (2012).

[46] Robischon, M., Du, J., Miura, E. & Groover, A. The Populus class III HD ZIP, popREVOLUTA, influences cambium initiation and patterning of woody stems. Plant Physiol. 155, 1214-25 (2011).

[47] Ko, J. H., Prassinos, C. & Han, K. H. Developmental and seasonal expression of PtaHB1, a Populus gene encoding a class III HD-Zip protein, is closely associated with secondary growth and inversely correlated with the level of microRNA (miR166). New Phytol. 169, 469-78 (2006).

[48] Wang, J. W., Czech, B. & Weigel, D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana. Cell 138, 738-49 (2009).

[49] Wang, J. W. et al. miRNA control of vegetative phase change in trees. PLoS Genet. 7, e1002012 (2011).

[50] Nilsson, J. et al. Dissecting the molecular basis of the regulation of wood formation by auxin in hybrid aspen. Plant Cell. 20, 843-55 (2008).

[51] Hayashi, K. The interaction and integration of auxin signaling components. Plant Cell Physiol. 53, 965-75 (2012).

[52] Baba, K. et al. Activity-dormancy transition in the cambial meristem involves stage-specific modulation of auxin response in hybrid aspen. Proc. Natl Acad. Sci. USA 108, 3418-23 (2011).

[53] Shuai, P., Liang, D., Zhang, Z., Yin, W. & Xia, X. Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis. BMC Genomics 14, 233 (2013).

[54] Lu, S., Sun, Y. H. & Chiang, V. L. Stress-responsive microRNAs in Populus. Plant J. 55, 131-51 (2008).

[55] Chen, M. & Cao, Z. M. Expression profiling of microRNAs in poplar upon infection with the foliar rust fungus Melampsora larici-populina. BMC Genomics 16, 696 (2015).

[56] Ramachandran, P., Carlsbecker, A. & Etchells, J. P. Class III HD-ZIPs govern vascular cell fate: an HD view on patterning and differentiation. J. Exp. Bot. 68, 55-69 (2017).

[57] Schob, H. & Grossniklaus, U. The first high-resolution DNA “methylome”. Cell 126, 1025-1028 (2006).

[58] Wang, X. F. et al. Genome-wide and organ-specific landscapes of epigenetic modifications and their relationships to mRNA and small RNA transcriptomes in maize. Plant Cell 21, 1053-1069 (2009).

[59] Lister, R. et al. Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell 133, 523-36 (2008).

[60] Cokus, S. J. et al. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature 452, 215-219 (2008).

[61] Kumar, G., Rattan, U. K. & Singh, A. K. Chilling-mediated DNA methylation changes during dormancy and its release reveal the importance of epigenetic regulation during winter dormancy in Apple (Malus x domestica Borkh.). PLoS ONE 11, e0149934 (2016).

[62] Turco, G. M. et al. DNA methylation and gene expression regulation associated with vascularization in Sorghum bicolor. New Phytol. 214, 1213-1229 (2017).

[63] Feng, S. et al. Conservation and divergence of methylation patterning in plants and animals. Proc. Natl Acad. Sci. USA 107, 8689-94 (2010).

[64] Li, J. et al. Multi-omics analyses reveal epigenomics basis for cotton somatic embryogenesis through successive regeneration acclimation process. Plant Biotechnol. J. 17, 435-450 (2019).

[65] Guo, H. et al. Somatic embryogenesis critical initiation stage-specific^mCHH hypomethylation reveals epigenetic basis underlying embryogenic redifferentiation in cotton . Plant Biotechnol. J. 18, 1648-1650 (2020).

[66] Chen, X. et al. Genome-wide investigation of DNA methylation dynamics reveals a critical role of DNA demethylation during the early somatic embryogenesis of Dimocarpus longan Lour. Tree Physiol. 40, 1807-1826 (2020).

[67] Ji, L. et al. Genome-wide reinforcement of DNA methylation occurs during somatic embryogenesis in soybean. Plant Cell 31, 2315-2331 (2019).

[68] Xu, J. et al. Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple. Plant Biotechnol. J. 16, 672-687 (2018).

[69] Chan, S. W., Henderson, I. R. & Jacobsen, S. E. Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat. Rev. Genet. 6, 351-60 (2005).

[70] Tran, R. K. et al. DNA methylation profiling identifies CG methylation clusters in Arabidopsis genes. Curr. Biol. 15, 154-9 (2005).

[71] Li, X. et al. Single-base resolution maps of cultivated and wild rice methylomes and regulatory roles of DNA methylation in plant gene expression. BMC Genomics 13, 300 (2012).

[72] Lindroth, A. M. et al. Requirement of CHROMOMETHYLASE3 for maintenance of CpXpG methylation. Science 292, 2077-80 (2001).

[73] Gent, J. I. et al. CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize. Genome Res. 23, 628-37 (2013).

[74] Schultz, M. D., Schmitz, R. J. & Ecker, J. R. ‘Leveling’ the playing field for analyses of single-base resolution DNA methylomes. Trends Genet. 28, 583-585 (2012).