Heterochromatin is a quantitative trait associated with spontaneous epiallele formation

Mapping parental DMRs predictive of local and distal methylome remodeling in epigenetic F1 hybrids

F1 hybrids derived from a cross between two inbred parental lines often display widespread changes in DNA methylation and gene expression patterns relative to their parents. An emerging challenge is to understand how parental epigenomic differences contribute to these events. Here, we generated a large mapping panel of F1 epigenetic hybrids, whose parents are isogenic but variable in their DNA methylation patterns. Using a combination of multi-omic profiling and epigenetic mapping strategies we show that differentially methylated regions in parental pericentromeres act as major reorganizers of hybrid methylomes and transcriptomes, even in the absence of genetic variation. These parental differentially methylated regions are associated with hybrid methylation remodeling events at thousands of target regions throughout the genome, both locally (in cis) and distally (in trans). Many of these distally-induced methylation changes lead to nonadditive expression of nearby genes and associate with phenotypic heterosis. Our study highlights the pleiotropic potential of parental pericentromeres in the functional remodeling of hybrid genomes and phenotypes.

BSXplorer: analytical framework for exploratory analysis of BS-seq data

BACKGROUND

Bisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily.

RESULTS

BSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication.

CONCLUSIONS

BSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

Heritable epigenetic variation facilitates long-term maintenance of epigenetic and genetic variation

How genetic and phenotypic variation are maintained has long been one of the fundamental questions in population and quantitative genetics. A variety of factors have been implicated to explain the maintenance of genetic variation in some contexts (e.g. balancing selection), but the potential role of epigenetic regulation to influence population dynamics has been understudied. It is well recognized that epigenetic regulation, including histone methylation, small RNA expression, and DNA methylation, helps to define differences between cell types and facilitate phenotypic plasticity. In recent years, empirical studies have shown the potential for epigenetic regulation to also be heritable for at least a few generations without selection, raising the possibility that differences in epigenetic regulation can act alongside genetic variation to shape evolutionary trajectories. Heritable differences in epigenetic regulation that arise spontaneously are termed "epimutations." Epimutations differ from genetic mutations in 2 key ways-they occur at a higher rate and the loci at which they occur often revert back to their original state within a few generations. Here, we present an extension of the standard population genetic model with selection to incorporate epigenetic variation arising via epimutation. Our model assumes a diploid, sexually reproducing population with random mating. In addition to spontaneous genetic mutation, we included parameters for spontaneous epimutation and back-epimutation, allowing for 4 potential epialleles at a single locus (2 genetic alleles, each with 2 epigenetic states), each of which affect fitness. We then analyzed the conditions under which stable epialleles were maintained. Our results show that highly reversible epialleles can be maintained in long-term equilibrium under neutral conditions in a manner that depends on the epimutation and back-epimutation rates, which we term epimutation-back-epimutation equilibrium. On the other hand, epialleles that compensate for deleterious mutations cause deviations from the expectations of mutation-selection balance by a simple factor that depends on the epimutation and back-epimutation rates. We also numerically analyze several sets of fitness parameters for which large deviations from mutation-selection balance occur. Together, these results demonstrate that transient epigenetic regulation may be an important factor in the maintenance of both epigenetic and genetic variation in populations.

Natural methylation epialleles correlate with gene expression in maize

DNA methylation in plants is depleted from cis-regulatory elements in and near genes but is present in some gene bodies, including exons. Methylation in exons solely in the CG context is called gene body methylation (gbM). Methylation in exons in both CG and non-CG contexts is called TE-like methylation (teM). Assigning functions to both forms of methylation in genes has proven to be challenging. Toward that end, we utilized recent genome assemblies, gene annotations, transcription data, and methylome data to quantify common patterns of gene methylation and their relations to gene expression in maize. We found that gbM genes exist in a continuum of CG methylation levels without a clear demarcation between unmethylated genes and gbM genes. Analysis of expression levels across diverse maize stocks and tissues revealed a weak but highly significant positive correlation between gbM and gene expression except in endosperm. gbM epialleles were associated with an approximately 3% increase in steady-state expression level relative to unmethylated epialleles. In contrast to gbM genes, which were conserved and were broadly expressed across tissues, we found that teM genes, which make up about 12% of genes, are mainly silent, are poorly conserved, and exhibit evidence of annotation errors. We used these data to flag teM genes in the 26 NAM founder genome assemblies. While some teM genes are likely functional, these data suggest that the majority are not, and their inclusion can confound the interpretation of whole-genome studies.

Stochasticity in gene body methylation

Gene body methylation (gbM) is a widely conserved epigenetic feature of plant genomes. Efforts to delineate the mechanisms by which gbM contributes to transcriptional regulation remain largely inconclusive, and its evolutionary significance continues to be debated. Curiously, although steady-state gbM levels are remarkably stable across mitotic and meiotic cell divisions, the methylation status of individual CG dinucleotides in gbM genes is highly stochastic. How can these two seemingly contradictory observations be reconciled? Here, we discuss how stochastic processes relate to gbM maintenance dynamics. We show that a quantitative understanding of these processes can shed deeper insights into the molecular and evolutionary biology of this enigmatic epigenetic trait.

An evolutionary epigenetic clock in plants

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

On the causes of gene-body methylation variation in Arabidopsis thaliana

Gene-body methylation (gbM) refers to sparse CG methylation of coding regions, which is especially prominent in evolutionarily conserved house-keeping genes. It is found in both plants and animals, but is directly and stably (epigenetically) inherited over multiple generations in the former. Studies in Arabidopsis thaliana have demonstrated that plants originating from different parts of the world exhibit genome-wide differences in gbM, which could reflect direct selection on gbM, but which could also reflect an epigenetic memory of ancestral genetic and/or environmental factors. Here we look for evidence of such factors in F2 plants resulting from a cross between a southern Swedish line with low gbM and a northern Swedish line with high gbM, grown at two different temperatures. Using bisulfite-sequencing data with nucleotide-level resolution on hundreds of individuals, we confirm that CG sites are either methylated (nearly 100% methylation across sampled cells) or unmethylated (approximately 0% methylation across sampled cells), and show that the higher level of gbM in the northern line is due to more sites being methylated. Furthermore, methylation variants almost always show Mendelian segregation, consistent with their being directly and stably inherited through meiosis. To explore how the differences between the parental lines could have arisen, we focused on somatic deviations from the inherited state, distinguishing between gains (relative to the inherited 0% methylation) and losses (relative to the inherited 100% methylation) at each site in the F2 generation. We demonstrate that deviations predominantly affect sites that differ between the parental lines, consistent with these sites being more mutable. Gains and losses behave very differently in terms of the genomic distribution, and are influenced by the local chromatin state. We find clear evidence for different trans-acting genetic polymorphism affecting gains and losses, with those affecting gains showing strong environmental interactions (G×E). Direct effects of the environment were minimal. In conclusion, we show that genetic and environmental factors can change gbM at a cellular level, and hypothesize that these factors can also lead to transgenerational differences between individuals via the inclusion of such changes in the zygote. If true, this could explain genographic pattern of gbM with selection, and would cast doubt on estimates of epimutation rates from inbred lines in constant environments.

Epigenetic Changes Occurring in Plant Inbreeding

Inbreeding is the crossing of closely related individuals in nature or a plantation or self-pollinating plants, which produces plants with high homozygosity. This process can reduce genetic diversity in the offspring and decrease heterozygosity, whereas inbred depression (ID) can often reduce viability. Inbred depression is common in plants and animals and has played a significant role in evolution. In the review, we aim to show that inbreeding can, through the action of epigenetic mechanisms, affect gene expression, resulting in changes in the metabolism and phenotype of organisms. This is particularly important in plant breeding because epigenetic profiles can be linked to the deterioration or improvement of agriculturally important characteristics.

Transgenerational epigenetic control of constitutive heterochromatin, transposons, and centromeres

Epigenetic mechanisms are important not only for development but also for genome stability and chromosome dynamics. The latter types of epigenetic controls can often be transgenerational. Here, we review recent progress in two examples of transgenerational epigenetic control: i) the control of constitutive heterochromatin and transposable elements and ii) epigenetic mechanisms that regulate centromere specification and functions. We also discuss the biological significance of enigmatic associations among centromeres, transposons, and constitutive heterochromatin.

Climate-responsive DNA methylation is involved in the biosynthesis of lignin in birch

Lignin is one of the most important secondary metabolites and essential to the formation of cell walls. Changes in lignin biosynthesis have been reported to be associated with environmental variations and can influence plant fitness and their adaptation to abiotic stresses. However, the molecular mechanisms underlying this association remain unclear. In this study, we evaluated the relations between the lignin biosynthesis and environmental factors and explored the role of epigenetic modification (DNA methylation) in contributing to these relations if any in natural birch. Significantly negative correlations were observed between the lignin content and temperature ranges. Analyzing the transcriptomes of birches in two habitats with different temperature ranges showed that the expressions of genes and transcription factors (TFs) involving lignin biosynthesis were significantly reduced at higher temperature ranges. Whole-genome bisulfite sequencing revealed that promoter DNA methylation of two NAC-domain TFs, BpNST1/2 and BpSND1, may be involved in the inhibition of these gene expressions, and thereby reduced the content of lignin. Based on these results we proposed a DNA methylation-mediated lignin biosynthesis model which responds to environmental factors. Overall, this study suggests the possibility of environmental signals to induce epigenetic variations that result in changes in lignin content, which can aid to develop resilient plants to combat ongoing climate changes or to manipulate secondary metabolite biosynthesis for agricultural, medicinal, or industrial values.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture—once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

DNA methylation is an epigenetic mark common across eukaryotes. It is associated with transcriptional silencing, in particular of transposable elements. Multiple elements, including epigenetic inheritance, shape DNA methylation patterns, and the complexity makes it challenging to dissect the regulation. Our work on the 1001 Arabidopsis Epigenomes project led to the unexpected discovery that much of the natural variation for CHH methylation (mCHH; H is A, T, or C) on transposable elements could be attributed to allelic variation in three genes known to be involved in epigenetic regulation. However, our analysis of methylation in other sequence contexts (mCHG and mCG) revealed no genetic regulator.

Here we show that if mCHG variation is analyzed while taking mCHH into account, mCHG is also strongly influenced by allelic variation in a small number of genes with known or highly plausible direct roles in epigenetic regulation. The presence of common allelic variation of large effect is suggestive of some form of local adaptation. The nature of this adaptation remains obscure, but we present further evidence that allelic variation regulating DNA methylation influences transposon mobilization.

H3K9 demethylases IBM1 and JMJ27 are required for male meiosis in Arabidopsis thaliana

Dimethylation of histone H3 lysine 9 (H3K9me2), a crucial modification for heterochromatin formation and transcriptional silencing, is essential for proper meiotic prophase progression in mammals. We analyzed meiotic defects and generated genome-wide profiles of H3K9me2 and transcriptomes for the mutants of H3K9 demethylases. Moreover, we also identified proteins interacting with H3K9 demethylases. H3K9me2 is usually found at transposable elements and repetitive sequences but is absent from the bodies of protein-coding genes. In this study, we show that the Arabidopsis thaliana H3K9 demethylases IBM1 and JMJ27 cooperatively regulate crossover formation and chromosome segregation. They protect thousands of protein-coding genes from ectopic H3K9me2, including genes essential for meiotic prophase progression. In addition to removing H3K9me2, IBM1 and JMJ27 interact with the Precocious Dissociation of Sisters 5 (PDS5) cohesin complex cofactors. The pds5 mutant shared similar transcriptional alterations with ibm1 jmj27, including meiosis-essential genes, yet without affecting H3K9me2 levels. Hence, PDS5s, together with IBM1 and JMJ27, regulate male meiosis and gene expression independently of H3K9 demethylation. These findings uncover a novel role of H3K9me2 removal in meiosis and a new function of H3K9 demethylases and cohesin cofactors in meiotic transcriptional regulation.

Crosstalk among pathways to generate DNA methylome

Cytosine is methylated in both CpG and non-CpG contexts (mCG and mCH, respectively) in plant genomes. Although mCG and mCH are almost independent in regard to their "maintenance," recent studies uncovered crosstalk between them during their "establishment," which unexpectedly functions in both RNAi-dependent and -independent pathways. In addition, the importance of linker histone H1 and variants of histone H2A to DNA methylation dynamics is starting to be understood. We summarize these new aspects of mechanisms to generate DNA methylomes and discuss future prospects.

Molecular properties of epimutation hotspots

Mistakes in the maintenance of CG methylation are a source of heritable epimutations in plants. Multigenerational surveys indicate that the rate of these stochastic events varies substantially across the genome, with some regions harbouring localized 'epimutation hotspots'. Using Arabidopsis as a model, we show that epimutation hotspots are indexed by a specific set of chromatin states that map to subregions of gene body methylation genes. Although these regions comprise only ~12% of all CGs in the genome, they account for ~63% of all epimutation events per unit time. Molecular profiling revealed that these regions contain unique sequence features, harbour steady-state intermediate methylation levels and act as putative targets of antagonistic DNA methylation pathways. We further demonstrate that experimentally induced shifts in steady-state methylation in these hotspot regions are sufficient to significantly alter local epimutation intensities. Our work lays the foundation for dissecting the molecular mechanisms and evolutionary consequences of epimutation hotspots in plants.

MethylScore, a pipeline for accurate and context-aware identification of differentially methylated regions from population-scale plant whole-genome bisulfite sequencing data

Whole-genome bisulfite sequencing (WGBS) is the standard method for profiling DNA methylation at single-nucleotide resolution. Different tools have been developed to extract differentially methylated regions (DMRs), often built upon assumptions from mammalian data. Here, we present MethylScore, a pipeline to analyse WGBS data and to account for the substantially more complex and variable nature of plant DNA methylation. MethylScore uses an unsupervised machine learning approach to segment the genome by classification into states of high and low methylation. It processes data from genomic alignments to DMR output and is designed to be usable by novice and expert users alike. We show how MethylScore can identify DMRs from hundreds of samples and how its data-driven approach can stratify associated samples without prior information. We identify DMRs in the A. thaliana 1,001 Genomes dataset to unveil known and unknown genotype-epigenotype associations .