Impacts of in-vitro zebularine treatment on genome-wide DNA methylation and transcriptomic profiles in Salix purpurea L.

Renewable energy resources such as biomass are crucial for a sustainable global society. Trees are a major source of lignocellulosic biomass, which can vary in response to different environmental factors owing to epigenetic regulation, such as DNA C-methylation. To investigate the effects of DNA methylation on plant development and wood formation, and its impacts on gene expression, with a focus on secondary cell wall (SCW)-associated genes, Salix purpurea plantlets were cloned from buds derived from a single hybrid tree for both treatment and control conditions. For the treatment condition, buds were exposed to 50 µM zebularine in vitro and a combined strategy of whole-genome bisulfite sequencing (WGBS) and RNA-seq was employed to examine the methylome and transcriptome profiles of different tissues collected at various time points under both conditions. Transcriptomic and methylome data revealed that most of the promoter and gene body demethylation had no marked effects on the expression profiles of genes. Nevertheless, gene expression tended to decrease with the increased methylation levels of genes with highly methylated promoters. Results indicated that demethylation is less evident in centromeric regions and sex chromosomes. Promoters of secondary cell wall-associated genes, such as 4-coumarate-CoA ligase-like and Rac-like GTP-binding protein RHO, were differentially methylated in the secondary xylem samples collected from two-month potted treated plants compared to control samples. Our results provide novel insights into DNA methylation and gene expression landscapes and a basis for investigating the epigenetic regulation of wood formation in S. purpurea as a model plant for bioenergy species.

Revealing the transitory and local effect of zebularine on development and on proteome dynamics of Salix purpurea.

Introduction: DNA methylation plays major roles in the epigenetic regulation of gene expression, transposon and transcriptional silencing, and DNA repair, with implications in developmental processes and phenotypic plasticity. Relevantly for woody species, DNA methylation constitutes a regulative layer in cell wall dynamics associated with xylogenesis. The use of methyltransferase and/or demethylase inhibitors has been proven informative to shed light on the methylome dynamics behind the regulation of these processes. Methods: The present work employs the cytidine analog zebularine to inhibit DNA methyltransferases and induce DNA hypomethylation in Salix purpurea plantlets grown in vitro and in soil. An integrative approach was adopted to highlight the effects of zebularine on proteomic dynamics, revealing age-specific (3 weeks of in vitro culture and 1 month of growth in soil) and tissue-specific (stem and root) effects. Results and discussion: After 3 weeks of recovery from zebularine treatment, a decrease of 5-mC levels was observed in different genomic contexts in the roots of explants that were exposed to zebularine, whereas a functionally heterogeneous subset of protein entries was differentially accumulated in stem samples, including entries related to cell wall biosynthesis, tissue morphogenesis, and hormonal regulation. Significant proteomic remodeling was revealed in the development from in vitro to in-soil culture, but no significant changes in 5-mC levels were observed. The identification of tissue-specific proteomic hallmarks in combination with hypomethylating agents provides new insights into the role of DNA methylation and proteome in early plant development in willow species. Proteomic data are available via ProteomeXchange with identifier PXD045653. WGBS data are available under BioProject accession PRJNA889596.

A comparative transcriptomic approach to understanding the formation of cork.

KEY MESSAGE: The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.