RNAseq based variant dataset in a black poplar association panel.

OBJECTIVE: Black poplar (Populus nigra L.) is a species native to Eurasia with a wide distribution area. It is an ecologically important species from riparian ecosystems, that is used as a parent of interspecific (P. deltoides x P. nigra) cultivated poplar hybrids. Variant detection from transcriptomics sequences of 241 P. nigra individuals, sampled in natural populations from 11 river catchments (in four European countries) is described here. These data provide new valuable resources for population structure analysis, population genomics and genome-wide association studies. DATA DESCRIPTION: We generated transcriptomics data from a mixture of young differentiating xylem and cambium tissues of 480 Populus nigra trees sampled in a common garden experiment located at Orléans (France), corresponding to 241 genotypes (2 clonal replicates per genotype, at maximum) by using RNAseq technology. We launched on the resulting sequences an in-silico pipeline that allowed us to obtain 878,957 biallelic polymorphisms without missing data. More than 99% of these positions are annotated and 98.8% are located on the 19 chromosomes of the P. trichocarpa reference genome. The raw RNAseq sequences are available at the NCBI Sequence Read Archive SPR188754 and the variant dataset at the Recherche Data Gouv repository under https://doi.org/10.15454/8DQXK5 .

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus.

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

Zebularine, a DNA Methylation Inhibitor, Activates Anthocyanin Accumulation in Grapevine Cells.

Through its role in the regulation of gene expression, DNA methylation can participate in the control of specialized metabolite production. We have investigated the link between DNA methylation and anthocyanin accumulation in grapevine using the hypomethylating drug, zebularine and Gamay Teinturier cell suspensions. In this model, zebularine increased anthocyanin accumulation in the light, and induced its production in the dark. To unravel the underlying mechanisms, cell transcriptome, metabolic content, and DNA methylation were analyzed. The up-regulation of stress-related genes, as well as a decrease in cell viability, revealed that zebularine affected cell integrity. Concomitantly, the global DNA methylation level was only slightly decreased in the light and not modified in the dark. However, locus-specific analyses demonstrated a decrease in DNA methylation at a few selected loci, including a CACTA DNA transposon and a small region upstream from the UFGT gene, coding for the UDP glucose:flavonoid-3-O-glucosyltransferase, known to be critical for anthocyanin biosynthesis. Moreover, this decrease was correlated with an increase in UFGT expression and in anthocyanin content. In conclusion, our data suggest that UFGT expression could be regulated through DNA methylation in Gamay Teinturier, although the functional link between changes in DNA methylation and UFGT transcription still needs to be demonstrated.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Developmental Methylome of the Medicinal Plant Catharanthus roseus Unravels the Tissue-Specific Control of the Monoterpene Indole Alkaloid Pathway by DNA Methylation.

Catharanthus roseus produces a wide spectrum of monoterpene indole alkaloids (MIAs). MIA biosynthesis requires a tightly coordinated pathway involving more than 30 enzymatic steps that are spatio-temporally and environmentally regulated so that some MIAs specifically accumulate in restricted plant parts. The first regulatory layer involves a complex network of transcription factors from the basic Helix Loop Helix (bHLH) or AP2 families. In the present manuscript, we investigated whether an additional epigenetic layer could control the organ-, developmental- and environmental-specificity of MIA accumulation. We used Whole-Genome Bisulfite Sequencing (WGBS) together with RNA-seq to identify differentially methylated and expressed genes among nine samples reflecting different plant organs and experimental conditions. Tissue specific gene expression was associated with specific methylation signatures depending on cytosine contexts and gene parts. Some genes encoding key enzymatic steps from the MIA pathway were found to be simultaneously differentially expressed and methylated in agreement with the corresponding MIA accumulation. In addition, we found that transcription factors were strikingly concerned by DNA methylation variations. Altogether, our integrative analysis supports an epigenetic regulation of specialized metabolisms in plants and more likely targeting transcription factors which in turn may control the expression of enzyme-encoding genes.

Hypomethylation of the aquatic invasive plant, Ludwigia grandiflora subsp. hexapetala mimics the adaptive transition into the terrestrial morphotype.

Ongoing global changes affect ecosystems and open up new opportunities for biological invasion. The ability of invasive species to rapidly adapt to new environments represents a relevant model for studying short-term adaptation mechanisms. The aquatic invasive plant, Ludwigia grandiflora subsp. hexapetala, is classified as harmful in European rivers. In French wet meadows, this species has shown a rapid transition from aquatic to terrestrial environments with emergence of two distinct morphotypes in 5 years. To understand the heritable mechanisms involved in adjustment to such a new environment, we investigate both genetic and epigenetic as possible sources of flexibility involved in this fast terrestrial transition. We found a low overall genetic differentiation between the two morphotypes arguing against the possibility that terrestrial morphotype emerged from a new adaptive genetic capacity. Artificial hypomethylation was induced on both morphotypes to assess the epigenetic hypothesis. We analyzed global DNA methylation, morphological changes, phytohormones and metabolite profiles of both morphotype responses in both aquatic and terrestrial conditions in shoot and root tissues. Hypomethylation significantly affected morphological variables, phytohormone levels and the amount of some metabolites. The effects of hypomethylation depended on morphotypes, conditions and plant tissues, which highlighted differences among the morphotypes and their plasticity. Using a correlative integrative approach, we showed that hypomethylation of the aquatic morphotype mimicked the characteristics of the terrestrial morphotype. Our data suggest that DNA methylation rather than a new adaptive genetic capacity is playing a key role in L. grandiflora subsp. hexapetala plasticity during its rapid aquatic to terrestrial transition.

Assessing Global DNA Methylation Changes Associated with Plasticity in Seven Highly Inbred Lines of Snapdragon Plants (Antirrhinum majus).

Genetic and epigenetic variations are commonly known to underlie phenotypic plastic responses to environmental cues. However, the role of epigenetic variation in plastic responses harboring ecological significance in nature remains to be assessed. The shade avoidance response (SAR) of plants is one of the most prevalent examples of phenotypic plasticity. It is a phenotypic syndrome including stem elongation and multiple other traits. Its ecological significance is widely acknowledged, and it can be adaptive in the presence of competition for light. Underlying genes and pathways were identified, but evidence for its epigenetic basis remains scarce. We used a proven and accessible approach at the population level and compared global DNA methylation between plants exposed to regular light and three different magnitudes of shade in seven highly inbred lines of snapdragon plants (Antirrhinum majus) grown in a greenhouse. Our results brought evidence of a strong SAR syndrome for which magnitude did not vary between lines. They also brought evidence that its magnitude was not associated with the global DNA methylation percentage for five of the six traits under study. The magnitude of stem elongation was significantly associated with global DNA demethylation. We discuss the limits of this approach and why caution must be taken with such results. In-depth approaches at the DNA sequence level will be necessary to better understand the molecular basis of the SAR syndrome.

Winter-dormant shoot apical meristem in poplar trees shows environmental epigenetic memory.

Trees have a long lifespan and must continually adapt to environmental pressures, notably in the context of climate change. Epigenetic mechanisms are doubtless involved in phenotypic plasticity and in stress memory; however, little evidence of the role of epigenetic processes is available for trees growing in fields. Here, we analyzed the possible involvement of epigenetic mechanisms in the winter-dormant shoot apical meristem of Populus × euramericana clones in memory of the growing conditions faced during the vegetative period. We aimed to estimate the range of genetic and environmentally induced variations in global DNA methylation and to evaluate their correlation with changes in biomass production, identify differentially methylated regions (DMRs), and characterize common DMRs between experiments. We showed that the variations in global DNA methylation between conditions were genotype dependent and correlated with biomass production capacity. Microarray chip analysis allowed detection of DMRs 6 months after the stressful summer period. The 161 DMRs identified as common to three independent experiments most notably targeted abiotic stress and developmental response genes. Results are consistent with a winter-dormant shoot apical meristem epigenetic memory of stressful environmental conditions that occurred during the preceding summer period. This memory may facilitate tree acclimation.

Changes in the epigenome and transcriptome of the poplar shoot apical meristem in response to water availability affect preferentially hormone pathways.

The adaptive capacity of long-lived organisms such as trees to the predicted climate changes, including severe and successive drought episodes, will depend on the presence of genetic diversity and phenotypic plasticity. Here, the involvement of epigenetic mechanisms in phenotypic plasticity toward soil water availability was examined in Populus×euramericana. This work aimed at characterizing (i) the transcriptome plasticity, (ii) the genome-wide plasticity of DNA methylation, and (iii) the function of genes affected by a drought-rewatering cycle in the shoot apical meristem. Using microarray chips, differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were identified for each water regime. The rewatering condition was associated with the highest variations of both gene expression and DNA methylation. Changes in methylation were observed particularly in the body of expressed genes and to a lesser extent in transposable elements. Together, DEGs and DMRs were significantly enriched in genes related to phytohormone metabolism or signaling pathways. Altogether, shoot apical meristem responses to changes in water availability involved coordinated variations in DNA methylation, as well as in gene expression, with a specific targeting of genes involved in hormone pathways, a factor that may enable phenotypic plasticity.

Epigenomics and bolting tolerance in sugar beet genotypes.

In sugar beet (Beta vulgaris altissima), bolting tolerance is an essential agronomic trait reflecting the bolting response of genotypes after vernalization. Genes involved in induction of sugar beet bolting have now been identified, and evidence suggests that epigenetic factors are involved in their control. Indeed, the time course and amplitude of DNA methylation variations in the shoot apical meristem have been shown to be critical in inducing sugar beet bolting, and a few functional targets of DNA methylation during vernalization have been identified. However, molecular mechanisms controlling bolting tolerance levels among genotypes are still poorly understood. Here, gene expression and DNA methylation profiles were compared in shoot apical meristems of three bolting-resistant and three bolting-sensitive genotypes after vernalization. Using Cot fractionation followed by 454 sequencing of the isolated low-copy DNA, 6231 contigs were obtained that were used along with public sugar beet DNA sequences to design custom Agilent microarrays for expression (56k) and methylation (244k) analyses. A total of 169 differentially expressed genes and 111 differentially methylated regions were identified between resistant and sensitive vernalized genotypes. Fourteen sequences were both differentially expressed and differentially methylated, with a negative correlation between their methylation and expression levels. Genes involved in cold perception, phytohormone signalling, and flowering induction were over-represented and collectively represent an integrative gene network from environmental perception to bolting induction. Altogether, the data suggest that the genotype-dependent control of DNA methylation and expression of an integrative gene network participate in bolting tolerance in sugar beet, opening up perspectives for crop improvement.

In search of markers for somatic embryo maturation in hybrid larch (Larix × eurolepis): global DNA methylation and proteomic analyses.

A global DNA methylation and proteomics approach was used to investigate somatic embryo maturation in hybrid larch. Each developmental step during somatic embryogenesis was associated with a distinct and significantly different global DNA methylation level: from 45.8% mC for undifferentiated somatic embryos (1-week proliferation) to 61.5% mC for immature somatic embryos (1-week maturation), while maturation was associated with a decrease in DNA methylation to 53.4% for mature cotyledonary somatic embryos (8-weeks maturation). The presence of 5-azacytidine (hypo-methylating agent) or hydroxyurea (hyper-methylating agent) in the maturation medium altered the global DNA methylation status of the embryogenic cultures, and significantly reduced both their relative growth rate and embryogenic potential, suggesting an important role for DNA methylation in embryogenesis. Maturation was also assessed by examining changes in the total protein profile. Storage proteins, identified as legumin- and vicilin-like, appeared at the precotyledonary stage. In the proteomic study, total soluble proteins were extracted from embryos after 1 and 8 weeks of maturation, and separated by two-dimensional gel electrophoresis. There were 147 spots which showed significant differences between the stages of maturation; they were found to be involved mainly in primary metabolism and the stabilization of the resulting metabolites. This indicated that the somatic embryo was still metabolically active at 8 weeks of maturation. This is the first report of analyses of global DNA methylation (including the effects of hyper- and hypo-treatments) and proteome during somatic embryogenesis in hybrid larch, and thus provides novel insights into maturation of conifer somatic embryos.

Effects of polysaccharide elicitors on secondary metabolite production and antioxidant response in Hypericum perforatum L. shoot cultures.

The effects of polysaccharide elicitors such as chitin, pectin, and dextran on the production of phenylpropanoids (phenolics and flavonoids) and naphtodianthrones (hypericin and pseudohypericin) in Hypericum perforatum shoot cultures were studied. Nonenzymatic antioxidant properties (NEAOP) and peroxidase (POD) activity were also observed in shoot extracts. The activities of phenylalanine ammonia lyase (PAL) and chalcone-flavanone isomerase (CHFI) were monitored to estimate channeling in phenylpropanoid/flavonoid pathways of elicited shoot cultures. A significant suppression of the production of total phenolics and flavonoids was observed in elicited shoots from day 14 to day 21 of postelicitation. This inhibition of phenylpropanoid production was probably due to the decrease in CHFI activity in elicited shoots. Pectin and dextran promoted accumulation of naphtodianthrones, particularly pseudohypericin, within 21 days of postelicitation. The enhanced accumulation of naphtodianthrones was positively correlated with an increase of PAL activity in elicited shoots. All tested elicitors induced NEAOP at day 7, while chitin and pectin showed increase in POD activity within the entire period of postelicitation. The POD activity was in significantly positive correlation with flavonoid and hypericin contents, suggesting a strong perturbation of the cell redox system and activation of defense responses in polysaccharide-elicited H. perforatum shoot cultures.

Identification of differentially methylated regions during vernalization revealed a role for RNA methyltransferases in bolting.

Sugar beet (Beta vulgaris altissima) is a biennial root crop with an absolute requirement for cold exposure to bolt and flower, a process called vernalization. Global DNA methylation variations have been reported during vernalization in several plants. However, few genes targeted by DNA methylation during vernalization have been described. The objectives of this study were to identify differentially methylated regions and to study their involvement in bolting induction and tolerance. Restriction landmark genome scanning was applied to DNA from shoot apical meristems of sugar beet genotypes, providing a direct quantitative epigenetic assessment of several CG methylated genes without prior knowledge of gene sequence. Several differentially methylated regions exhibiting variations of gene-body DNA methylation and expression during cold exposure and/or between genotypes were identified, including an AROGENATE DEHYDRATASE and two RNA METHYLCYTOSINE TRANSFERASE sequences. One RNA METHYLCYTOSINE TRANSFERASE sequence displayed gene-body hypermethylation and activation of expression, while the other was hypomethylated and inhibited by cold exposure. Global RNA methylation and phenolic compound levels changed during cold exposure in a genotype-dependent way. The use of methyl RNA immunoprecipitation of total RNA and reverse transcription-PCR analysis revealed mRNA methylation in a vernalized bolting-resistant genotype for the FLOWERING LOCUS 1 gene, a repressor of flowering. Finally, Arabidopsis mutants for RNA METHYLCYTOSINE TRANSFERASE and AROGENATE DEHYDRATASE were shown to exhibit, under different environmental conditions, early or late bolting phenotypes, respectively. Overall, the data identified functional targets of DNA methylation during vernalization in sugar beet, and it is proposed that RNA methylation and phenolic compounds play a role in the floral transition.

Methylome of DNase I sensitive chromatin in Populus trichocarpa shoot apical meristematic cells: a simplified approach revealing characteristics of gene-body DNA methylation in open chromatin state.

DNA methylation is involved in the control of plant development and adaptation to the environment through modifications of chromatin compaction and gene expression. In poplar (Populus trichocarpa), a perennial plant, variations in DNA methylation have been reported between genotypes and tissues or in response to drought. Nevertheless, the relationships between gene-body DNA methylation, gene expression and chromatin compaction still need clarification. Here, DNA methylation was mapped in the noncondensed chromatin fraction from P. trichocarpa shoot apical meristematic cells, the center of plant morphogenesis, where DNA methylation variations could influence the developmental trajectory. DNase I was used to isolate the noncondensed chromatin fraction. Methylated sequences were immunoprecipitated, sequenced using Illumina/Solexa technology and mapped on the v2.0 poplar genome. Bisulfite sequencing of candidate sequences was used to confirm mapping data and to assess cytosine contexts and methylation levels. While the methylated DNase I hypersensitive site fraction covered 1.9% of the poplar genome, it contained sequences corresponding to 74% of poplar gene models, mostly exons. The level and cytosine context of gene-body DNA methylation varied with the structural characteristics of the genes. Taken together, our data show that DNA methylation is widespread and variable among genes in open chromatin of meristematic cells, in agreement with a role in their developmental trajectory.

Genic DNA methylation changes during in vitro organogenesis: organ specificity and conservation between parental lines of epialleles.

During differentiation, in vitro organogenesis calls for the adjustment of the gene expression program toward a new fate. The role of epigenetic mechanisms including DNA methylation is suggested but little is known about the loci affected by DNA methylation changes, particularly in agronomic plants for witch in vitro technologies are useful such as sugar beet. Here, three pairs of organogenic and non-organogenic in vitro cell lines originating from different sugar beet (Beta vulgaris altissima) cultivars were used to assess the dynamics of DNA methylation at the global or genic levels during shoot or root regeneration. The restriction landmark genome scanning for methylation approach was applied to provide a direct quantitative epigenetic assessment of several CG methylated genes without prior knowledge of gene sequence that is particularly adapted for studies on crop plants without a fully sequenced genome. The cloned sequences had putative roles in cell proliferation, differentiation or unknown functions and displayed organ-specific DNA polymorphism for methylation and changes in expression during in vitro organogenesis. Among them, a potential ubiquitin extension protein 6 (UBI6) was shown, in different cultivars, to exhibit repeatable variations of DNA methylation and gene expression during shoot regeneration. In addition, abnormal development and callogenesis were observed in a T-DNA insertion mutant (ubi6) for a homologous sequence in Arabidopsis. Our data showed that DNA methylation is changed in an organ-specific way for genes exhibiting variations of expression and playing potential role during organogenesis. These epialleles could be conserved between parental lines opening perspectives for molecular markers.

Time course and amplitude of DNA methylation in the shoot apical meristem are critical points for bolting induction in sugar beet and bolting tolerance between genotypes.

An epigenetic control of vernalization has been demonstrated in annual plants such as Arabidopsis and cereals, but the situation remains unclear in biennial plants such as sugar beet that has an absolute requirement for vernalization. The role of DNA methylation in flowering induction and the identification of corresponding target loci also need to be clarified. In this context, sugar beet (Beta vulgaris altissima) genotypes differing in bolting tolerance were submitted to various bolting conditions such as different temperatures and/or methylating drugs. DNA hypomethylating treatment was not sufficient to induce bolting while DNA hypermethylation treatment inhibits and delays bolting. Vernalizing and devernalizing temperatures were shown to affect bolting as well as DNA methylation levels in the shoot apical meristem. In addition, a negative correlation was established between bolting and DNA methylation. Genotypes considered as resistant or sensitive to bolting could also be distinguished by their DNA methylation levels. Finally, sugar beet homologues of the Arabidopsis vernalization genes FLC and VIN3 exhibited distinct DNA methylation marks during vernalization independently to the variations of global DNA methylation. These vernalization genes also displayed differences in mRNA accumulation and methylation profiles between genotypes resistant or sensitive to bolting. Taken together, the data suggest that the time course and amplitude of DNA methylation variations are critical points for the induction of sugar beet bolting and represent an epigenetic component of the genotypic bolting tolerance, opening up new perspectives for sugar beet breeding.

Relationship between DNA methylation and histone acetylation levels, cell redox and cell differentiation states in sugarbeet lines.

In order to evaluate the permanent chromatin remodeling in plant allowing their high developmental plasticity, three sugarbeet cell lines (Beta vulgaris L. altissima) originating from the same mother plant and exhibiting graduate states of differentiation were analyzed. Cell differentiation has been estimated by the cell redox state characterized by 36 biochemical parameters as reactive oxygen species steady-state levels, peroxidation product contents and enzymatic or non-enzymatic protective systems. Chromatin remodeling has been estimated by the measurement of levels of DNA methylation, histone acetylation and corresponding enzyme activities that were shown to differ between cell lines. Furthermore, distinct loci related to proteins involved in cell cycle, gene expression regulation and cell redox state were shown by restriction landmark genome scanning or bisulfite sequencing to display differential methylation states in relation to the morphogenic capacity of the lines. DNA methylating, demethylating and/or histone acetylating treatments allowed to generate a collection of sugarbeet cell lines differing by their phenotypes (from organogenic to dedifferentiated), methylcytosine percentages (from 15.0 to 43.5%) and acetylated histone ratios (from 0.37 to 0.52). Correlations between methylcytosine or acetylated histone contents and levels of various parameters (23 or 7, respectively, out of 36) of the cell redox state could be established. These data lead to the identification of biomarkers of sugarbeet morphogenesis in vitro under epigenetic regulation and provide evidence for a connection between plant morphogenesis in vitro, cell redox state and epigenetic mechanisms.

DNA methylating and demethylating treatments modify phenotype and cell wall differentiation state in sugarbeet cell lines.

In plants organogenesis, cell differentiation and dedifferentiation are fundamental processes allowing high developmental plasticity. Such plasticity involved epigenetic mechanisms but limited knowledge is available concerning quantitative aspects. Three sugarbeet (Beta vulgaris L. altissima) cell lines originating from the same mother plant and exhibiting graduate states of morphogenesis were used to assess whether these differences could be related or not to changes in DNA methylation levels. Methylcytosine percentages from 18.3 to 28.8% and distinct levels of DNA methyltransferase (EC 2.1.1.37) activities were shown in the three cell lines. The lowest methylcytosine percentage was associated to organogenesis. In order to test the plasticity of these cell lines, various treatments causing DNA hypo or hypermethylation were performed at different times and concentrations. In this collection of treated lines with+/-10% of methylcytosine percentages, loss of organogenic properties and cell dedifferentiation were observed. As cell wall formation fits well with cell differentiation state, the lignification process was further investigated in treated and untreated lines as a biochemical marker of the phenotypic changes. For example, peroxidase specific activities (EC 1.11.1.7) varied from 0.7 to 0.02 pkat mg(-1) of protein in organogenic and dedifferentiated lines, respectively. A negative relationship between peroxidase activities, incorporation of cell wall-bound phenolic compounds as ferulate and sinapate derivatives and methylcytosine percentages was obtained. This is the first biochemical evidence that phenotypic changes in plant cells induced by DNA hypo- or hypermethylating treatments are correlated in a linear relationship to modifications of the cell wall differentiation state.