CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize

A transposable element in a NAC gene is associated with drought tolerance in maize seedlings

Drought represents a major constraint on maize production worldwide. Understanding the genetic basis for natural variation in drought tolerance of maize may facilitate efforts to improve this trait in cultivated germplasm. Here, using a genome-wide association study, we show that a miniature inverted-repeat transposable element (MITE) inserted in the promoter of a NAC gene (ZmNAC111) is significantly associated with natural variation in maize drought tolerance. The 82-bp MITE represses ZmNAC111 expression via RNA-directed DNA methylation and H3K9 dimethylation when heterologously expressed in Arabidopsis. Increasing ZmNAC111 expression in transgenic maize enhances drought tolerance at the seedling stage, improves water-use efficiency and induces upregulation of drought-responsive genes under water stress. The MITE insertion in the ZmNAC111 promoter appears to have occurred after maize domestication and spread among temperate germplasm. The identification of this MITE insertion provides insight into the genetic basis for natural variation in maize drought tolerance.

Drought is a major cause of yield loss in maize and understanding the genetic determinants of natural variation in drought tolerance may aid breeding programs produce more tolerant varieties. Here, Mao et al. identify a MITE transposon insertion in a NAC transcription factor, which is associated with natural variation in drought tolerance.

The Methylome of Soybean Roots during the Compatible Interaction with the Soybean Cyst Nematode

The soybean cyst nematode (SCN; Heterodera glycines) induces the formation of a multinucleated feeding site, or syncytium, whose etiology includes massive gene expression changes. Nevertheless, the genetic networks underlying gene expression control in the syncytium are poorly understood. DNA methylation is a critical epigenetic mark that plays a key role in regulating gene expression. To determine the extent to which DNA methylation is altered in soybean (Glycine max) roots during the susceptible interaction with SCN, we generated whole-genome cytosine methylation maps at single-nucleotide resolution. The methylome analysis revealed that SCN induces hypomethylation to a much higher extent than hypermethylation. We identified 2,465 differentially hypermethylated regions and 4,692 hypomethylated regions in the infected roots compared with the noninfected control. In addition, 703 and 1,346 unique genes were identified as overlapping with hyper- or hypomethylated regions, respectively. The differential methylation in genes apparently occurs independently of gene size and GC content but exhibits strong preference for recently duplicated paralogs. Furthermore, a set of 278 genes was identified as specifically syncytium differentially methylated genes. Of these, we found genes associated with epigenetic regulation, phytohormone signaling, cell wall architecture, signal transduction, and ubiquitination. This study provides, to our knowledge, new evidence that differential methylation is part of the regulatory mechanisms controlling gene expression changes in the nematode-induced syncytium.

Examining the Causes and Consequences of Context-Specific Differential DNA Methylation in Maize

DNA methylation is a stable modification of chromatin that can contribute to epigenetic variation through the regulation of genes or transposons. Profiling of DNA methylation in five maize (Zea mays) inbred lines found that while DNA methylation levels for more than 99% of the analyzed genomic regions are similar, there are still 5,000 to 20,000 context-specific differentially methylated regions (DMRs) between any two genotypes. The analysis of identical-by-state genomic regions that have limited genetic variation provided evidence that DMRs can occur without local sequence variation, but they are less common than in regions with genetic variation. Characterization of the sequence specificity of DMRs, location of DMRs relative to genes and transposons, and patterns of DNA methylation in regions flanking DMRs reveals a distinct subset of DMRs. Transcriptome profiling of the same tissue revealed that only approximately 20% of genes with qualitative (on-off) differences in gene expression are associated with DMRs, and there is little evidence for association of DMRs with genes that show quantitative differences in gene expression. We also identify a set of genes that may represent cryptic information that is silenced by DNA methylation in the reference B73 genome. Many of these genes exhibit natural variation in other genotypes, suggesting the potential for selection to act upon existing epigenetic natural variation. This study provides insights into the origin and influences of DMRs in a crop species with a complex genome organization.

Trans-Homolog Interactions Facilitating Paramutation in Maize

Paramutations represent locus-specific trans-homolog interactions affecting the heritable silencing properties of endogenous alleles. Although examples of paramutation are well studied in maize (Zea mays), the responsible mechanisms remain unclear. Genetic analyses indicate roles for plant-specific DNA-dependent RNA polymerases that generate small RNAs, and current working models hypothesize that these small RNAs direct heritable changes at sequences often acting as transcriptional enhancers. Several studies have defined specific sequences that mediate paramutation behaviors, and recent results identify a diversity of DNA-dependent RNA polymerase complexes operating in maize. Other reports ascribe broader roles for some of these complexes in normal genome function. This review highlights recent research to understand the molecular mechanisms of paramutation and examines evidence relevant to small RNA-based modes of transgenerational epigenetic inheritance.

A Comparative Epigenomic Analysis of Polyploidy-Derived Genes in Soybean and Common Bean

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a paleopolyploidy (whole-genome duplication [WGD]) event, approximately 56.5 million years ago, followed by a genus Glycine-specific polyploidy, approximately 10 million years ago. Cytosine methylation is an epigenetic mark that plays an important role in the regulation of genes and transposable elements (TEs); however, the role of DNA methylation in the fate/evolution of genes following polyploidy and speciation has not been fully explored. Whole-genome bisulfite sequencing was used to produce nucleotide resolution methylomes for soybean and common bean. We found that, in soybean, CG body-methylated genes were abundant in WGD genes, which were, on average, more highly expressed than single-copy genes and had slower evolutionary rates than unmethylated genes, suggesting that WGD genes evolve more slowly than single-copy genes. CG body-methylated genes were also enriched in shared single-copy genes (single copy in both species) that may be responsible for the broad and high expression patterns of this class of genes. In addition, diverged methylation patterns in non-CG contexts between paralogs were due mostly to TEs in or near genes, suggesting a role for TEs and non-CG methylation in regulating gene expression post polyploidy. Reference methylomes for both soybean and common bean were constructed, providing resources for investigating epigenetic variation in legume crops. Also, the analysis of methylation patterns of duplicated and single-copy genes has provided insights into the functional consequences of polyploidy and epigenetic regulation in plant genomes.

DNA methylation and gene expression in Mimulus guttatus

Background

The presence of methyl groups on cytosine nucleotides across an organism’s genome (methylation) is a major regulator of genome stability, crossing over, and gene regulation. The capacity for DNA methylation to be altered by environmental conditions, and potentially passed between generations, makes it a prime candidate for transgenerational epigenetic inheritance. Here we conduct the first analysis of the Mimulus guttatus methylome, with a focus on the relationship between DNA methylation and gene expression.

Results

We present a whole genome methylome for the inbred line Iron Mountain 62 (IM62). DNA methylation varies across chromosomes, genomic regions, and genes. We develop a model that predicts gene expression based on DNA methylation (R² = 0.2). Post hoc analysis of this model confirms prior relationships, and identifies novel relationships between methylation and gene expression. Additionally, we find that DNA methylation is significantly depleted near gene transcriptional start sites, which may explain the recently discovered elevated rate of recombination in these same regions.

Conclusions

The establishment here of a reference methylome will be a useful resource for the continued advancement of M. guttatus as a model system. Using a model-based approach, we demonstrate that methylation patterns are an important predictor of variation in gene expression. This model provides a novel approach for differential methylation analysis that generates distinct and testable hypotheses regarding gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1668-0) contains supplementary material, which is available to authorized users.

Transposable elements, mRNA expression level and strand-specificity of small RNAs are associated with non-additive inheritance of gene expression in hybrid plants

BACKGROUND

Gene expression inheritance patterns in Arabidopsis hybrid plants were investigated for correlation with the presence of transposable elements (TEs) and small RNA profile.

RESULTS

The presence of TEs in a gene and the expression of small RNA matching a gene were both found to be associated with non-additive mRNA inheritance patterns in hybrids. Expression levels below mid-parent values in the hybrids were associated with low mRNA expression in parents, with the presence of small RNA from both strands, and with the presence of TEs. High-parent dominance of mRNA levels was found to be associated with high parental mRNA expression levels, the absence of TEs, and for some genes, with small RNA fragments that are predominantly from the sense strand. These small RNAs exhibit a broader size distribution than siRNA and reduced nucleotide end bias, which are consistent with an origin from degraded mRNA. Thus, increased as well as decreased gene expression in hybrids relative to the parental mean is associated with gene expression levels, TE presence and small RNA fragments with differing characteristics.

CONCLUSIONS

The data presented here is consistent with a role for differential mRNA decay kinetics as one mechanism contributing to high-parent dominance in gene expression. Our evidence is also consistent with trans repression by siRNA and TEs as the cause of low-parent dominance.

Crop Epigenomics: Identifying, Unlocking, and Harnessing Cryptic Variation in Crop Genomes

DNA methylation is a key chromatin modification in plant genomes that is meiotically and mitotically heritable, and at times is associated with gene expression and morphological variation. Benefiting from the increased availability of high-quality reference genome assemblies and methods to profile single-base resolution DNA methylation states, DNA methylomes for many crop species are available. These efforts are making it possible to begin answering crucial questions, including understanding the role of DNA methylation in developmental processes, its role in crop species evolution, and whether DNA methylation is dynamically altered and heritable in response to changes in the environment. These genome-wide maps provide evidence for the existence of silent epialleles in plant genomes which, once identified, can be targeted for reactivation leading to phenotypic variation.

Controlling for conservation in genome-wide DNA methylation studies

BACKGROUND

A commonplace analysis in high-throughput DNA methylation studies is the comparison of methylation extent between different functional regions, computed by averaging methylation states within region types and then comparing averages between regions. For example, it has been reported that methylation is more prevalent in coding regions as compared to their neighboring introns or UTRs, leading to hypotheses about novel forms of epigenetic regulation.

RESULTS

We have identified and characterized a bias present in these seemingly straightforward comparisons that results in the false detection of differences in methylation intensities across region types. This bias arises due to differences in conservation rates, rather than methylation rates, and is broadly present in the published literature. When controlling for conservation at coding start sites the differences in DNA methylation rates disappear. Moreover, a re-evaluation of methylation rates at intronexon junctions reveals that the magnitude of previously reported differences is greatly exaggerated. We introduce two correction methods to address this bias, an inferencebased matrix completion algorithm and an averaging approach, tailored to address different underlying biological questions. We evaluate how analysis using these corrections affects the detection of differences in DNA methylation across functional boundaries.

CONCLUSIONS

We report here on a bias in DNA methylation comparative studies that originates in conservation rate differences and manifests itself in the false discovery of differences in DNA methylation intensities and their extents. We have characterized this bias and its broad implications, and show how to control for it so as to enable the study of a variety of biological questions.

MSAP markers and global cytosine methylation in plants: a literature survey and comparative analysis for a wild-growing species

Methylation of DNA cytosines affects whether transposons are silenced and genes are expressed, and is a major epigenetic mechanism whereby plants respond to environmental change. Analyses of methylation-sensitive amplification polymorphism (MS-AFLP or MSAP) have been often used to assess methyl-cytosine changes in response to stress treatments and, more recently, in ecological studies of wild plant populations. MSAP technique does not require a sequenced reference genome and provides many anonymous loci randomly distributed over the genome for which the methylation status can be ascertained. Scoring of MSAP data, however, is not straightforward, and efforts are still required to standardize this step to make use of the potential to distinguish between methylation at different nucleotide contexts. Furthermore, it is not known how accurately MSAP infers genome-wide cytosine methylation levels in plants. Here, we analyse the relationship between MSAP results and the percentage of global cytosine methylation in genomic DNA obtained by HPLC analysis. A screening of literature revealed that methylation of cytosines at cleavage sites assayed by MSAP was greater than genome-wide estimates obtained by HPLC, and percentages of methylation at different nucleotide contexts varied within and across species. Concurrent HPLC and MSAP analyses of DNA from 200 individuals of the perennial herb Helleborus foetidus confirmed that methyl-cytosine was more frequent in CCGG contexts than in the genome as a whole. In this species, global methylation was unrelated to methylation at the inner CG site. We suggest that global HPLC and context-specific MSAP methylation estimates provide complementary information whose combination can improve our current understanding of methylation-based epigenetic processes in nonmodel plants.

FEATnotator: A tool for integrated annotation of sequence features and variation, facilitating interpretation in genomics experiments

As approaches are sought for more efficient and democratized uses of non-model and expanded model genomics references, ease of integration of genomic feature datasets is especially desirable in multidisciplinary research communities. Valuable conclusions are often missed or slowed when researchers refer experimental results to a single reference sequence that lacks integrated pan-genomic and multi-experiment data in accessible formats. Association of genomic positional information, such as results from an expansive variety of next-generation sequencing experiments, with annotated reference features such as genes or predicted protein binding sites, provides the context essential for conclusions and ongoing research. When the experimental system includes polymorphic genomic inputs, rapid calculation of gene structural and protein translational effects of sequence variation from the reference can be invaluable. Here we present FEATnotator, a lightweight, fast and easy to use open source software program that integrates and reports overlap and proximity in genomic information from any user-defined datasets including those from next generation sequencing applications. We illustrate use of the tool by summarizing whole genome sequence variation of a widely used natural isolate of Arabidopsis thaliana in the context of gene models of the reference accession. Previous discovery of a protein coding deletion influencing root development is replicated rapidly. Appropriate even in investigations of a single gene or genic regions such as QTL, comprehensive reports provided by FEATnotator better prepare researchers for interpretation of their experimental results. The tool is available for download at http://featnotator.sourceforge.net.

Recombination in diverse maize is stable, predictable, and associated with genetic load

Among the fundamental evolutionary forces, recombination arguably has the largest impact on the practical work of plant breeders. Varying over 1,000-fold across the maize genome, the local meiotic recombination rate limits the resolving power of quantitative trait mapping and the precision of favorable allele introgression. The consequences of low recombination also theoretically extend to the species-wide scale by decreasing the power of selection relative to genetic drift, and thereby hindering the purging of deleterious mutations. In this study, we used genotyping-by-sequencing (GBS) to identify 136,000 recombination breakpoints at high resolution within US and Chinese maize nested association mapping populations. We find that the pattern of cross-overs is highly predictable on the broad scale, following the distribution of gene density and CpG methylation. Several large inversions also suppress recombination in distinct regions of several families. We also identify recombination hotspots ranging in size from 1 kb to 30 kb. We find these hotspots to be historically stable and, compared with similar regions with low recombination, to have strongly differentiated patterns of DNA methylation and GC content. We also provide evidence for the historical action of GC-biased gene conversion in recombination hotspots. Finally, using genomic evolutionary rate profiling (GERP) to identify putative deleterious polymorphisms, we find evidence for reduced genetic load in hotspot regions, a phenomenon that may have considerable practical importance for breeding programs worldwide.

MethylC-seq library preparation for base-resolution whole-genome bisulfite sequencing

Current high-throughput DNA sequencing technologies enable acquisition of billions of data points through which myriad biological processes can be interrogated, including genetic variation, chromatin structure, gene expression patterns, small RNAs and protein-DNA interactions. Here we describe the MethylC-sequencing (MethylC-seq) library preparation method, a 2-d protocol that enables the genome-wide identification of cytosine DNA methylation states at single-base resolution. The technique involves fragmentation of genomic DNA followed by adapter ligation, bisulfite conversion and limited amplification using adapter-specific PCR primers in preparation for sequencing. To date, this protocol has been successfully applied to genomic DNA isolated from primary cell culture, sorted cells and fresh tissue from over a thousand plant and animal samples.

Nascent transcription affected by RNA polymerase IV in Zea mays

All eukaryotes use three DNA-dependent RNA polymerases (RNAPs) to create cellular RNAs from DNA templates. Plants have additional RNAPs related to Pol II, but their evolutionary role(s) remain largely unknown. Zea mays (maize) RNA polymerase D1 (RPD1), the largest subunit of RNA polymerase IV (Pol IV), is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs), and transcriptional regulation of specific alleles. Here, we define the nascent transcriptomes of rpd1 mutant and wild-type (WT) seedlings using global run-on sequencing (GRO-seq) to identify the broader targets of RPD1-based regulation. Comparisons of WT and rpd1 mutant GRO-seq profiles indicate that Pol IV globally affects transcription at both transcriptional start sites and immediately downstream of polyadenylation addition sites. We found no evidence of divergent transcription from gene promoters as seen in mammalian GRO-seq profiles. Statistical comparisons identify genes and TEs whose transcription is affected by RPD1. Most examples of significant increases in genic antisense transcription appear to be initiated by 3'-proximal long terminal repeat retrotransposons. These results indicate that maize Pol IV specifies Pol II-based transcriptional regulation for specific regions of the maize genome including genes having developmental significance.

Whole-genome DNA methylation patterns and complex associations with gene structure and expression during flower development in Arabidopsis

Flower development is a complex process requiring proper spatiotemporal expression of numerous genes. Accumulating evidence indicates that epigenetic mechanisms, including DNA methylation, play essential roles in modulating gene expression. However, few studies have examined the relationship between DNA methylation and floral gene expression on a genomic scale. Here we present detailed analyses of DNA methylomes at single-base resolution for three Arabidopsis floral periods: meristems, early flowers and late flowers. We detected 1.5 million methylcytosines, and estimated the methylation levels for 24 035 genes. We found that many cytosine sites were methylated de novo from the meristem to the early flower stage, and many sites were demethylated from early to late flowers. A comparison of the transcriptome data of the same three periods revealed that the methylation and demethylation processes were correlated with expression changes of >3000 genes, many of which are important for normal flower development. We also found different methylation patterns for three sequence contexts ((m) CG, (m) CHG and (m) CHH) and in different genic regions, potentially with different roles in gene expression.

Genetic perturbation of the maize methylome

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Accessible DNA and relative depletion of H3K9me2 at maize loci undergoing RNA-directed DNA methylation

RNA-directed DNA methylation (RdDM) in plants is a well-characterized example of RNA interference-related transcriptional gene silencing. To determine the relationships between RdDM and heterochromatin in the repeat-rich maize (Zea mays) genome, we performed whole-genome analyses of several heterochromatic features: dimethylation of lysine 9 and lysine 27 (H3K9me2 and H3K27me2), chromatin accessibility, DNA methylation, and small RNAs; we also analyzed two mutants that affect these processes, mediator of paramutation1 and zea methyltransferase2. The data revealed that the majority of the genome exists in a heterochromatic state defined by inaccessible chromatin that is marked by H3K9me2 and H3K27me2 but that lacks RdDM. The minority of the genome marked by RdDM was predominantly near genes, and its overall chromatin structure appeared more similar to euchromatin than to heterochromatin. These and other data indicate that the densely staining chromatin defined as heterochromatin differs fundamentally from RdDM-targeted chromatin. We propose that small interfering RNAs perform a specialized role in repressing transposons in accessible chromatin environments and that the bulk of heterochromatin is incompatible with small RNA production.

Methylated DNA is over-represented in whole-genome bisulfite sequencing data

The development of whole-genome bisulfite sequencing (WGBS) has resulted in a number of exciting discoveries about the role of DNA methylation leading to a plethora of novel testable hypotheses. Methods for constructing sodium bisulfite-converted and amplified libraries have recently advanced to the point that the bottleneck for experiments that use WGBS has shifted to data analysis and interpretation. Here we present empirical evidence for an over-representation of reads from methylated DNA in WGBS. This enrichment for methylated DNA is exacerbated by higher cycles of PCR and is influenced by the type of uracil-insensitive DNA polymerase used for amplifying the sequencing library. Future efforts to computationally correct for this enrichment bias will be essential to increasing the accuracy of determining methylation levels for individual cytosines. It is especially critical for studies that seek to accurately quantify DNA methylation levels in populations that may segregate for allelic DNA methylation states.

Epigenetic regulation and functional exaptation of transposable elements in higher plants

Transposable elements (TEs) are mobile genetic elements that can proliferate in their host genomes. Because of their robust amplification, TEs have long been considered 'selfish DNA', harmful insertions that can threaten host genome integrity. The idea of TEs as junk DNA comes from analysis of epigenetic silencing of their mobility in plants and animals. This idea contrasts with McClintock's characterization of TEs as 'controlling elements'. Emerging studies on the regulatory functions of TEs in plant genomes have updated McClintock's characterization, indicating exaptation of TEs for genetic regulation. In this review, we summarize recent progress in TE silencing, particularly in Arabidopsis and rice, and show that TEs provide an abundant, natural source of regulation for the host genome.

DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis

Background

DNA demethylases regulate DNA methylation levels in eukaryotes. Arabidopsis encodes four DNA demethylases, DEMETER (DME), REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 (DML2), and DML3. While DME is involved in maternal specific gene expression during seed development, the biological function of the remaining DNA demethylases remains unclear.

Results

We show that ROS1, DML2, and DML3 play a role in fungal disease resistance in Arabidopsis. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3), shows increased susceptibility to the fungal pathogen Fusarium oxysporum. We identify 348 genes differentially expressed in rdd relative to wild type, and a significant proportion of these genes are downregulated in rdd and have functions in stress response, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. The rdd-downregulated stress response genes are enriched for short transposable element sequences in their promoters. Many of these transposable elements and their surrounding sequences show localized DNA methylation changes in rdd, and a general reduction in CHH methylation, suggesting that RNA-directed DNA methylation (RdDM), responsible for CHH methylation, may participate in DNA demethylase-mediated regulation of stress response genes. Many of the rdd-downregulated stress response genes are downregulated in the RdDM mutants nrpd1 and nrpe1, and the RdDM mutants nrpe1 and ago4 show enhanced susceptibility to F. oxysporum infection.

Conclusions

Our results suggest that a primary function of DNA demethylases in plants is to regulate the expression of stress response genes by targeting promoter transposable element sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0458-3) contains supplementary material, which is available to authorized users.

pENCODE: a plant encyclopedia of DNA elements

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

Consistent and heritable alterations of DNA methylation are induced by tissue culture in maize

Plants regenerated from tissue culture and their progenies are expected to be identical clones, but often display heritable molecular and phenotypic variation. We characterized DNA methylation patterns in callus, primary regenerants, and regenerant-derived progenies of maize using immunoprecipitation of methylated DNA (meDIP) to assess the genome-wide frequency, pattern, and heritability of DNA methylation changes. Although genome-wide DNA methylation levels remained similar following tissue culture, numerous regions exhibited altered DNA methylation levels. Hypomethylation events were observed more frequently than hypermethylation following tissue culture. Many of the hypomethylation events occur at the same genomic sites across independent regenerants and cell lines. The DNA methylation changes were often heritable in progenies produced from self-pollination of primary regenerants. Methylation changes were enriched in regions upstream of genes and loss of DNA methylation at promoters was associated with altered expression at a subset of loci. Differentially methylated regions (DMRs) found in tissue culture regenerants overlap with the position of naturally occurring DMRs more often than expected by chance with 8% of tissue culture hypomethylated DMRs overlapping with DMRs identified by profiling natural variation, consistent with the hypotheses that genomic stresses similar to those causing somaclonal variation may also occur in nature, and that certain loci are particularly susceptible to epigenetic change in response to these stresses. The consistency of methylation changes across regenerants from independent cultures suggests a mechanistic response to the culture environment as opposed to an overall loss of fidelity in the maintenance of epigenetic states.

Cytosine methylation of an ancient satellite family in the wild beet Beta procumbens

DNA methylation is an essential epigenetic feature for the regulation and maintenance of heterochromatin. Satellite DNA is a repetitive sequence component that often occurs in large arrays in heterochromatin of subtelomeric, intercalary and centromeric regions. Knowledge about the methylation status of satellite DNA is important for understanding the role of repetitive DNA in heterochromatization. In this study, we investigated the cytosine methylation of the ancient satellite family pEV in the wild beet Beta procumbens. The pEV satellite is widespread in species-specific pEV subfamilies in the genus Beta and most likely originated before the radiation of the Betoideae and Chenopodioideae. In B. procumbens, the pEV subfamily occurs abundantly and spans intercalary and centromeric regions. To uncover its cytosine methylation, we performed chromosome-wide immunostaining and bisulfite sequencing of pEV satellite repeats. We found that CG and CHG sites are highly methylated while CHH sites show only low levels of methylation. As a consequence of the low frequency of CG and CHG sites and the preferential occurrence of most cytosines in the CHH motif in pEV monomers, this satellite family displays only low levels of total cytosine methylation.

The CHH motif in sugar beet satellite DNA: a modulator for cytosine methylation

Methylation of DNA is important for the epigenetic silencing of repetitive DNA in plant genomes. Knowledge about the cytosine methylation status of satellite DNAs, a major class of repetitive DNA, is scarce. One reason for this is that arrays of tandemly arranged sequences are usually collapsed in next-generation sequencing assemblies. We applied strategies to overcome this limitation and quantified the level of cytosine methylation and its pattern in three satellite families of sugar beet (Beta vulgaris) which differ in their abundance, chromosomal localization and monomer size. We visualized methylation levels along pachytene chromosomes with respect to small satellite loci at maximum resolution using chromosome-wide fluorescent in situ hybridization complemented with immunostaining and super-resolution microscopy. Only reduced methylation of many satellite arrays was obtained. To investigate methylation at the nucleotide level we performed bisulfite sequencing of 1569 satellite sequences. We found that the level of methylation of cytosine strongly depends on the sequence context: cytosines in the CHH motif show lower methylation (44-52%), while CG and CHG motifs are more strongly methylated. This affects the overall methylation of satellite sequences because CHH occurs frequently while CG and CHG are rare or even absent in the satellite arrays investigated. Evidently, CHH is the major target for modulation of the cytosine methylation level of adjacent monomers within individual arrays and contributes to their epigenetic function. This strongly indicates that asymmetric cytosine methylation plays a role in the epigenetic modification of satellite repeats in plant genomes.

Transposable elements, a treasure trove to decipher epigenetic variation: insights from Arabidopsis and crop epigenomes

In the past decade, plant biologists and breeders have developed a growing interest in the field of epigenetics, which is defined as the study of heritable changes in gene expression that cannot be explained by changes in the DNA sequence. Epigenetic marks can be responsive to the environment, and evolve faster than genetic changes. Therefore, epigenetic diversity may represent an unexplored resource of natural variation that could be used in plant breeding programmes. On the other hand, crop genomes are largely populated with transposable elements (TEs) that are efficiently targeted by epigenetic marks, and part of the epigenetic diversity observed might be explained by TE polymorphisms. Characterizing the degree to which TEs influence epigenetic variation in crops is therefore a major goal to better use epigenetic variation. To date, epigenetic analyses have been mainly focused on the model plant Arabidopsis thaliana, and have provided clues on epigenome features, components that silence pathways, and effects of silencing impairment. But to what extent can Arabidopsis be used as a model for the epigenomics of crops? In this review, we discuss the similarities and differences between the epigenomes of Arabidopsis and crops. We explore the relationship between TEs and epigenomes, focusing on TE silencing control and escape, and the impact of TE mobility on epigenomic variation. Finally, we provide insights into challenges to tackle, and future directions to take in the route towards using epigenetic diversity in plant breeding programmes.

Chromatin structure and gene expression changes associated with loss of MOP1 activity in Zea mays

Though the mechanisms governing nuclear organization are not well understood, it is apparent that epigenetic modifications coordinately modulate chromatin organization as well as transcription. In maize, MEDIATOR OF PARAMUTATION1 (MOP1) is required for 24 nt siRNA-mediated epigenetic regulation and transcriptional gene silencing via a putative Pol IV- RdDM pathway. To elucidate the mechanisms of nuclear chromatin organization, we investigated the relationship between chromatin structure and transcription in response to loss of MOP1 function. We used a microarray based micrococcal nuclease sensitivity assay to identify genome-wide changes in chromatin structure in mop1-1 immature ears and observed an increase in chromatin accessibility at chromosome arms associated with loss of MOP1 function. Within the many genes misregulated in mop1 mutants, we identified one subset likely to be direct targets of epigenetic transcriptional silencing via Pol-IV RdDM. We found that target specificity for MOP1-mediated RdDM activity is governed by multiple signals that include accumulation of 24 nt siRNAs and the presence of specific classes of gene-proximal transposons, but neither of these attributes alone is sufficient to predict transcriptional misregulation in mop1-1 homozygous mutants. Our results suggest a role for MOP1 in regulation of higher-order chromatin organization where loss of MOP1 activity at a subset of loci triggers a broader cascade of transcriptional consequences and genome-wide changes in chromatin structure.

Epigenetics and plant genome evolution

Epigenetics was envisioned as a topic to inform evolutionary theory, but the interplay between epigenetics and evolution has received little attention. With the advent of high-throughput methods, it is now routine to measure the genome-wide distribution of epigenetic marks, and these genome-wide patterns are providing insights into evolutionary processes. For example, DNA methylation is associated with transposable element silencing but also with repression of the expression of nearby genes, perhaps caused by the spread of methylation into regulatory regions. This repressive effect, which is typically deleterious, is acted upon by purifying selection. These dynamics may also govern the outcome of hybridization and polyploid events by affecting homoeolog expression. Finally, genes are also often methylated, but the implications of genic methylation for plant gene and genome evolution are not yet characterized fully.

A MITE transposon insertion is associated with differential methylation at the maize flowering time QTL Vgt1

One of the major quantitative trait loci for flowering time in maize, the Vegetative to generative transition 1 (Vgt1) locus, corresponds to an upstream (70 kb) noncoding regulatory element of ZmRap2.7, a repressor of flowering. At Vgt1, a miniature transposon (MITE) insertion into a conserved noncoding sequence was previously found to be highly associated with early flowering in independent studies. Because cytosine methylation is known to be associated with transposons and to influence gene expression, we aimed to investigate how DNA methylation patterns in wild-type and mutant Vgt1 correlate with ZmRap2.7 expression. The methylation state at Vgt1 was assayed in leaf samples of maize inbred and F1 hybrid samples, and at the syntenic region in sorghum. The Vgt1-linked conserved noncoding sequence was very scarcely methylated both in maize and sorghum. However, in the early maize Vgt1 allele, the region immediately flanking the highly methylated MITE insertion was significantly more methylated and showed features of methylation spreading. Allele-specific expression assays revealed that the presence of the MITE and its heavy methylation appear to be linked to altered ZmRap2.7 transcription. Although not providing proof of causative connection, our results associate transposon-linked differential methylation with allelic state and gene expression at a major flowering time quantitative trait locus in maize.

Differential methylation during maize leaf growth targets developmentally regulated genes

DNA methylation is an important and widespread epigenetic modification in plant genomes, mediated by DNA methyltransferases (DMTs). DNA methylation is known to play a role in genome protection, regulation of gene expression, and splicing and was previously associated with major developmental reprogramming in plants, such as vernalization and transition to flowering. Here, we show that DNA methylation also controls the growth processes of cell division and cell expansion within a growing organ. The maize (Zea mays) leaf offers a great tool to study growth processes, as the cells progressively move through the spatial gradient encompassing the division zone, transition zone, elongation zone, and mature zone. Opposite to de novo DMTs, the maintenance DMTs were transcriptionally regulated throughout the growth zone of the maize leaf, concomitant with differential CCGG methylation levels in the four zones. Surprisingly, the majority of differentially methylated sequences mapped on or close to gene bodies and not to repeat-rich loci. Moreover, especially the 5' and 3' regions of genes, which show overall low methylation levels, underwent differential methylation in a developmental context. Genes involved in processes such as chromatin remodeling, cell cycle progression, and growth regulation, were differentially methylated. The presence of differential methylation located upstream of the gene anticorrelated with transcript expression, while gene body differential methylation was unrelated to the expression level. These data indicate that DNA methylation is correlated with the decision to exit mitotic cell division and to enter cell expansion, which adds a new epigenetic level to the regulation of growth processes.

Accessing epigenetic variation in the plant methylome

Cytosine DNA methylation is the addition of a methyl group to the 5' position of a cytosine, which plays a crucial role in plant development and gene silencing. Genome-wide profiling of DNA methylation is now possible using various techniques and strategies. Using these technologies, we are beginning to elucidate the extent and impact of variation in DNA methylation between individuals and/or tissues. Here, we review the different techniques used to analyze the methylomes at the whole-genome level and their applications to better understand epigenetic variations in plants.

Sequencing-based large-scale genomics approaches with small numbers of isolated maize meiocytes

High-throughput sequencing has become the large-scale approach of choice to study global gene expression and the distribution of specific chromatin marks and features. However, the limited availability of large amounts of purified cells made it very challenging to apply sequencing-based techniques in plant meiosis research in the past. In this paper, we describe a method to isolate meiocytes from maize anthers and detailed protocols to successfully perform RNA-seq, smRNA-seq, H3K4me3-ChIP-seq, and DNA bisulfite conversion sequencing with 5000-30,000 isolated maize male meiotic cells. These methods can be adjusted for other flowering plant species as well.

Vernalization treatment induces site-specific DNA hypermethylation at the VERNALIZATION-A1 (VRN-A1) locus in hexaploid winter wheat

BACKGROUND

Certain temperate species require prolonged exposure to low temperature to initiate transition from vegetative growth to flowering, a process known as vernalization. In wheat, winter cultivars require vernalization to initiate flowering, making vernalization requirement a trait of key importance in wheat agronomy. The genetic bases of vernalization response have been largely studied in wheat, leading to the characterization of a regulation pathway that involves the key gene VERNALIZATION1 (VRN1). While previous studies in wheat and barley have revealed the functional role of histone modification in setting VRN1 expression, other mechanisms might also be involved. Here, we were interested in determining whether the cold-induced expression of the wheat VRN-A1 gene is associated with a change in DNA methylation.

RESULTS

We provide the first DNA methylation analysis of the VRN-A1 gene, and describe the existence of methylation at CG but also at non CG sites. While CG sites show a bell-shape profile typical of gene-body methylation, non CG methylation is restricted to the large (8.5 kb) intron 1, in a region harboring fragments of transposable elements (TEs). Interestingly, cold induces a site-specific hypermethylation at these non CG sites. This increase in DNA methylation is transmitted through mitosis, and is reset to its original level after sexual reproduction.

CONCLUSIONS

These results demonstrate that VRN-A1 has a particular DNA methylation pattern, exhibiting rapid shift within the life cycle of a winter wheat plant following exposure to particular environmental conditions. The finding that this shift occurs at non CG sites in a TE-rich region opens interesting questions onto the possible consequences of this type of methylation in gene expression.

Epigenetic and genetic influences on DNA methylation variation in maize populations

DNA methylation is a chromatin modification that is frequently associated with epigenetic regulation in plants and mammals. However, genetic changes such as transposon insertions can also lead to changes in DNA methylation. Genome-wide profiles of DNA methylation for 20 maize (Zea mays) inbred lines were used to discover differentially methylated regions (DMRs). The methylation level for each of these DMRs was also assayed in 31 additional maize or teosinte genotypes, resulting in the discovery of 1966 common DMRs and 1754 rare DMRs. Analysis of recombinant inbred lines provides evidence that the majority of DMRs are heritable. A local association scan found that nearly half of the DMRs with common variation are significantly associated with single nucleotide polymorphisms found within or near the DMR. Many of the DMRs that are significantly associated with local genetic variation are found near transposable elements that may contribute to the variation in DNA methylation. Analysis of gene expression in the same samples used for DNA methylation profiling identified over 300 genes with expression patterns that are significantly associated with DNA methylation variation. Collectively, our results suggest that DNA methylation variation is influenced by genetic and epigenetic changes that are often stably inherited and can influence the expression of nearby genes.

Gene body methylation is conserved between plant orthologs and is of evolutionary consequence

DNA methylation is a common feature of eukaryotic genomes and is especially common in noncoding regions of plants. Protein coding regions of plants are often methylated also, but the extent, function, and evolutionary consequences of gene body methylation remain unclear. Here we investigate gene body methylation using an explicit comparative evolutionary approach. We generated bisulfite sequencing data from two tissues of Brachypodium distachyon and compared genic methylation patterns to those of rice (Oryza sativa ssp. japonica). Gene body methylation was strongly conserved between orthologs of the two species and affected a biased subset of long, slowly evolving genes. Because gene body methylation is conserved over evolutionary time, it shapes important features of plant genome evolution, such as the bimodality of G+C content among grass genes. Our results superficially contradict previous observations of high cytosine methylation polymorphism within Arabidopsis thaliana genes, but reanalyses of these data are consistent with conservation of methylation within gene regions. Overall, our results indicate that the methylation level is a long-term property of individual genes and therefore of evolutionary consequence.