Loss of Gene Body Methylation in Eutrema salsugineum Is Associated with Reduced Gene Expression

Patterns of methylation and transcriptional plasticity during thermal acclimation in a reef-building coral

Phenotypic plasticity can buffer organisms against short-term environmental fluctuations. For example, previous exposure to increased temperatures can increase thermal tolerance in many species. Prior studies have found that acclimation to higher temperature can influence the magnitude of transcriptional response to subsequent acute thermal stress (hereafter, "transcriptional response modulation"). However, mechanisms mediating this gene expression response and, ultimately, phenotypic plasticity remain largely unknown. Epigenetic modifications are good candidates for modulating transcriptional response, as they broadly correlate with gene expression. Here, we investigate changes in DNA methylation as a possible mechanism controlling shifts in gene expression plasticity and thermal acclimation in the reef-building coral Acropora nana. We find that gene expression response to acute stress is altered in corals acclimated to different temperatures, with many genes exhibiting a dampened response to heat stress in corals pre-conditioned to higher temperatures. At the same time, we observe shifts in methylation during both acclimation (11 days) and acute heat stress (24 h). We observed that the acute heat stress results in shifts in gene-level methylation and elicits an acute transcriptional response in distinct gene sets. Further, acclimation-induced shifts in gene expression plasticity and differential methylation also largely occur in separate sets of genes. Counter to our initial hypothesis no overall correlation between the magnitude of differential methylation and the change in gene expression plasticity. We do find a small but statistically significant overlap in genes exhibiting both dampened expression response and shifts in methylation (14 genes), which could be candidates for further inquiry. Overall, our results suggest transcriptional response modulation occurs independently from methylation changes induced by thermal acclimation.

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.

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.

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.

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

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

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

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

Asymmetric subgenomic chromatin architecture impacts on gene expression in resynthesized and natural allopolyploid Brassica napus

Although asymmetric subgenomic epigenetic modification and gene expression have been revealed in the successful establishment of allopolyploids, the changes in chromatin accessibility and their relationship with epigenetic modifications and gene expression are poorly understood. Here, we synthetically analyzed chromatin accessibility, four epigenetic modifications and gene expression in natural allopolyploid Brassica napus, resynthesized allopolyploid B. napus, and diploid progenitors (B. rapa and B. oleracea). “Chromatin accessibility shock” occurred in both allopolyploidization and natural evolutionary processes, and genic accessible chromatin regions (ACRs) increased after allopolyploidization. ACRs associated with H3K27me3 modifications were more accessible than those with H3K27ac or H3K4me3. Although overall chromatin accessibility may be defined by H3K27me3, the enrichment of H3K4me3 and H3K27ac and depletion of DNA methylation around transcriptional start sites up-regulated gene expression. Moreover, we found that subgenome C_n exhibited higher chromatin accessibility than A_n, which depended on the higher chromatin accessibility of C_n-unique genes but not homologous genes.

Changes in chromatin accessibility occuring during the process of allopolyploidization of Brassica napus are analysed using ATAC and ChIPseq, with differences in asymmetric chromatin accessibility between subgenomes of B. napus investigated.

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.

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.

Asymmetric epigenetic modification and homoeolog expression bias in the establishment and evolution of allopolyploid Brassica napus

This study explores how allopolyploidization reshapes the biased expression and asymmetric epigenetic modification of homoeologous gene pairs, and examines the regulation types and epigenetic basis of expression bias. We analyzed the gene expression and four epigenetic modifications (DNA methylation, H3K4me3, H3K27me3 and H3K27ac) of 29 976 homoeologous gene pairs in resynthesized, natural allopolyploid Brassica napus and an in silico 'hybrid'. We comprehensively elucidated the biased gene expression, asymmetric epigenetic modifications and the generational transmission characteristics of these homoeologous gene pairs in B. napus. We analyzed cis/trans effects and the epigenetic basis of homoeolog expression bias. There was a significant positive correlation between two active histone modifications and biased gene expression. We revealed that parental legacy was the dominant principle in the remodeling of homoeolog expression bias and asymmetric epigenetic modifications in B. napus, and further clarified that this depends on whether there were differences in the expression/epigenetic modifications of gene pairs in parents/progenitors. The maternal genome was dominant in the homoeolog expression bias of resynthesized B. napus, and this phenomenon was attenuated in natural B. napus. Furthermore, cis rather than trans effects were dominant when epigenetic modifications potentially affected biased expression of gene pairs in B. napus.

Repression of CHROMOMETHYLASE 3 prevents epigenetic collateral damage in Arabidopsis

DNA methylation has evolved to silence mutagenic transposable elements (TEs) while typically avoiding the targeting of endogenous genes. Mechanisms that prevent DNA methyltransferases from ectopically methylating genes are expected to be of prime importance during periods of dynamic cell cycle activities including plant embryogenesis. However, virtually nothing is known regarding how DNA methyltransferase activities are precisely regulated during embryogenesis to prevent the induction of potentially deleterious and mitotically stable genic epimutations. Here, we report that microRNA-mediated repression of CHROMOMETHYLASE 3 (CMT3) and the chromatin features that CMT3 prefers help prevent ectopic methylation of thousands of genes during embryogenesis that can persist for weeks afterwards. Our results are also consistent with CMT3-induced ectopic methylation of promoters or bodies of genes undergoing transcriptional activation reducing their expression. Therefore, the repression of CMT3 prevents epigenetic collateral damage on endogenous genes. We also provide a model that may help reconcile conflicting viewpoints regarding the functions of gene-body methylation that occurs in nearly all flowering plants.

Gene body DNA methylation in seagrasses: inter- and intraspecific differences and interaction with transcriptome plasticity under heat stress

The role of DNA methylation and its interaction with gene expression and transcriptome plasticity is poorly understood, and current insight comes mainly from studies in very few model plant species. Here, we study gene body DNA methylation (gbM) and gene expression patterns in ecotypes from contrasting thermal environments of two marine plants with contrasting life history strategies in order to explore the potential role epigenetic mechanisms could play in gene plasticity and responsiveness to heat stress. In silico transcriptome analysis of CpGO/E ratios suggested that the bulk of Posidonia oceanica and Cymodocea nodosa genes possess high levels of intragenic methylation. We also observed a correlation between gbM and gene expression flexibility: genes with low DNA methylation tend to show flexible gene expression and plasticity under changing conditions. Furthermore, the empirical determination of global DNA methylation (5-mC) showed patterns of intra and inter-specific divergence that suggests a link between methylation level and the plants' latitude of origin and life history. Although we cannot discern whether gbM regulates gene expression or vice versa, or if other molecular mechanisms play a role in facilitating transcriptome responsiveness, our findings point to the existence of a relationship between gene responsiveness and gbM patterns in marine plants.

Cross species multi-omics reveals cell wall sequestration and elevated global transcript abundance as mechanisms of boron tolerance in plants

Boron toxicity is a world-wide problem for crops, yet we have a limited understanding of the genetic responses and adaptive mechanisms to this stress in plants. We employed a cross-species comparison between boron stress-sensitive Arabidopsis thaliana and its boron stress-tolerant extremophyte relative Schrenkiella parvula, and a multi-omics approach integrating genomics, transcriptomics, metabolomics and ionomics to assess plant responses and adaptations to boron stress. Schrenkiella parvula maintains lower concentrations of total boron and free boric acid than Arabidopsis when grown with excess boron. Schrenkiella parvula excludes excess boron more efficiently than Arabidopsis, which we propose is partly driven by SpBOR5, a boron transporter that we functionally characterize in this study. Both species use cell walls as a partial sink for excess boron. When accumulated in the cytoplasm, excess boron appears to interrupt RNA metabolism. The extremophyte S. parvula facilitates critical cellular processes while maintaining the pool of ribose-containing compounds that can bind with boric acid. The S. parvula transcriptome is pre-adapted to boron toxicity. It exhibits substantial overlaps with the Arabidopsis boron-stress responsive transcriptome. Cell wall sequestration and increases in global transcript levels under excess boron conditions emerge as key mechanisms for sustaining plant growth under boron toxicity.

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.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha

Marchantia polymorpha has recently become a prime model for cellular, evo-devo, synthetic biological, and evolutionary investigations. We present a pseudomolecule-scale assembly of the M. polymorpha genome, making comparative genome structure analysis and classical genetic mapping approaches feasible. We anchored 88% of the M. polymorpha draft genome to a high-density linkage map resulting in eight pseudomolecules. We found that the overall genome structure of M. polymorpha is in some respects different from that of the model moss Physcomitrella patens. Specifically, genome collinearity between the two bryophyte genomes and vascular plants is limited, suggesting extensive rearrangements since divergence. Furthermore, recombination rates are greatest in the middle of the chromosome arms in M. polymorpha like in most vascular plant genomes, which is in contrast with P. patens where recombination rates are evenly distributed along the chromosomes. Nevertheless, some other properties of the genome are shared with P. patens. As in P. patens, DNA methylation in M. polymorpha is spread evenly along the chromosomes, which is in stark contrast with the angiosperm model Arabidopsis thaliana, where DNA methylation is strongly enriched at the centromeres. Nevertheless, DNA methylation and recombination rate are anticorrelated in all three species. Finally, M. polymorpha and P. patens centromeres are of similar structure and marked by high abundance of retroelements unlike in vascular plants. Taken together, the highly contiguous genome assembly we present opens unexplored avenues for M. polymorpha research by linking the physical and genetic maps, making novel genomic and genetic analyses, including map-based cloning, feasible.

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.

Excess Light Priming in Arabidopsis thaliana Genotypes with Altered DNA Methylomes

Plants must continuously react to the ever-fluctuating nature of their environment. Repeated exposure to stressful conditions can lead to priming, whereby prior encounters heighten a plant's ability to respond to future events. A clear example of priming is provided by the model plant Arabidopsis thaliana (Arabidopsis), in which photosynthetic and photoprotective responses are enhanced following recurring light stress. While there are various post-translational mechanisms underpinning photoprotection, an unresolved question is the relative importance of transcriptional changes toward stress priming and, consequently, the potential contribution from DNA methylation - a heritable chemical modification of DNA capable of influencing gene expression. Here, we systematically investigate the potential molecular underpinnings of physiological priming against recurring excess-light (EL), specifically DNA methylation and transcriptional regulation: the latter having not been examined with respect to EL priming. The capacity for physiological priming of photosynthetic and photoprotective parameters following a recurring EL treatment was not impaired in Arabidopsis mutants with perturbed establishment, maintenance, or removal of DNA methylation. Importantly, no differences in development or basal photoprotective capacity were identified in the mutants that may confound the above result. Little evidence for a causal transcriptional component of physiological priming was identified; in fact, most alterations in primed plants presented as a transcriptional 'dampening' in response to an additional EL exposure, likely a consequence of physiological priming. However, a set of transcripts uniquely regulated in primed plants provide preliminary evidence for a novel transcriptional component of recurring EL priming, independent of physiological changes. Thus, we propose that physiological priming of recurring EL in Arabidopsis occurs independently of DNA methylation; and that the majority of the associated transcriptional alterations are a consequence, not cause, of this physiological priming.

EvoChromo: towards a synthesis of chromatin biology and evolution

Over the past few years, interest in chromatin and its evolution has grown. To further advance these interests, we organized a workshop with the support of The Company of Biologists to debate the current state of knowledge regarding the origin and evolution of chromatin. This workshop led to prospective views on the development of a new field of research that we term 'EvoChromo'. In this short Spotlight article, we define the breadth and expected impact of this new area of scientific inquiry on our understanding of both chromatin and evolution.

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.

Evolutionary and Experimental Loss of Gene Body Methylation and Its Consequence to Gene Expression

In flowering plants, gene body methylation (gbM) is associated with a subset of constitutively expressed genes. It has been proposed that gbM modulates gene expression. Here, we show that there are no consistent and direct differences to expression following the loss of gbM. By comparing expression of gbM genes in Arabidopsis thaliana accessions to orthologous genes in two Eutrema salsugineum genotypes, we identified both positive and negative expression differences associated with gbM loss. However, expression is largely unaffected by gbM loss in E. salsugineum Expression differences between species were within the variation of expression observed within A. thaliana accessions that displayed variation in gbM. Furthermore, experimentally induced loss of gbM did not consistently lead to differences in expression compared to wild type. To date, there is no convincing data to support a direct causal link between the presence/absence of gbM and the modulation of expression in flowering plants.

Single-cell expression noise and gene-body methylation in Arabidopsis thaliana

Gene-body methylation (gbM) refers to an increased level of methylated cytosines specifically in a CG sequence context within genes. gbM is found in plant genes with intermediate expression level, which evolve slowly, and is often broadly conserved across millions of years of evolution. Intriguingly however, some plants lack gbM, and thus it remains unclear whether gbM has a function. In animals, there is support for a role of gbM in reducing erroneous transcription and transcription noise, but so far most studies in plants have tested for an effect of gbM on expression level, not noise. Here, we therefore tested whether gbM was associated with reduced expression noise in Arabidopsis thaliana, using single-cell transcriptome sequencing data from root quiescent centre cells. We find that gbM genes have lower expression noise levels than unmethylated genes. However, an analysis of covariance revealed that, if other genomic features are taken into account, this association disappears. Nonetheless, gbM genes were more consistently expressed across single-cell samples, supporting previous inference that gbM genes are constitutively expressed. Finally, we observed that fewer RNAseq reads map to introns of gbM genes than to introns of unmethylated genes, which indicates that gbM might be involved in reducing erroneous transcription by reducing intron retention.