Gene body DNA methylation in plants

Dynamic evolution of the heterochromatin sensing histone demethylase IBM1

Heterochromatin is critical for maintaining genome stability, especially in flowering plants, where it relies on a feedback loop involving the H3K9 methyltransferase, KRYPTONITE (KYP), and the DNA methyltransferase CHROMOMETHYLASE3 (CMT3). The H3K9 demethylase INCREASED IN BONSAI METHYLATION 1 (IBM1) counteracts the detrimental consequences of KYP-CMT3 activity in transcribed genes. IBM1 expression in Arabidopsis is uniquely regulated by methylation of the 7th intron, allowing it to monitor global H3K9me2 levels. We show the methylated intron is prevalent across flowering plants and its underlying sequence exhibits dynamic evolution. We also find extensive genetic and expression variations in KYP, CMT3, and IBM1 across flowering plants. We identify Arabidopsis accessions resembling weak ibm1 mutants and Brassicaceae species with reduced IBM1 expression or deletions. Evolution towards reduced IBM1 activity in some flowering plants could explain the frequent natural occurrence of diminished or lost CMT3 activity and loss of gene body DNA methylation, as cmt3 mutants in A. thaliana mitigate the deleterious effects of IBM1.

Targeted bisulfite sequencing of Scots pine adaptation-related genes

During environmental changes, epigenetic processes can enable adaptive responses faster than natural selection. In plants, very little is known about the role of DNA methylation during long-term adaptation. Scots pine is a widely distributed coniferous species which must adapt to different environmental conditions throughout its long lifespan. Thus, epigenetic modifications may contribute towards this direction. We provide bisulfite next-generation sequencing data from the putative promoters and exons of eight adaptation-related genes (A3IP2, CCA1, COL1, COL2, FTL2, MFT1, PHYO, and ZTL) in three Scots pine populations located in northern and southern parts of Finland. DNA methylation levels were studied in the two seed tissues: the maternal megagametophyte which contributes to embryo viability, and the biparental embryo which represents the next generation. In most genes, differentially methylated cytosines (DMCs) were in line with our previously demonstrated gene expression differences found in the same Scots pine populations. In addition, we found a strong correlation of total methylation levels between the embryo and megagametophyte tissues of a given individual tree, which indicates that DNA methylation can be inherited from the maternal parent. In conclusion, our results imply that DNA methylation differences may contribute to the adaptation of Scots pine populations in different climatic conditions.

Differentially methylated genes involved in reproduction and ploidy levels in recent diploidized and tetraploidized Eragrostis curvula genotypes

Epigenetics studies changes in gene activity without changes in the DNA sequence. Methylation is an epigenetic mechanism important in many pathways, such as biotic and abiotic stresses, cell division, and reproduction. Eragrostis curvula is a grass species reproducing by apomixis, a clonal reproduction by seeds. This work employed the MCSeEd technique to identify deferentially methylated positions, regions, and genes in the CG, CHG, and CHH contexts in E. curvula genotypes with similar genomic backgrounds but with different reproductive modes and ploidy levels. In this way, we focused the analysis on the cvs. Tanganyika INTA (4x, apomictic), Victoria (2x, sexual), and Bahiense (4x, apomictic). Victoria was obtained from the diploidization of Tanganyika INTA, while Bahiense was produced from the tetraploidization of Victoria. This study showed that polyploid/apomictic genotypes had more differentially methylated positions and regions than the diploid sexual ones. Interestingly, it was possible to observe fewer differentially methylated positions and regions in CG than in the other contexts, meaning CG methylation is conserved across the genotypes regardless of the ploidy level and reproductive mode. In the comparisons between sexual and apomictic genotypes, we identified differentially methylated genes involved in the reproductive pathways, specifically in meiosis, cell division, and fertilization. Another interesting observation was that several differentially methylated genes between the diploid and the original tetraploid genotype recovered their methylation status after tetraploidization, suggesting that methylation is an important mechanism involved in reproduction and ploidy changes.

Recent Advances in Studies of Genomic DNA Methylation and Its Involvement in Regulating Drought Stress Response in Crops

As global arid conditions worsen and groundwater resources diminish, drought stress has emerged as a critical impediment to plant growth and development globally, notably causing declines in crop yields and even the extinction of certain cultivated species. Numerous studies on drought resistance have demonstrated that DNA methylation dynamically interacts with plant responses to drought stress by modulating gene expression and developmental processes. However, the precise mechanisms underlying these interactions remain elusive. This article consolidates the latest research on the role of DNA methylation in plant responses to drought stress across various species, focusing on methods of methylation detection, mechanisms of methylation pattern alteration (including DNA de novo methylation, DNA maintenance methylation, and DNA demethylation), and overall responses to drought conditions. While many studies have observed significant shifts in genome-wide or gene promoter methylation levels in drought-stressed plants, the identification of specific genes and pathways involved remains limited. This review aims to furnish a reference for detailed research into plant responses to drought stress through epigenetic approaches, striving to identify drought resistance genes regulated by DNA methylation, specific signaling pathways, and their molecular mechanisms of action.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

The fungicide pyraclostrobin affects gene expression by altering the DNA methylation pattern in Magnaporthe oryzae

Introduction

Rice blast disease caused by Magnaporthe oryzae has long been the main cause of rice (Oryza sativa L.) yield reduction worldwide. The quinone external inhibitor pyraclostrobin is widely used as a fungicide to effectively control the spread of pathogenic fungi, including M. oryzae. However, M. oryzae can develop resistance through multiple levels of mutation, such as target protein cytb mutation G143A/S, leading to a decrease in the effectiveness of the biocide after a period of application. Therefore, uncovering the possible mutational mechanisms from multiple perspectives will further provide feasible targets for drug development.

Methods

In this work, we determined the gene expression changes in M. oryzae in response to pyraclostrobin stress and their relationship with DNA methylation by transcriptome and methylome.

Results

The results showed that under pyraclostrobin treatment, endoplasmic reticulum (ER)-associated and ubiquitin-mediated proteolysis were enhanced, suggesting that more aberrant proteins may be generated that need to be cleared. DNA replication and repair processes were inhibited. Glutathione metabolism was enhanced, while lipid metabolism was impaired. The number of alternative splicing events increased. These changes may be related to the elevated methylation levels of cytosine and adenine in gene bodies. Both hypermethylation and hypomethylation of differentially methylated genes (DMGs) mainly occurred in exons and promoters. Some DMGs and differentially expressed genes (DEGs) were annotated to the same pathways by GO and KEGG, including protein processing in the ER, ubiquitin-mediated proteolysis, RNA transport and glutathione metabolism, suggesting that pyraclostrobin may affect gene expression by altering the methylation patterns of cytosine and adenine.

Discussion

Our results revealed that 5mC and 6mA in the gene body are associated with gene expression and contribute to adversity adaptation in M. oryzae. This enriched the understanding for potential mechanism of quinone inhibitor resistance, which will facilitate the development of feasible strategies for maintaining the high efficacy of this kind of fungicide.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages

BACKGROUND

Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored.

METHODS

In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages.

RESULTS

Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles.

CONCLUSION

These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

DNA methylation: an emerging paradigm of gene regulation under drought stress in plants

Drought is an enormous threat to global crop production. In order to ensure food security for the burgeoning population, we must develop drought tolerant crop varieties. This necessitates the identification of drought-responsive genes and understanding the mechanisms involved in their regulation. DNA methylation is a widely studied mechanism of epigenetic regulation of gene expression, which is known to play vital role in conferring tolerance to various biotic and abiotic stress factors. The recent advances in next-generation sequencing (NGS) technologies, has allowed unprecedented access to genome-wide methylation marks, with single base resolution. The most important roles of DNA methylation have been studied in terms of gene body methylation (gbM), which is associated with regulation of both transcript abundance and its stability. The availability of mutants for the various genes encoding enzymes involved in methylation of DNA has allowed ascertainment of the biological significance of methylation. Even though a vast number of reports have emerged in the recent past, where both genome-wide methylation landscape and locus specific changes in DNA methylation have been studied, a conclusive picture with regards to the biological role of DNA methylation is still lacking. Compounding this, is the lack of sufficient evidence supporting the heritability of these epigenetic changes. Amongst the various epigenetic variations, the DNA methylation changes are observed to be the most stable. This review describes the drought-induced changes in DNA methylation identified across different plant species. We also briefly describe the stress memory contributed by these changes. The identification of heritable, drought-induced methylation marks would broaden the scope of crop improvement in the future.

Genome-wide methylome stability and parental effects in the worldwide distributed Lombardy poplar

BACKGROUND

Despite the increasing number of epigenomic studies in plants, little is known about the forces that shape the methylome in long-lived woody perennials. The Lombardy poplar offers an ideal opportunity to investigate the impact of the individual environmental history of trees on the methylome.

RESULTS

We present the results of three interconnected experiments on Lombardy poplar. In the first experiment, we investigated methylome variability during a growing season and across vegetatively reproduced generations. We found that ramets collected over Europe and raised in common conditions have stable methylomes in symmetrical CG-contexts. In contrast, seasonal dynamics occurred in methylation patterns in CHH context. In the second experiment, we investigated whether methylome patterns of plants grown in a non-parental environment correlate with the parental climate. We did not observe a biological relevant pattern that significantly correlates with the parental climate. Finally, we investigated whether the parental environment has persistent carry-over effects on the vegetative offspring's phenotype. We combined new bud set observations of three consecutive growing seasons with former published bud set data. Using a linear mixed effects analysis, we found a statistically significant but weak short-term, parental carry-over effect on the timing of bud set. However, this effect was negligible compared to the direct effects of the offspring environment.

CONCLUSIONS

Genome-wide cytosine methylation patterns in symmetrical CG-context are stable in Lombardy poplar and appear to be mainly the result of random processes. In this widespread poplar clone, methylation patterns in CG-context can be used as biomarkers to infer a common ancestor and thus to investigate the recent environmental history of a specific Lombardy poplar. The Lombardy poplar shows high phenotypic plasticity in a novel environment which enabled this clonal tree to adapt and survive all over the temperate regions of the world.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

Loops, crosstalk, and compartmentalization: it takes many layers to regulate DNA methylation

DNA methylation is a conserved epigenetic modification associated with transposon silencing and gene regulation. The stability of this modification relies on intimate connections between DNA and histone modifications that generate self-reinforcing loops wherein the presence of one mark promotes the other. However, it is becoming increasingly clear that the efficiency of these loops is affected by cross-talk between pathways and by chromatin accessibility, which is heavily influenced by histone variants. Focusing primarily on plants, this review provides an update on the aforementioned self-reinforcing loops, highlights recent advances in understanding how DNA methylation pathways are restricted to prevent encroachment on genes, and discusses the roles of histone variants in compartmentalizing epigenetic pathways within the genome. This multilayered approach facilitates two essential, yet opposing functions, the ability to maintain heritable DNA methylation patterns while retaining the flexibility to modify these patterns during development.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

Arabidopsis thaliana: a powerful model organism to explore histone modifications and their upstream regulations

Histones are subjected to extensive covalent modifications that affect inter-nucleosomal interactions as well as alter chromatin structure and DNA accessibility. Through switching the corresponding histone modifications, the level of transcription and diverse downstream biological processes can be regulated. Although animal systems are widely used in studying histone modifications, the signalling processes that occur outside the nucleus prior to histone modifications have not been well understood due to the limitations including non viable mutants, partial lethality, and infertility of survivors. Here, we review the benefits of using Arabidopsis thaliana as the model organism to study histone modifications and their upstream regulations. Similarities among histones and key histone modifiers such as the Polycomb group (PcG) and Trithorax group (TrxG) in Drosophila, Human, and Arabidopsis are examined. Furthermore, prolonged cold-induced vernalization system has been well-studied and revealed the relationship between the controllable environment input (duration of vernalization), its chromatin modifications of FLOWERING LOCUS C (FLC), following gene expression, and the corresponding phenotypes. Such evidence suggests that research on Arabidopsis can bring insights into incomplete signalling pathways outside of the histone box, which can be achieved through viable reverse genetic screenings based on the phenotypes instead of direct monitoring of histone modifications among individual mutants. The potential upstream regulators in Arabidopsis can provide cues or directions for animal research based on the similarities between them.

DNA methylome analysis provides insights into gene regulatory mechanism for better performance of rice under fluctuating environmental conditions: epigenomics of adaptive plasticity

MAIN CONCLUSION

Roots play an important role in adaptive plasticity of rice under dry/direct-sown conditions. However, hypomethylation of genes in leaves (resulting in up-regulated expression) complements the adaptive plasticity of Nagina-22 under DSR conditions. Rice is generally cultivated by transplanting which requires plenty of water for irrigation. Such a practice makes rice cultivation a challenging task under global climate change and reducing water availability. However, dry-seeded/direct-sown rice (DSR) has emerged as a resource-saving alternative to transplanted rice (TPR). Though some of the well-adapted local cultivars are used for DSR, only limited success has been achieved in developing DSR varieties mainly because of a limited knowledge of adaptability of rice under fluctuating environmental conditions. Based on better morpho-physiological and agronomic performance of Nagina-22 (N-22) under DSR conditions, N-22 and IR-64 were grown by transplanting and direct-sowing and used for whole genome methylome analysis to unravel the epigenetic basis of adaptive plasticity of rice. Comparative methylome and transcriptome analyses indicated a large number (4078) of genes regulated through DNA methylation/demethylation in N-22 under DSR conditions. Gene × environment interactions play important roles in adaptive plasticity of rice under direct-sown conditions. While genes for pectinesterase, LRK10, C2H2 zinc-finger protein, splicing factor, transposable elements, and some of the unannotated proteins were hypermethylated, the genes for regulation of transcription, protein phosphorylation, etc. were hypomethylated in CG context in the root of N-22, which played important roles in providing adaptive plasticity to N-22 under DSR conditions. Hypomethylation leading to up-regulation of gene expression in the leaf complements the adaptive plasticity of N-22 under DSR conditions. Moreover, differential post-translational modification of proteins and chromatin assembly/disassembly through DNA methylation in CHG context modulate adaptive plasticity of N-22. These findings would help developing DSR cultivars for increased water-productivity and ecological efficiency.

Millennia-long epigenetic fluctuations generate intragenic DNA methylation variance in Arabidopsis populations

Methylation of CG dinucleotides (mCGs), which regulates eukaryotic genome functions, is epigenetically propagated by Dnmt1/MET1 methyltransferases. How mCG is established and transmitted across generations despite imperfect enzyme fidelity is unclear. Whether mCG variation in natural populations is governed by genetic or epigenetic inheritance also remains mysterious. Here, we show that MET1 de novo activity, which is enhanced by existing proximate methylation, seeds and stabilizes mCG in Arabidopsis thaliana genes. MET1 activity is restricted by active demethylation and suppressed by histone variant H2A.Z, producing localized mCG patterns. Based on these observations, we develop a stochastic mathematical model that precisely recapitulates mCG inheritance dynamics and predicts intragenic mCG patterns and their population-scale variation given only CG site spacing. Our results demonstrate that intragenic mCG establishment, inheritance, and variance constitute a unified epigenetic process, revealing that intragenic mCG undergoes large, millennia-long epigenetic fluctuations and can therefore mediate evolution on this timescale.

Insights into the Evolution of Ohnologous Sequences and Their Epigenetic Marks Post-WGD in Malus Domestica

A Whole Genome Duplication (WGD) event occurred several Ma in a Rosaceae ancestor, giving rise to the Maloideae subfamily which includes today many pome fruits such as pear (Pyrus communis) and apple (Malus domestica). This complete and well-conserved genome duplication makes the apple an organism of choice to study the early evolutionary events occurring to ohnologous chromosome fragments. In this study, we investigated gene sequence evolution and expression, transposable elements (TE) density, and DNA methylation level. Overall, we identified 16,779 ohnologous gene pairs in the apple genome, confirming the relatively recent WGD. We identified several imbalances in QTL localization among duplicated chromosomal fragments and characterized various biases in genome fractionation, gene transcription, TE densities, and DNA methylation. Our results suggest a particular chromosome dominance in this autopolyploid species, a phenomenon that displays similarities with subgenome dominance that has only been described so far in allopolyploids.

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.

Warmer temperature during asexual reproduction induce methylome, transcriptomic, and lasting phenotypic changes in Fragaria vesca ecotypes

Plants must adapt with increasing speed to global warming to maintain their fitness. One rapid adaptation mechanism is epigenetic memory, which may provide organisms sufficient time to adapt to climate change. We studied how the perennial Fragaria vesca adapted to warmer temperatures (28°C vs. 18°C) over three asexual generations. Differences in flowering time, stolon number, and petiole length were induced by warmer temperature in one or more ecotypes after three asexual generations and persisted in a common garden environment. Induced methylome changes differed between the four ecotypes from Norway, Iceland, Italy, and Spain, but shared methylome responses were also identified. Most differentially methylated regions (DMRs) occurred in the CHG context, and most CHG and CHH DMRs were hypermethylated at the warmer temperature. In eight CHG DMR peaks, a highly similar methylation pattern could be observed between ecotypes. On average, 13% of the differentially methylated genes between ecotypes also showed a temperature-induced change in gene expression. We observed ecotype-specific methylation and expression patterns for genes related to gibberellin metabolism, flowering time, and epigenetic mechanisms. Furthermore, we observed a negative correlation with gene expression when repetitive elements were found near (±2 kb) or inside genes. In conclusion, lasting phenotypic changes indicative of an epigenetic memory were induced by warmer temperature and were accompanied by changes in DNA methylation patterns. Both shared methylation patterns and transcriptome differences between F. vesca accessions were observed, indicating that DNA methylation may be involved in both general and ecotype-specific phenotypic variation.

Epitype-inducing temperatures drive DNA methylation changes during somatic embryogenesis in the long-lived gymnosperm Norway spruce

An epigenetic memory of the temperature sum experienced during embryogenesis is part of the climatic adaptation strategy of the long-lived gymnosperm Norway spruce. This memory has a lasting effect on the timing of bud phenology and frost tolerance in the resulting epitype trees. The epigenetic memory is well characterized phenotypically and at the transcriptome level, but to what extent DNA methylation changes are involved have not previously been determined. To address this, we analyzed somatic epitype embryos of Norway spruce clones produced at contrasting epitype-inducing conditions (18 and 28°C). We screened for differential DNA methylation in 2744 genes related mainly to the epigenetic machinery, circadian clock, and phenology. Of these genes, 68% displayed differential DNA methylation patterns between contrasting epitype embryos in at least one methylation context (CpG, CHG, CHH). Several genes related to the epigenetic machinery (e.g., DNA methyltransferases, ARGONAUTE) and the control of bud phenology (FTL genes) were differentially methylated. This indicates that the epitype-inducing temperature conditions induce an epigenetic memory involving specific DNA methylation changes in Norway spruce.

Depositing centromere repeats induces heritable intragenic heterochromatin establishment and spreading in Arabidopsis

Stable transmission of non-DNA-sequence-based epigenetic information contributes to heritable phenotypic variants and thus to biological diversity. While studies on spontaneous natural epigenome variants have revealed an association of epialleles with a wide range of biological traits in both plants and animals, the function, transmission mechanism, and stability of an epiallele over generations in a locus-specific manner remain poorly investigated. Here, we invented a DNA sequence deposition strategy to generate a locus-specific epiallele by depositing CEN180 satellite repeats into a euchromatic target locus in Arabidopsis. Using CRISPR/Cas9-mediated knock-in system, we demonstrated that depositing CEN180 repeats can induce heterochromatin nucleation accompanied by DNA methylation, H3K9me2, and changes in the nucleosome occupancy at the insertion sites. Interestingly, both DNA methylation and H3K9me2 are restricted within the depositing sites and depletion of an H3K9me2 demethylase IBM1 enables the outward heterochromatin propagation into the neighboring regions, leading to inheritable target gene silencing to persist for at least five generations. Together, these results demonstrate the promise of employing a cis-engineering system for the creation of stable and site-specific epialleles and provide important insights into functional epigenome studies and locus-specific transgenerational epigenetic inheritance.

Transcriptomic and methylomic analyses provide insights into the molecular mechanism and prediction of heterosis in rice

Heterosis has been widely used in multiple crops. However, the molecular mechanism and prediction of heterosis remains elusive. We generated five F1 hybrids [four showing better-parent heterosis (BPH) and one showing mid-parent heterosis], and performed the transcriptomic and methylomic analyses to identify the candidate genes for BPH and explore the molecular mechanism of heterosis and the potential predictors for heterosis. Transcriptomic results showed that most of the differentially expressed genes shared in the four better-parent hybrids were significantly enriched into the terms of molecular function, and the additive and dominant effects played crucial roles for BPH. DNA methylation level, especially in CG context, significantly and positively correlated with grain yield per plant. The ratios of differentially methylated regions in CG context in exons to transcription start sites between the parents exhibited significantly negative correlation with the heterosis levels of their hybrids, as was further confirmed in 24 pairwise comparisons of other rice lines, implying that this ratio could be a feasible predictor for heterosis level, and this ratio of less than 5 between parents in early growth stages might be a critical index for judging that their F1 hybrids would show BPH. Additionally, we identified some important genes showing differential expression and methylation, such as OsDCL2, Pi5, DTH2, DTH8, Hd1 and GLW7 in the four better-parent hybrids as the candidate genes for BPH. Our findings helped shed more light on the molecular mechanism and heterosis prediction.

Distinct regulatory pathways contribute to dynamic CHH methylation patterns in transposable elements throughout Arabidopsis embryogenesis

CHH methylation (mCHH) increases gradually during embryogenesis across dicotyledonous plants, indicating conserved mechanisms of targeting and conferral. Although it is suggested that methylation increase during embryogenesis enhances transposable element silencing, the detailed epigenetic pathways underlying this process remain unclear. In Arabidopsis, mCHH is regulated by both small RNA-dependent DNA methylation (RdDM) and RNA-independent Chromomethylase 2 (CMT2) pathways. Here, we conducted DNA methylome profiling at five stages of Arabidopsis embryogenesis, and classified mCHH regions into groups based on their dependency on different methylation pathways. Our analysis revealed that the gradual increase in mCHH in embryos coincided with the expansion of small RNA expression and regional mCHH spreading to nearby sites at numerous loci. We identified distinct methylation dynamics in different groups of mCHH targets, which vary according to transposon length, location, and cytosine frequency. Finally, we highlight the characteristics of transposable element loci that are targeted by different mCHH machinery, showing that short, heterochromatic TEs with lower mCHG levels are enriched in loci that switch from CMT2 regulation in leaves, to RdDM regulation during embryogenesis. Our findings highlight the interplay between the length, location, and cytosine frequency of transposons and the mCHH machinery in modulating mCHH dynamics during embryogenesis.

On the prediction of non-CG DNA methylation using machine learning

DNA methylation can be detected and measured using sequencing instruments after sodium bisulfite conversion, but experiments can be expensive for large eukaryotic genomes. Sequencing nonuniformity and mapping biases can leave parts of the genome with low or no coverage, thus hampering the ability of obtaining DNA methylation levels for all cytosines. To address these limitations, several computational methods have been proposed that can predict DNA methylation from the DNA sequence around the cytosine or from the methylation level of nearby cytosines. However, most of these methods are entirely focused on CG methylation in humans and other mammals. In this work, we study, for the first time, the problem of predicting cytosine methylation for CG, CHG and CHH contexts on six plant species, either from the DNA primary sequence around the cytosine or from the methylation levels of neighboring cytosines. In this framework, we also study the cross-species prediction problem and the cross-context prediction problem (within the same species). Finally, we show that providing gene and repeat annotations allows existing classifiers to significantly improve their prediction accuracy. We introduce a new classifier called AMPS (annotation-based methylation prediction from sequence) that takes advantage of genomic annotations to achieve higher accuracy.

The plant response to high CO2 levels is heritable and orchestrated by DNA methylation

Plant responses to abiotic environmental challenges are known to have lasting effects on the plant beyond the initial stress exposure. Some of these lasting effects are transgenerational, affecting the next generation. The plant response to elevated carbon dioxide (CO₂ ) levels has been well studied. However, these investigations are typically limited to plants grown for a single generation in a high CO₂ environment while transgenerational studies are rare. We aimed to determine transgenerational growth responses in plants after exposure to high CO₂ by investigating the direct progeny when returned to baseline CO₂ levels. We found that both the flowering plant Arabidopsis thaliana and seedless nonvascular plant Physcomitrium patens continue to display accelerated growth rates in the progeny of plants exposed to high CO₂ . We used the model species Arabidopsis to dissect the molecular mechanism and found that DNA methylation pathways are necessary for heritability of this growth response. More specifically, the pathway of RNA-directed DNA methylation is required to initiate methylation and the proteins CMT2 and CMT3 are needed for the transgenerational propagation of this DNA methylation to the progeny plants. Together, these two DNA methylation pathways establish and then maintain a cellular memory to high CO₂ exposure.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

DDM1-mediated gene body DNA methylation is associated with inducible activation of defense-related genes in Arabidopsis

BACKGROUND

Plants memorize previous pathogen attacks and are "primed" to produce a faster and stronger defense response, which is critical for defense against pathogens. In plants, cytosines in transposons and gene bodies are reported to be frequently methylated. Demethylation of transposons can affect disease resistance by regulating the transcription of nearby genes during defense response, but the role of gene body methylation (GBM) in defense responses remains unclear.

RESULTS

Here, we find that loss of the chromatin remodeler decrease in DNA methylation 1 (ddm1) synergistically enhances resistance to a biotrophic pathogen under mild chemical priming. DDM1 mediates gene body methylation at a subset of stress-responsive genes with distinct chromatin properties from conventional gene body methylated genes. Decreased gene body methylation in loss of ddm1 mutant is associated with hyperactivation of these gene body methylated genes. Knockout of glyoxysomal protein kinase 1 (gpk1), a hypomethylated gene in ddm1 loss-of-function mutant, impairs priming of defense response to pathogen infection in Arabidopsis. We also find that DDM1-mediated gene body methylation is prone to epigenetic variation among natural Arabidopsis populations, and GPK1 expression is hyperactivated in natural variants with demethylated GPK1.

CONCLUSIONS

Based on our collective results, we propose that DDM1-mediated GBM provides a possible regulatory axis for plants to modulate the inducibility of the immune response.

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.

The SPL transcription factor genes are potential targets for epigenetic regulation in response to drought stress in chickpea (C. arietinum L.)

BACKGROUND

Crop improvement for tolerance to various biotic and abiotic stress factors necessitates understanding the key gene regulatory mechanisms. One such mechanism of gene regulation involves changes in cytosine methylation at the gene body and flanking regulatory sequences. The present study was undertaken to identify genes which might be potential targets of drought-induced DNA methylation in chickpea.

METHODS AND RESULTS

Two chickpea genotypes, which contrast for drought tolerance, were subjected to drought stress conditions and their differential response was studied by analysing different morpho-physiological traits. Utilizing the in-house, high throughput sequencing data, the SQUAMOSA promoter-binding (SBP) protein-like (SPL) transcription factor genes were identified to be differentially methylated and expressed amongst the two genotypes, in response to drought stress. The methylation status of one of these genes was examined and validated through bisulfite PCR (BS-PCR). The identified genes could be possible homologs to known epialleles and can therefore serve as potential epialleles which can be utilized for crop improvement in chickpea.

CONCLUSION

The SPL TF genes are potential targets of epigenetic regulation in response to drought stress in chickpea. Since these are TFs, they might play important roles in controlling the expression of other genes, thus contributing to differential drought response of the two genotypes.

Transcriptome and DNA methylome divergence of inflorescence development between two ecotypes in Panicum hallii

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the two major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMRs-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

The Chromatin Remodeling Factor BrCHR39 Targets DNA Methylation to Positively Regulate Apical Dominance in Brassica rapa

The SHPRH (SNF2, histone linker, PHD, RING, helicase) subfamily belonging to ATP-dependent chromatin remodeling factor is the effective tumor-suppressor, which can polyubiquitinate PCNA (proliferating cell nuclear antigen) and participate in post-replication repair in human. However, little is known about the functions of SHPRH proteins in plants. In this study, we identified a novel SHPRH member BrCHR39 and obtained BrCHR39-silenced transgenic Brassica rapa. In contrast to wild-type plants, transgenic Brassica plants exhibited a released apical dominance phenotype with semi-dwarfism and multiple lateral branches. Furthermore, a global alteration of DNA methylation in the main stem and bud appeared after silencing of BrCHR39. Based on the GO (gene ontology) functional annotation and KEGG (Kyoto encyclopedia of genes and genomes) pathway analysis, the plant hormone signal transduction pathway was clearly enriched. In particular, we found a significant increase in the methylation level of auxin-related genes in the stem, whereas auxin- and cytokinin-related genes were hypomethylated in the bud of transgenic plants. In addition, further qRT-PCR (quantitative real-time PCR) analysis revealed that DNA methylation level always had an opposite trend with gene expression level. Considered together, our findings indicated that suppression of BrCHR39 expression triggered the methylation divergence of hormone-related genes and subsequently affected transcription levels to regulate the apical dominance in Brassica rapa.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

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.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C₁₈ Δ⁹-elongase and a C₂₀ Δ⁵-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ⁵-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Changes in the rearing environment cause reorganization of molecular networks associated with DNA methylation

Disentangling the interaction between the genetic basis and environmental context underlying phenotypic variation is critical for understanding organismal evolution. Environmental change, such as increased rates of urbanization, can induce shifts in phenotypic plasticity with some individuals adapting to city life while others are displaced. A key trait that can facilitate adaptation is the degree at which animals respond to stressors. This stress response, which includes elevation of baseline circulating concentrations of glucocorticoids, has a heritable component and exhibits intra- and inter-individual variation. However, the mechanisms behind this variability and whether they might be responsible for adaptation to different environments are not known. Variation in DNA methylation can be a potential mechanism that mediates environmental effects on the stress response, as early-life stressors increase glucocorticoid concentrations and change adult phenotype. We used an inter- and intra-environmental cross-foster experiment to analyse the contribution of DNA methylation to early-life phenotypic variation. We found that at hatching, urban house wren (Troglodytes aedon) offspring had higher methylation frequencies compared with their rural counterparts. We also observed age-related patterns in offspring methylation, indicating the developmental effects of the rearing environment on methylation. At fledgling, differential methylation analyses showed that cellular respiration genes were differentially methylated in broods of different origins and behavioural and metabolism genes were differentially methylated in broods of different rearing environments. Lastly, hyper-methylation of a single gene (CNTNAP2) is associated with decreased glucocorticoid levels and the rearing environment. These differential methylation patterns linked to a specific physiological phenotype suggest that DNA methylation may be a mechanism by which individuals adjust to novel environments during their lifespan. Characterizing genetic and environmental influences on methylation is critical for understanding the role of epigenetic mechanisms in evolutionary adaptation.

The giant diploid faba genome unlocks variation in a global protein crop

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

Re-analysis of publicly available methylomes using signal detection yields new information

Cytosine methylation is an epigenetic mark that participates in regulation of gene expression and chromatin stability in plants. Advancements in whole genome sequencing technologies have enabled investigation of methylome dynamics under different conditions. However, the computational methods for analyzing bisulfite sequence data have not been unified. Contention remains in the correlation of differentially methylated positions with the investigated treatment and exclusion of noise, inherent to these stochastic datasets. The prevalent approaches apply Fisher's exact test, logistic, or beta regression, followed by an arbitrary cut-off for differences in methylation levels. A different strategy, the MethylIT pipeline, utilizes signal detection to determine cut-off based on a fitted generalized gamma probability distribution of methylation divergence. Re-analysis of publicly available BS-seq data from two epigenetic studies in Arabidopsis and applying MethylIT revealed additional, previously unreported results. Methylome repatterning in response to phosphate starvation was confirmed to be tissue-specific and included phosphate assimilation genes in addition to sulfate metabolism genes not implicated in the original study. During seed germination plants undergo major methylome reprogramming and use of MethylIT allowed us to identify stage-specific gene networks. We surmise from these comparative studies that robust methylome experiments must account for data stochasticity to achieve meaningful functional analyses.

Transcriptional and post-transcriptional controls for tuning gene expression in plants

Plant biotechnologists seek to modify plants through genetic reprogramming, but our ability to precisely control gene expression in plants is still limited. Here, we review transcription and translation in the model plants Arabidopsis thaliana and Nicotiana benthamiana with an eye toward control points that may be used to predictably modify gene expression. We highlight differences in gene expression requirements between these plants and other species, and discuss the ways in which our understanding of gene expression has been used to engineer plants. This review is intended to serve as a resource for plant scientists looking to achieve precise control over gene expression.

Contribution of RdDM to the ecotype-specific differential methylation on conserved as well as highly variable regions between Arabidopsis ecotypes

BACKGROUND

Several studies showed genome-wide DNA methylation during Arabidopsis embryogenesis and germination. Although it has been known that the change of DNA methylation mainly occurs at CHH context mediated by small RNA-directed DNA methylation pathway during seed ripening and germination, the causality of the methylation difference exhibited in natural Arabidopsis ecotypes has not been thoroughly studied.

RESULTS

In this study we compared DNA methylation difference using comparative pairwise multi-omics dynamics in Columbia-0 (Col) and Cape Verde Island (Cvi) ecotypes. Arabidopsis genome was divided into two regions, common regions in both ecotypes and Col-specific regions, depending on the reads mapping of whole genome bisulfite sequencing libraries from both ecotypes. Ecotype comparison was conducted within common regions and the levels of DNA methylation on common regions and Col-specific regions were also compared. we confirmed transcriptome were relatively dynamic in stage-wise whereas the DNA methylome and small RNAome were more ecotype-dependent. While the global CG methylation remains steady during maturation and germination, we found genic CG methylation differs the most between the two accessions. We also found that ecotype-specific differentially methylated regions (eDMR) are positively correlated with ecotype-specifically expressed 24-nt small RNA clusters. In addition, we discovered that Col-specific regions enriched with transposable elements (TEs) and structural variants that tend to become hypermethylated, and TEs in Col-specific regions were longer in size, more pericentromeric, and more hypermethylated than those in the common regions. Through the analysis of RdDM machinery mutants, we confirmed methylation on Col-specific region as well as on eDMRs in common region are contributed by RdDM pathway. Lastly, we demonstrated that highly variable sequences between ecotypes (HOT regions) were also affected by RdDM-mediated regulation.

CONCLUSIONS

Through ecotype comparison, we revealed differences and similarities of their transcriptome, methylome and small RNAome both in global and local regions. We validated the contribution of RdDM causing differential methylation of common regions. Hypermethylated ecotype-specific regions contributed by RNA-directed DNA methylation pathway largely depend on the presence of TEs and copy-gain structural variations. These ecotype-specific regions are frequently associated with HOT regions, providing evolutionary insights into the epigenome dynamics within a species.

Unravelling Differential DNA Methylation Patterns in Genotype Dependent Manner under Salinity Stress Response in Chickpea

DNA methylation is one of the epigenetic mechanisms that govern gene regulation in response to abiotic stress in plants. Here, we analyzed the role of epigenetic variations by exploring global DNA methylation and integrating it with differential gene expression in response to salinity stress in tolerant and sensitive chickpea genotypes. Genome-wide DNA methylation profiles showed higher CG methylation in the gene body regions and higher CHH methylation in the TE body regions. The analysis of differentially methylated regions (DMRs) suggested more hyper-methylation in response to stress in the tolerant genotype compared to the sensitive genotype. We observed higher enrichment of CG DMRs in genes and CHH DMRs in transposable elements (TEs). A positive correlation of gene expression with CG gene body methylation was observed. The enrichment analysis of DMR-associated differentially expressed genes revealed they are involved in biological processes, such as lateral root development, transmembrane transporter activity, GTPase activity, and regulation of gene expression. Further, a high correlation of CG methylation with CHG and CHH methylation under salinity stress was revealed, suggesting crosstalk among the methylation contexts. Further, we observed small RNA-mediated CHH hypermethylation in TEs. Overall, the interplay between DNA methylation, small RNAs, and gene expression provides new insights into the regulatory mechanism underlying salinity stress response in chickpeas.

The creation of autotetraploid provides insights into critical features of DNA methylome changes after genome doubling in water spinach (Ipomoea aquatica Forsk)

Water spinach (Ipomoea aquatica Forsk) is an essential green leafy vegetable in Asia. In this study, we induced autotetraploid water spinach by colchicine. Furthermore, DNA methylation and transcriptome of tetraploid and diploid were compared using Whole Genome Bisulfite Sequencing (WGBS) and RNA-sequencing techniques. Autotetraploid water spinach was created for the first time. Compared with the diploid parent, autotetraploid water spinach had wider leaves, thicker petioles and stems, thicker and shorter adventitious roots, longer stomas, and larger parenchyma cells. The whole genome methylation level of the autotetraploid was slightly higher than that of the diploid. Compared with the diploid, 12281 Differentially Methylated Regions (DMRs)were found in the autotetraploid, including 2356 hypermethylated and 1310 hypomethylated genes, mainly enriched in 'Arginine and Proline metabolism', 'beta - Alanine metabolism', 'Plant homone signal translation', 'Ribome', and 'Plant - pathgen interaction' pathways. Correlation analysis of transcriptome and DNA methylation data showed that 121 differentially expressed genes undergone differential methylation, related to four pathways 'Other types of O-glycan biosynthesis', 'Terpenoid backbone biosynthesis', 'Biosynthesis of secondary metabolites', and 'Metabolic paths'. This work obtained important autotetraploid resources of water spinach and revealed the genomic DNA methylation changes after genome doubling, being helpful for further studying the molecular mechanism of variations caused by polyploids of the Ipomoea genus.

DNA methylation in clonal duckweed (Lemna minor L.) lineages reflects current and historical environmental exposures

Environmentally induced DNA methylation variants may mediate gene expression responses to environmental changes. If such induced variants are transgenerationally stable, there is potential for expression responses to persist over multiple generations. Our current knowledge in plants, however, is almost exclusively based on studies conducted in sexually reproducing species where the majority of DNA methylation changes are subject to resetting in germlines, limiting the potential for transgenerational epigenetics stress memory. Asexual reproduction circumvents germlines, and may therefore be more conducive to long-term inheritance of epigenetic marks. Taking advantage of the rapid clonal reproduction of the common duckweed Lemna minor, we hypothesize that long-term, transgenerational stress memory from exposure to high temperature can be detected in DNA methylation profiles. Using a reduced representation bisulphite sequencing approach (epiGBS), we show that temperature stress induces DNA hypermethylation at many CG and CHG cytosine contexts but not CHH. Additionally, differential methylation in CHG context that was observed was still detected in a subset of cytosines, even after 3-12 generations of culturing in a common environment. This demonstrates a memory effect of stress reflected in the methylome and that persists over multiple clonal generations. Structural annotation revealed that this memory effect in CHG methylation was enriched in transposable elements. The observed epigenetic stress memory is probably caused by stable transgenerational persistence of temperature-induced DNA methylation variants across clonal generations. To the extent that such epigenetic memory has functional consequences for gene expression and phenotypes, this result suggests potential for long-term modulation of stress responses in asexual plants.

Integrated transcriptome and methylome analyses reveal the molecular regulation of drought stress in wild strawberry (Fragaria nilgerrensis)

BACKGROUND

Fragaria nilgerrensis, which is a diploid wild strawberry with excellent drought-resistance, would provide useful candidate genes for improving drought resistance of cultivated strawberry. So far, its molecular regulatory networks involved in drought stress are unclear. We therefore investigated the drought response regulatory networks of F. nilgerrensis based on the integrated analysis of DNA methylation, transcriptome and physiological traits during four time points under drought stress.  RESULTS: The most differentially expressed genes and the physiological changes were found at 8 days (T8) compared with 0 day (T0, control). Methylome analysis revealed slight dynamic changes in genome-wide mC levels under drought conditions, while the most hypomethylated and hypermethylated regions were identified at T4 and T8. Association analysis of the methylome and transcriptome revealed that unexpressed genes exhibited expected hypermethylation levels in mCHG and mCHH contexts, and highly expressed genes exhibited corresponding hypomethylation levels in the gene body, but mCG contexts showed the opposite trend. Then, 835 differentially methylated and expressed genes were identified and grouped into four clustering patterns to characterize their functions. The genes with either negative or positive correlation between methylation and gene expression were mainly associated with kinases, Reactive Oxygen Species (ROS) synthesis, scavenging, and the abscisic acid (ABA) signal pathway. Consistently, weighted gene co-expression network analysis (WGCNA) revealed Hub genes including NCED, CYP707A2, PP2Cs and others that play important roles in the ABA signaling pathway.

CONCLUSION

F. nilgerrensis drought is dominated by ABA-dependent pathways, possibly accompanied by ABA-independent crosstalk. DNA methylation may affect gene expression, but their correlation was more subtle and multiple types of association exist. Maintaining the balance between ROS regeneration and scavenging is an important factor in drought resistance in F. nilgerrensis. These results deepen our understanding of drought resistance and its application in breeding in strawberry plants.

Gene expression and DNA methylation altering lead to the high oil content in wild allotetraploid peanut (A. monticola)

INTRODUCTION

The wild allotetraploid peanut Arachis monticola contains a higher oil content than the cultivated allotetraploid Arachis hypogaea. Besides the fact that increasing oil content is the most important peanut breeding objective, a proper understanding of its molecular mechanism controlling oil accumulation is still lacking.

METHODS

We investigated this aspect by performing comparative transcriptomics from developing seeds between three wild and five cultivated peanut varieties.

RESULTS

The analyses not only showed species-specific grouping transcriptional profiles but also detected two gene clusters with divergent expression patterns between two species enriched in lipid metabolism. Further analysis revealed that expression alteration of lipid metabolic genes with co-expressed transcription factors in wild peanut led to enhanced activity of oil biogenesis and retarded the rate of lipid degradation. In addition, bisulfite sequencing was conducted to characterize the variation of DNA methylation between wild allotetraploid (245, WH 10025) and cultivated allotetraploid (Z16, Zhh 7720) genotypes. CG and CHG context methylation was found to antagonistically correlate with gene expression during seed development. Differentially methylated region analysis and transgenic assay further illustrated that variations of DNA methylation between wild and cultivated peanuts could affect the oil content via altering the expression of peroxisomal acyl transporter protein (Araip.H6S1B).

DISCUSSION

From the results, we deduced that DNA methylation may negatively regulate lipid metabolic genes and transcription factors to subtly affect oil accumulation divergence between wild and cultivated peanuts. Our work provided the first glimpse on the regulatory mechanism of gene expression altering for oil accumulation in wild peanut and gene resources for future breeding applications.

Systematic genome-wide and expression analysis of RNA-directed DNA methylation pathway genes in grapes predicts their involvement in multiple biological processes

RNA-directed DNA methylation (RdDM) is an important epigenetic pathway in plants and mediates transcriptional silencing by siRNAs. Different gene families have role in the regulation of the RdDM pathway and there is a lack of information about these gene families in the grapes (Vitis vinifera L.). Here, we mentioned the genome-wide identification, bioinformatics analysis, evolutionary history, and expression profiling of VvRdDM pathway genes against various stresses, hormonal treatments as well as in different organs. Sixty VvRdDM genes belonging to fourteen different families were identified. All the genes were unevenly distributed and chromosome 4 contained the highest number of genes (7). Most of the genes showed similar exon-intron and motif distribution patterns within the same subfamilies. Out of 14 families, only members of 4 families underwent duplication events during the evolutionary process and 50% of members of the AGO family are the result of duplication events. Based on Ka/Ks ratio all duplicated gene pairs have a negative mode of selection. VvRdDM pathway genes showed differential spatiotemporal expression patterns against different hormone and stress treatments. Further, with multiple transcriptome analysis, some VvRdDM genes showed a broad spectrum of high expression in different organs at various stages, and VvRdDM genes also displayed different expression in seeded and seedless cultivars during different phases of seed development. This proposed that VvRdDM genes may play multiple roles in grape growth and development, especially in seed development. qRT-PCR analysis of selected genes further verified the critical roles of RdDM genes in multiple biological processes, especially in seed development/ovule abortion i.e., VvIDN2a, VvDRD1a, VvRDR1a, and VvRDR6. Our study provides detailed information about VvRdDM genes in perspective of gene structure and evolution, as well as expression pattern against different stress, hormones and in different plants parts. It provides new candidate gene resources for further functional characterization and molecular breeding of grapes.

Genome-wide DNA methylation reveals potential epigenetic mechanism of age-dependent viral susceptibility in grass carp

Background

Grass carp are an important farmed fish in China that are infected by many pathogens, especially grass carp reovirus (GCRV). Notably, grass carp showed age-dependent susceptibility to GCRV; that is, grass carp not older than one year were sensitive to GCRV, while those over three years old were resistant to this virus. However, the underlying mechanism remains unclear. Herein, whole genome-wide DNA methylation and gene expression variations between susceptible five-month-old (FMO) and resistant three-year-old (TYO) grass carp were investigated aiming to uncover potential epigenetic mechanisms.

Results

Colorimetric quantification revealed that the global methylation level in TYO fish was higher than that in FMO fish. Whole-genome bisulfite sequencing (WGBS) of the two groups revealed 6214 differentially methylated regions (DMRs) and 4052 differentially methylated genes (DMGs), with most DMRs and DMGs showing hypermethylation patterns in TYO fish. Correlation analysis revealed that DNA hypomethylation in promoter regions and DNA hypermethylation in gene body regions were associated with gene expression. Enrichment analysis revealed that promoter hypo-DMGs in TYO fish were significantly enriched in typical immune response pathways, whereas gene body hyper-DMGs in TYO fish were significantly enriched in terms related to RNA transcription, biosynthesis, and energy production. RNA-seq analysis of the corresponding samples indicated that most of the genes in the above terms were upregulated in TYO fish. Moreover, gene function analysis revealed that the two genes involved in energy metabolism displayed antiviral effects.

Conclusions

Collectively, these results revealed genome-wide variations in DNA methylation between grass carp of different ages. DNA methylation and gene expression variations in genes involved in immune response, biosynthesis, and energy production may contribute to age-dependent susceptibility to GCRV in grass carp. Our results provide important information for disease-resistant breeding programs for grass carp and may also benefit research on age-dependent diseases in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12979-022-00285-w.

Gene body methylation in cancer: molecular mechanisms and clinical applications

DNA methylation is an important epigenetic mechanism that regulates gene expression. To date, most DNA methylation studies have focussed on CpG islands in the gene promoter region, and the mechanism of methylation and the regulation of gene expression after methylation have been clearly elucidated. However, genome-wide methylation studies have shown that DNA methylation is widespread not only in promoters but also in gene bodies. Gene body methylation is widely involved in the expression regulation of many genes and is closely related to the occurrence and progression of malignant tumours. This review focusses on the formation of gene body methylation patterns, its regulation of transcription, and its relationship with tumours, providing clues to explore the mechanism of gene body methylation in regulating gene transcription and its significance and application in the field of oncology.

Evidence of chromatin and transcriptional dynamics for cold development in peach flower bud

In temperate zones, fruit trees regulate their annual growth cycle to seasonal environmental changes. During the cold season, growth is limited by both environmental and genetic factors. After the exposure to low temperature and fulfillment of chilling requirements, mild temperatures promote the growth and flowering. However, an insufficient chilling exposure may lead to nonuniform blooming, with a negative impact on fruit set. To gain insights into flower development in the fruit tree buds, peach is an interesting model, the flower and vegetative bud being distinct organs. To understand how flower bud development is regulated, we integrated cytological observations and epigenetic and chromatin genome-wide data with transcriptional changes to identify the main regulatory factors involved in flower development during chilling accumulation. We demonstrated that growth cessation does not occur in peach flower buds during chilling accumulation, but that there are changes in transcript abundance of key genes of hormone metabolism and flower bud development, distribution of histone modifications (H3K4me3 and H3K27me3) and DNA methylation. Altogether, our findings indicate that during the cold season the flower bud is in a nondormant state and that the chilling experience allows flower differentiation to be completed.