Plant stress biology in epigenomic era

Plant stress biology in epigenomic era

Journal Pre-proof Plant Stress Biology in Epigenomic Era Anna Perrone PII: S0168-9452(19)31549-3 DOI: https://doi.org/10.1016/j.plantsci.2019.1103...

5MB Sizes 0 Downloads 64 Views

Journal Pre-proof Plant Stress Biology in Epigenomic Era Anna Perrone

PII:

S0168-9452(19)31549-3

DOI:

https://doi.org/10.1016/j.plantsci.2019.110376

Reference:

PSL 110376

To appear in:

Plant Science

Federico Martinelli

PII:

S0168-9452(19)31549-3

DOI:

https://doi.org/10.1016/j.plantsci.2019.110376

Reference:

PSL 110376

To appear in:

Plant Science

Received Date:

12 April 2019

Please cite this article as: Perrone A, Martinelli F, Plant Stress Biology in Epigenomic Era, Plant Science (2019), doi: https://doi.org/10.1016/j.plantsci.2019.110376

This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2019 Published by Elsevier.

Plant Stress Biology in Epigenomic Era Anna Perronea, Federico Martinellib* a

Department of Biological, Chemical and Pharmaceutical Sciences and Technologies (STEBICEF), University of Palermo, Viale delle Scienze, Palermo, 90128, Italy b Department of Biology, University of Firenze, Sesto Fiorentino, Florence, 50019, Italy

ro of

*Corresponding authors: Federico Martinelli Email: [email protected] Anna Perrone Email: [email protected]

Highlights

na

ABSTRACT

lP

re

-p

 Different DNA methylation enzymes are involved in different plant processes  An integrated systems biology approach is needed in crop stress biology  Priming drives a wide range of epigenetic mechanisms  An epigenetic code is present in response to environmental stresses in plants

Recent progress in "omics" methodologies allow us to gain insight into the complex molecular

ur

regulatory networks underlying plant responses to environmental stresses. Among the different

Jo

genome-wide analysis, epigenomics is the most under-investigated “omic” approach requiring more critical and speculative discussion about approaches, methods and experimental designs. Epigenomics allows us to gain insight into the molecular adaptation of plants in response to environmental stresses. The identification of epigenetic marks transmitted during filial generations enables new theories to be developed on the evolution of living organisms in relation to environmental changes. The molecular mechanisms driving the capacity of plants to memorize a stress and to generate stress-resistant progenies are still unclear and scarcely investigated. The 1

elucidation of these cryptic molecular switches will assist breeders in designing crops characterized by minimally compromised productivity in relation to stresses caused by climate change. The aim of this review is to briefly describe the most uptodate epigenomic approaches, update recent progresses in crop epigenomics in plant stress biology, and to stimulate the discussion of new epigenomic methods and approaches in the new era of “omic” sciences.

Keywords: epigenomics, omics, plant stress biology, priming, stress memory, transgenerational

1. Epigenetic mechanisms and techniques in plants

ro of

modifications.

-p

Epigenetic mechanisms refer to all the modifications of gene expression that occur without altering the DNA sequence. The main epigenetic-driven modifications belong to three types of mechanisms:

re

1) chromatin restructuring and histone post-translational modifications (PTMs), 2) DNA

lP

methylation changes, and 3) post-transcriptional modifications due to non-coding RNAs (microRNA or long non–coding RNAs).

na

1.1 Chromatin restructuring and histone post-translational modifications Histone post-transcriptional modifications are key modifications significantly altering gene

ur

activities due to modifications at the N-terminal tail of core histones of the nucleosome. These

Jo

chemical modifications that provoke chromatin restructuring consists of acetylation, methylation, phosphorylation and ubiquitination [1]. The first two modifications, histone acetylation and methylation, are generally due to the activity of large (HATs) and histone lysine methyltransferases (HKMTs) enzyme families. These epigenetic marks can be reversed by other enzymes such as histone deacetylases (HDACs) and histone demethylases (HDMs), respectively [2]. While phosphorylation of histone H3 Ser-10 has been linked with responses to salinity, cold stress and abscissic acid treatment [3], ubiquitination has been mainly associated with increases in 2

transcription [4]. Ubiquitination conjugation is a post-translational modification that is a key player in the signal transduction of oxidative stress responses inducing DNA damage repair, protein degradation, and trafficking [5]. Particularly changes in K63 ubiquitination in response to oxidative stresses have been shown to affect polysome stability and protein expression. These results implied that ubiquitination has a regulatory role in redox processes [5]. Moreover, the tails of the N-terminal histones modified through, for example, acetylation processes, have been shown to be crucial events both for defense mechanisms and for plant development [6].

ro of

Studies on the pluripotency of plant cells demonstrate the existence of regulatory mechanisms that depend on important epigenetic factors in Arabidopsis. In this context, epigenetic modifications could be key players in reprogramming lineage-specified plant cells to pluripotency. Specifically,

-p

regeneration mechanisms [7] involve the histone acetyltransferase HAG1 / AtGCN5. Callus development causes the catalyzing of histone acetylation processes in different root-meristem gene

re

loci, including WOX5, WOX14, SCR, PLT1 and PLT2, providing an epigenetic platform for their

1.1.1 Main techniques

lP

transcriptional activation.

na

Chromatin Immuno-Precipitation combined with chip (DNA microarray) analysis (ChIP-CHIP) has been the most common technique used for the genome-wide analysis of chromatin epigenetic

ur

modifications before the arrival of the next-generation of sequencing technologies. While the ChIP-

Jo

CHIP technique is effective in easily detecting histone positions with specific modifications of their N-terminal tails [8], ChIP-sequencing combines chromatin immunoprecipitation (ChIP) with genome-wide DNA sequencing providing a deeper level of identification of the DNA sequences bound with transcription factors and other regulating proteins. ChIP-seq consists of the selection of DNA fragments linked with the studied proteins and the analysis of the DNA fragments with deep sequencing platforms (Illumina, PacBio, Roche etc..). The combination of chromatin precipitation and DNA next-generation sequencing has highly advanced our understanding of the mechanism of 3

histone modifications in regulating gene expression [9]. ChIP-sequencing requires a known genomic sequence reference that might be represented by the studied species or a phylogenetically closely related one. Chip-seq is currently the most-desirable approach for studying histone modifications because it allows a better depth/costs ratio and the identification of new unknown SNPs responsible for chromatin restructuring, absent in arrays used in the ChIP-CHIP approach.

1.2 DNA methylation analysis

ro of

One of the most commonly analyzed epigenetic modification is the methylation of cytosine occurring at the fifth carbon position. This epigenetic modification is frequently transient in response to the environment although DNA methylation can be inherited in the next generation

-p

[10]. DNA methylation plays a key role in transcriptional regulation activity during embryonic development, in key biological phenomenon such as genome stability (transposon silencing) and

re

genomic imprinting during cell differentiation. Although DNA methylation has been frequently

lP

linked with the repression of gene expression, in some instances it has been responsible for transcription activation [11]. For example, in Arabidopsis thaliana Harris et al. found that DNA methylation promoted the recruitment of a transcriptional activation protein complex to chromatin

na

sites inducing transposon silencing in neighbouring genes. It has been debated that DNA methylation is a mechanism that has occurred during evolution in

ur

order to silence transposons considered to be genomic parasites in the host genomes. This is

Jo

corroborated by the fact that transposons are generally highly methylated and represent the main targets for epigenetic silencing. Unfortunately, the study of epigenetic modifications in genomic regions rich in transposon elements is not easy due to the difficulties in determining the link between DNA methylation and transposon polymorphism for their repetitive nature and their capacity to create large insertion/deletion polymorphisms. Methylation is site specific [12] and frequently occurs at pericentric heterochromatin regions driven by histone H3 lysine 9 (H3-K9) methylation [13]. Several findings support the idea that DNA methylation also plays a key role in 4

modifying genome size, structure and genetic variations. Frequently, cytosine methylation occurs when it is present in three sequences: CG, CHG and CHH (where H = A, T or C). DNA methylation has been clearly linked to abiotic stress responses in plants (Figure 1). Each of these methylation marks are catalyzed by different enzymes. Methyltransferase1 (MET1) is responsible for the CG methylation while Chromomethylase3 (CMT3) provokes CHG methylation [14]. Domains Rearranged Methyltransferases (DRM1 and DRM2) that are guided by small interfering RNAs (siRNAs) drive the methylation of CHH [15]. CHH methylation is due to RNA-dependent DNA

ro of

methylation (RdDM) and requires plant-specific RNA polymerases IV and V. Pol IV produces precursor transcripts of 24-nt small RNAs (sRNAs) that are able to target scaffold transcripts using the complementary sequence. The maintaining of the cytosine methylation during replication

-p

depends on proteins such as Decrease in DNA Methylation 1 [DDM1; 14] while DNA methylations are generally deleted by specific classes of DNA glycosylases such as DEMETER (DME),

re

REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 and DML3 [15]. In plant stress

lP

responses, DNA methylation is a key modulator of both pathogen and symbiont interactions with plants [16] and previous works have linked DNA methylation with host susceptibility [17]. Methylome analyses permitted a better clarification of epigenomic plasticity mechanisms. For

na

example, studies of epigenetic analysis on the bean [18]

have facilitated the drawing up of

preliminary epigenomic maps based on chromatin accessibility, histone modifications and DNA

Jo

ur

methylation linked with important agronomic traits.

1.2.1 Main techniques

Next-generation sequencing technologies represent the most powerful method of identifying epigenetic modifications occurring at the DNA level. It has permitted new discoveries in the mechanisms of epigenetic regulation, through the understanding of the biological meaning of DNA methylation, chromatin modeling and the role of small RNA. NGS technologies are represented by various popular platforms, including the 454 FLX (Roche) [19], the Genome Analyzer / Hiseq 5

(Illumina Solexa) [20,21] and the SOLiD (Life Technologies) and new platforms such as Heliscope (Helicos) [22] and PacBio RS (Pacific Biosciences) [23] for single molecular sequencing, and Ion Torrent (Life Technologies), based on semiconductor chips [24]. The most common technique for genome-wide methylation analysis is the Whole Genome Bisulphite sequencing (WGBS) that consists in the treatment of DNA with sodium bisulphite followed by next-generation sequencing. Bisulphite treatment of DNA provokes the conversion of unmodified cytosines to uracil maintaining 5mC unchanged. The use of restriction enzymes sensitive to DNA methylation permits the

ro of

enriching of DNA in methylated or un-methylated sequences. WGBS allows DNA methylation sites to be identified in a genome-wide manner. This approach has been used to generate different plant methylomes [25]. Bisulphite treatment of DNA sequences precedes the hybridization of DNA

-p

fragments to tilling arrays or their direct sequencing. In addition, antibodies specific for m5C can be used to analyze them before the analysis using arrays or sequencing platforms (MeDip-ChIP or

re

MeDip-Seq). The bisulphite-treated DNA can be analyzed with PCR for target epigenetic analysis

lP

followed by sequencing [26]. However, the PCR amplification that is also included in the genomewide analysis technique is not exempt from bias. False positive and false negative results can also be produced due the incomplete cytosine conversion by sodium bisulphite or the over-conversion of

na

5mC. Indeed, PacBio sequencing that does not require bisulphite treatment is able to eliminate this bias. Although methylome analysis through sequencing might be relatively expensive for large

ur

plant genomes [27], the reduction of the cost of sequencing is doubling every 18 months, rendering

Jo

these analyses affordable for small projects (< 10000-20000 $). A method for DNA methylation is the identification of target regions with padlock-probe based targeting coupled with promoter tiling arrays [28]. An alternative is a reduced representation bisulphite sequencing (RRBS), which consists of the fractionation of the genome and the enrichment of fragments with high CpG content regions before sequencing [29]. In addition, pyrosequencing has been used without the use of methylation specific PCR [30]. This is employed to determine the specific CpG sites in the amplified product. Based on the amount of incorporated C and T, it is possible to determine the 6

ratio of C and T and individual sites of methylation. However, the costs for this technique are high. An application of this technique is represented by the analysis of genomic imprinting through the use of allele specific primers with single nucleotide polymorphisms. Another technique for the analysis of DNA methylation is the methylation-sensitive single-strand conformation analysis (MSSSCA). This method employs the single-strand conformation analysis to resolve sequence differences between methylated and unmethylated samples due to the conversion of bisulphite of cytosine to thiamine while the methylated cytosine remains unchanged. The efficiency of SSCA is

ro of

70-95% [31]. High-resolution melting analysis (HRM) also permits the analysis of DNA methylation when quantitative real-time PCR analysis is completed. It consists of applying melting cycles that allow converted to be distinguished from unconverted bisulphite-treated DNA. The PCR

-p

products are analyzed by temperature ramping that will result in the liberation of intercalating fluorescent dye during the melting process. The rapidity of the melting and consequent release of

re

the dye is due to the C to T content that represents the DNA methylation amount [32]. The

lP

weakness of this method is the ability to analyze only one or a few amplicons each time [34-35]. Methylation-sensitive single-nucleotide primer extension (Ms-SNuPE) is a technique suitable for a quick quantitation of methylation occurrence at individual CpG sites [35]. It uses the primer

na

extension method that has been optimized for detecting single-nucleotide polymorphisms. In this technique, bisulphite treatment is used together with PCR amplifications with bisulphite-specific

ur

primers. They are designed to anneal to the sequence up to the base pair immediately before the site

Jo

of DNA methylation. This is a simple method for screening DNA methylation in a large number of samples. This method also promotes genomic sequencing since the electrophoretic bands may be isolated and directly sequenced. Other methodologies employ matrix-assisted laser desorption ionization/time-of-flight (MALDI-TOF) mass spectrometry analysis to differentiate between the two polymorphic primer extension products based on the GOOD assay designed for SNP genotyping. Base-specific cleavage/MALDI-TOF consists in the in vitro transcription of the region under study using an RNA polymerase promoter site. RNAse A treatment follows retro7

transcription and the base specificity is obtained by adding incorporating cleavage resistant dTTP or dCTP if cytosine or uracil specific cleaves are desired, respectively. The fragments obtained by the cleavage

are

analyzed

with

MALDI-TOF.

The

bisulphite

treatment

provokes

the

introduction/removal of cleavage sites by C to U conversions or a shift in fragment mass by G to A conversions in the amplified reverse strand. C-specific cleavage will only work at the methylated CpG sites. This technique has a great efficiency, allowing identification of CpG sites in multiple tissues with low costs. The EpiTYPER™ assay allows the detection and quantification of DNA

ro of

methylation in a specific region through a high-resolution scan of selected regions. This method enables us to perform an epigenomic analysis of a region surrounding a single CpGs and can be employed to validate the CpGs detected by genome-wide approaches. This technique is considered

-p

a rapid method for the targeted quantification of individual CpGs in a high throughput manner through the base-specific cleavage after bisulphite-conversion of DNA coupled with MALDI-TOF

lP

re

MS [36].

1.3 Non-coding RNA-mediated epigenetic regulation It is known that 90% of the eukariotic genome is transcribed but only 1-2% of the genome is

na

translated to proteins. This means that there is a large number of RNA molecules that are circulating in the cells whose functions are unknown but most probably having a key role in the post-

ur

transcriptional regulation [37]. These RNA molecules are noncoding RNAs (ncRNAs) that are the

Jo

third main actors in epigenetic modifications. Noncoding RNA molecules can be divided in different groups based on their size, function and regulation mechanisms: microRNAs (miRNAs), short interfering RNAs (siRNAs), piwi-interacting RNAs (piRNAs), and long non-coding (lncRNAs). Another classification has divided them into two categories: short ncRNAs (<30 nts; including miRNAs, siRNAs, and piRNAs) and long ncRNAs (>200 nts). Small RNAs are important activators/inhibitors of plant environmental stresses. For example, it has been shown that the plant pathogen Botrytis cinerea delivers small RNAs into the host cells to repress plant defense genes and 8

thus deactivate host immunity [38]. Besides lncRNAs that clearly function in chromatin modifications, a few other ncRNAs also play roles in other biological processes in plants [38]. For example, the soybean early nodulin gene, Enod40, is involved in the regulation of symbiotic interactions between leguminous plants and soil bacteria [38]. Several stress-responsive miRNAs have been found in different crops leading to the possibility of using miRNA-driven RNAinterference (RNAi) to improve key traits and provide phenotypic plasticity in a climate change scenario [39]. RNA-mediated chromatin-level silencing is strongly involved in stress responses and

ro of

natural epigenetic variation and consequently in phenotypic diversity, physiological plasticity and evolutionary change [40]. Epigenomic regulation of chromatin structure and genome stability is essential for genetic interpretation information and key roles have emerged for small RNAs,

-p

proteins with domains that bind methylated DNA and DNA glycosylases in these processes [41,42]. The analysis of cytosine methylome (methylC-seq), transcriptome (mRNA-seq) and small RNA

re

transcriptome (small RNA-seq) in Arabidopsis inflorescences showed methylation patterns on a

lP

genomic scale and a direct relationship between the sRNA position and DNA methylation [43]. In addition, genomic analysis in Arabidopsisis showed that various methylated cytosines have been detected by bisulphite sequencing using NGS technologies, precisely through the so-called BS-seq

na

[44-46].

ur

1.3.1 Main techniques

Jo

The most potent method for analyzing ncRNA (including miRNA) is the use of next-generation sequencing because it permits us not only to conduct a quantitative analysis of their expression but also the identification of single nuclear polymorphisms in RNA sequences [47]. Novaseq 6000 (Illumina) is one of the most-desired platforms for studying rare transcripts and splice variants due to the higher number of reads with longer reads than previous platforms (i.e. Hiseq 2500) (250 bp instead of 75 bp). In silico methods are used since they are able to distinguish lncRNAs from protein-coding transcripts based on the absence of discernible open-reading frames. Full-length 9

targeted cDNA sequencing may be the best method to determine the coding potential of transcripts, although this approach is time-consuming and expensive. Alternatively, tiling DNA microarray for genome-wide transcriptome analysis can be employed to detect and determine the expression of lncRNAs in plants. Tiling DNA microarray analyses have been used to identify new stress-induced lncRNAs [48]. High-throughput deep RNA sequencing is able to identify missing or incomplete transcripts of lncRNA and to determine the levels of lncRNA expression. A similar approach was used to identify rare alternative splicing variants of a lncRNA, HOTAIR [49]. An approach with a

ro of

combination of different kinds of techniques (genome-wide histone modification analysis and deep sequencing) may permit the identification of the number of long intergenic ncRNAs (lincRNAs) transcribed from a K4-K36 domain, which marks active promoters with trimethylation [49].

-p

A summary of all the mentioned techniques divided by the three main epigenetic mechanisms is

re

provided in Table 1.

lP

1. What is the role of epigenomics in system biology?

An integrated systems biology approach is required to gain insight into the complex molecular regulatory networks that modulate the expression of plant traits. Epigenomics is one of the most

na

unexplored “omic” approaches to understand the adaptation of living organisms to rapid local environment change, weather or other stresses. However, the analysis of epigenomic changes in

ur

relation to plant environmental stresses requires a strong integration of the different “omic”

Jo

analyses (Figure 2). Transcriptome analysis is essential to validate the gene silencing that is predicted by DNA methylation data. The use of the other analyses (proteomics and metabolomics) are commonly used to gain insight into the post-transcriptional and post-translational regulation mechanisms linked to specific epigenetic changes. The epigenomic modification has a real effect at the phenotypic level if the consequent change in transcription is not reversed at the posttranscriptional level. Indeed, proteomic and metabolomic platforms are essential to determine if the epigenetic modification has effects on composition at protein and metabolite levels. The power of 10

epigenomic analysis can be better exploited for plant species whose genomes have been sequenced. Transcriptomic approaches permit analysis of the expression of many thousands of genes, identifying splice variants, rare genes, single nucleotide polymorphisms. However, transcriptomic results are usually inconclusive and not strictly correlated to phenotype changes since cells frequently activate multiple regulatory mechanisms for their adaptation to stress responses. Proteomic analysis permits us to understand how protein changes are modulated by rapid or chronic environmental stresses in crops [50]. A recent integrated approach has been proposed correlating

ro of

genes and protein expression, named proteogenomics [51]. Proteomics can be used to validate the expression of the genes while transcriptomic data permits the construction of customized protein databases, avoiding the need for homology-based approaches [52-54]. Metabolomics is the

-p

untargeted high-throughput technique that can be linked to epigenetic modifications. Metabolomic approaches allow a comprehensive analysis of the cell metabolism and could validate epigenomics.

re

This “omic” approach has been less employed compared with the other “omic” analyses especially

lP

for trait mapping and plant selections [55]. It allows for deep characterization of how plant metabolism is affected in response to environmental stimuli and changing climatic conditions. The recent progresses in metabolomics allow quantification of a wide range of metabolites from a small

na

amount of analyzed samples (> 10 mg). Metabolomics has been used to understand cell dehydration, water stress responses and plant-microbe interactions [55], providing a great

ur

contribution to shedding light on the link between genotype and metabolic pathways leading to the

Jo

phenotype. They help identify key points of molecular network regulation. Plants biosynthesize more than a hundred thousand metabolites with different chemical structures and quantity, each one with key roles in physiological, development, and stress response processes [56-57]. Integrated multi-omic approaches are highly effective in discovering epigenomic mechanisms of regulation of multi-factorial quantitative traits. The real challenge in analyzing both “omic” data is understanding the complex gene/protein regulatory networks identification due to gene-gene and protein-protein interactions. Nodes, molecular players (genes or proteins), are generally indicated in different colors 11

depending on their trend of expression (up, down or un-regulated between control and stress condition), different forms based on their functions, and different sizes proportional to the number of connections (Figure S1). Through the construction of these networks it is possible to identify not only the highly interactive nodes (hubs), but also low interactive proteins that play a key role in the networks because they connect hubs or cores of highly interconnected proteins/genes. In this regard, the interaction between key highly interactive proteins and non-coding RNA often plays an important role in plant stress responses [56]. For example, WRKY transcription factors are typically

ro of

key players in modulating environmental stress responses in plants at the interactome level. Several targets of miRNAs encode WRKY transcription factors [57] implying that miRNAs are involved

in the regulation of the expression of WRKYs. On the other hand, some miRNA promoters are

-p

abundant in W box sequences, the cis-elements recognized by WRKYs [58]. For example, the response to powdery mildew has been shown to be modulated by WRKY-smRNA interactome

re

[59]

lP

2. Approaches for investigating epigenomic transient and transgenerational mechanisms We propose a general workflow to investigate plant transgenerational epigenomic modifications in

na

response to environmental stresses. First, the stress-epigenetic response should be analyzed in controlled conditions to reduce undesired environmental variability. Basing ourselves on the

ur

evidence that DNA methylation generally causes the silencing of the methylated gene, the decrease in transcript abundance of the mRNA should be observed. In fact, the most desired approach is to

Jo

integrate epigenomics with transcriptomic and metabolomic approaches to observe the inhibition of expression of the methylated genes [60-63]. The use of proteomic approaches would allow the confirmation of the epigenetic-driven changes in transcript abundance although the magnitude of analysis is generally 10-fold less than RNA-seq. ITRAQ proteomic approaches could partially contribute to fill this gap through the analysis of more than 6000 proteins using small biological samples (> 50-100 mg). The gap that still exists between transcriptomics and proteomics (even with iTRAQ) precludes the conclusion that the absence of the protein targeted by the DNA methylation 12

is a consequence of a variation in its expression. It might be due that the protein was not present in the protein extract. The heterogenous chemical nature of the proteins renders the extraction complex and difficult, and does not allow the recovery of all types of proteins during the extraction processes. In this regard, Phenol-based protocols are commonly used to extract proteins from recalcitrant plant tissues [64]. Once DNA methylation analysis, transcriptomics and proteomics are conducted, a list of potential epigenetic marks can be delivered and associated with the considered stress response. Quantitative RT-PCRs may be conducted for the targeted list of methylated genes

ro of

to confirm the trend of altered expression in relation to the detected DNA methylation. Seeds can be collected and DNA may be extracted from zygotic tissues to perform a genome-wide DNA methylation analysis in order to see if the methylation signal is transmitted in the filial generation.

-p

The epigenetic patterns of the seeds can be compared with the one related to the same tissue (leaf, fruit, shoots etc.) coming from the first generation to see if the DNA methylation is a transient

re

stimulus that occurrs only when the stress is present, or if it is provided by inheritable stress

lP

memory. In addition, the epigenetic modification should be analyzed after the occurrence of the stress if the unstressed condition could be restored after the treatment (irrigation after drought etc.). The presence of the same epigenetic mark in parental and filial generations in absence of stress

na

should prove the transgenerational inheritance excluding its transient nature. Natural epigenetic variation may originate from polymorphisms in transposon insertions and repeats. Epigenetic works

ur

in model plant species such as Arabidopsis thaliana have shown that epigenetic regulation is closely

Jo

associated with functional aspects. However, epigenetic studies on Arabidopsis thaliana have limitations, such as the low value methylated cytosine in the genome (5%) while many crops contain more than 20% [63]. In addition, this plant contains very few transposable elements compared to other crops (review by Lee and Kim, 2014) [65]. Mutations in epigenetic regulators appear to have a greater impact in crops than Arabidopsis. Furthermore, while the distribution of genes along the Arabidopsis chromosomes is fairly homogeneous, this situation may differ in crops. For example, the genome of Vitis vinifera (L.) is characterized by the alternation of large regions 13

with high and low genetic density [66-67]. Several epigenetic studies have been performed in crops such as tomato fruits which constitute an important model for studies on grapes where DNA methylation and histone modifications should play a pivotal role in fruit veraison [68]. In woody crops, epigenetic modifications can be investigated using clonal propagation. The stress applied in the first vegetative generation can be analyzed in the second upcoming vegetation derived from the meristematic apex previously subjected to the studied stress. The propagation through cuttings is a desirable way to identify how the epigenetic modifications is

ro of

stable or transient. Grape has recently been proposed as a model plant to study epigenomics in perennial plants [68]. This is due for the following reasons: 1) it is a model for non-climateric fruits and flower development, 2) it presents flower buds that occur one year after they open, and 3) it is

-p

clonally propagated and commonly grafted.

It is well-known that environmental conditions of the previous year have a significant impact on the

re

flower and subsequently on the development of the fruit. Other specific requirements that affect

lP

flowering and fructification are represented by clonal grafting and propagation. In this context, the epigenetic variability added to the genetic diversity of the grape could shape the phenotypic variations of the plant. The clonal diversity within the grape varieties was identified using the

na

methylation-sensitive amplified polymorphism technique, allowing the understanding of the utility of epigenetic markers in intra-varietal diversity studies [69]. Another important impact on the

ur

epigenetic state of a plant is represented by grafting. The understanding of these mechanisms is

Jo

useful to improve grape tolerance to highly threatening environmental stresses enhancing the product quality. Interestingly, systemic acquired resistance to downy mildew pathogen has been shown to be due to DNA methylation inherited in offspring [70]. Another crop that might be considered as a model for epigenomic studies is tomato due to the easy availability of the mutants. Chemical treatments inducing DNA methylations may be useful for generating epimutations although they might not be stable as genetic mutants. Without the need of using genetic recombination, epimutagenesis may represent an important means of exploring allelic 14

variations and the new combinations of alleles. The analysis of DNA methylation in grape in relation to fruit qualitative traits could represent an important tool of identifying interesting epialleles usable in crop breeding. Specifically, DNA methylation can generate epialleles with various levels of expression, leading to a continuous quantitative variation of a crop trait. Herbaceous auto-gamous species are great candidates to understand if epigenetic modifications have a transient or transgenerational nature. Another model plant for cereals is Brachypodium since it has a rapid life cycle, belongs to the Graminaceae family, is easy to cultivate, and has a small

3. Epigenetic regulation of plant stress responses

ro of

diploid high-quality genome that has been recently sequenced.

-p

The recent progresses in the genome-wide analysis of epigenetic marks have allowed gaining insights into the dynamics and plasticity of plant adaptation to altered environmental conditions.

re

Plants are sessile living organisms and they are particularly characterized by a finely-tuned

lP

adaptation to changing environmental conditions through the trigger or repression of epigenetic mechanisms. Most of the epigenetic features are transiently activated in response to the stress occurrence and reverted when the environmental stimuli are no longer present. However, a certain

na

part of these epigenetic marks is “memorized” by the plant determining stress priming. This phenomenon consists in an improved resistance/response to an environmental stress due to the

ur

activation of responses (priming fingerprinting) which remain latent but prompted to be activated

Jo

again in the presence of a subsequent stress. This induced stress resistance could be inherited in offspring resulting in an improved response to a particular stress that has been previously occurred in parents. Priming is directly linked to the concept of plant stress memory. It is known that plants do not have the same group of germline cells for all their life span as animals do. New germ cells come from meristematic cells during their entire life so that any epigenetic modification that occurrs during vegetative growth may be transmitted to the offspring. Any stable epigenetic information acquired by the genetic structures of the somatic cells might be transmitted to the next generation. 15

The presence of a stress may induce epigenetic heritable traits that represents a ‘memory' inherited in the next generations which helps the colonization of different environments. This imprint, or stress memory, can be defined as the structural, genetic, and biochemical modifications occurred as a consequence of stress exposure rendering a plant more resistant (although it might also be more sensitive in some cases) to future exposure to the same stress factor. The most common histone modification associated with stress memory is H3K4me3 although also H3K4me2 and H3K27me3 have also been found in systemic acquired response and bacterial stress

ro of

(Figure 3). Chromatin modifications, nucleosome positioning, and DNA methylation are playing a key role in plant stress memory and several epigenetic mechanisms have been shown to be responsible for priming of plants [71] (Figure 4). Although priming may compromise plant

-p

productivity in the short-term, for example through a reduction of photosynthesis, it usually induces tolerance to subsequent stress and therefore favors productivity in the long-term perspective [72].

re

Abscisic acid is known to be involved in plant responses to reiterated abiotic stresses and might be

lP

involved in priming [73. In some plant species, it was observed that ABA levels are higher under drought conditions if the plants have previously been subjected to water stress, implying that ABA may be linked to drought stress memory [74]. In this regard, Virlouvet and Fromm (2015) [75]

na

showed that previously stressed plants have the capacity to maintain stomatal apertures that remain partially closed during a recovery period, which reduces transpiration during a subsequent drought

ur

episode. An increased expression of 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (NCED3) and

Jo

ALDEHYDE OXIDASE 3 (AAO3), key modulators of ABA biosynthesis, was linked with drought memory.

3.1 Abiotic stresses Recent studies have demonstrated that chromatin and DNA methylation changes play an important role in resistance and tolerance mechanisms to different abiotic stresses such as cold, salt stress, drought, osmolality, or mineral changes (Figure 5), thereby pointing out the importance of 16

epigenetic regulations in these aspects [77-75]. In addition, genome-wide analyses of histone PTMs have revealed global epigenomic reprogramming in plants under abiotic stresses. Recently, a specialized histone H1 variant was shown to be required for DNA methylation changes linked with stress responses in Arabidopsis [75]. Two DEAD-box RNA helicases are responsible for epigenetic silencing of gene expression leading to inactivation of abiotic stress responses in Arabidopsis [76]. Arabidopsis mutants deficient in genes involved in the RdDM pathway or with defects in maintenance of CHG methylation showed a variable sensitivity of regulation of the stomatal index

ro of

in conditions of low relative humidity [77]. In addition, they showed hypersensitivity to heat stress [78] or high sensitivity to phosphate starvation [79]. These results are consistent with an important

-p

function of the DNA methylation dynamics in the regulation of abiotic stress-responsive genes.

3.1.1 Drought and salt stress

re

Drought causes extensive remodeling of DNA methylation patterns as it was observed in different

lP

plants species 4-81]. Modulation of DNA methylation at repetitive elements seems to be important to modulate the expression of adjacent genes. H3K4me3 modifications under drought stress have been observed [82]. It has been confirmed that H3K4me3 might be a marker for highly expressed

na

genes in stress conditions. In populus, epigenetic modifications have been linked to drought responses, affecting the transcriptomic changes in response to stress [83]. In apple, a variation of

ur

approximately 54%, 38% and 8.5% of methylation respectively at CG, CHG and CHH sequences

Jo

has been linked to changes in drought-sensitivity/tolerance. A part of these epigenetic changes occurred at genes encoding transcription factors (TFs) and transposable elements (TEs). It is worth noting that epigenetic effects of drought stress vary depending on many factors, such as the intensity and duration of the stress, plant genotype or growth phase, and also genomic imprinting left by previous stress episodes. In barley, drought stress-induced genome-wide DNA methylation changes were associated with water deficit [84]. Drought-stressed poplar has showed genome-wide changes in DNA methylation implying that epigenetic mechanisms are involved in adaption of trees 17

to environmental changes [85]. This has been suggested in corn where transposon elements could function as local enhancers for stress responsive genes [86]. Epigenetic responses to salinity are known to play an important role in plant resilience to this type of stress. Salt stress affected DNA methylation in key stress-related genes such as WRKYs in Arabidopsis thaliana [87], poplar [88], rice [89] forest species[90] and wheat [91]. Acetylation, methylation and phosphorylation of histones play an important role in plant responses to salt stress. In addition, salinity triggers the activation of cell wall related genes such as expansin β2 and

ro of

xyloglucan endotransglucosylase in roots and these gene expression changes have been linked with enhanced acetylation of histones in proximity of the promoter and the open reading frames of these genes [91]. An interaction between salt and drought stress has been observed. Salinity promoted

-p

drought tolerance in Arabidopsis through the modification of H3K27me3. A prolonged absence of H3K27me3 induced a transient and quick transcription of high affinity K+transporter 1 [68].

re

Another epigenetic mechanism involved in plant stress memory modulates the accumulation of high

lP

levels of proline, an aminoacid with antioxidant activities during abiotic stresses [50]. The intraspecific variability in salt tolerance in wheat has been associated with changes in differentially methylated regions (DMRs) containing genes closely associated with transposons [92]. The picture

na

of the molecular regulatory mechanisms is complex since hypermethylated regions responsible for silencing of heterochromatin has also been linked to the activity of small RNAs [93]. Other

ur

evidences in different crops have confrimed the role of epigenomic modifications in drought and

Jo

salt responses [94-97].

3.1.2 Metal-related stresses

Manganese oxide nanoparticles provoke DNA hypomethylation in the moss Physcomitrella patens [98]. Epigenetic studies have been conducted in Populus, a model plant for investigating molecular long-term responses under metal excess conditions. Chip-Seq and RNA-seq data were integrated in poplar showing that genes with a H3K4me3 modification were usually upregulated, while genes 18

with a H3K27me3 modification on the 5-UTR were repressed under excess of Zn [99]. The epigenome of leaves of poplar plants in response to AMF inoculation on multi-metal polluted or unpolluted soil was investigated using the MSAP approach. Few changes in cytosine methylation patterns were observed after 4 months after from planting, whereas increased hypomethylation was observed after 6 months from mycorrhizal treatment in the presence of heavy metals [100]. Transgenerational inheritance of changes in DNA methylation patterns in response to heavy metal stress was shown in rice (Oryza sativa L.) [101]. The DNA methylation changes were observed in

ro of

leaves immediately after heavy metal treatment and they consist in CHG hypomethylation; CHGdemethylated states that were heritable to offsprings. MSAP technique was useful to determine the effect of metal contamination on cytosine residues in CCGG motifs in the genome of Red maple

-p

(Acer rubum), a common deciduous tree in northern America. Zinc stress has been associated with enhanced methylation linked to reduced levels of key hormones such as gibberellic acid, zeatin and

re

indole acetic acid and increased levels of abscisic acid [102]. Finally, Aluminum stress altered DNA

retrotransposons [103].

na

3.1.3 Cold stress

lP

methylation in corn and it was associated with polymorphisms in long terminal repeat

Non-coding RNAs are involved in epigenetic regulation of vernalization through the silencing of

ur

FLC even after winter cold. Molecular studies showed that chromatin-remodeling processes are

Jo

responsible for the regulation of FLC expression. A high expression level of FLC has been associated with delayed flowering, whereas flowering is promoted when FLC is repressed by vernalization mediated by non-coding RNA regulation. VERNALIZATION INSENSITIVE 3 (VIN3) was identified as an essential gene to trigger repression of FLC by vernalization. VIN3 encodes a plant homeodomain (PHD) finger protein that is induced only during a cold period. The PHD finger motif in VIN3 is often found in various components of chromatin-remodeling complexes [104]. Vernalization provokes an increased association of PHD-PRC2 with FLC chromatin and an 19

enhanced deposition of H3K27me3 marks at the FLC chromatin [104]. Increased enrichments of PRC2 and H3K27me3 are hallmarks for the stable repression of FLC by vernalization. Another lncRNA involved in vernalization is a cold-inducible intronic lncRNA, COLDAIR [105]. Transcriptional silencing of MET1 gene in corn roots in response to cold causes demethylation of the Ac/Ds transposon system [106]. The Tam3 transposon methylation at CHH sites was observed in snapdragon during cold conditions [107]. Plants have developed mechanisms to survive in an ever-changing environment through the active regulation of chromatin. Particularly, it was shown

ro of

that the interaction of HOS15 with histone deacetylase 2C (HD2C) modulates the action of promoters of cold-responsive COR genes. Epigenetic mechanisms have been observed to play a key role in modulating cold tolerance in several plants such as brassica rapa and sugar beet [108-

-p

109].

re

3.2 Biotic stresses

lP

Plants are sessile organisms that have evolved different mechanisms compared to animals to escape the exposure of adverse environmental factors, including pathogens. Plants have evolved complex mechanisms of protection, among which the ability to establish symbiotic interactions to prevent

na

infections is particularly relevant. Recent researches have evidenced several regulatory mechanisms controlling the plant-microbe interactions. In this context, epigenetic mechanisms are important key

ur

modulators of both pathogenic and symbiotic interactions between plants and microbes (Figure 6).

Jo

DNA methylation is an epigenetic modification of DNA that prevents the loss of gene expression in plant defenses against pathogens and plays a crucial role in microbial-plant interactions. For example, several research reports confirmed the role played by recurrent DNA methylations in the induction of enhanced resistance to pathogens [110-111]. Responses to biotic and abiotic stresses are stored by intracellular molecules that allow plants to respond more efficiently to future environmental stresses. Systemic acquired resistance (SAR) or defense priming is a classic example of memory response to future stresses able to promote hypersensitivity and rapid expression of 20

defensive genes through the action of the NPR1 gene as it was found in Arabidopsis [112]. Regulation of the epigenetic expression of NPR1/NIM1 genes is an important example of H3 histone modifications affecting stress memory in response to biotic stresses [113]. Regulation of gene expression by DNA methylation is crucial to modulate the transcription of genes involved in the biosynthesis of the SA affecting plant defense mechanisms. DNA methylation could mediate the response to biotic stress through the increase of SA and the modulation of the downstream genes [1114]. These data suggested that methylation levels within TEs may dynamically control the

ro of

expressions of transposon and the proximal gene in response to stress in Arabidopsis [115]. Chromatin remodeling plays a crucial role in plant immunity against both necrotrophic and biotrophic pathogens [110]. As previously reported, mutant bon1-1 with loss of- function leads to

-p

constitutive defense response mechanisms to the virulent pathogen Pst DC3000 through the activation of HUB1 and HUB2 and consequently the induction of autoimmune responses [116]. It

re

was found that dimethylated or trimethylated histone H3 Lys 27 (H3K27me2/3) is involved in the

lP

expression of responsive genes to biotic stresses in rice. Jumonji C (jmjC) protein gene encodes a histone Lys demethylase that specifically reverses the activity of H3K27me2/3. During pathogenic infections, the expression of JMJ705 is induced. The overexpression of JMJ705 in transgenic plants

na

removes levels of H3K27me2/3, inducing their expression and enhancing plant resistance to biotic stress through the involvement of jasmonic acid. JMJ704 positively regulates rice tolerance to

ur

bacterial blight through the epigenetic suppression of defense regulators, that is a distinct way from

Jo

its homolog JMJ705 [116]. Furthermore, during plant infection processes, chromatin-based transcriptional regulation can also affect the expression of the effector gene in fungi [117]. Different studies have focused on the concept that pathogen infections induce histone modifications affecting key response genes. Changes in DNA methylation has been shown to affect plant resistance to diseases. For example, significant differences in cytosine methylation have been observed between resistant and sensitive genotypes to Fusarium wilt in chickpea [118]. DNA hypomethylation at specific pericentromeric regions affects priming and biotic stress resistance [119]. Responses to 21

geminivirus infection are modulated by cytosine DNA methylation. In this regard, the interaction between V2 protein of Tomato Yellow Leaf Curl Virus (TYLCV) and the host HDA6 affects the recruitment of MET1 by HDA6, causing hypomethylation of the viral DNA genome and the consequent enhancement in plant susceptibility to TYLCV infection [120]. DNA methylation and histone modifications are involved in priming since the modification of chromatin structure modifies the genetic responses to a subsequent attack. It has been shown that MeJA promotes histone modifications in the promotion of defense-related genes, as well as DNA methylation

ro of

modifications in the genome. This supports the hypothesis that epigenetic regulation plays a key role in priming onset. Transposable elements are involved in plant immunity as shown by the transposon ATCOPIA93 long-terminal repeat (LTR) that is negatively regulated by DNA

-p

methylation and modifies immune-responsive genes during pathogen attacks [121]. RNA-directed DNA-methylation (RdDM) inhibits transcription through post-transcriptional gene silencing driven

re

by Dicer-like (DCL) proteins that are responsible for the formation of small-interfering RNAs

expression,

while

lP

(siRNAs). Indeed, hypomethylation is linked to biotic stress, which is linked with enhanced DCL4 hypermethylation

is

associated

with

DCL4

down-regulation

[122].

Hypomethylation of centromeric regions and terminal heterochromatic blocks have been observed

na

in response to infections of virulent and avirulent strains of the fungus Colletotrichum lindemuthianum [123]. Finally, it has been demonstrated that histone methyltransferase plays an

ur

important role in plant defenses against fungal pathogens through the modulation of genes involved

Jo

in the jasmonic acid (JA) and/or ethylene signaling pathway. Methylation of histone H3 lysine 36 mediated by the histone methyltransferase SET DOMAIN GROUP8 (SDG8) may work as a memory of permissive expression of several defense genes promoting priming [124].

4

A proposed experimental plan for epigenomics in crop stress biology

A typical experimental plan to identify epigenetic marks linked with plant abiotic stress responses might be the one shown in Figure S2. Thanks to the agronomic evaluations, we should be able to 22

select a resistant and a susceptible genotype to a particular stress (i.e. drought, salinity etc..) in a well-known studied crop such as bean. Both DNA methylations and chromatin remodeling marks might be identified at genome-wide levels using techniques such as whole genome bisulphite sequencing (WGBS) and ChIP-seq. A typical experimental design might provide the use of two genotypes for three stress condition (control, drought, salinity). If we analyze two tissues (roots and leaves), one time point (early stress) and three biological replicates we would have a total of 36 samples. Transcriptomic analysis will provide the list of the differentially regulated genes and will

ro of

guide the epigenomic analysis focusing only on key genes affecting resistance/susceptibility to the studied stress. The integration of epigenomic and transcriptomic allows filtering of this huge preliminary list of epigenetic patterns focusing on those genes that have been previously known to

-p

play a key role in resistance to the studied stress. A reasonable number of genes (i.e. 25-50) might be analyzed with a targeted epigenetic approach. RNA-seq using the same samples analysed at

re

epigenomic level (one resistant and one susceptible for each stress) should be employed. The

lP

analysis of more than one time point is preferred for transcriptomics due to the low cost of analysis compared to epigenomics. After read trimming, mapping, annotation, and differential expression analysis, the list of the genes differentially regulated by each stress will be determined for each

na

analyzed genotype. A possible targeted approach might to perform PCRs coupled with bisulphite sequencing on the filial generation comparing progenies of stressed and unstressed parents. It is

ur

expected that most of the epigenetic marks should be linked to environmental variability. The

Jo

identification of those epigenetic marks strictly linked to stress resistance is hard. An intensive bibliographic research is required in order to identify DNA methylation patterns involved in the modulation of stress responses. Particular attention should be paid to several gene categories typically affected by abiotic stresses such as transcription factors (WRKYs, MYBs etc.), biosynthesis of water soluble molecules (proline, sugars, etc.), hormone biosynthesis, signaling and response, secondary metabolism, protein folding, stress signaling receptors, and antioxidant pathways. The use of the same samples for both epigenomic and transcriptomic should be allowed 23

to reduce environmental variability external to the studied stress. Proteomics analysis would be another “omic” approach that is desirable to include when trying to identifying genome-wide epigenetic marks. Isobaric tags for relative and absolute quantification systems (iTRAQ) would allow analysis of at least 4000 proteins/sample. Proteins should be extracted by powdered tissues using procedures that allow extraction of a high amount of high-quality proteins. The digested peptides are analyzed using a QExactive mass spectrometer coupled with an Easy-LC and a nanospray ionization source. Data are acquired using a data-dependent ms/ms method. Raw data

ro of

might be analyzed using X!Tandem peptide sequences, which are annotated when searching for them against Uniprot databases. The list of differentially abundant proteins between genotypic and control/stress comparisons are compared with the list of genes identified from transcriptomic

-p

analysis. Functional data mining tools for proteomic data analysis might be the same used for RNAseq data (i.e. Mapman, Pageman, DAVID, STRING etc..). Considering that the resolution power of

re

proteomics is lower than RNA-seq, some genes found in RNA-seq analysis will not show

lP

correspondent proteins in the proteomic data. In addition, protein extraction procedures do not allow us to obtain sufficient amounts of proteins to be analyzed by iTRAQ. Indeed, it is advisable to not exclude epigenetic marks related to those genes that do not show a correspondent change in

na

protein abundance. This is due to the inability of proteomic technologies to analyze the entire

Crop epigenetic adaptation to environment

Jo

5

ur

proteome.

Epigenetic systems should play an important role in sensing any environmental perturbation modulating gene expression. Rapid epigenetic changes (transient or stably inherited) should be the key for the high flexibility of plant responses to the environmental changes [125]. Epigenetic mechanisms control genomic structure, particularly the heterochromatic state of centromeric telomeric regions. Epigenetics together with genetic variation are probably the most involved events in phenotypic diversity and in the plasticity of plants. Epigenetic variation depends strictly 24

on environmental inputs. Dynamic environmental or experimental conditions induce the generation of epialleles that may a have a different phenotype from the one subjected to alleles generated by genetic variation. The epigenetic variants can be inherited and lead to better adaptation to the environment. Examples of experimentally induced epialleles were observed in Arabidopsis by generating Epigenetic Recombinant Inbred Line populations (EpiRIL) derived from the decrease in DNA methylation 1-2 (ddm1-2) [126]. These EpiRILs have been shown to exhibit changes in growth capacity and are markedly susceptible to salinity stress compared to the Col0 parent line,

ro of

demonstrating that the epigenotypes derived from decrease in DNA methylation influence the ability to adapt to salt stress [127]. Methylation processes have also been observed in natural plants and may strongly depend on specific environmental traits. Therefore, the induced or natural

-p

epigenetic diversity could represent an unexplored resource of phenotypic variation and constitute a useful means in plant breeding programs. Polyploidy plays a prominent role in plant evolution,

re

especially in flowering plants. It is an important evolutionary process in plants where hybridization

lP

and chromosome doubling produce epigenetic modifications related to retrotransposons. In allopolyploids, interspecific hybridization may promote epigenetic modifications that may be stabilized and transmitted as epialleles into the filial generation progeny, which are subject to

na

natural selection and adaptation. LTR retrotransposons are the most abundant type of transposable elements (TEs) in plants. The presence of cis-regulatory motifs in the LTRs plays a pivotal role in

ur

stress responses and in the complex evolutionary mechanisms occurring in plants. The insertions of

Jo

LTR retrotransposons close to genes have been linked to alternative splicing and epigenetic control. These characteristics can become an active part of the evolution of gene families actively involved in the evolution such as resistance genes. Epigenetic modifications are induced by stresses causing genome modification through the action of retrotransposons that are moving in the genome generating rearrangements that can be stably inherited. Climate change due to increased emission of greenhouse gases should negatively affect plant diversity and productivity. This adverse effect might be reduced through plant adaptive mechanisms to abiotic stresses such as changes in the 25

production of reactive oxygen species triggering signal transduction networks that provide protection against environmental stresses [127].

6

Conclusions

A greater elucidation of the epigenetic modifications in both perennial and autogamous species would allow to decipher the complex interactions between genotype and environment, identifying the real impact of global climate changes in the adaptation and consequently in the quality and

ro of

productivity of commercial crop species. In addition, the construction of crop methylomes would help to develop new molecular markers associated with epigenetic changes that will assist breeding activities. The next-generation markers might be integrated in current methods of genetic

-p

characterization such as those using DNA barcoding [128]. It will allow the identification of transient stress conditions to manage plant stresses in a cautious and sustainable way. Indeed, we

re

propose that an “epigenetic code” should exist in addition to the well-known genetic code

lP

composed by 4 bases. The existence of an “epigenetic alphabeth” might explain the low correlation observed between genetic and phenotypic variance. Key players in the epigenetic code should be represented by DNA methylation in CG, CH, and CHH sequences. Other epigenetic mechanisms

na

driven by chromatin restructuring should be driven by histone methylations such as H3K27me2/3 occurring in response to bacterial stress. The deciphering of this epigenetic alphabet is one of the

ur

most fascinating fields of modern genetics. This may greatly help to gain insight into complex

Jo

mechanisms of plant epigenomic responses to environmental stresses in a climate change scenario.

AUTHOR CONTRIBUTIONS Both FM and AP contributed in writing the manuscript. Both authors have made substantial, direct and intellectual contribution to the work.

26

Conflict of Interest None

References 1. A. Berr, S. Shafiq, W.H. Shen, Histone modifications in transcriptional activation during plant development, Biochim. Biophys. Acta 1809 (2011) 567–576.

ro of

2. C.S. Pikaard, O.M. Scheid, Epigenetic regulation in plants, Cold Spring Harb. Perspect. Biol. 6 (2015) a019315. doi: 10.1101/cshperspect. a019315.

3. A. Sokol, A. Kwiatkowska, A. Jerzmanowski, M. Prymakowska Bosak, Up-regulation of

-p

stress-inducible genes in tobacco and Arabidopsis cells in response to abiotic stresses and ABA treatment correlates with dynamic changes in histone H3 and H4 modifications,

re

Planta. 227 (2007) 245-54.

lP

4. V.V. Sridhar, A. Kapoor, K. Zhang, J. Zhu, T. Zhou, P.M. Hasegawa, R.A. Bressan, J.K. Zhu, Control of DNA methylation and heterochromatic silencing by histone H2B deubiquitination, Nature 447 (2007) 735–73.

na

5. G.M. Silva, D. Finley, C. Vogel, K63 polyubiquitination is a new modulator of the oxidative stress response, Nat Struct Mol Biol. 22 (2015) 116-23. doi: 10.1038/nsmb.2955.

ur

6. A. Sokol, A. Kwiatkowska, A. Jerzmanowski, M. Prymakowska-Bosak, Up-regulation of

Jo

stress-inducible genes in tobacco and Arabidopsis cells in response to abiotic stresses and ABA treatment correlates with dynamic changes in histone H3 and H4 modifications, Planta. 227(2007) 245-54.

7. J.Y. Kim, W. Yang, J. Forner, J.U. Lohmann, B. Noh, Y.S. Noh, Epigenetic reprogramming by histone acetyltransferase HAG1/AtGCN5 is required for pluripotency acquisition in Arabidopsis, EMBO J. 37 (2018) (20). pii: e98726. doi: 10.15252/embj.201798726.

27

8. M. Pellegrini, R. Ferrari, Epigenetic analysis: ChIP-chip and ChIP-seq, Methods Mol. Biol. 802 (2012) 377-87. doi: 10.1007/978-1-61779-400-1_25. 9. X. Chen, V. Bhadauria, B. Ma, ChIP-Seq: A Powerful Tool for Studying Protein-DNA Interactions

in

Plants,

Curr.

Issues

Mol.

Biol.

27

(2018)

171-180.

doi:

10.21775/cimb.027.171. Review. 10. X. Wang, D. Moazed, DNA sequence-dependent epigenetic inheritance of gene silencing and histone H3K9 methylation, Science. 356 (2017) 88-91. doi: 10.1126/science.aaj2114.

ro of

11. C.J. Harris, M. Scheibe, S.P. Wongpalee, W. Liu, E.M. Cornett, R.M. Vaughan, X. Li, W. Chen, Y. Xue, Z. Zhong, L. Yen, W.D. Barshop, S. Rayatpisheh, J. Gallego-Bartolome, M. Groth, Z. Wang, J.A. Wohlschlegel, J. Du, S.B. Rothbart, F. Butter, S.E. Jacobsen, A DNA

-p

methylation reader complex that enhances gene transcription, Science. 362 (2018) 11821186. doi: 10.1126/science.aar7854.

re

12. C.E. Niederhuth, A.J. Bewick, L. Ji, M.S. Alabady, K.D. Kim, Q. Li, N.A. Rohr, A.

Springer, R.J.

lP

Rambani, J.M. Burke, J.A. Udall, C. Egesi, J. Schmutz, J., Grimwood, S.A. Jackson, N.M. Schmitz, Widespread natural variation of DNA methylation within

angiosperms, Genome Biol. 17 (2016) 194.

na

13. J.P. Jackson, A.M. Lindroth, X. Cao, S.E. Jacobsen, Control of CpNpG DNA methylation by the KRYPTONITE histone H3 methyltransferase, Nature 416 (2002) 556–560.

ur

14. A.M. Lindroth, D. Shultis, Z. Jasencakova, J. Fuchs, L. Johnson, D. Schubert, D. Patnaik,

Jo

S. Pradhan, J. Goodrich, I. Schubert, T. Jenuwein, S. Khorasanizadeh, S.E. Jacobsen, Dual histone H3 methylation marks at lysines 9 and 27 required for interaction with CHROMOMETHYLASE3, EMBO. J. 23 (2004) 4286-96. Epub 2004 Sep 30.

15. M.A. Matzke, R.A. Mosher, RNA-directed DNA methylation: an epigenetic pathway of increasing complexity, Nat. Rev. Genet. 15 (2014) 394-408. doi: 10.1038/nrg3683. 16. J.A. Jeddeloh, J. Bender, E.J. Richards, The DNA methylation locus DDM1 is required for maintenance of gene silencing in Arabidopsis, Genes Dev. 12 (1998) 1714-25. 28

17. Z. Gong, T. Morales-Ruiz, R.R. Ariza, T. Roldán-Arjona, L. David, J.K. Zhu, ROS1, a repressor

of

transcriptional

gene

silencing

in

Arabidopsis,

encodes

a

DNA

glycosylase/lyase, Cell. 111 (2002) 803-14. 18. M. Crampton, V. R. Sripathi, K.Hossain, V. Kalavacharla, Analyses of Methylomes Derived from Meso-American Common Bean (Phaseolus vulgaris L.) Using MeDIP-Seq and Whole Genome Sodium Bisulphite-Sequencing, Front Plant Sci. 7 (2016) 447. Published online 2016 Apr 26. doi: 10.3389/fpls.2016.00447.

ro of

19. M. Margulies, M Egholm, W.E. Altman, S. Attiya, J.S. Bader, L.A Bemben, J. Berka, M.S. Braverman, Y.J. Chen, Z. Chen, S.B. Dewell, L. Du, J.M. Fierro, X.V. Gomes, B.C. Godwin, W. He, S. Helgesen, C.H. Ho, G.P. Irzyk, S.C. Jando, M.L. Alenquer, T.P. Jarvie,

-p

K.B. Jirage, J.B. Kim, J.R. Knight, J.R. Lanza, J.H. Leamon, S.M. Lefkowitz, M. Lei, J. Li, K.L. Lohman, H. Lu, V.B. Makhijani, K.E. McDade, M.P. McKenna, E.W. Myers, E.

re

Nickerson, J.R. Nobile, R. Plant, B.P. Puc, M.T. Ronan, G.T. Roth, G.J. Sarkis, J.F. Simons,

lP

J.W. Simpson, M. Srinivasan, K.R. Tartaro, A. Tomasz, K.A. Vogt, G.A. Volkmer, S.H. Wang, Y. Wang, M.P. Weiner, P. Yu, R.F. Begley, J.M. Rothberg, Genome sequencing in microfabricated high-density picolitre reactors, Nature. 437 (2005) 376–380.

na

20. S. Bennett, Solexa Ltd, Pharmacogenomics 5 (2004) 433–438. 21. S. Bennett, C. Barnes, A. Cox, L. Davies, C. Brown, Toward the 1,000 dollars human

ur

genome, Pharmacogenomics. 6 (2005) 373–382.

Jo

22. P. Milos, Helicos BioSciences, Pharmacogenomics. 9 (2008) 477–480. 23. J. Eid, A. Fehr, J. Gray, K. Luong, J. Lyle, G. Otto, P. Peluso, D. Rank, P. Baybayan, B. Bettman, A. Bibillo, K. Bjornson, B. Chaudhuri, F. Christians, R. Cicero, S. Clark, R. Dalal, A. Dewinter, J. Dixon, M. Foquet, A. Gaertner, P. Hardenbol, C. Heiner, K. Hester, D. Holden, G. Kearns, X. Kong, R. Kuse, Y. Lacroix, S. Lin, P. Lundquist, C. Ma, P. Marks, M. Maxham, D. Murphy, I. Park, T. Pham, M. Phillips, J. Roy, R. Sebra, G. Shen, J. Sorenson, A. Tomaney, K. Travers, M. Trulson, J. Vieceli, J. Wegener, D. Wu, A. Yang, D. 29

Zaccarin, P. Zhao, F. Zhong, J. Korlach, S. Turner, Real-time DNA sequencing from single polymerase molecules, Science 323 (2009) 133–138. 24. J.M. Rothberg, W. Hinz, T.M. Rearick, J. Schultz, W. Mileski, M. Davey, . et al. (2011) J.H. Leamon, K. Johnson, M.J. Milgrew, M. Edwards, J. Hoon, J.F. Simons, D. Marran, J.W. Myers, J.F. Davidson, A. Branting, J.R. Nobile, B.P. Puc, D. Light, T.A. Clark, M. Huber, J.T. Branciforte, I.B. Stoner, S.E. Cawley, M. Lyons, Y. Fu, N. Homer, M. Sedova , X. Miao, B. Reed, J. Sabina, E. Feierstein, M. Schorn, M. Alanjary, E. Dimalanta, D.

ro of

Dressman, R. Kasinskas, T. Sokolsky, J.A. Fidanza, E. Namsaraev, K.J. McKernan, A. Williams, G.T. Roth, J. Bustillo, An integrated semiconductor device enabling non-optical genome sequencing, Nature 475 (2011) 348–352.

-p

25. A.G. Davis-Richardson, J.T. Russell, R. Dias, A.J. McKinlay, R. Canepa, J.R. Fagen, K.T. Rusoff, J.C. Drew, B. Kolaczkowski, D.W. Emerich, E.W. Triplett, Integrating DNA

re

methylation and gene expression data in the development of the soybean-Bradyrhizobium

lP

N2-fixing symbiosis, Front. Microbiol. 7 (2016) 518. doi: 10.3389/fmicb.2016.00518. 26. A. Yu, G. Lepère, F. Jay, J. Wang, L. Bapaume, Y. Wang, A.L. Abraham, J. Penterman, R.L. Fischer, O. Voinnet, L. Navarro, Dynamics and biological relevance of DNA

na

demethylation in Arabidopsis antibacterial defense, Proc. Natl. Acad. Sci. U S A. 110 (2013) 2389-94. doi: 10.1073/pnas.1211757110.

ur

27. H. Stroud, M.V. Greenberg, S. Feng, Y.V. Bernatavichute, S.E. Jacobsen, Comprehensive

Jo

analysis of silencing mutants reveals complex regulation of the Arabidopsis methylome, Cell. 152 (2013) 352-64. doi: 10.1016/j.cell.2012.10.054.

28. Ö. Çelik, S. Meriç, A. Ayan, Ç. Atak, Epigenetic analysis of WRKY transcription factor genes in salt stressed rice (Oryza sativa L.) plants, Environ. Exp. Bot. 159 (2019) 121–131. doi: 10.1016/j.envexpbot.2018.12.015. 29. N. Olova, F. Krueger, S. Andrews, D. Oxley, R.V. Berrens, M.R. Branco, W. Reik, Comparison of whole-genome bisulphite sequencing library preparation strategies identifies 30

sources of biases affecting DNA methylation data, Genome Biol. 19 (2018) 33. doi: 10.1186/s13059-018-1408-2. 30. T.C. Mockler, S. Chan, A. Sundaresan, H. Chen, S.E. Jacobsen, J.R. Ecker, Applications of DNA

tiling

arrays

for

whole-genome

analysis,

Genomics.

85

(2005)

1–15.

doi:10.1016/j.ygeno.2004.10.005. 31. A. Meissner, A. Gnirke, G.W. Bell, B. Ramsahoye, E.S. Lander, R. Jaenisch, Reduced representation bisulphite sequencing for comparative high-resolution DNA methylation

ro of

analysis, Nucleic Acids Res. 33 (2005) 5868-77. Print 2005. 32. A. Agorio, S. Durand, E. Fiume, C. Brousse, I. Gy, M. Simon, S. Anava, O. Rechavi, O. Loudet, C. Camilleri, N. Bouché, An Arabidopsis Natural Epiallele Maintained by a Feed-

-p

Forward Silencing Loop between Histone and DNA, PLoS Genet.13 (2017) e1006551. doi: 10.1371/journal.pgen.1006551.

re

33. M. Zhang, J.L. Chen, S.K. Liang, G.H. Li, F.H. Wang, I. Ahmad, Differentially methylated

lP

genomic fragments related with sexual dimorphism of rice pests, Sogatella furcifera, Insect. Sci. 22 (2015) 731-8. doi: 10.1111/1744-7917.12179. 34. C.M. Rodríguez López, B. Guzmán Asenjo, A.J. Lloyd, M.J. Wilkinson, Direct detection

na

and quantification of methylation in nucleic acid sequences using high-resolution melting analysis, Anal. Chem. 82 (2010) 9100-8. doi: 10.1021/ac1024057.

ur

35. M.L. Gonzalgo, G. Liang, Methylation-sensitive single-nucleotide primer extension (Ms-

Jo

SNuPE) for quantitative measurement of DNA methylation, Nat. Protoc. 2 (2007) 1931-6. 36. S. Kunze, Quantitative Region-Specific DNA Methylation Analysis by the EpiTYPER™ Technology, Methods Mol. Biol. 1708 (2018) 515-535. doi: 10.1007/978-1-4939-74818_26. 37. B.T. Wilhelm, S. Marguerat, S. Watt, F. Schubert, V. Wood, I. Goodhead, C.J. Penkett, J. Rogers, J. Bahler, Dynamic repertoire of a eukaryotic transcriptome surveyed at singlenucleotide resolution, Nature. 453 (2008) 1239–1243. 31

38. A.Weiberg, M. Wang, F.M. Lin, H. Zhao, Z. Zhang, I. Kaloshian, H.D. Huang, H. Jin, Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways, Science. 342 (2013) 118-23. doi: 10.1126/science.1239705. 39. X. Sun, H. Zheng, N. Sui, Regulation mechanism of long non-coding RNA in plant response to

stress,

Biochem.

Biophys.

Res.

Commun.

503

(2018)

402-407.

doi:

10.1016/j.bbrc.2018.07.072. 40. M. Matzke, T. Kanno, L. Daxinger, B. Huettel, A.J. Matzke, RNA-mediated chromatin-

ro of

based silencing in plants, Curr. Opin. Cell Biol. 21 (2009) 367–376. 41. J.A. Law, S.E. Jacobsen, Establishing, maintaining and modifying DNA methylation patterns in plants and animals, Nat. Rev. Genet. 11 (2010) 204–220.

-p

42. G. He, A.A. Elling, X.W. Deng, The epigenome and plant development. Annu. Rev. Plant Biol. 62 (2011) 411–435.

re

43. R. Lister, B.D. Gregory, J.R. Ecker, Next is now: new technologies for sequencing of

lP

genomes, transcriptomes, and beyond. Curr. Opin. Plant Biol. 12 (2009) 107–118. 44. S.J. Cokus, S. Feng, X. Zhang, Z. Chen, B. Merriman, C.D. Haudenschild, S. Pradhan, S.F. Nelson, M. Pellegrini, S.E. Jacobsen, Shotgun bisulphite sequencing of the Arabidopsis

na

genome reveals DNA methylation patterning, Nature. 452 (2008) 215–219. 45. K. Fujita, K. Matsumoto, G.K. Ofosu-Budu, S. Ogata, Effect of Shading on Growth and

ur

Dinitrogen Fixation of Kudzu and Tropical Pasture Legumes, Soil Sci, Plant Nutr., 39

Jo

(1993) 43-54.

46. P. Xu, W. Cai, Function of Brassica napus BnABI3 in Arabidopsis gs1, an Allele of AtABI3, in Seed Development and Stress Response, Front. Plant Sci. 10 (2019) 67. doi: 10.3389/fpls.2019.00067. 47. F. Martinelli, G. Cannarozzi, B. Balan, F. Siegrist, A. Weichert, R. Blösch, Z. Tadele, Identification of miRNAs linked with the drought response of tef [Eragrostis tef (Zucc.) Trotter], J. Plant Physiol. 224-225 (2018) 163-172. doi: 10.1016/j.jplph.2018.02.011. 32

48. A. Matsui, J. Ishida, T. Morosawa, M. Okamoto, J.M. Kim, Y. Kurihara, M. Kawashima, M. Tanaka, T.K. To, K. Nakaminami, E. Kaminuma, T.A. Endo, Y. Mochizuki, S. Kawaguchi, N. Kobayashi, K. Shinozaki, T. Toyoda, M. Seki, Arabidopsis tiling array analysis to identify the stress-responsive genes, Methods Mol. Biol. 639 (2010) 141-55. doi: 10.1007/978-1-60761-702-0_8. 49. X. Huang, S. Li, P. Ni, Y. Gao, B. Jiang, Z. Zhou, A. Zhan, Rapid response to changing environments during biological invasions: DNA methylation perspectives, Mol Ecol. 26

ro of

(2017) 6621-6633. doi: 10.1111/mec.14382. Epub 2017 Nov 15. 50. G.H. Salekdeh, S. Komatsu, Crop proteomics: aim at sustainable agriculture of tomorrow, Proteomics. 7(2007) 2976-96. Review.

-p

51. T.R. Mercer, D.J. Gerhardt, M.E. Dinger, J. Crawford, C. Trapnell, J.A. Jeddeloh, J.S. Mattick, J.L. Rinn, Targeted RNA sequencing reveals the deep complexity of the human

re

transcriptome, Nat. Biotechnol. 30 (2011) 99–104.

lP

52. X. Zhang, Y.V. Bernatavichute, S. Cokus, M. Pellegrini, S.E Jacobsen, Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana, Genome Biol. 10( 2009) R62. doi: 10.1186/gb-2009-10-6-r62.

na

53. A. Nesvizhskii, Proteogenomics: Concepts, applications and computational strategies. Nature Methods 11 (2014) 1114-1125. doi: 10.1038/NMETH.3144.

ur

54. C. Voelckel, N. Gruenheit, P. Lockhart, Evolutionary transcriptomics and proteomics: into

plant

adaptation,

Trends

Plant

Sci.

22

(2017)

462-471.

doi:

Jo

insight

10.1016/j.tplants.2017.03.001.

55. R. Kumar, A. Bohra, A.K. Pandey, M.K. Pandey, Metabolomics for plant improvement: Status and prospects, Frontiers Plant Sci. 1302 (2017) doi: 10.3389/fpls.2017.01302. 56. R. Sunkar, V. Chinnusamy, J. Zhu, J.K. Zhu, Small RNAs as big players in plant abiotic stress responses and nutrient deprivation, Trends Plant Sci. 12(2007) 301-9.

33

57. U.J Phukan, G. S. Jeena, R. K. Shukla, WRKY Transcription Factors: Molecular Regulation and Stress Responses in Plants, Front Plant Sci. 7 (2016) 760. doi: 10.3389/fpls.2016.00760. 58. A.F.A. Samad, M. Sajad, N. Nazaruddin, I.A. Fauzi, A.M.A. Murad, Z. Zainal, I. Ismail, MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network, Front Plant Sci. 8 (2017) 565. doi: 10.3389/fpls.2017.00565. 59. S.P. Pandey, I.E. Somssich, The role of WRKY transcription factors in plant immunity, Plant Physiol. 150( 2009) 1648-55. doi: 10.1104/pp.109.138990. Review.

ro of

60. R. Tosetti, F. Martinelli, P. Tonutti, Metabolomics Approach to Studying Minimally Processed Peach (Prunus persica) Fruit, Acta Hortic. 934 (2012) 1017-1022.

61. S. Saia, P. Ruisi, V. Fileccia, G. Di Miceli, G. Amato, F. Martinelli, Metabolomics suggests

-p

that soil inoculation with arbuscular mycorrhizal fungi decreased free amino acid content in roots of durum wheat grown under N-limited, P-rich field conditions, Plos One 10 (6)

re

(2015) e0129591. doi: 10.1371/journal.pone.0129591.

lP

62. F. Martinelli, P. Tonutti, Flavonoid metabolism and gene expression in developing olive (Olea europaea L.) fruit, Plant Biosyst. 146 (2011) 164-170. 63. J. Hong, L. Yang, D. Zhang, J. Shi, Plant metabolomics: An indispensable system biology

na

tool for plant science, Int. J. Mol. Sci. 17 (2016) 767. doi: 10.3390/ijms17060767 64. M. Faurobert, E. Pelpoir, J. Chaïb, Phenol extraction of proteins for proteomic studies of

ur

recalcitrant plant tissues, Methods Mol. Biol. 355 (2007) 9-14.

Jo

65. S.I. Lee, N.S. Kim, Transposable elements and genome size variations in plants, Genomics Inform. 12 (2014) 87–97. 10.5808/GI.2014.12.3.87.

66. O. Jaillon, J.M. Aury, B. Noel, A. Policriti, C. Clepet, A. Casagrande, N. Choisne, S. Aubourg, N. Vitulo, C. Jubin, A. Vezzi, F. Legeai, P. Hugueney, C. Dasilva, D. Horner, E. Mica, D. Jublot, J. Poulain, C. Bruyère, A. Billault, B. Segurens, M. Gouyvenoux, E. Ugarte, F. Cattonaro, V. Anthouard, V. Vico, C. Del Fabbro, M. Alaux, G. Di Gaspero, V. Dumas, N. Felice, S. Paillard, I. Juman, M. Moroldo, S. Scalabrin, A. Canaguier, I. Le 34

Clainche, G. Malacrida, E. Durand, G. Pesole, V. Laucou, P. Chatelet, D. Merdinoglu, M. Delledonne, M. Pezzotti, A. Lecharny, C. Scarpelli, F. Artiguenave, M.E. Pè, G. Valle, M. Morgante, M. Caboche, A.F. Adam-Blondon, J. Weissenbach, F. Quétier, P. Wincker P; French-Italian Public Consortium for Grapevine Genome Characterization, The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla, Nature. 449 (2007) 463-7. 67. A.M. Fortes, P. Agudelo-Romero, M.S. Silva, K. Ali, L. Sousa, F. Maltese, Y.H. Choi, J.

ro of

Grimplet, J.M. Martinez-Zapater, R. Verpoorte, Transcript and metabolite analysis in Trincadeira cultivar reveals novel information regarding the dynamics of grape ripening, BMC Plant Biol. 11 (2011) 149. doi: 10.1186/1471-2229-11-149.

-p

68. A.M. Fortes, P. Gallusci, Plant stress responses and phenotypic plasticity in the epigenomics era: Perspectives on the grapevine scenario, a model for perennial crop plants, Frontiers

re

Plant Sci. 82 (2014) 82. doi: 10.3389/fpls.2017.00082.

lP

69. S. Amoah, S. Kurup, C.M. Rodriguez Lopez, S.J. Welham, S.J. Powers, C.J. Hopkins, A hypomethylated population of Brassica rapa for forward and reverse Epi-genetics, BMC Plant Biol. (2012) 12:193. doi: 10.1186/1471-2229- 12-193.

na

70. Ocaña, B. Walter, P. Schellenbaum, Stable MSAP markers for the distinction of Vitis vinifera cv Pinot noir clones, Mol. Biotechnol. 55 (2013) 236-48. doi: 10.1007/s12033-013-

ur

9675-3.

Jo

71. E. Luna, J. Ton, The epigenetic machinery controlling transgenerational systemic acquired resistance, Plant Signal Behav. 7 (2012) 615-618. doi: 10.4161/psb.20155.

72. J. Lämke, I. Bäurle, Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants, Genome Biol. 18 (2018) (1):124. doi: 10.1186/s13059-017-1263-6. Review. 73. T.J.A. Bruce, M.C. Matthes, J.A. Napier, J.A.Pickett, Stressful “memories” of plants: evidence and possible mechanisms, Plant Sci. 173 (2007) 603–608. 35

74. E. Fleta-Soriano, M. Pintó-Marijuan, S. Munné-Bosch, Evidence of Drought Stress Memory in the Facultative CAM, Aptenia cordifolia: Possible Role of Phytohormones, PLoS One. 10 (2015) e0135391. doi: 10.1371/journal.pone.0135391. 75. L. Virlouvet, M. Fromm, Physiological and transcriptional memory in guard cells during repetitive dehydration stress, New Phytol. 205 (2015) 596-607. doi: 10.1111/nph.13080. 76. E. Sani,

P. Herzyk, G. Perrella, V. Colot, A. Amtmann,

Hyperosmotic priming of

Arabidopsis seedlings establishes a long-term somatic memory accompanied by specific

ro of

changes of the epigenome, Genome Biol. 14 (2013) R59. doi: 10.1186/gb-2013-14-6-r59. 77. J.M. Kim, T. Sasaki, M. Ueda, K. Sako, M. Seki, Chromatin changes in response to drought, salinity, heat, and cold stresses in plants, Front. Plant Sci. 6 (2015) 114.

-p

doi:10.3389/fpls.2015.00114.

78. T. Liu, Y. Li, W. Duan, F. Huang, X. Hou, Cold acclimation alters DNA methylation

re

patterns and confers tolerance to heat and increases growth rate in Brassica rapa, J. Exp. Bot. 68 (2017) 1213-1224. doi: 10.1093/jxb/erw496.

lP

79. K. Rutowicz, M. Puzio, J. Halibart-Puzio, M. Lirski, M. Kotliński, M.A. Kroteń, L. Knizewski, B. Lange, A. Muszewska, K. Śniegowska-Świerk, J. Kościelniak, R. Iwanicka-

na

Nowicka, K. Buza, F. Janowiak, K. Żmuda, I. Jõesaar, K. Laskowska-Kaszub, A. Fogtman, H. Kollist, P. Zielenkiewicz, J. Tiuryn, P. Siedlecki, S. Swiezewski, K. Ginalski, M.

ur

Koblowska, R. Archacki, B. Wilczynski, M. Rapacz, A. Jerzmanowski, A Specialized

Jo

Histone H1 Variant Is Required for Adaptive Responses to Complex Abiotic Stress and Related DNA Methylation in Arabidopsis, Plant Physiol. 169 (2015) 2080-101. doi: 10.1104/pp.15.00493.

80. A. Khan, A. Garbelli, S. Grossi, A. Florentin, G. Batelli, T. Acuna, G. Zolla, Y. Kaye, L.K. Paul, J.K. Zhu, G. Maga, G. Grafi, S. Barak, The Arabidopsis STRESS RESPONSE SUPPRESSOR DEAD-box RNA helicases are nucleolar- and chromocenter-localized

36

proteins that undergo stress-mediated relocalization and are involved in epigenetic gene silencing, Plant J. 79 (2014) 28-43. doi: 10.1111/tpj.12533. 81. P.J. Tricker, J.G.

Gibbings, C.M. Rodríguez López, P. Hadley, M.J. Wilkinson, Low

relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development, J Exp. Bot. 63 (2012) 3799-813. doi: 10.1093/jxb/ers076. 82. O.V. Popova, H.Q. Dinh, W. Aufsatz, C. Jonak, The RdDM pathway is required for basal

ro of

heat tolerance in Arabidopsis, Mol. Plant. 6 (2013) 396-410. doi: 10.1093/mp/sst023. Epub 2013 Feb 1.

83. L. Yong-Villalobos, S.I. González-Morales, K. Wrobel, D. Gutiérrez-Alanis, S.A.

-p

Cervantes-Peréz, C. Hayano-Kanashiro, A. Oropeza-Aburto, A. Cruz-Ramírez, O. Martínez, L. Herrera-Estrella, Methylome analysis reveals an important role for epigenetic changes in

re

the regulation of the Arabidopsis response to phosphate starvation, Proc. Natl. Acad. Sci. U

lP

S A. 112 (2015) E7293-302. doi: 10.1073/pnas.1522301112. 84. K. Chwialkowska, U. Nowakowska , A. Mroziewicz, I. Szarejko, M. Kwasniewski, Waterdeficiency conditions differently modulate the methylome of roots and leaves in barley

na

(Hordeum vulgare L.), J. Exp. Bot. 67 (2016) 1109-21. doi: 10.1093/jxb/erv552. 85. D. Liang, Z. Zhang, H. Wu, C. Huang, P. Shuai, C.Y. Ye, S. Tang, Y. Wang, L. Yang, J.

ur

Wang, W. Yin, X. Xia, Single-base-resolution methylomes of Populus trichocarpa reveal the

Jo

association between DNA methylation and drought stress, BMC Genet. 15 (2014) Suppl 1:S9. doi: 10.1186/1471-2156-15-S1-S9.

86. I. Makarevitch, A.J. Waters, P.T. West, M. Stitzer, C.N. Hirsch, J. Ross-Ibarra, N.M, Springer, Transposable elements contribute to activation of maize genes in response to abiotic stress, PLoS Genet. 11 (2015) e1004915. doi: 10.1371/journal.pgen.1004915. 87. K. van Dijk, Y. Ding, S. Malkaram, J.J. Riethoven, R. Liu, J. Yang, P. Laczko, H. Chen, Y. Xia, I. Ladunga, Z. Avramova, M. Fromm, Dynamic changes in genome-wide histone H3 37

lysine 4 methylation patterns in response to dehydration stress in Arabidopsis thaliana, BMC Plant Biol. 10 (2010) 238. doi: 10.1186/1471-2229-10-238. 88. D. Gourcilleau, M.B. Bogeat-Triboulot, D. Le Thiec, C. Lafon-Placette, A. Delaunay, W. Abu El-Soud, F. Brignolas, S. Maury, DNA methylation and histone acetylation: Genotypic variations in hybrid poplars, impact of water deficit and relationships with productivity, Ann. For. Sci. 67 (2010) 208. DOI: 10.1051/forest/2009101. 89. W. Wang, X. Zhao, Y. Pan, L. Zhu, B. Fu, Z. Li, DNA methylation changes detected by

ro of

methylation-sensitive amplified polymorphism in two contrasting rice genotypes under salt stress, J. Genet. Genomics. 38 (2011) 419-24. doi: 10.1016/j.jgg.2011.07.006.

90. K. Bräutigam, K.J. Vining, C. Lafon-Placette, C.G. Fossdal, M. Mirouze, J.G Marcos, S.

-p

Fluch, M.F. Fraga, M.A. Guevara, D. Abarca, O. Johnsen, S. Maury, S.H. Strauss, M.M, Campbell, A. Rohde, C. Díaz-Sala, Epigenetic regulation of adaptive responses of forest

re

tree species to the environment, Ecol. Evol. 3 (2013) 399-415. doi: 10.1002/ece3.461.

lP

91. Y. Fei, Y. Xue, P. Du, S. Yang, X. Deng, Expression analysis and promoter methylation under osmotic and salinity stress of TaGAPC1 in wheat (Triticum aestivum L), Protoplasma. 254 (2017) 987-996. doi: 10.1007/s00709-016-1008-5.

na

92. S. Kumar, A.S. Beena, M. Awana, A. Singh.

Salt-Induced Tissue-Specific Cytosine

Methylation Downregulates Expression of HKT Genes in Contrasting Wheat (Triticum

ur

aestivum L.) Genotypes, DNA and Cell Biology (2017) 36 (4): 283-294.

Jo

93. C. Liu, H. Li, J. Lin, Y. Wang, X. Xu, Z.M. Cheng, Y. Chang, Genome-Wide Characterization of DNA Demethylase Genes and Their Association with Salt Response in Pyrus, Genes (Basel). 9 (2018) pii: E398. doi: 10.3390/genes9080398.

94. M.W. Yaish, A. Al-Lawati, I. Al-Harrasi, H.V. Patankar, Genome-wide DNA Methylation analysis in response to salinity in the model plant caliph medic (Medicago truncatula), BMC Genomics. 19 (2018) 78. doi: 10.1186/s12864-018-4484-5.

38

95. H. Li, S. Yan, L. Zhao, J. Tan, Q. Zhang, F. Gao, P. Wang, H. Hou, L. Li, Histone acetylation associated up-regulation of the cell wall related genes is involved in salt stress induced maize root swelling, BMC Plant Biol. 14 (2014) 105. doi: 10.1186/1471-2229-14105. 96. R. Garg, V. Narayana Chevala, R. Shankar, M. Jain, Divergent DNA methylation patterns associated with gene expression in rice cultivars with contrasting drought and salinity stress response, Sci. Rep. 9 (2015) 14922. doi: 10.1038/srep14922.

ro of

97. A. Banerjee, A. Roychoudhury , Epigenetic regulation during salinity and drought stress in plants: histone modifications and DNA methylation, Plant gene 11 (2017) 199-204. DOI: 10.1016/j.plgene.2017.05.011.

-p

98. I. Ghosh, A. Sadhu, Y. Moriyasu, M. Bandyopadhyay, A. Mukherjee, Manganese oxide nanoparticles induce genotoxicity and DNA hypomethylation in the moss Physcomitrella Mutat.

Res.

Genet.

Environ.

Mutagen.

DOI:

lP

10.1016/j.mrgentox.2018.12.006.

Toxicol.

re

patens,

99. A. Ariani, S. Romeo, A.T. Groover, L. Sebastiani, Comparative epigenomic and transcriptomic analysis of Populus roots under excess Zn, Environ. Sci. Pollut. R. 132

100.

na

(2016) 16-27. https://doi.org/10.1016/j.envexpbot.2016.08.005. A. Cicatelli, V. Todeschini, G. Lingua, S. Biondi, P. Torrigiani, S. Castiglione,

ur

Epigenetic control of heavy metal stress response in mycorrhizal versus non-mycorrhizal

Jo

poplar plants, Environ. Sci. Pollut. Res. Int. 21 (2014) 1723-1737. doi: 10.1007/s11356-0132072-4.

101.

X. Ou, Y. Zhang, C. Xu, X. Lin, Q. Zang, T. Zhuang, L. Jiang, D. von Wettstein, B.

Liu, Transgenerational inheritance of modified DNA methylation patterns and enhanced tolerance induced by heavy metal stress in rice (Oryza sativa L.), PLoS One. 7 (2012) e41143. doi: 10.1371/journal.pone.0041143.

39

102.

F.A. Erturk, G. Agar, E. Arslan, G. Nardemir, Analysis of genetic and epigenetic

effects of maize seeds in response to heavy metal (Zn) stress, Environ. Sci. Pollut. Res. Int. 22 (2015) 10291-7. doi: 10.1007/s11356-014-3886-4. 103.

M.S. Taspinar, M. Aydin, B. Sigmaz, S. Yagci, E. Arslan, G. Agar, Aluminum-

Induced Changes on DNA Damage, DNA Methylation and LTR Retrotransposon Polymorphism in Maize, Arab. J. Sci. Eng. 43 (2018) 123-131. 104.

D.H. Kim, S. Sung, Coordination of the vernalization response through a VIN3 and

ro of

FLC gene family regulatory network in Arabidopsis, Plant Cell. 25 (2013) 454-69. doi: 10.1105/tpc.112.104760. 105.

J.B. Heo, S. Sung, Vernalization-mediated epigenetic silencing by a long intronic

106.

-p

noncoding RNA, Science. 331 (2011) 76-9. doi: 10.1126/science.1197349.

N. Steward, T. Kusano, H. Sano, Expression of ZmMET1, a gene encoding a DNA

re

methyltransferase from maize, is associated not only with DNA replication in actively

lP

proliferating cells, but also with altered DNA methylation status in cold-stressed quiescent cells, Nucleic Acids Res. 28(2000) 3250-9. 107.

S.N. Hashida, T. Uchiyama, C. Martin, Y. Kishima, Y. Sano, T. Mikami, The

na

temperature-dependent change in methylation of the Antirrhinum transposon Tam3 is controlled by the activity of its transposase, Plant Cell. 18 (2006) 104-18. T. Liu, Y. Li, W. Duan, F. Huang, X.

ur

108.

Hou, Cold acclimation alters DNA

Jo

methylation patterns and confers tolerance to heat and increases growth rate in Brassica rapa, J. Exp. Bot. 68 (2017) 1213-1224. doi: 10.1093/jxb/erw496.

109.

C. Hébrard, D.G. Peterson, G. Willems, A. Delaunay, B. Jesson, M. Lefèbvre, S.

Barnes, S. Maury, Epigenomics and bolting tolerance in sugar beet genotypes, J. Exp. Bot. 67 (2016) 207-25. doi: 10.1093/jxb/erv449. 110.

C. Alonso, D. Ramos-Cruz, C. Becker, The role of plant epigenetics in biotic

interactions, New Phytol. 221 (2019) 731-737. doi: 10.1111/nph.15408. 40

111.

M.S. Alves, S.P. Dadalto, A.B. Gonçalves, G.B. De Souza,V.A. Barros, L.G. Fietto,

Plant bZIP transcription factors responsive to pathogens: a review, Int. J. Mol. Sci. 14 (2013) 7815-28. doi: 10.3390/ijms14047815. 112.

H. Jin, S.M. Choi, M.J. Kang, S.H. Yun, D.J. Kwon, Y.S. Noh, B. Noh, Salicylic

acid-induced

transcriptional

reprogramming

by

the

HAC-NPR1-TGA

histone

acetyltransferase complex in Arabidopsis, Nucleic Acids Res. 46 (2018) 11712-11725. doi: 10.1093/nar/gky847. R.H. Dowen, M. Pelizzola, R.J. Schmitz, R. Lister, J.M. Dowen, J.R. Nery, J.E.

ro of

113.

Dixon, J.R. Ecker, Widespread dynamic DNA methylation in response to biotic stress, Proc. Natl. Acad. Sci. U S A. 109 (2012) E2183-91. doi: 10.1073/pnas.1209329109.

B. Ding, G.L. Wang, Chromatin versus pathogens: the function of epigenetics in

-p

114.

plant immunity, Front. Plant Sci. 6 (2015) 675. doi: 10.3389/fpls.2015.00675. B. Zou, D.L. Yang, Z. Shi, H. Dong, J. Hua, Monoubiquitination of histone 2B at the

re

115.

lP

disease resistance gene locus regulates its expression and impacts immune responses in Arabidopsis, Plant Physiol. 165 (2014) 309-18. doi: 10.1104/pp.113.227801. 116.

Y. Hou, L. Wang, L. Wang, L. Liu, L. Li, L. Sun, Q. Rao, J. Zhang, S. Huang,

na

JMJ704 positively regulates rice defense response against Xanthomonas oryzae pv. oryzae infection via reducing H3K4me2/3 associated with negative disease resistance regulators,

J.L. Soyer, T. Rouxel, I. Fudal, Chromatin-based control of effector gene expression

Jo

117.

ur

BMC Plant Biol. 9 (2015) 15:286. doi: 10.1186/s12870-015-0674-3.

in

plant-associated

fungi,

Curr.

Opin.

Plant

Biol.

26

(2015)

51-6.

doi:

10.1016/j.pbi.2015.05.025.

118.

P. Mohammadi, B. Bahramnejad, H. Badakhshan, H. Kanouni, DNA methylation

changes in fusarium wilt resistant and sensitive chickpea genotypes (Cicer arietinum L.), Physiol. Mol. Plant Pathol. 91 (2015) 72-80. https://doi.org/10.1016/j.pmpp.2015.06.001.

41

119.

L. Furci, R. Jain, J. Stassen, O. Berkowitz, J. Whelan, D. Roquis, V. Baillet, V.

Colot, F. Johannes, J. Ton, Identification and characterisation of hypomethylated DNA loci controlling

quantitative

resistance

in

Arabidopsis,

(2019)

pii:

e40655.

doi:

10.7554/eLife.40655. 120.

B. Wang, X. Yang, Y. Wang, Y. Xie, X. Zhou, Tomato Yellow Leaf Curl Virus V2

Interacts with Host Histone Deacetylase 6 To Suppress Methylation-Mediated Transcriptional Gene Silencing in Plants, J. Virol.

92 (2018)

pii: e00036-18. doi:

ro of

10.1128/JVI.00036-18. 121.

J. Zervudacki, A.Yu, D. Amesefe, J. Wang, J. Drouaud, L. Navarro, A. Deleris,

Transcriptional control and exploitation of an immune-responsive family of plant

-p

retrotransposons, Embo. J. 37 (2018). pii: e98482. doi: 10.15252/embj.201798482. 122.

N. Pumplin, A. Sarazin, P.E. Jullien, N.G. Bologna, S. Oberlin, O. Voinnet, DNA

re

Methylation Influences the Expression of DICER-LIKE4 Isoforms, Which Encode Proteins

lP

of Alternative Localization and Function, Plant Cell. 28 (2016) 2786-2804. 123.

A. Fonsêca, M.M. Richard,V. Geffroy, A. Pedrosa-Harand, Epigenetic analyses and

the distribution of repetitive DNA and resistance genes reveal the complexity of common

na

bean (Phaseolus vulgaris L., Fabaceae) heterochromatin, Cytogenet. Genome Res. 143 (2014) 168-78. doi: 10.1159/000360572. A. Berr, E.J. McCallum, A. Alioua, D. Heintz, T. Heitz, W.H. Shen, Arabidopsis

ur

124.

methyltransferase

SET DOMAIN

GROUP8

mediates

induction of the

Jo

histone

jasmonate/ethylene pathway genes in plant defense response to necrotrophic fungi, Plant Physiol. 154 (2010) 1403-14. doi: 10.1104/pp.110.161497.

125.

L.J. Ferreira, M.T.A. Donoghue, F. Borges F, P. Barros, N.J. Saibo, A.P. Santos,

M.M. Oliveira, Differentially Methylated Regions (DMRs): uncovering new points for regulation of salt stress responsive genes in rice. Epigenomes, 3 (2019) 4. doi:10.3390/epigenomes3010004. 42

126.

R.T. Grant-Downton, H.G. Dickinson, Epigenetics and its implications for plant

biology. 1. The epigenetic network in plants, Ann. Bot. 96 (2005) 1143-64. 127.

J. Reinders, J. Paszkowski, Unlocking the Arabidopsis epigenome, Epigenetics 4

(2009) 557-63. 128.

A. Giovino, P. Marino, G. Domina, A. Scialabba, R. Schicchi, G. Diliberto, C. Rizza,

S. Scibetta, Evaluation of the DNA barcoding approach to develop a reference data-set for

Jo

ur

na

lP

re

-p

ro of

the threatened flora of Sicily, Plant Biosyst. 150 (2016) 631-640.

43

Figure legends Figure 1. Main DNA methylation mechanisms at cytosine position in CG, CHH, CHG sequences

lP

re

-p

ro of

and their role in plants.

na

Figure 2. Integration of different “omic” approaches to develop a systems-based approaches to

Jo

ur

investigate plant responses to environmental stresses.

44

ro of -p re

Jo

ur

na

environmental stresses.

lP

Figure 3. Main epigenetic modifications linked to stress memory in relation to different

45

ro of -p re lP

Jo

ur

na

Figure 4. Epigenetic mechanisms linked to priming in plants.

46

-p

ro of

Figure 5. Main epigenetic mechanisms involved in abiotic stress responses in plants.

Jo

ur

na

lP

re

Figure 6. Main epigenetic mechanisms involved in plant biotic stress responses.

47

re

-p

ro of

Figure 7

Jo

ur

na

lP

Figure 8

48

Tables

Jo

ur

na

lP

re

-p

ro of

Table 1. A selected list of references dealing with epigenetic/epigenomic studies in plants, the type of mechanism, major findings or article’s subject, and technique were indicated.

49

Epigenetic Mechanisms

Main Findings/ Subjects

Used/Studied Techniques

Histone modification

ChIP-CHIP

Chen et al., 2018

DNA methylation

Zhong et al., 2013

DNA methylation

Zilberman and Henikoff, 2007

DNA methylation me-wide methylation analysis

The technique can be used to identify target sequences of transcription factors, and positions of histones with modified N-terminal tails Description of the ChIP-Seq workflow and its various applications in plants Plant epigenome is dynamic during fruit development and play a key role in processes such as ripening Combining different techniques with DNA microarrays and highthroughput sequencing is possible to map DNA methylation at genomewide scale Wholegenome arrays can be used for epigenetic analysis such as chromatinimmunoprecipitationchip studies, methylome characterization Opportunities and limitations of RRBS iwere described for nonmodel plant species

-p

re

lP DNA methylation

ur

na

Mockler et al., 2005

Paun et al., 2019

DNA methylation

Zhang et al., 2015

DNA methylation

Kunze, 2018

DNA methylation

Jo

ro of

References dealing with Epigenetic/Epigenomic Technique Pellegrini and Ferrari, 2012

ChIP-Seq

Whole Genome Bisulphite sequencing

Different techniques

Padlock-probe combined with promoter tiling arrays

Reduced representation bisulphite sequencing (RRBS) Differentially Methylationmethylated genomic sensitive singlefragments were linked strand with sexual dimorphism conformation of pests in rice analysis The EpiTYPER™ assay EpiTYPER™ allows a fast and assay reproducible targeted quantification of 50 individual CpGs in a high throughput manner

DNA methylation

Agorio et al., 2017

DNA methylation

Rodriguez-Lopez et al., 2010

DNA methylation

Gonzalgo and Liang, 2007

DNA methylation

Van den Boom and Ehrich, 2009

Methylation levels in WRKY genes were negatively correlated with gene expression in response to NaCl treatments Data revealed a main role of this loop in maintaining a natural epiallele

Bisulphite coupled with PCR

High-resolution melting analysis allows adjusting thermal parameters for DNA hybridization and PCRbased techniques analyzing the impact of DNA methylation marks on the thermostability of regulatory sequences Methylation-sensitive single-nucleotide primer extension can be used to quantify DNA methylation

High-resolution melting analysis

Bisulphite treatment coupled with pyrosequencing

Methylationsensitive singlenucleotide primer extension ((MsSNuPE) Matrix-assisted laser desorption ionization/timeof-flight

Martinelli et al., 2018

Non-coding RNA

Next-generation sequencing

Matsui et al. 2010

Non-coding RNA

re lP

-p

ro of

Celik et al., 2019

Jo

ur

na

DNA methylation

Matrix-assisted laser desorption ionization/time-of-flight combines candidate gene amplification with base-specific cleavage or primer extension methods and MALDITOF mass spectrometric analysis to overcome high dynamicity of DNA methylation patterns MiRNAs may be included in nextgeneration breeding programs Tiling array can identify the stress-responsive genes in Arabidopsis thaliana

Tiling DNA microarray

51

Table 2. A selected list of main epigenetic modifications in response to different environmental stresses. Levels of analysis (target or genome-wide analysis) and references were indicated.

DNA methylation Chromatin restructuring DNA methylation Histone modification

Level of Environmental Analysis Stresses (Genome-wide or target epigenetic specific) Abiotic Stress Responses Genome-wide Osmotic Genome-wide Cold

References

Sani et al., 2013 Kim et al., 2015

ro of

Epigenetic modifications

Genome-wide Target epigenetic specific Target epigenetic specific Target epigenetic specific

Cold Drought, cold,

CHG methylation

Genome-wide

Phosphate starvation

DNA methylation

Genome-wide

Histone modification

Genome-wide

Drought

H3K4me3

Genome-wide

Drought

Target epigenetic specific

Salinity

re

-p

Low humidity Heat

Drought

Target epigenetic specific

Salinity

DNA methylation

Genome-wide

Salinity

DNA methylation

Genome-wide

Salinity

H3K27me3 H3K27me3 Differentially methylated regions Hypermethylated regions; heterochromatin

Genome-wide Genome-wide Genome-wide

Salinity Drought Salinity

Genome-wide

Salinity

Jo

DNA methylation

General abiotic stresses

lP

ur

DNA methylation (WRKY, GAPC1)

na

Non coding RNA (RNA helicases) CHG methylation

Liu et al., 2017 Rutowicz et al., 2015 Khan et al., 2014 Tricker et al., 2012; Popova et al., 2013; Yong-Villalobos et al., 2015 Liang et al.,
 2014; Bräutigam et al., 2013 Gourcilleau et al., 2010; van Dijk et al., 2010; Çelik et al., 2019; Fei et al., 2017; Liu et al., 2018; Ferreira et al., 2019 Yaish et al., 2018 Al-Lawati et al., 2016 Sani et al., 2013 Garg et al., 2015 Banerjee and Roychoudhury, 2017 52

silencing DNA methylation

Target epigenetic specific

Aluminum stress

H3K27me3

Genome-wide

Zinc stress

CHG hypomethylation

Genome-wide

Heavy metals

Chromatin-remodeling processes; Non-coding RNAs MET1 transcriptional silencing Tam3 transposon methylation DNA methylation

Genome-wide

Cold (vernalization)

Agnieszka, 2018; Ghosh et al., 2018 Ariani et al., 2016; Erturk et al., 2015 Cicatelli et al., 2014 Kim and Sung 2013 Steward et al., 2000 Hashida et al., 2006 Liu et al., 2017

Cytosine methylation

Target epigenetic specific Genome-wide

Responses to pathogen infections Responses to Fusarium wilt

DNA hypomethylation

Genome-wide

General biotic stress

Cytosine DNA methylation DNA methylation

Genome-wide

Responses to geminivirus infection Responses to pathogen infections General biotic stress

Jaskiewicz et al., 2011 Mohammadi, et al., 2015 Furci et al., 2019 Wang et al., 2018 Zervudacki et al., 2018 Pumplin et al., 2016 Fonsêca et al.,

DNA methylation

-p

DNA methylation

ro of

Target epigenetic Cold specific Target epigenetic Cold specific Target epigenetic Cold specific Biotic Stress Responses Target epigenetic General biotic stress specific Target epigenetic Stress memory specific Target epigenetic Stress memory specific Target epigenetic General biotic stress specific Target epigenetic Responses to necrotrophic specific and biotrophic pathogens Target epigenetic Response to pathogen Pst specific DC3000 Target epigenetic Response to pathogen Pst specific DC3000 Target epigenetic General biotic stress specific Target epigenetic Responses to fungal attacks specific

Chromatin remodeling

Jo

ur

na

Activation of HUB1 and HUB2 Histone deacetylase (HDA19) H3K27me2/3; Jumonji C Chromatin-based transcriptional regulation Histone modifications

lP

re

H3 histone modifications DNA methylation

Target epigenetic specific DNA methylation and Target epigenetic small-interfering RNAs specific Hypomethylation of Genome-wide

Responses to fungus C.

Alonso et al., 2019 Alves et al., 2013 Jin et al., 2018 Dowen et al., 2012 Ding and Wang, 2015 Zou et al., 2014 Choi et al., 2012 Hou et al., 2015 Soyer et al., 2014

53

Target epigenetic specific

lindemuthianum Responses to fungal pathogens

2014 Berr et al., 2010

Jo

ur

na

lP

re

-p

ro of

centromeric regions Methylation of histone H3 lysine 36

54