Limited Contribution of DNA Methylation Variation to Expression Regulation in Arabidopsis thaliana

Divergent Flowering Time Responses to Increasing Temperatures Are Associated With Transcriptome Plasticity and Epigenetic Modification Differences at FLC Promoter Region of Arabidopsis thaliana

Understanding the genetic, and transcriptomic changes that drive the phenotypic plasticity of fitness traits is a central question in evolutionary biology. In this study, we utilised 152 natural Swedish Arabidopsis thaliana accessions with re-sequenced genomes, transcriptomes and methylomes and measured flowering times (FTs) under two temperature conditions (10°C and 16°C) to address this question. We revealed that the northern accessions exhibited advanced flowering in response to decreased temperature, whereas the southern accessions delayed their flowering, indicating a divergent flowering response. This contrast in flowering responses was associated with the isothermality of their native ranges, which potentially enables the northern accessions to complete their life cycle more rapidly in years with shorter growth seasons. At the transcriptome level, we observed extensive rewiring of gene co-expression networks, with the expression of 25 core genes being associated with the mean FT and its plastic variation. Notably, variations in FLC expression sensitivity between northern and southern accessions were found to be associated with the divergence FT response. Further analysis suggests that FLC expression sensitivity is associated with differences in CG, CHG and CHH methylation at the promoter region. Overall, our study revealed the association between transcriptome plasticity and flowering time plasticity among different accessions, providing evidence for its relevance in ecological adaptation. These findings offer deeper insights into the genetics of rapid responses to environmental changes and ecological adaptation.

Targeted bisulfite sequencing of Scots pine adaptation-related genes

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

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

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

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

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.

Concepts, mechanisms and implications of long-term epigenetic inheritance

Many modes and mechanisms of epigenetic inheritance have been elucidated in eukaryotes. Most of them are relatively short-term, generally not exceeding one or a few organismal generations. However, emerging evidence indicates that one mechanism, cytosine DNA methylation, can mediate epigenetic inheritance over much longer timescales, which are mostly or completely inaccessible in the laboratory. Here we discuss the evidence for, and mechanisms and implications of, such long-term epigenetic inheritance. We argue that compelling evidence supports the long-term epigenetic inheritance of gene body methylation, at least in the model angiosperm Arabidopsis thaliana, and that variation in such methylation can therefore serve as an epigenetic basis for phenotypic variation in natural populations.

CG hypermethylation of the bHLH39 promoter regulates its expression and Fe deficiency responses in tomato roots

Iron (Fe) is an essential micronutrient for all organisms, including plants, whose limited bioavailability restricts plant growth, yield, and nutritional quality. While the transcriptional regulation of plant responses to Fe deficiency have been extensively studied, the contribution of epigenetic modulations, such as DNA methylation, remains poorly understood. Here, we report that treatment with a DNA methylase inhibitor repressed Fe deficiency-induced responses in tomato (Solanum lycopersicum) roots, suggesting the importance of DNA methylation in regulating Fe deficiency responses. Dynamic changes in the DNA methylome in tomato roots responding to short-term (12 hours) and long-term (72 hours) Fe deficiency identified many differentially methylated regions (DMRs) and DMR-associated genes. Most DMRs occurred at CHH sites under short-term Fe deficiency, whereas they were predominant at CG sites following long-term Fe deficiency. Furthermore, no correlation was detected between the changes in DNA methylation levels and the changes in transcript levels of the affected genes under either short-term or long-term treatments. Notably, one exception was CG hypermethylation at the bHLH39 promoter, which was positively correlated with its transcriptional induction. In agreement, we detected lower CG methylation at the bHLH39 promoter and lower bHLH39 expression in MET1-RNA interference lines compared with wild-type seedlings. Virus-induced gene silencing of bHLH39 and luciferase reporter assays revealed that bHLH39 is positively involved in the modulation of Fe homeostasis. Altogether, we propose that dynamic epigenetic DNA methylation in the CG context at the bHLH39 promoter is involved in its transcriptional regulation, thus contributing to the Fe deficiency response of tomato.

Transcriptome analysis indicates the involvement of herbicide-responsive and plant-pathogen interaction pathways in the development of resistance to ACCase inhibitors in Apera spica-venti

BACKGROUND

The continuous use of the herbicides contributes to the emergence of the resistant populations of numerous weed species that are tolerant to multiple herbicides with different modes of action (multiple resistance) which is provided by non-target-site resistance mechanisms. In this study, we addressed the question of rapid acquisition of herbicide resistance to pinoxaden (acetyl CoA carboxylase inhibitor) in Apera spica-venti, which endangers winter cereal crops and has high adaptation capabilities to inhabit many rural locations. To this end, de novo transcriptome of Apera spica-venti was assembled and RNA-sequencing analysis of plants resistant and susceptible to pinoxaden treated with this herbicide was performed.

RESULTS

The obtained data showed that the prime candidate genes responsible for herbicide resistance were those encoding 3-ketoacyl-CoA synthase 12-like, UDP-glycosyltransferases (UGT) including UGT75K6, UGT75E2, UGT83A1-like, and glutathione S-transferases (GSTs) such as GSTU1 and GSTU6. Also, such highly accelerated herbicide resistance emergence may result from the enhanced constitutive expression of a wide range of genes involved in detoxification already before herbicide treatment and may also influence response to biotic stresses, which was assumed by the detection of expression changes in genes encoding defence-related proteins, including receptor kinase-like Xa21. Moreover, alterations in the expression of genes associated with methylation in non-treated herbicide-resistant populations were identified.

CONCLUSION

The obtained results indicated genes that may be involved in herbicide resistance. Moreover, they provide valuable insight into the possible effect of resistance on the weed interaction with the other stresses by indicating pathways associated with both abiotic and biotic stresses. © 2023 Society of Chemical Industry.

Local diversity of drought resistance and resilience in Populus tomentosa correlates with the variation of DNA methylation

Little information is known about DNA methylation variation in shaping environment-specific drought resistance and resilience for tree adaptation. In this study, we leveraged RNA sequencing and whole-genome bisulfite sequencing data to dissect the distinction of epigenetic regulation under drought stress and rewater condition of Populus tomentosa accessions from three geographical regions. We demonstrated low resistance and high resilience for accessions from South. Non-CG methylation levels in promoter regions of Southern accessions were lower than accessions from higher latitudes and negatively regulated gene expression. CHH context methylation was more sensitive to drought stress, and the geographical-specific differentially methylated regions were scarcely changed by environmental fluctuation. We identified 60 conserved hub genes within the co-expression networks that correlate with photosynthetic and stomatal morphological traits. Epigenome-wide association studies and genome-wide association studies of these 60 hub genes revealed the interdependency between genetic and epigenetic variation in GATA9 and LECRK-VIII.2, which was associated with stomatal morphology and chlorophyll content. The natural epigenetic variation in GATA9 was also faithfully transmitted to progenies in two family-based F1 populations. This study indicates a functional relationship of DNA methylation diversity with drought resistance and resilience which offers new insights into plants' local adaptation to a stressful environment.

Canalization of genome-wide transcriptional activity in Arabidopsis thaliana accessions by MET1-dependent CG methylation

BACKGROUND

Despite its conserved role on gene expression and transposable element (TE) silencing, genome-wide CG methylation differs substantially between wild Arabidopsis thaliana accessions.

RESULTS

To test our hypothesis that global reduction of CG methylation would reduce epigenomic, transcriptomic, and phenotypic diversity in A. thaliana accessions, we knock out MET1, which is required for CG methylation, in 18 early-flowering accessions. Homozygous met1 mutants in all accessions suffer from common developmental defects such as dwarfism and delayed flowering, in addition to accession-specific abnormalities in rosette leaf architecture, silique morphology, and fertility. Integrated analysis of genome-wide methylation, chromatin accessibility, and transcriptomes confirms that MET1 inactivation greatly reduces CG methylation and alters chromatin accessibility at thousands of loci. While the effects on TE activation are similarly drastic in all accessions, the quantitative effects on non-TE genes vary greatly. The global expression profiles of accessions become considerably more divergent from each other after genome-wide removal of CG methylation, although a few genes with diverse expression profiles across wild-type accessions tend to become more similar in mutants. Most differentially expressed genes do not exhibit altered chromatin accessibility or CG methylation in cis, suggesting that absence of MET1 can have profound indirect effects on gene expression and that these effects vary substantially between accessions.

CONCLUSIONS

Systematic analysis of MET1 requirement in different A. thaliana accessions reveals a dual role for CG methylation: for many genes, CG methylation appears to canalize expression levels, with methylation masking regulatory divergence. However, for a smaller subset of genes, CG methylation increases expression diversity beyond genetically encoded differences.

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.

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.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

Genome-Wide Profiling of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Potato Response to DON Stress

Potato is an important food crop that occupies lesser area but has greater production than rice and wheat. However, potato production is affected by numerous biotic and abiotic stresses, among which Fusarium dry rot is a disease that has significant effect on potato production, storage, and processing. However, the role of DNA methylation in regulating potato response to Fusarium toxin deoxynivalenol (DON) stress is still not fully understood. In this study, we performed DNA methylome and transcriptome analyses of potato tubers treated with five concentrations of DON. The global DNA methylation levels in potato tubers treated with different concentrations of DON showed significant changes relative to those in the control. In particular, the 20 ng/ml treatment showed the largest decrease in all three contexts of methylation levels, especially CHH contexts in transposon regions. The differentially methylated region (DMR)-associated differentially expressed genes (DEGs) were significantly enriched in resistance-related metabolic pathways, indicating that DNA methylation plays an essential role in potato response to DON stress. Furthermore, we examined lesions on potato tubers infested with Fusarium after treatment. Furthermore, the potato tubers treated with 5 and 35 ng/ml DON had lesions of significantly smaller diameters than those of the control, indicating that DON stress may induce resistance. We speculate that this may be related to epigenetic memory created after DNA methylation changes. The detailed DNA methylome and transcriptome profiles suggest that DNA methylation plays a vital role in potato disease resistance and has great potential for enhancing potato dry rot resistance.

DNA comethylation analysis reveals a functional association between BRCA1 and sperm DNA fragmentation

OBJECTIVES

To identify the DNA comethylation patterns associated with sperm DNA fragmentation (SDF) and to explore the potential associations of hub genes with SDF.

DESIGN

Prospective study.

SETTING

University-affiliated reproductive medicine center.

PATIENT(S)

A total of 300 male patients consulting for couple infertility were recruited from the First Affiliated Hospital of Wenzhou Medical University.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Comethylation network analysis based on the genome-wide methylation profile of spermatozoal DNA from 20 men was performed to identify hub modules and genes involved in SDF. Human spermatozoa were used for targeted bisulfite amplicon sequencing (267 men) or droplet digital polymerase chain reaction (45 men). The potential role of Brca1 in DNA damage was explored in mouse GC2 spermatocyte cells. Oxidative damage to spermatocytes was modeled by incubating GC2 cells with H₂O₂ (25 mM) for 90 minutes.

RESULT(S)

BRCA1 was identified as a hub gene in SDF. Promoter hypermethylation of BRCA1 was observed in those samples with a high DNA fragmentation index (DFI) compared to those with a low DFI. Concomitantly, BRCA1 mRNA expression was lower in samples with a high DFI than with a low DFI. In the GC2 cell model, Brca1 knockdown reduced cell proliferation and increased sensitivity to oxidative stress. Moreover, it increased double-strand breaks and decreased the protein levels of the DNA repair genes MRE11 and RAD51.

CONCLUSION(S)

A prominent cluster of comethylated patterns associated with SDF was identified. BRCA1 may be the hub gene involved in sperm DNA damage.

Methylome and transcriptome analyses of three different degrees of albinism in apple seedlings

BACKGROUND

Leaf colour mutations are universally expressed at the seedling stage and are ideal materials for exploring the chlorophyll biosynthesis pathway, carotenoid metabolism and the flavonoid biosynthesis pathway in plants.

RESULTS

In this research, we analysed the different degrees of albinism in apple (Malus domestica) seedlings, including white-leaf mutants (WM), piebald leaf mutants (PM), light-green leaf mutants (LM) and normal leaves (NL) using bisulfite sequencing (BS-seq) and RNA sequencing (RNA-seq). There were 61,755, 79,824, and 74,899 differentially methylated regions (DMRs) and 7566, 3660, and 3546 differentially expressed genes (DEGs) identified in the WM/NL, PM/NL and LM/NL comparisons, respectively.

CONCLUSION

The analysis of the methylome and transcriptome showed that 9 DMR-associated DEGs were involved in the carotenoid metabolism and flavonoid biosynthesis pathway. The expression of different transcription factors (TFs) may also influence the chlorophyll biosynthesis pathway, carotenoid metabolism and the flavonoid biosynthesis pathway in apple leaf mutants. This study provides a new method for understanding the differences in the formation of apple seedlings with different degrees of albinism.

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.

Phenotypic and Methylome Responses to Salt Stress in Arabidopsis thaliana Natural Accessions

Salt stress threatens plant growth, development and crop yields, and has become a critical global environmental issue. Increasing evidence has suggested that the epigenetic mechanism such as DNA methylation can mediate plant response to salt stress through transcriptional regulation and transposable element (TE) silencing. However, studies exploring genome-wide methylation dynamics under salt stress remain limited, in particular, for studies on multiple genotypes. Here, we adopted four natural accessions of the model species Arabidopsis thaliana and investigated the phenotypic and genome-wide methylation responses to salt stress through whole-genome bisulfite sequencing (WGBS). We found that salt stress significantly changed plant phenotypes, including plant height, rosette diameter, fruit number, and aboveground biomass, and the change in biomass tended to depend on accessions. Methylation analysis revealed that genome-wide methylation patterns depended primarily on accessions, and salt stress caused significant methylation changes in ∼ 0.1% cytosines over the genomes. About 33.5% of these salt-induced differential methylated cytosines (DMCs) were located to transposable elements (TEs). These salt-induced DMCs were mainly hypermethylated and accession-specific. TEs annotated to have DMCs (DMC-TEs) across accessions were found mostly belonged to the superfamily of Gypsy, a type II transposon, indicating a convergent DMC dynamic on TEs across different genetic backgrounds. Moreover, 8.0% of salt-induced DMCs were located in gene bodies and their proximal regulatory regions. These DMCs were also accession-specific, and genes annotated to have DMCs (DMC-genes) appeared to be more accession-specific than DMC-TEs. Intriguingly, both accession-specific DMC-genes and DMC-genes shared by multiple accessions were enriched in similar functions, including methylation, gene silencing, chemical homeostasis, polysaccharide catabolic process, and pathways relating to shifts between vegetative growth and reproduction. These results indicate that, across different genetic backgrounds, methylation changes may have convergent functions in post-transcriptional, physiological, and phenotypic modulation under salt stress. These convergent methylation dynamics across accession may be autonomous from genetic variation or due to convergent genetic changes, which requires further exploration. Our study provides a more comprehensive picture of genome-wide methylation dynamics under salt stress, and highlights the importance of exploring stress response mechanisms from diverse genetic backgrounds.

Potential Role of Domains Rearranged Methyltransferase7 in Starch and Chlorophyll Metabolism to Regulate Leaf Senescence in Tomato

Deoxyribonucleic acid (DNA) methylation is an important epigenetic mark involved in diverse biological processes. Here, we report the critical function of tomato (Solanum lycopersicum) Domains Rearranged Methyltransferase7 (SlDRM7) in plant growth and development, especially in leaf interveinal chlorosis and senescence. Using a hairpin RNA-mediated RNA interference (RNAi), we generated SlDRM7-RNAi lines and observed pleiotropic developmental defects including small and interveinal chlorosis leaves. Combined analyses of whole genome bisulfite sequence (WGBS) and RNA-seq revealed that silencing of SlDRM7 caused alterations in both methylation levels and transcript levels of 289 genes, which are involved in chlorophyll synthesis, photosynthesis, and starch degradation. Furthermore, the photosynthetic capacity decreased in SlDRM7-RNAi lines, consistent with the reduced chlorophyll content and repression of genes involved in chlorophyll biosynthesis, photosystem, and photosynthesis. In contrast, starch granules were highly accumulated in chloroplasts of SlDRM7-RNAi lines and associated with lowered expression of genes in the starch degradation pathway. In addition, SlDRM7 was activated by aging- and dark-induced senescence. Collectively, these results demonstrate that SlDRM7 acts as an epi-regulator to modulate the expression of genes related to starch and chlorophyll metabolism, thereby affecting leaf chlorosis and senescence in tomatoes.

Multi-Omics Analysis Reveals the Mechanism Underlying the Edaphic Adaptation in Wild Barley at Evolution Slope (Tabigha)

At the microsite "Evolution Slope", Tabigha, Israel, wild barley (Hordeum spontaneum) populations adapted to dry Terra Rossa soil, and its derivative abutting wild barley population adapted to moist and fungi-rich Basalt soil. However, the mechanisms underlying the edaphic adaptation remain elusive. Accordingly, whole genome bisulfite sequencing, RNA-sequencing, and metabolome analysis are performed on ten wild barley accessions inhabiting Terra Rossa and Basalt soil. A total of 121 433 differentially methylated regions (DMRs) and 10 478 DMR-genes are identified between the two wild barley populations. DMR-genes in CG context (CG-DMR-genes) are enriched in the pathways related with the fundamental processes, and DMR-genes in CHH context (CHH-DMR-genes) are mainly associated with defense response. Transcriptome and metabolome analysis reveal that the primary and secondary metabolisms are more active in Terra Rossa and Basalt wild barley populations, respectively. Multi-omics analysis indicate that sugar metabolism facilitates the adaptation of wild barley to dry Terra Rossa soil, whereas the enhancement of phenylpropanoid/phenolamide biosynthesis is beneficial for wild barley to inhabit moist and fungi pathogen-rich Basalt soil. The current results make a deep insight into edaphic adaptation of wild barley and provide elite genetic and epigenetic resources for developing barley with high abiotic stress tolerance.

Gene Body Methylation Confers Transcription Robustness in Mangroves During Long-Term Stress Adaptation

Whether induced epigenetic changes contribute to long-term adaptation remains controversial. Recent studies indicate that environmentally cued changes in gene body methylation (gbM) can facilitate acclimatization. However, such changes are often associated with genetic variation and their contribution to long-term stress adaptation remains unclear. Using whole-genome bisulfite sequencing, we examined evolutionary gains and losses of gbM in mangroves that adapted to extreme intertidal environments. We treated mangrove seedlings with salt stress, and investigated expression changes in relation with stress-induced or evolutionarily-acquired gbM changes. Evolution and function of gbM was compared with that of genetic variation. Mangroves gained much more gbM than their terrestrial relatives, mainly through convergent evolution. Genes that convergently gained gbM during evolution are more likely to become methylated in response to salt stress in species where they are normally not marked. Stress-induced and evolutionarily convergent gains of gbM both correlate with reduction in expression variation, conferring genome-wide expression robustness under salt stress. Moreover, convergent gbM evolution is uncoupled with convergent sequence evolution. Our findings suggest that transgenerational inheritance of acquired gbM helps environmental canalization of gene expression, facilitating long-term stress adaptation of mangroves in the face of a severe reduction in genetic diversity.

Beyond the gene: epigenetic and cis-regulatory targets offer new breeding potential for the future

For millennia, natural and artificial selection has combined favourable alleles for desirable traits in crop species. While modern plant breeding has achieved steady increases in crop yields over the last century, on the current trajectory we will simply not meet demand by 2045. Novel breeding strategies and sources of genetic variation will be required to sustainably fill predicted yield gaps and meet new consumer preferences. Here, we highlight that stepping up to meet this grand challenge will increasingly require thinking 'beyond the gene'. Significant progress has been made in understanding the contributions of both epigenetic variation and cis-regulatory variation to plant traits. This non-genic variation has great potential in future breeding, synthetic biology and biotechnology applications.

Genome-wide DNA methylation analysis and biochemical responses provide insights into the initial domestication of halophyte Puccinellia tenuiflora

Puccinellia tenuiflora was domesticated for two years by growing it under non-saline conditions, providing epigenetic and biochemical insights into the initial domestication of extreme halophytes. Some halophytes have economic value as crop species. The domestication of halophytes may offer hope in solving the problem of soil salinization. We domesticated a wild halophyte, Puccinellia tenuiflora, for two years by growing it under non-saline conditions in a greenhouse and used re-sequencing, genome-wide DNA methylation, biochemical, and transcriptome analyses to uncover the mechanisms underlying alterations in the halophyte's tolerance to saline following domestication. Our results showed that non-saline domestication altered the methylation status for a number of genes and transposable elements, resulting in a much higher frequency of hypomethylation than hypermethylation. These modifications to DNA methylation were observed in many critical salinity-tolerance genes, particularly their promoter regions or transcriptional start sites. Twenty-nine potassium channel genes were hypomethylated and three were hypermethylated, suggesting that the DNA methylation status of potassium channel genes was influenced by domestication. The accelerated uptake of potassium is a major salinity tolerance characteristic of P. tenuiflora. We propose that modifications to the DNA methylation of potassium channel genes may be associated with the development of salinity tolerance in this species. By assessing whether non-saline domestication could change the salinity tolerance of P. tenuiflora, we demonstrated that non-saline domesticated plants are less tolerant to saline, which may be attributable to altered sucrose metabolism. DNA methylation and transposable elements may, therefore, be integrated into an environment-sensitive molecular engine that promotes the rapid domestication of P. tenuiflora to enable its use as a crop plant.

Epigenetics and the success of invasive plants

Biological invasions impose ecological and economic problems on a global scale, but also provide extraordinary opportunities for studying contemporary evolution. It is critical to understand the evolutionary processes that underly invasion success in order to successfully manage existing invaders, and to prevent future invasions. As successful invasive species sometimes are suspected to rapidly adjust to their new environments in spite of very low genetic diversity, we are obliged to re-evaluate genomic-level processes that translate into phenotypic diversity. In this paper, we review work that supports the idea that trait variation, within and among invasive populations, can be created through epigenetic or other non-genetic processes, particularly in clonal invaders where somatic changes can persist indefinitely. We consider several processes that have been implicated as adaptive in invasion success, focusing on various forms of 'genomic shock' resulting from exposure to environmental stress, hybridization and whole-genome duplication (polyploidy), and leading to various patterns of gene expression re-programming and epigenetic changes that contribute to phenotypic variation or even novelty. These mechanisms can contribute to transgressive phenotypes, including hybrid vigour and novel traits, and may thus help to understand the huge successes of some plant invaders, especially those that are genetically impoverished. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Salt stress and salt shock differently affect DNA methylation in salt-responsive genes in sugar beet and its wild, halophytic ancestor

Here we determined the impact of salt shock and salt stress on the level of DNA methylation in selected CpG islands localized in promoters or first exons of sixteen salt-responsive genes in beets. Two subspecies differing in salt tolerance were subjected for analysis, a moderately salt-tolerant sugar beet Beta vulgaris ssp. vulgaris cv. Huzar and a halophytic beet, Beta vulgaris ssp. maritima. The CpG island methylation status was determined. All target sequences were hyper- or hypomethylated under salt shock and/or salt stress in one or both beet subspecies. It was revealed that the genomic regions analyzed were highly methylated in both, the salt treated plants and untreated controls. Methylation of the target sequences changed in a salt-dependent manner, being affected by either one or both treatments. Under both shock and stress, the hypomethylation was a predominant response in sugar beet. In Beta vulgaris ssp. maritima, the hypermethylation occurred with higher frequency than hypomethylation, especially under salt stress and in the promoter-located CpG sites. Conversely, the hypomethylation of the promoter-located CpG sites predominated in sugar beet plants subjected to salt stress. This findings suggest that DNA methylation may be involved in salt-tolerance and transcriptomic response to salinity in beets.

Positive selective sweeps of epigenetic mutations regulating specialized metabolites in plants

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana. In comparison with single nucleotide polymorphisms, three types of methylation—methylated CGs (mCGs), mCHGs, and mCHHs—contributed highly to variable gene expression levels among an A. thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A. thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

Herbicide stress-induced DNA methylation changes in two Zea mays inbred lines differing in Roundup® resistance

DNA methylation plays a crucial role in the regulation of gene expression, activity of transposable elements, defense against foreign DNA, and inheritance of specific gene expression patterns. The link between stress exposure and sequence-specific changes in DNA methylation was hypothetical until it was shown that stresses can induce changes in the gene expression through hypomethylation or hypermethylation of DNA. To detect changes in DNA methylation under herbicide stress in two local Zea mays inbred lines exhibiting differential susceptibility to Roundup®, the methylation-sensitive amplified polymorphism (MSAP) technique was used. The overall DNA methylation levels were determined at approximately 60% for both tested lines. The most significant changes were observed for the more sensitive Z. mays line, where 6 h after the herbicide application, a large increase in the level of DNA methylation (attributed to the increase in fully methylated bands (18.65%)) was noted. DNA sequencing revealed that changes in DNA methylation profiles occurred in genes encoding heat shock proteins, membrane proteins, transporters, kinases, lipases, methyltransferases, zinc-finger proteins, cytochromes, and transposons. Herbicide stress-induced changes depended on the Z. mays variety, and the large increase in DNA methylation level in the sensitive line resulted in a lower ability to cope with stress conditions.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

The methylation landscape and its role in domestication and gene regulation in the chicken

Domestication is one of the strongest examples of artificial selection and has produced some of the most extreme within-species phenotypic variation known. In the case of the chicken, it has been hypothesized that DNA methylation may play a mechanistic role in the domestication response. By inter-crossing wild-derived red junglefowl with domestic chickens, we mapped quantitative trait loci for hypothalamic methylation (methQTL), gene expression (eQTL) and behaviour. We find large, stable methylation differences, with 6,179 cis and 2,973 trans methQTL identified. Over 46% of the trans effects were genotypically controlled by five loci, mainly associated with increased methylation in the junglefowl genotype. In a third of eQTL, we find that there is a correlation between gene expression and methylation, while statistical causality analysis reveals multiple instances where methylation is driving gene expression, as well as the reverse. We also show that methylation is correlated with some aspects of behavioural variation in the inter-cross. In conclusion, our data suggest a role for methylation in the regulation of gene expression underlying the domesticated phenotype of the chicken.

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2-year common-garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non-native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.

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.

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.

Double-edged sword: The evolutionary consequences of the epigenetic silencing of transposable elements

Transposable elements (TEs) are genomic parasites that selfishly replicate at the expense of host fitness. Fifty years of evolutionary studies of TEs have concentrated on the deleterious genetic effects of TEs, such as their effects on disrupting genes and regulatory sequences. However, a flurry of recent work suggests that there is another important source of TEs' harmful effects-epigenetic silencing. Host genomes typically silence TEs by the deposition of repressive epigenetic marks. While this silencing reduces the selfish replication of TEs and should benefit hosts, a picture is emerging that the epigenetic silencing of TEs triggers inadvertent spreading of repressive marks to otherwise expressed neighboring genes, ultimately jeopardizing host fitness. In this Review, we provide a long-overdue overview of the recent genome-wide evidence for the presence and prevalence of TEs' epigenetic effects, highlighting both the similarities and differences across mammals, insects, and plants. We lay out the current understanding of the functional and fitness consequences of TEs' epigenetic effects, and propose possible influences of such effects on the evolution of both hosts and TEs themselves. These unique evolutionary consequences indicate that TEs' epigenetic effect is not only a crucial component of TE biology but could also be a significant contributor to genome function and evolution.

Machine learning analyses of methylation profiles uncovers tissue-specific gene expression patterns in wheat

DNA methylation is a mechanism of epigenetic modification in eukaryotic organisms. Generally, methylation within genes promoter inhibits regulatory protein binding and represses transcription, whereas gene body methylation is associated with actively transcribed genes. However, it remains unclear whether there is interaction between methylation levels across genic regions and which site has the biggest impact on gene regulation. We investigated and used the methylation patterns of the bread wheat cultivar Chinese Spring to uncover differentially expressed genes (DEGs) between roots and leaves, using six machine learning algorithms and a deep neural network. As anticipated, genes with higher expression in leaves were mainly involved in photosynthesis and pigment biosynthesis processes whereas genes that were not differentially expressed between roots and leaves were involved in protein processes and membrane structures. Methylation occurred preponderantly (60%) in the CG context, whereas 35 and 5% of methylation occurred in CHG and CHH contexts, respectively. Methylation levels were highly correlated (r = 0.7 to 0.9) between all genic regions, except within the promoter (r = 0.4 to 0.5). Machine learning models gave a high (0.81) prediction accuracy of DEGs. There was a strong correlation (p-value = 9.20×10^-10 ) between all features and gene expression, suggesting that methylation across all genic regions contribute to gene regulation. However, the methylation of the promoter, the CDS and the exon in CG context was the most impactful. Our study provides more insights into the interplay between DNA methylation and gene expression and paves the way for identifying tissue-specific genes using methylation profiles.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Shading Effects on Leaf Color Conversion and Biosynthesis of the Major Secondary Metabolites in the Albino Tea Cultivar "Yujinxiang"

Albino became a novel kind of tea cultivar in China recently. In this study, transcriptome and whole-genome bisulfite sequencing (WGBS) were employed to investigate the shading effects on leaf color conversion and biosynthesis of three major secondary metabolites in the albino tea cultivar "Yujinxiang". The increased leaf chlorophyll level was likely the major cause for shaded leaf greening from young pale or yellow leaf. In comparison with the control, the total catechin level of the shading group was significantly decreased and the abundance of caffeine was markedly increased, while the theanine level was nearly not influenced. Meanwhile, differentially expressed genes (DEGs) enriched in some biological processes and pathways were identified by transcriptome analysis. Furthermore, whole-genome DNA methylation analysis revealed that the global genomic DNA methylation patterns of the shading period were remarkably altered in comparison with the control. In addition, differentially methylated regions (DMRs) and the DMR-related DEG analysis indicated that the DMR-related DEGs were the critical participants in biosynthesis of the major secondary metabolites. These findings suggest that DNA methylation is probably responsible for changes in the contents of the major secondary metabolites in Yujinxiang.

Inhibition of DNA demethylation enhances plant tolerance to cadmium toxicity by improving iron nutrition

Although the alteration of DNA methylation due to abiotic stresses, such as exposure to the toxic metal cadmium (Cd), has been often observed in plants, little is known about whether such epigenetic changes are linked to the ability of plants to adapt to stress. Herein, we report a close linkage between DNA methylation and the adaptational responses in Arabidopsis plants under Cd stress. Exposure to Cd significantly inhibited the expression of three DNA demethylase genes ROS1/DML2/DML3 (RDD) and elevated DNA methylation at the genome-wide level in Col-0 roots. Furthermore, the profile of DNA methylation in Cd-exposed Col-0 roots was similar to that in the roots of rdd triple mutants, which lack RDD, indicating that Cd-induced DNA methylation is associated with the inhibition of RDD. Interestingly, the elevation in DNA methylation in rdd conferred a higher tolerance against Cd stress and improved cellular Fe nutrition in the root tissues. In addition, lowering the Fe supply abolished improved Cd tolerance due to the lack of RDD in rdd. Together, these data suggest that the inhibition of RDD-mediated DNA demethylation in the roots by Cd would in turn enhance plant tolerance to Cd stress by improving Fe nutrition through a feedback mechanism.

Population-level analysis reveals the widespread occurrence and phenotypic consequence of DNA methylation variation not tagged by genetic variation in maize

BACKGROUND

DNA methylation can provide a source of heritable information that is sometimes entirely uncoupled from genetic variation. However, the extent of this uncoupling and the roles of DNA methylation in shaping diversity of both gene expression and phenotypes are hotly debated. Here, we investigate the genetic basis and biological functions of DNA methylation at a population scale in maize.

RESULTS

We perform targeted DNA methylation profiling for a diverse panel of 263 maize inbred genotypes. All genotypes show similar levels of DNA methylation globally, highlighting the importance of DNA methylation in maize development. Nevertheless, we identify more than 16,000 differentially methylated regions (DMRs) that are distributed across the 10 maize chromosomes. Genome-wide association analysis with high-density genetic markers reveals that over 60% of the DMRs are not tagged by SNPs, suggesting the presence of unique information in DMRs. Strong associations between DMRs and the expression of many genes are identified in both the leaf and kernel tissues, pointing to the biological significance of methylation variation. Association analysis with 986 metabolic traits suggests that DNA methylation is associated with phenotypic variation of 156 traits. There are some traits that only show significant associations with DMRs and not with SNPs.

CONCLUSIONS

These results suggest that DNA methylation can provide unique information to explain phenotypic variation in maize.

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.

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.

Natural depletion of histone H1 in sex cells causes DNA demethylation, heterochromatin decondensation and transposon activation

Transposable elements (TEs), the movement of which can damage the genome, are epigenetically silenced in eukaryotes. Intriguingly, TEs are activated in the sperm companion cell - vegetative cell (VC) - of the flowering plant Arabidopsis thaliana. However, the extent and mechanism of this activation are unknown. Here we show that about 100 heterochromatic TEs are activated in VCs, mostly by DEMETER-catalyzed DNA demethylation. We further demonstrate that DEMETER access to some of these TEs is permitted by the natural depletion of linker histone H1 in VCs. Ectopically expressed H1 suppresses TEs in VCs by reducing DNA demethylation and via a methylation-independent mechanism. We demonstrate that H1 is required for heterochromatin condensation in plant cells and show that H1 overexpression creates heterochromatic foci in the VC progenitor cell. Taken together, our results demonstrate that the natural depletion of H1 during male gametogenesis facilitates DEMETER-directed DNA demethylation, heterochromatin relaxation, and TE activation.

Genome-Wide DNA Methylation Profiling in the Lotus (Nelumbo nucifera) Flower Showing its Contribution to the Stamen Petaloid

DNA methylation is a vital epigenetic modification. Methylation has a significant effect on the gene expression influencing the regulation of different physiological processes. Current studies on DNA methylation have been conducted on model plants. Lotus (Nelumbo nucifera) is a basic eudicot exhibiting variations during development, especially in flower formation. DNA methylation profiling was conducted on different flower tissues of lotuses through whole genome bisulfite sequencing (WGBS) to investigate the effects of DNA methylation on its stamen petaloid. A map of methylated cytosines at the single base pair resolution for the lotus was constructed. When the stamen was compared with the stamen petaloid, the DNA methylation exhibited a global decrease. Genome-wide relationship analysis between DNA methylation and gene expression identified 31 different methylation region (DMR)-associated genes, which might play crucial roles in floral organ formation, especially in the stamen petaloid. One out of 31 DMR-associated genes, NNU_05638 was homolog with Plant U-box 33 (PUB33). The DNA methylation status of NNU_05638 promoter was distinct in three floral organs, which was confirmed by traditional bisulfite sequencing. These results provide further insights about the regulation of stamen petaloids at the epigenetic level in lotus.

Dynamic Cytosine DNA Methylation Patterns Associated with mRNA and siRNA Expression Profiles in Alternate Bearing Apple Trees

Cytosine DNA methylation plays an important role in plants: it can mediate gene expression to affect plant growth and development. However, little is known about the potential involvement of cytosine DNA methylation in apple trees as well as in response to alternate bearing. Here, we performed whole-genome bisulfate sequencing to investigate genomic CG, CHG, and CHH methylation patterns, together with their global mRNA accumulation and small RNA expression in "Fuji" apple trees. Results showed that "Fuji" apple trees have a higher CHH methylation than Arabidopsis. Moreover, genomic methylation analysis revealed that CG and CHG methylation were robustly maintained at the early stage of flower induction. Additionally, differentially methylated regions (DMRs), including hypermethylated and hypomethylated DMRs, were also characterized in alternate bearing (AB) apple trees. Intriguingly, the DMRs were enriched in hormones, redox state, and starch and sucrose metabolism, which affected flowering. Further global gene expression evaluation based on methylome analysis revealed a negative correlation between gene body methylation and gene expression. Subsequent small RNA analyses showed that 24-nucleotide small interfering RNAs were activated and maintained in non-CG methylated apple trees. Our whole-genome DNA methylation analysis and RNA and small RNA expression profile construction provide valuable information for future studies.

DNA Methylation Analysis of the Citrullus lanatus Response to Cucumber Green Mottle Mosaic Virus Infection by Whole-Genome Bisulfite Sequencing

DNA methylation is an important epigenetic mark associated with plant immunity, butlittle is known about its roles in viral infection of watermelon. We carried out whole-genomebisulfite sequencing of watermelon leaves at 0 h (ck), 48 h, and 25 days post-inoculation withCucumber green mottle mosaic virus (CGMMV). The number of differentially methylated regions(DMRs) increased during CGMMV infection and 2788 DMR-associated genes (DMGs) werescreened out among three libraries. Most DMRs and DMGs were obtained under the CHH context.These DMGs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG)pathways of secondary biosynthesis and metabolism, plant-pathogen interactions, Toll-likereceptor signaling, and ABC transporters. Additionally, DMGs encoding PR1a, CaMs, calciumbindingprotein, RIN4, BAK1, WRKYs, RBOHs, STKs, and RLPs/RLKs were involved in thewatermelon-CGMMV interaction and signaling. The association between DNA methylation andgene expression was analyzed by RNA-seq and no clear relationship was detected. Moreover,downregulation of genes in the RdDM pathway suggested the reduced RdDM-directed CHHmethylation plays an important role in antiviral defense in watermelon. Our findings providegenome-wide DNA methylation profiles of watermelon and will aid in revealing the molecularmechanism in response to CGMMV infection at the methylation level.

Adaptation and Phenotypic Diversification in Arabidopsis through Loss-of-Function Mutations in Protein-Coding Genes

According to the less-is-more hypothesis, gene loss is an engine for evolutionary change. Loss-of-function (LoF) mutations resulting in the natural knockout of protein-coding genes not only provide information about gene function but also play important roles in adaptation and phenotypic diversification. Although the less-is-more hypothesis was proposed two decades ago, it remains to be explored on a large scale. In this study, we identified 60,819 LoF variants in 1071 Arabidopsis (Arabidopsis thaliana) genomes and found that 34% of Arabidopsis protein-coding genes annotated in the Columbia-0 genome do not have any LoF variants. We found that nucleotide diversity, transposable element density, and gene family size are strongly correlated with the presence of LoF variants. Intriguingly, 0.9% of LoF variants with minor allele frequency larger than 0.5% are associated with climate change. In addition, in the Yangtze River basin population, 1% of genes with LoF mutations were under positive selection, providing important insights into the contribution of LoF mutations to adaptation. In particular, our results demonstrate that LoF mutations shape diverse phenotypic traits. Overall, our results highlight the importance of the LoF variants for the adaptation and phenotypic diversification of plants.

Global Methylomic and Transcriptomic Analyses Reveal the Broad Participation of DNA Methylation in Daily Gene Expression Regulation of Populus trichocarpa

Changes in DNA methylation patterns in different tissues, at various developmental stages, and under environmental stimuli have been investigated in plants. However, the involvement of DNA methylation in daily gene expression regulation and the plant circadian clock have not been reported. Here, we investigated DNA methylomes and mRNA transcriptomes from leaves of P. trichocarpa over 24 h by high-throughput sequencing. We found that approximately 15.63-19.50% of the genomic cytosine positions were methylated in mature poplar leaves, with approximately half being in the form of asymmetric CHH sites. Repetitive sequences and transposable elements (TEs) were heavily methylated, and the hAT and CMC-EnSpm transposons were more heavily methylated than other TEs. High methylation levels were observed upstream and downstream of the transcribed region, medium in exon and intron, low in untranslated region (5'-UTR and 3'-UTR) of genic regions. In total, about 53,689 differentially methylated regions (DMRs) were identified and CHH context was the most abundant type among daily DNA methylation changes. The DMRs overlapped with over one third of the total poplar genes, including plant defense genes. In addition, a positive correlation between expression levels and DNA methylation levels in the gene body region were observed in DMR overlapping genes. About 1,895 circadian regulated genes overlapped with DMRs, with 871 hypermethylated genes with down-regulated expression levels and 881 hypomethylated genes with up-regulated expression levels, indicating the possible regulation of DNA methylation on the daily rhythmic expression of these genes. But rhythmic DNA methylation changes were not detected in any oscillator component genes controlling the plant circadian clock. Our results suggest that DNA methylation participates widely in daily gene expression regulation, but is not the main mechanism modulating the plant circadian clock.

Methylome and transcriptome analyses of apple fruit somatic mutations reveal the difference of red phenotype

Background

Fruit peel colour is an important agronomic trait for fruit quality. Cytosine methylation plays an important role in gene regulation. Although the DNA methylation level of a single gene is important to affect the phenotype of mutation, there are large unknown of difference of the DNA methylation in plant and its mutants.

Results

Using bisulfite sequencing (BS-Seq) and RNA-sequencing (RNA-Seq), we analysed three deep-red-skinned apple (Malus × domestica) mutants (Yanfu 3, YF3; Yanfu 8, YF8; Shannonghong, SNH) and their lighter-skinned parents (Nagafu 2, NF2; Yanfu 3, YF3; Ralls, RL) to explore the different changes in methylation patterns associated with anthocyanin concentrations. We identified 13,405, 13,384, and 10,925 differentially methylated regions (DMRs) and 1987, 956, and 1180 differentially expressed genes (DEGs) in the NF2/YF3, YF3/YF8, and RL/SNH comparisons, respectively. And we found two DMR-associated DEGs involved in the anthocyanin pathway: ANS (MD06G1071600) and F3H (MD05G1074200). These genes exhibited upregulated expression in apple mutants, and differences were observed in the methylation patterns of their promoters. These results suggested that both the regulatory and structural genes may be modified by DNA methylation in the anthocyanin pathway. However, the methylation of structural genes was not the primary reason for expression-level changes. The expression of structural genes may be synergistically regulated by transcription factors and methylation changes. Additionally, the expression of the transcription factor gene MYB114 (MD17G1261100) was upregulated in the deep-red-skinned apple.

Conclusion

Through the analysis of global methylation and transcription, we did not find the correlation between gene expression and the DNA methylation. However, we observed that the upregulated expression of ANS (MD06G1071600) and F3H (MD05G1074200) in apple mutants results in increased anthocyanin contents. Moreover, MYB114 (MD17G1261100) is likely another regulatory gene involved in apple coloration. Our data provided a new understanding about the differences in formation of apple colour mutants.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5499-2) contains supplementary material, which is available to authorized users.

A DNA methylation reader complex that enhances gene transcription

DNA methylation generally functions as a repressive transcriptional signal, but it is also known to activate gene expression. In either case, the downstream factors remain largely unknown. By using comparative interactomics, we isolated proteins in Arabidopsis thaliana that associate with methylated DNA. Two SU(VAR)3-9 homologs, the transcriptional antisilencing factor SUVH1, and SUVH3, were among the methyl reader candidates. SUVH1 and SUVH3 bound methylated DNA in vitro, were associated with euchromatic methylation in vivo, and formed a complex with two DNAJ domain-containing homologs, DNAJ1 and DNAJ2. Ectopic recruitment of DNAJ1 enhanced gene transcription in plants, yeast, and mammals. Thus, the SUVH proteins bind to methylated DNA and recruit the DNAJ proteins to enhance proximal gene expression, thereby counteracting the repressive effects of transposon insertion near genes.

The contrasting response to drought and waterlogging is underpinned by divergent DNA methylation programs associated with transcript accumulation in sesame

DNA methylation is a heritable epigenetic mechanism that participates in gene regulation under abiotic stresses in plants. Sesame (Sesamum indicum) is typically considered a drought-tolerant crop but highly susceptible to waterlogging, probably because of its origin in Africa or India. Understanding DNA methylation patterns under drought and waterlogging conditions can provide insights into the regulatory mechanisms underlying sesame contrasting responses to these abiotic stresses. We combined Methylation-Sensitive Amplified Polymorphism and transcriptome analyses to profile cytosine methylation patterns, transcript accumulation, and their interplay in drought-tolerant and waterlogging-tolerant sesame genotypes. Drought stress strongly induced de novo methylation (DNM) whereas most of the loci were demethylated (DM) during the recovery phase. In contrast, waterlogging stress decreased the level of methylation but during the recovery phase, both DM and DNM were concomitantly deployed. In both stresses, the levels of the differentially accumulated transcripts (DATs) highly correlated with the methylation patterns. We observed that DM was associated with an increase of DAT levels while DNM was correlated with a decrease of DAT levels. Altogether, sesame has divergent epigenetic programs that respond to drought and waterlogging stresses and an interplay among DNA methylation and transcript accumulation may partly modulate the contrasting responses to these stresses.