Evolution of DNA methylation patterns in the Brassicaceae is driven by differences in genome organization

Almond Grafting for Plum Pox Virus Resistance Triggers Significant Transcriptomic and Epigenetic Shifts in Peaches

Sharka disease, caused by the plum pox virus (PPV), negatively impacts stone fruit production, resulting in economic losses. It has been demonstrated that grafting the almond (Prunus dulcis (Miller) D.A. Webb) variety 'Garrigues' into susceptible peach (Prunus persica (L.) Batsch) rootstocks can result in PPV resistance. The molecular circuits related to grafting in Prunus species, however, have not been fully investigated. In this study, susceptible peach rootstocks 'GF305' were either heterografted with 'Garrigues' almond or homografted with the same cultivar. Peach samples were collected at two stages of scion development, with ungrafted plants utilized as controls. Profiles of transcripts, small RNAs (sRNAs), and DNA methylation were obtained and analyzed on a genome-wide scale. Homografting and heterografting significantly altered the transcriptome and methylome of peach rootstocks, with these modifications being more pronounced during the early stages of scion development. The profiles of sRNAs were significantly more impacted when almonds were used as a scion as opposed to peaches, likely due to the transmission of PPV-unrelated viral sequences. Gene expression differences resulting from DNA methylation alterations are more thoroughly documented at the promoter sequences of genes than within their bodies. This study suggests that the 'Garrigues' almond variety triggers a complex defense response in the peach rootstock, potentially involving the interplay of epigenetic modifications and small RNA-mediated priming of antiviral defenses, which ultimately may contribute to PPV resistance.

Genome-wide DNA methylation and their transgenerational pattern differ in Arabidopsis thaliana populations originated along the elevation of West Himalaya

Methylation at 5' cytosine of DNA molecule is an important epigenetic mark. It is known to play critical role in adaptation of organisms under different biotic and abiotic stressors via modulating gene expression and/or chromatin architecture. Plant populations evolved under variable climatic conditions may have evolved different epigenetic marks including DNA methylation. Here we, describe the genome-wide DNA methylation pattern under native field, F1 and F6 generation followed by their association with phenotypes, climate and global gene expression in the three Arabidopsis thaliana populations originated at different elevation ranges of Indian West Himalaya. We show that the global methyl cytosine (mC) content is more or less similar in the three populations but differ in their distribution across genome. There was an increase in differential methylation between the populations as elevation increased. The methylation divergence was the highest between the low and the high elevation populations. The high elevation populations were hypo-methylated than the low elevation population. The methylation in the genes was associated with population specific phenotypes and climate of the region. The genes which were differentially methylated as well as differentially expressed between the low and high elevation populations were mostly related to abiotic stresses. When grown under controlled condition, there was gain of differential methylation over native condition and the maximum percent changes was observed in CHH-sequence context. Further ~ 99.8% methylated cytosines were stably passed on from F1 to F6 generation. Overall, our data suggest that high elevation population is epigenetically more plastic under changing environmental condition.Background Arabidopsis thaliana is the model plant species and has been extensively studied to understand plants life processes. There are numerous reports on its origin, demography, evolution, epigenomes and adaptation etc. however, Indian populations of Arabidopsis thaliana evolved along wide elevation ranging from ~ 700 m amsl to ~ 3400 m amsl not explored yet. Here we, describe the genome-wide DNA methylation pattern under native field, F1 and F6 generation followed by their association with phenotypes, climate and global gene expression in the three Arabidopsis thaliana populations originated at different elevation ranges of Indian West Himalaya.Results In our study we found that total mCs percent was more or less similar in the three populations but differ in their distribution across genome. The proportion of CG-mCs was the highest, followed by CHH-mCs and CHG-mCs in all the three populations. Under native field condition the methylation divergence was more prominent between low and high elevation populations and the high elevation populations were hypo-methylated than the low elevation population. The methylation in the genes was linked to population-specific phenotypes and the regional climate. The genes that showed differential methylation and expression between low and high elevation populations were primarily associated with abiotic stress responses. When grown under controlled condition, there was gain of differential methylation compared to the native condition and the maximum percent changes was observed in CHH-sequence context. Further 99.8% methylated cytosines were stably passed on from F1 to F6 generation.Conclusions The populations of A. thaliana adapted at different climatic conditions were significantly differentially methylated both under native and controlled condition. However, the magnitude and extent of gain or loss of methylation were most significant between the low and the high elevation populations. Overall, our data suggest that high elevation population is epigenetically more plastic under changing environmental condition.

The evolutionary consequences of interactions between the epigenome, the genome and the environment

The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.

DNA Methylation Analysis of Growth Differences between Upright and Inverted Cuttings of Populus yunnanensis

DNA methylation is an important mechanism for epigenetic modifications that have been shown to be associated with responses to plant development. Previous studies found that inverted Populus yunnanensis cuttings were still viable and could develop into complete plants. However, the growth status of inverted cuttings was weaker than that of upright cuttings, and the sprouting time of inverted cuttings was later than that of upright cuttings. There is currently no research on DNA methylation patterns in inverted cuttings of Populus yunnanensis. In this study, we detected genome-wide methylation patterns of stem tips of Populus yunnanensis at the early growth stage and the rapid growth stage by Oxford Nanopore Technologies (ONT) methylation sequencing. We found that the methylation levels of CpG, CHG, CHH, and 6mA were 41.34%, 33.79%, 17.27%, and 12.90%, respectively, in the genome of inverted poplar cuttings, while the methylation levels of the four methylation types were higher in the genome of upright poplar cuttings than in inverted cuttings, 41.90%, 34.57%, 18.09%, and 14.11%, suggesting important roles for DNA methylation in poplar cells. In all comparison groups, CpG-type methylation genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were annotated to pathways associated with carbon metabolism, ribosome biogenesis in eukaryotes, glycolysis/gluconeogenesis, pyruvate metabolism, and mRNA detection pathways, suggesting that different biological processes are activated in upright and inverted cuttings. The results show that methylation genes are commonly present in the poplar genome, but only a few of them are involved in the regulation of expression in the growth and development of inverted cuttings. From this, we screened the DET2 gene for significant differences in methylation levels in upright or inverted cuttings. The DET2 gene is a key gene in the Brassinolide (BRs) synthesis pathway, and BRs have an important influence on the growth and development process of plants. These results provide important clues for studying DNA methylation patterns in P. yunnanensis.

High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Differentiation of Morphological Traits and Genome-Wide Expression Patterns between Rice Subspecies Indica and Japonica

Changes in gene expression patterns can lead to the variation of morphological traits. This phenomenon is particularly evident in recent evolution events such as crop domestication and responses to environmental stress, where alterations in expression levels can efficiently give rise to domesticated syndromes and adaptive phenotypes. Rice (Oryza sativa L.), one of the world's most crucial cereal crops, comprises two morphologically distinct subspecies, Indica and Japonica. To investigate the morphological divergence between these two rice subspecies, this study planted a total of 315 landrace individuals of both Indica and Japonica under identical cultivation conditions. Out of the 16 quantitative traits measured in this study, 12 exhibited significant differences between the subspecies. To determine the genetic divergence between Indica and Japonica at the whole-genome sequence level, we constructed a phylogenetic tree using a resequencing dataset encompassing 95 rice landrace accessions. The samples formed two major groups that neatly corresponded to the two subspecies, Indica and Japonica. Furthermore, neighbor-joining (NJ) trees based on the expression quantity of effectively expressed genes (EEGs) across five different tissues categorized 12 representative samples into two major clades aligning with the two subspecies. These results imply that divergence in genome-wide expression levels undergoes stabilizing selection under non-stressful conditions, with evolutionary trends in expression levels mirroring sequence variation levels. This study further supports the pivotal role of changes in genome-wide expression regulation in the divergence of the two rice subspecies, Indica and Japonica.

Dynamic DNA methylation turnover in gene bodies is associated with enhanced gene expression plasticity in plants

BACKGROUND

In several eukaryotes, DNA methylation occurs within the coding regions of many genes, termed gene body methylation (GbM). Whereas the role of DNA methylation on the silencing of transposons and repetitive DNA is well understood, gene body methylation is not associated with transcriptional repression, and its biological importance remains unclear.

RESULTS

We report a newly discovered type of GbM in plants, which is under constitutive addition and removal by dynamic methylation modifiers in all cells, including the germline. Methylation at Dynamic GbM genes is removed by the DRDD demethylation pathway and added by an unknown source of de novo methylation, most likely the maintenance methyltransferase MET1. We show that the Dynamic GbM state is present at homologous genes across divergent lineages spanning over 100 million years, indicating evolutionary conservation. We demonstrate that Dynamic GbM is tightly associated with the presence of a promoter or regulatory chromatin state within the gene body, in contrast to other gene body methylated genes. We find Dynamic GbM is associated with enhanced gene expression plasticity across development and diverse physiological conditions, whereas stably methylated GbM genes exhibit reduced plasticity. Dynamic GbM genes exhibit reduced dynamic range in drdd mutants, indicating a causal link between DNA demethylation and enhanced gene expression plasticity.

CONCLUSIONS

We propose a new model for GbM in regulating gene expression plasticity, including a novel type of GbM in which increased gene expression plasticity is associated with the activity of DNA methylation writers and erasers and the enrichment of a regulatory chromatin state.

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.

Multi-omic analyses reveal the unique properties of chia (Salvia hispanica) seed metabolism

Chia (Salvia hispanica) is an emerging crop considered a functional food containing important substances with multiple potential applications. However, the molecular basis of some relevant chia traits, such as seed mucilage and polyphenol content, remains to be discovered. This study generates an improved chromosome-level reference of the chia genome, resolving some highly repetitive regions, describing methylation patterns, and refining genome annotation. Transcriptomic analysis shows that seeds exhibit a unique expression pattern compared to other organs and tissues. Thus, a metabolic and proteomic approach is implemented to study seed composition and seed-produced mucilage. The chia genome exhibits a significant expansion in mucilage synthesis genes (compared to Arabidopsis), and gene network analysis reveals potential regulators controlling seed mucilage production. Rosmarinic acid, a compound with enormous therapeutic potential, was classified as the most abundant polyphenol in seeds, and candidate genes for its complex pathway are described. Overall, this study provides important insights into the molecular basis for the unique characteristics of chia seeds.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Genome-wide DNA methylation predicts environmentally driven life history variation in a marine fish

Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation.

The pattern of DNA methylation alteration, and its association with the expression changes of non-coding RNAs and mRNAs in Moso bamboo under abiotic stress

Epigenetic changes play an important role in plant growth and development and in stress response. However, DNA methylation pattern and its relationship with the expression changes of non-coding RNAs and mRNAs of Moso bamboo in response to abiotic stress is still largely unknown. In this work, we used whole-genome bisulfite sequencing in combination with whole-transcriptome sequencing to analyze the DNA methylation and transcription patterns of mRNAs and non-coding RNAs in Moso bamboo under abiotic stresses such as cold, heat, ultraviolet (UV) and salinity. We found that CHH methylation in the promoter region was positively correlated with gene expression, while CHG and CHH methylations in the gene body regions were negatively associated with gene expression. Moreover, CG and CHG methylations in the promoter regions were negatively correlated with the transcript abundance of long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs). Similarly, the methylation levels of three contexts in the genic regions were negatively correlated with the transcript abundance of lncRNAs and miRNAs but positively correlated with that of circRNAs. In addition, we suggested that the reduction of 21-nt and 24-nt small interfering RNA (siRNA) expression tended to increase methylation levels in the genic regions. We found that stress-responsive genes such as CRPK1, HSFB2A and CIPK were differentially methylated and expressed. Our results also proposed that DNA methylation may regulate the expression of the transcription factors (TFs) and plant hormone signalling genes such as IAA9, MYC2 and ERF110 in response to abiotic stress. This study firstly reports the abiotic stress-responsive DNA methylation pattern and its involvement of expression of coding RNAs and non-coding RNAs in Moso bamboo. The results expand the knowledge of epigenetic mechanisms in Moso bamboo under abiotic stress and support in-depth deciphering of the function of specific non-coding RNAs in future studies.

Predictable and stable epimutations induced during clonal plant propagation with embryonic transcription factor

Clonal propagation is frequently used in commercial plant breeding and biotechnology programs because it minimizes genetic variation, yet it is not uncommon to observe clonal plants with stable phenotypic changes, a phenomenon known as somaclonal variation. Several studies have linked epigenetic modifications induced during regeneration with this newly acquired phenotypic variation. However, the factors that determine the extent of somaclonal variation and the molecular changes underpinning this process remain poorly understood. To address this gap in our knowledge, we compared clonally propagated Arabidopsis thaliana plants derived from somatic embryogenesis using two different embryonic transcription factors- RWP-RK DOMAIN-CONTAINING 4 (RKD4) or LEAFY COTYLEDON2 (LEC2) and from two epigenetically distinct founder tissues. We found that both the epi(genetic) status of the explant and the regeneration protocol employed play critical roles in shaping the molecular and phenotypic landscape of clonal plants. Phenotypic variation in regenerated plants can be largely explained by the inheritance of tissue-specific DNA methylation imprints, which are associated with specific transcriptional and metabolic changes in sexual progeny of clonal plants. For instance, regenerants were particularly affected by the inheritance of root-specific epigenetic imprints, which were associated with an increased accumulation of salicylic acid in leaves and accelerated plant senescence. Collectively, our data reveal specific pathways underpinning the phenotypic and molecular variation that arise and accumulate in clonal plant populations.

Genetic and environmental drivers of large-scale epigenetic variation in Thlaspi arvense

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

DNA hypermethylation promotes the flowering of orchardgrass during vernalization

Vernalization, influenced by environmental factors, is an essential process associated with the productivity of temperate crops, during which epigenetic regulation of gene expression plays an important role. Although DNA methylation is one of the major epigenetic mechanisms associated with the control of gene expression, global changes in DNA methylation in the regulation of gene expression during vernalization-induced flowering of temperate plants remain largely undetermined. To characterize vernalization-associated DNA methylation dynamics, we performed whole-genome bisulfite-treated sequencing and transcriptome sequencing in orchardgrass (Dactylis glomerata) during vernalization. The results revealed that increased levels of genome DNA methylation during the early vernalization of orchardgrass were associated with transcriptional changes in DNA methyltransferase and demethylase genes. Upregulated expression of vernalization-related genes during early vernalization was attributable to an increase in mCHH in the promoter regions of these genes. Application of an exogenous DNA methylation accelerator or overexpression of orchardgrass NUCLEAR POLY(A) POLYMERASE (DgPAPS4) promoted earlier flowering, indicating that DNA hypermethylation plays an important role in vernalization-induced flowering. Collectively, our findings revealed that vernalization-induced hypermethylation is responsible for floral primordium initiation and development. These observations provide a theoretical foundation for further studies on the molecular mechanisms underlying the control of vernalization in temperate grasses.

Solving the grand challenge of phenotypic integration: allometry across scales

Phenotypic integration is a concept related to the cascade of trait relationships from the lowest organizational levels, i.e. genes, to the highest, i.e. whole-organism traits. However, the cause-and-effect linkages between traits are notoriously difficult to determine. In particular, we still lack a mathematical framework to model the relationships involved in the integration of phenotypic traits. Here, we argue that allometric models developed in ecology offer testable mathematical equations of trait relationships across scales. We first show that allometric relationships are pervasive in biology at different organizational scales and in different taxa. We then present mechanistic models that explain the origin of allometric relationships. In addition, we emphasized that recent studies showed that natural variation does exist for allometric parameters, suggesting a role for genetic variability, selection and evolution. Consequently, we advocate that it is time to examine the genetic determinism of allometries, as well as to question in more detail the role of genome size in subsequent scaling relationships. More broadly, a possible-but so far neglected-solution to understand phenotypic integration is to examine allometric relationships at different organizational levels (cell, tissue, organ, organism) and in contrasted species.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

Epigenetic reprogramming of H3K27me3 and DNA methylation during leaf-to-callus transition in peach

Plant tissues are capable of developing unorganized cell masses termed calluses in response to the appropriate combination of auxin and cytokinin. Revealing the potential epigenetic mechanisms involved in callus development can improve our understanding of the regeneration process of plant cells, which will be beneficial for overcoming regeneration recalcitrance in peach. In this study, we report on single-base resolution mapping of DNA methylation and reprogramming of the pattern of trimethylation of histone H3 at lysine 27 (H3K27me3) at the genome-wide level during the leaf-to-callus transition in peach. Overall, mCG and mCHH were predominant at the genome-wide level and mCG was predominant in genic regions. H3K27me3 deposition was mainly detected in the gene body and at the TSS site, and GAGA repetitive sequences were prone to recruit H3K27me3 modification. H3K27me3 methylation was negatively correlated with gene expression. In vitro culture of leaf explants was accompanied by DNA hypomethylation and H3K27me3 demethylation, which could activate auxin- and cytokinin-related regulators to induce callus development. The DNA methylation inhibitor 5-azacytidine could significantly increase callus development, while the H3K27me3 demethylase inhibitor GSK-J4 dramatically reduced callus development. These results demonstrate the roles of DNA methylation and H3K27me3 modification in mediating chromatin status during callus development. Our study provides new insights into the epigenetic mechanisms through which differentiated cells acquire proliferative competence to induce callus development in plants.

Gene-body methylation in plants: mechanisms, functions and important implications for understanding evolutionary processes

Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants.

Genome-wide high-resolution mapping of DNA methylation reveals epigenetic variation in the offspring of sexual and asexual propagation in Robinia pseudoacacia

We detected the genome-wide pattern of DNA methylation and its association with gene expression in sexual and asexual progenies of mature Robinia pseudoacacia trees. DNA methylation plays an important role in plant reproduction and development. Although some studies on sexual reproduction have been carried out in model plants, little is known about the dynamic changes in DNA methylation and their effect on gene expression in sexual and asexual progeny of woody plants. Here, through whole-genome bisulfite sequencing, we revealed DNA methylation patterns in the sexual and asexual progenies of mature Robinia pseudoacacia to understand the regulation of gene expression by DNA methylation in juvenile seedlings. An average of 53% CG, 34% CHG and 5% CHH contexts was methylated in the leaves of mature and juvenile individuals. The CHH methylation level of asexually propagated seedlings was significantly lower than that of seed-derived seedlings and mature trees. The intergenic regions had the highest methylation level. Analysis of differentially methylated regions (DMRs) showed that most of them were hypermethylated and located in the gene upstream and introns. A total of 24, 108 and 162 differentially expressed genes containing DMRs were identified in root sprouts (RSs), root cuttings (RCs) and seed-derived seedlings (SSs), respectively, and a large proportion of them showed hypermethylation. In addition, DMRs were enriched within GO subcategories including catalytic activity, metabolic process and cellular process. The results reveal widespread DNA methylation changes between mature plants and their progenies through sexual/asexual reproduction, which provides novel insights into DNA methylation reprogramming and the regulation of gene expression in woody plants.

Insight Between the Epigenetics and Transcription Responding of Cotton Hypocotyl Cellular Elongation Under Salt-Alkaline Stress

Gossypium barbadense is a cultivated cotton not only known for producing superior fiber but also for its salt and alkaline resistance. Here, we used Whole Genome Bisulfite Sequencing (WGBS) technology to map the cytosine methylation of the whole genome of the G. barbadense hypocotyl at single base resolution. The methylation sequencing results showed that the mapping rates of the three samples were 75.32, 77.54, and 77.94%, respectively. In addition, the Bisulfite Sequence (BS) conversion rate was 99.78%. Approximately 71.03, 53.87, and 6.26% of the cytosine were methylated at CG, CHG, and CHH sequence contexts, respectively. A comprehensive analysis of DNA methylation and transcriptome data showed that the methylation level of the promoter region was a positive correlation in the CHH context. Saline-alkaline stress was related to the methylation changes of many genes, transcription factors (TFs) and transposable elements (TEs), respectively. We explored the regulatory mechanism of DNA methylation in response to salt and alkaline stress during cotton hypocotyl elongation. Our data shed light into the relationship of methylation regulation at the germination stage of G. barbadense hypocotyl cell elongation and salt-alkali treatment. The results of this research help understand the early growth regulation mechanism of G. barbadense in response to abiotic stress.

CHH Methylation Islands: A Nonconserved Feature of Grass Genomes That Is Positively Associated with Transposable Elements but Negatively Associated with Gene-Body Methylation

Methylated CHH (mCHH) islands are peaks of CHH methylation that occur primarily upstream to genes. These regions are actively targeted by the methylation machinery, occur at boundaries between heterochromatin and euchromatin, and tend to be near highly expressed genes. Here we took an evolutionary perspective by studying upstream mCHH islands across a sample of eight grass species. Using a statistical approach to define mCHH islands as regions that differ from genome-wide background CHH methylation levels, we demonstrated that mCHH islands are common and associate with 39% of genes, on average. We hypothesized that islands should be more frequent in genomes of large size, because they have more heterochromatin and hence more need for defined boundaries. We found, however, that smaller genomes tended to have a higher proportion of genes associated with 5' mCHH islands. Consistent with previous work suggesting that islands reflect the silencing of the edge of transposable elements (TEs), genes with nearby TEs were more likely to have mCHH islands. However, the presence of mCHH islands was not a function solely of TEs, both because the underlying sequences of islands were often not homologous to TEs and because genic properties also predicted the presence of 5' mCHH islands. These genic properties included length and gene-body methylation (gbM); in fact, in three of eight species, the absence of gbM was a stronger predictor of a 5' mCHH island than TE proximity. In contrast, gene expression level was a positive but weak predictor of the presence of an island. Finally, we assessed whether mCHH islands were evolutionarily conserved by focusing on a set of 2,720 orthologs across the eight species. They were generally not conserved across evolutionary time. Overall, our data establish additional genic properties that are associated with mCHH islands and suggest that they are not just a consequence of the TE silencing machinery.

Epimutations Define a Fast-Ticking Molecular Clock in Plants

Stochastic gains and losses of DNA methylation at CG dinucleotides are a frequent occurrence in plants. These spontaneous 'epimutations' occur at a rate that is 100 000 times higher than the genetic mutation rate, are effectively neutral at the genome-wide scale, and are stably inherited across mitotic and meiotic cell divisions. Mathematical models have been extraordinarily successful at describing how epimutations accumulate in plant genomes over time, making this process one of the most predictable epigenetic phenomena to date. Here, we propose that their high rate and effective neutrality make epimutations a powerful new molecular clock for timing evolutionary events of the recent past and for age dating of long-lived perennials such as trees.

Intertwined evolution of plant epigenomes and genomes

DNA methylation is found across eukaryotes; however, plants have evolved patterns and pathways of DNA methylation that are distinct from animals and fungi. DNA methylation shapes the evolution of genomes through its direct roles in transposon silencing, gene expression, genome stability, and its impact on mutation rates. In return the diversity of DNA methylation across species is shaped by genome sequence evolution. Extensive diversification of key DNA methylation pathways has continued in plants through gene duplication and loss. Meanwhile, frequent movement of transposons has altered local DNA methylation patterns and the genes affected. Only recently has the diversity and evolutionary history of plant DNA methylation become evident with the availability of increasing genomic and epigenomic data. However, much remains unresolved regarding the evolutionary forces that have shaped the dynamics of the complex and intertwined history of plant genome and epigenome evolution.

Gene body methylation is under selection in Arabidopsis thaliana

In plants, mammals and insects, some genes are methylated in the CG dinucleotide context, a phenomenon called gene body methylation (gbM). It has been controversial whether this phenomenon has any functional role. Here, we took advantage of the availability of 876 leaf methylomes in Arabidopsis thaliana to characterize the population frequency of methylation at the gene level and to estimate the site-frequency spectrum of allelic states. Using a population genetics model specifically designed for epigenetic data, we found that genes with ancestral gbM are under significant selection to remain methylated. Conversely, ancestrally unmethylated genes were under selection to remain unmethylated. Repeating the analyses at the level of individual cytosines confirmed these results. Estimated selection coefficients were small, on the order of 4 N_es = 1.4, which is similar to the magnitude of selection acting on codon usage. We also estimated that A. thaliana is losing gbM threefold more rapidly than gaining it, which could be due to a recent reduction in the efficacy of selection after a switch to selfing. Finally, we investigated the potential function of gbM through its link with gene expression. Across genes with polymorphic methylation states, the expression of gene body methylated alleles was consistently and significantly higher than unmethylated alleles. Although it is difficult to disentangle genetic from epigenetic effects, our work suggests that gbM has a small but measurable effect on fitness, perhaps due to its association to a phenotype-like gene expression.

Integrated analysis of DNA methylome and transcriptome reveals epigenetic regulation of CAM photosynthesis in pineapple

BACKGROUND

Crassulacean acid metabolism (CAM) photosynthesis is an important carbon fixation pathway especially in arid environments because it leads to higher water-use efficiency compared to C3 and C4 plants. However, the role of DNA methylation in regulation CAM photosynthesis is not fully understood.

RESULTS

Here, we performed temporal DNA methylome and transcriptome analysis of non-photosynthetic (white base) and photosynthetic (green tip) tissues of pineapple leaf. The DNA methylation patterns and levels in these two tissues were generally similar for the CG and CHG cytosine sequence contexts. However, CHH methylation was reduced in white base leaf tissue compared with green tip tissue across diel time course in both gene and transposon regions. We identified thousands of local differentially methylated regions (DMRs) between green tip and white base at different diel periods. We also showed that thousands of genes that overlapped with DMRs were differentially expressed between white base and green tip leaf tissue across diel time course, including several important CAM pathway-related genes, such as beta-CA, PEPC, PPCK, and MDH.

CONCLUSIONS

Together, these detailed DNA methylome and transcriptome maps provide insight into DNA methylation changes and enhance our understanding of the relationships between DNA methylation and CAM photosynthesis.

Dynamics of DNA Methylation and Its Functions in Plant Growth and Development

Epigenetic modifications in DNA bases and histone proteins play important roles in the regulation of gene expression and genome stability. Chemical modification of DNA base (e.g., addition of a methyl group at the fifth carbon of cytosine residue) switches on/off the gene expression during developmental process and environmental stresses. The dynamics of DNA base methylation depends mainly on the activities of the writer/eraser guided by non-coding RNA (ncRNA) and regulated by the developmental/environmental cues. De novo DNA methylation and active demethylation activities control the methylation level and regulate the gene expression. Identification of ncRNA involved in de novo DNA methylation, increased DNA methylation proteins guiding DNA demethylase, and methylation monitoring sequence that helps maintaining a balance between DNA methylation and demethylation is the recent developments that may resolve some of the enigmas. Such discoveries provide a better understanding of the dynamics/functions of DNA base methylation and epigenetic regulation of growth, development, and stress tolerance in crop plants. Identification of epigenetic pathways in animals, their existence/orthologs in plants, and functional validation might improve future strategies for epigenome editing toward climate-resilient, sustainable agriculture in this era of global climate change. The present review discusses the dynamics of DNA methylation (cytosine/adenine) in plants, its functions in regulating gene expression under abiotic/biotic stresses, developmental processes, and genome stability.

DNA Methylation Silences Exogenous Gene Expression in Transgenic Birch Progeny

The genetic stability of exogenous genes in the progeny of transgenic trees is extremely important in forest breeding; however, it remains largely unclear. We selected transgenic birch (Betula platyphylla) and its hybrid F1 progeny to investigate the expression stability and silencing mechanism of exogenous genes. We found that the exogenous genes of transgenic birch could be transmitted to their offspring through sexual reproduction. The exogenous genes were segregated during genetic transmission. The hybrid progeny of transgenic birch WT1×TP22 (184) and WT1×TP23 (212) showed higher Bgt expression and greater insect resistance than their parents. However, the hybrid progeny of transgenic birch TP23×TP49 (196) showed much lower Bgt expression, which was only 13.5% of the expression in its parents. To elucidate the mechanism underlying the variation in gene expression between the parents and progeny, we analyzed the methylation rates of Bgt in its promoter and coding regions. The hybrid progeny with normally expressed exogenous genes showed much lower methylation rates (0-29%) than the hybrid progeny with silenced exogenous genes (32.35-45.95%). These results suggest that transgene silencing in the progeny is mainly due to DNA methylation at cytosine residues. We further demonstrated that methylation in the promoter region, rather than in the coding region, leads to gene silencing. We also investigated the relative expression levels of three methyltransferase genes: BpCMT, BpDRM, and BpMET. The transgenic birch line 196 with a silenced Gus gene showed, respectively, 2.54, 9.92, and 4.54 times higher expression levels of BpCMT, BpDRM, and BpMET than its parents. These trends are consistent with and corroborate the high methylation levels of exogenous genes in the transgenic birch line 196. Therefore, our study suggests that DNA methylation in the promoter region leads to silencing of exogenous genes in transgenic progeny of birch.

Impact of Transposable Elements on Methylation and Gene Expression across Natural Accessions of Brachypodium distachyon

Transposable elements (TEs) constitute a large fraction of plant genomes and are mostly present in a transcriptionally silent state through repressive epigenetic modifications, such as DNA methylation. TE silencing is believed to influence the regulation of adjacent genes, possibly as DNA methylation spreads away from the TE. Whether this is a general principle or a context-dependent phenomenon is still under debate, pressing for studying the relationship between TEs, DNA methylation, and nearby gene expression in additional plant species. Here, we used the grass Brachypodium distachyon as a model and produced DNA methylation and transcriptome profiles for 11 natural accessions. In contrast to what is observed in Arabidopsis thaliana, we found that TEs have a limited impact on methylation spreading and that only few TE families are associated with a low expression of their adjacent genes. Interestingly, we found that a subset of TE insertion polymorphisms is associated with differential gene expression across accessions. Thus, although not having a global impact on gene expression, distinct TE insertions may contribute to specific gene expression patterns in B. distachyon.

Epigenetics and epigenomics: underlying mechanisms, relevance, and implications in crop improvement

Epigenetics is defined as changes in gene expression that are not associated with changes in DNA sequence but due to the result of methylation of DNA and post-translational modifications to the histones. These epigenetic modifications are known to regulate gene expression by bringing changes in the chromatin state, which underlies plant development and shapes phenotypic plasticity in responses to the environment and internal cues. This review articulates the role of histone modifications and DNA methylation in modulating biotic and abiotic stresses, as well as crop improvement. It also highlights the possibility of engineering epigenomes and epigenome-based predictive models for improving agronomic traits.

The distribution pattern of 5-methylcytosine in rye (Secale L.) chromosomes

The rye (Secale L.) genome is large, and it contains many classes of repetitive sequences. Secale species differ in terms of genome size, heterochromatin content, and global methylation level; however, the organization of individual types of sequences in chromosomes is relatively similar. The content of the abundant subtelomeric heterochromatin fraction in rye do not correlate with the global level of cytosine methylation, hence immunofluorescence detection of 5-methylcytosine (5-mC) distribution in metaphase chromosomes was performed. The distribution patterns of 5-methylcytosine in the chromosomes of Secale species/subspecies were generally similar. 5-methylcytosine signals were dispersed along the entire length of the chromosome arms of all chromosomes, indicating high levels of methylation, especially at retrotransposon sequences. 5-mC signals were absent in the centromeric and telomeric regions, as well as in subtelomeric blocks of constitutive heterochromatin, in each of the taxa studied. Pericentromeric domains were methylated, however, there was a certain level of polymorphism in these areas, as was the case with the nucleolus organizer region. Sequence methylation within the region of the heterochromatin intercalary bands were also demonstrated to be heterogenous. Unexpectedly, there was a lack of methylation in rye subtelomeres, indicating that heterochromatin is a very diverse fraction of chromatin, and its epigenetic regulation or potential influence on adjacent regions can be more complex than has conventionally been thought. Like telomeres and centromeres, subtelomeric heterochromatin can has a specific role, and the absence of 5-mC is required to maintain the heterochromatin state.

Genome-Wide Analysis Reveals Dynamic Epigenomic Differences in Soybean Response to Low-Phosphorus Stress

Low-phosphorus (low-P) stress has a significant limiting effect on crop yield and quality. Although the molecular mechanisms of the transcriptional level responsible for the low-P stress response have been studied in detail, the underlying epigenetic mechanisms in gene regulation remain largely unknown. In this study, we evaluated the changes in DNA methylation, gene expression and small interfering RNAs (siRNAs) abundance genome-wide in response to low-P stress in two representative soybean genotypes with different P-efficiencies. The DNA methylation levels were slightly higher under low-P stress in both genotypes. Integrative methylation and transcription analysis suggested a complex regulatory relationship between DNA methylation and gene expression that may be associated with the type, region, and extent of methylation. Association analysis of low-P-induced differential methylation and gene expression showed that transcriptional alterations of a small part of genes were associated with methylation changes. Dynamic methylation alterations in transposable element (TE) regions in the CHH methylation context correspond with changes in the amount of siRNA under low-P conditions, indicating an important role of siRNAs in modulating TE activity by guiding CHH methylation in TE regions. Together, these results could help to elucidate the epigenetic regulation mechanisms governing the responses of plants to abiotic stresses.

The Seagrass Methylome Is Associated With Variation in Photosynthetic Performance Among Clonal Shoots

Evolutionary theory predicts that clonal organisms are more susceptible to extinction than sexually reproducing organisms, due to low genetic variation and slow rates of evolution. In agreement, conservation management considers genetic variation as the ultimate measure of a population's ability to survive over time. However, clonal plants are among the oldest living organisms on our planet. Here, we test the hypothesis that clonal seagrass meadows display epigenetic variation that complements genetic variation as a source of phenotypic variation. In a clonal meadow of the seagrass Zostera marina, we characterized DNA methylation among 42 shoots. We also sequenced the whole genome of 10 shoots to correlate methylation patterns with photosynthetic performance under exposure to and recovery from 27°C, while controlling for somatic mutations. Here, we show for the first time that clonal seagrass shoots display DNA methylation variation that is independent from underlying genetic variation, and associated with variation in photosynthetic performance under experimental conditions. It remains unknown to what degree this association could be influenced by epigenetic responses to transplantation-related stress, given that the methylomes showed a strong shift under acclimation to laboratory conditions. The lack of untreated control samples in the heat stress experiment did not allow us to distinguish methylome shifts induced by acclimation from such induced by heat stress. Notwithstanding, the co-variation in DNA methylation and photosynthetic performance may be linked via gene expression because methylation patterns varied in functionally relevant genes involved in photosynthesis, and in the repair and prevention of heat-induced protein damage. While genotypic diversity has been shown to enhance stress resilience in seagrass meadows, we suggest that epigenetic variation plays a similar role in meadows dominated by a single genotype. Consequently, conservation management of clonal plants should consider epigenetic variation as indicator of resilience and stability.

Comparative analysis of genome-wide DNA methylation in Neurospora

DNA methylation is an epigenetic mark that plays an important role in genetic regulation in eukaryotes. Major progress has been made in dissecting the molecular pathways that regulate DNA methylation. Yet, little is known about DNA methylation variation over evolutionary time. Here we present an investigation of the variation of DNA methylation and transposable element (TE) content in species of the filamentous ascomycetes Neurospora. We generated genome-wide DNA methylation data at single-base resolution, together with genomic TE content and gene expression data, of 10 individuals representing five closely related Neurospora species. We found that the methylation levels were low (ranging from 1.3% to 2.5%) and varied among the genomes in a species-specific way. Furthermore, we found that the TEs over 400 bp long were targeted by DNA methylation, and in all genomes, high methylation correlated with low GC, confirming a conserved link between DNA methylation and Repeat Induced Point (RIP) mutations in this group of fungi. Both TE content and DNA methylation pattern showed phylogenetic signal, and the species with the highest TE load (N. crassa) also exhibited the highest methylation level per TE. Our results suggest that DNA methylation is an evolvable trait and indicate that the genomes of Neurospora are shaped by an evolutionary arms race between TEs and host defence.

Principles of Epigenetic Homeostasis Shared Between Flowering Plants and Mammals

In diverse eukaryotes, epigenetic information such as DNA methylation is stably propagated over many cell divisions and generations, and can remain the same over thousands or millions of years. However, this stability is the product of dynamic processes that add and remove DNA methylation by specialized enzymatic pathways. The activities of these dynamic pathways must therefore be finely orchestrated in order to ensure that the DNA methylation landscape is maintained with high fidelity - a concept we term epigenetic homeostasis. In this review, we summarize recent insights into epigenetic homeostasis mechanisms in flowering plants and mammals, highlighting analogous mechanisms that have independently evolved to achieve the same goal of stabilizing the epigenetic landscape.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Cell type-specific genome scans of DNA methylation divergence indicate an important role for transposable elements

In population genomics, genetic diversity measures play an important role in genome scans for divergent sites. In population epigenomics, comparable tools are rare although the epigenome can vary at several levels of organization. We propose a model-free, information-theoretic approach, the Jensen-Shannon divergence (JSD), as a flexible diversity index for epigenomic diversity. Here, we demonstrate how JSD uncovers the relationship between genomic features and cell type-specific methylome diversity in Arabidopsis thaliana. However, JSD is applicable to any epigenetic mark and any collection of individuals, tissues, or cells, for example to assess the heterogeneity in healthy organs and tumors.

A blind and independent benchmark study for detecting differeally methylated regions in plants

MOTIVATION

Bisulfite sequencing (BS-seq) is a state-of-the-art technique for investigating methylation of the DNA to gain insights into the epigenetic regulation. Several algorithms have been published for identification of differentially methylated regions (DMRs). However, the performances of the individual methods remain unclear and it is difficult to optimally select an algorithm in application settings.

RESULTS

We analyzed BS-seq data from four plants covering three taxonomic groups. We first characterized the data using multiple summary statistics describing methylation levels, coverage, noise, as well as frequencies, magnitudes and lengths of methylated regions. Then, simulated datasets with most similar characteristics to real experimental data were created. Seven different algorithms (metilene, methylKit, MOABS, DMRcate, Defiant, BSmooth, MethylSig) for DMR identification were applied and their performances were assessed. A blind and independent study design was chosen to reduce bias and to derive practical method selection guidelines. Overall, metilene had superior performance in most settings. Data attributes, such as coverage and spread of the DMR lengths, were found to be useful for selecting the best method for DMR detection. A decision tree to select the optimal approach based on these data attributes is provided. The presented procedure might serve as a general strategy for deriving algorithm selection rules tailored to demands in specific application settings.

AVAILABILITY AND IMPLEMENTATION

Scripts that were used for the analyses and that can be used for prediction of the optimal algorithm are provided at https://github.com/kreutz-lab/DMR-DecisionTree. Simulated and experimental data are available at https://doi.org/10.6084/m9.figshare.11619045.

CONTACT

ckreutz@imbi.uni-freiburg.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Viral Fitness Determines the Magnitude of Transcriptomic and Epigenomic Reprograming of Defense Responses in Plants

Although epigenetic factors may influence the expression of defense genes in plants, their role in antiviral responses and the impact of viral adaptation and evolution in shaping these interactions are still poorly explored. We used two isolates of turnip mosaic potyvirus with varying degrees of adaptation to Arabidopsis thaliana to address these issues. One of the isolates was experimentally evolved in the plant and presented increased load and virulence relative to the ancestral isolate. The magnitude of the transcriptomic responses was larger for the evolved isolate and indicated a role of innate immunity systems triggered by molecular patterns and effectors in the infection process. Several transposable elements located in different chromatin contexts and epigenetic-related genes were also affected. Correspondingly, mutant plants having loss or gain of repressive marks were, respectively, more tolerant and susceptible to turnip mosaic potyvirus, with a more efficient response against the ancestral isolate. In wild-type plants, both isolates induced similar levels of cytosine methylation changes, including in and around transposable elements and stress-related genes. Results collectively suggested that apart from RNA silencing and basal immunity systems, DNA methylation and histone modification pathways may also be required for mounting proper antiviral defenses and that the effectiveness of this type of regulation strongly depends on the degree of viral adaptation to the host.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

Natural variation in DNA methylation homeostasis and the emergence of epialleles

In plants and mammals, DNA methylation plays a critical role in transcriptional silencing by delineating heterochromatin from transcriptionally active euchromatin. A homeostatic balance between heterochromatin and euchromatin is essential to genomic stability. This is evident in many diseases and mutants for heterochromatin maintenance, which are characterized by global losses of DNA methylation coupled with localized ectopic gains of DNA methylation that alter transcription. Furthermore, we have shown that genome-wide methylation patterns in Arabidopsis thaliana are highly stable over generations, with the exception of rare epialleles. However, the extent to which natural variation in the robustness of targeting DNA methylation to heterochromatin exists, and the phenotypic consequences of such variation, remain to be fully explored. Here we describe the finding that heterochromatin and genic DNA methylation are highly variable among 725 A. thaliana accessions. We found that genic DNA methylation is inversely correlated with that in heterochromatin, suggesting that certain methylation pathway(s) may be redirected to genes upon the loss of heterochromatin. This redistribution likely involves a feedback loop involving the DNA methyltransferase, CHROMOMETHYLASE 3 (CMT3), H3K9me2, and histone turnover, as highly expressed, long genes with a high density of CMT3-preferred CWG sites are more likely to be methylated. Importantly, although the presence of CG methylation in genes alone may not affect transcription, genes containing CG methylation are more likely to become methylated at non-CG sites and silenced. These findings are consistent with the hypothesis that natural variation in DNA methylation homeostasis may underlie the evolution of epialleles that alter phenotypes.

Interspecific variation across angiosperms in global DNA methylation: phylogeny, ecology and plant features in tropical and Mediterranean communities

The interspecific range of epigenetic variation and the degree to which differences between angiosperm species are related to geography, evolutionary history, ecological settings or species-specific traits, remain essentially unexplored. Genome-wide global DNA cytosine methylation is a tractable 'epiphenotypic' feature suitable for exploring these relationships. Global cytosine methylation was estimated in 279 species from two distant, ecologically disparate geographical regions: Mediterranean Spain and tropical México. At each region, four distinct plant communities were analyzed. Global methylation spanned a 10-fold range among species (4.8-42.2%). Interspecific differences were related to evolutionary trajectories, as denoted by a strong phylogenetic signal. Genomes of tropical species were on average less methylated than those of Mediterranean ones. Woody plants have genomes with lower methylation than perennial herbs, and genomes of widespread species were less methylated than those of species with restricted geographical distribution. The eight communities studied exhibited broad and overlapping interspecific variances in global cytosine methylation and only two of them differed in average methylation. Altogether, our broad taxonomic survey supported global methylation as a plant 'epiphenotypic' trait largely associated with species evolutionary history, genome size, range size and woodiness. Additional studies are required for better understanding the environmental components underlying local and geographical variation.

Genome invasion by a hypomethylated satellite repeat in Australian crucifer Ballantinia antipoda

Repetitive sequences are ubiquitous components of all eukaryotic genomes. They contribute to genome evolution and the regulation of gene transcription. However, the uncontrolled activity of repetitive sequences can negatively affect genome functions and stability. Therefore, repetitive DNAs are embedded in a highly repressive heterochromatic environment in plant cell nuclei. Here, we analyzed the sequence, composition and the epigenetic makeup of peculiar non-pericentromeric heterochromatic segments in the genome of the Australian crucifer Ballantinia antipoda. By the combination of high throughput sequencing, graph-based clustering and cytogenetics, we found that the heterochromatic segments consist of a mixture of unique sequences and an A-T-rich 174 bp satellite repeat (BaSAT1). BaSAT1 occupies about 10% of the B. antipoda nuclear genome in >250 000 copies. Unlike many other highly repetitive sequences, BaSAT1 repeats are hypomethylated; this contrasts with the normal patterns of DNA methylation in the B. antipoda genome. Detailed analysis of several copies revealed that these non-methylated BaSAT1 repeats were also devoid of heterochromatic histone H3K9me2 methylation. However, the factors decisive for the methylation status of BaSAT1 repeats remain currently unknown. In summary, we show that even highly repetitive sequences can exist as hypomethylated in the plant nuclear genome.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.