CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize

Epigenome profiling reveals distinctive regulatory features and cis-regulatory elements in pepper

BACKGROUND

Pepper (Capsicum annuum) is one of the earliest and most widely cultivated vegetable crops worldwide. While the large and complex genome of pepper severely hampered the understanding of its functional genome, it also indicates a rich yet unexplored reservoir of regulatory elements (REs). In fact, variations in the REs represent a major driving force in evolution and domestication in both plants and animals. However, identification of the REs remains difficult especially for plants with complex genomes.

RESULTS

Here, we present a comprehensive epigenomic landscape of Capsicum annuum, Zhangshugang (ST-8), including chromatin accessibility, histone modifications, DNA methylation, and transcriptome. We also develop comparative crosslinked immunoprecipitation mass spectrometry to reveal the proteome associated with certain chromatin features. Through integrated analysis of these epigenetic features, we profile promoters and enhancers involved in development, heat stress and cucumber mosaic virus challenges. We generate stress responsive expression networks composed of potential transcription activators and their target genes. Through population genetics analysis, we demonstrate that some regulatory elements show lower nucleotide diversity compare to other genomic regions during evolution.

CONCLUSIONS

We demonstrate that variations in the REs may contribute to more diversified and agronomically desired phenotypes. Our study provides a foundation not only for studying gene regulation, but also for targeted genetic and epigenetic manipulation for pepper improvement.

Epigenetic variation in maize agronomical traits for breeding and trait improvement

Epigenetics-mediated breeding (epibreeding) involves engineering crop traits and stress responses through the targeted manipulation of key epigenetic features to enhance agricultural productivity. While conventional breeding methods raise concerns about reduced genetic diversity, epibreeding propels crop improvement through epigenetic variations that regulate gene expression, ultimately impacting crop yield. Epigenetic regulation in crops encompasses various modes, including histone modification, DNA modification, RNA modification, non-coding RNA, and chromatin remodeling. This review summarizes the epigenetic mechanisms underlying major agronomic traits in maize and identifies candidate epigenetic landmarks in the maize breeding process. We propose a valuable strategy for improving maize yield through epibreeding, combining CRISPR/Cas-based epigenome editing technology and Synthetic Epigenetics (SynEpi). Finally, we discuss the challenges and opportunities associated with maize trait improvement through epibreeding.

Epigenetic gene regulation in plants and its potential applications in crop improvement

DNA methylation, also known as 5-methylcytosine, is an epigenetic modification that has crucial functions in plant growth, development and adaptation. The cellular DNA methylation level is tightly regulated by the combined action of DNA methyltransferases and demethylases. Protein complexes involved in the targeting and interpretation of DNA methylation have been identified, revealing intriguing roles of methyl-DNA binding proteins and molecular chaperones. Structural studies and in vitro reconstituted enzymatic systems have provided mechanistic insights into RNA-directed DNA methylation, the main pathway catalysing de novo methylation in plants. A better understanding of the regulatory mechanisms will enable locus-specific manipulation of the DNA methylation status. CRISPR-dCas9-based epigenome editing tools are being developed for this goal. Given that DNA methylation patterns can be stably transmitted through meiosis, and that large phenotypic variations can be contributed by epimutations, epigenome editing holds great promise in crop breeding by creating additional phenotypic variability on the same genetic material.

Epigenetic regulation of organ-specific functions in Mikania micrantha and Mikania cordata: insights from DNA methylation and siRNA integration

BACKGROUND

DNA methylation is a crucial epigenetic mechanism that regulates gene expression during plant growth and development. However, the role of DNA methylation in regulating the organ-specific functions of the invasive weed Mikania micrantha remains unknown.

RESULTS

Here, we generated DNA methylation profiles for M. micrantha and a local congeneric species, Mikania cordata, in three vegetative organs (root, stem, and leaf) using whole-genome bisulfite sequencing. The results showed both differences and conservation in methylation levels and patterns between the two species. Combined with transcriptome data, we found that DNA methylation generally inhibited gene expression, with varying effects depending on the genomic region and sequence context (CG, CHG, and CHH). Genes overlapping with differentially methylated regions (DMRs) were more likely to be differentially expressed between organs, and DMR-associated upregulated differentially expressed genes (DEGs) were enriched in organ-specific pathways. A comparison between photosynthetic (leaf) and non-photosynthetic (root) organs of M. micrantha further confirmed the regulatory role of DNA methylation in leaf-specific photosynthesis. Integrating small RNA-Seq data revealed that 24-nt small interfering RNAs (siRNAs) were associated with CHH methylation in gene-rich regions and regulated CHH methylation in the flanking regions of photosynthesis-related genes.

CONCLUSION

This study provides insights into the complex regulatory role of DNA methylation and siRNAs in organ-specific functions and offers valuable information for exploring the invasive characteristics of M. micrantha from an epigenetic perspective.

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.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Can DNA methylation shape climate response in trees?

Woody plants create the ecosystems they occupy and shape their biodiversity. Today's rapidly warming climate threatens these long-lived species by creating new environments in which existing populations become maladapted. Plants show enormous phenotypic diversity in response to environmental change, which can be caused by genotype or epigenetic mechanisms that influence the expression of the underlying DNA sequence. Whether epigenetics can affect ecologically important traits in trees is an important and controversial question. We explore the evidence that DNA methylation can affect gene expression, both directly and indirectly via its interaction with transposable elements (TEs), and subsequently shapes phenotypic variation in natural tree populations. Furthermore, we consider the potential of epigenetic approaches to assist in their conservation management strategies.

PICKLE and HISTONE DEACETYLASE6 coordinately regulate genes and transposable elements in Arabidopsis

Chromatin dynamics play essential roles in transcriptional regulation. The chromodomain helicase DNA-binding domain 3 chromatin remodeler PICKLE (PKL) and HISTONE DEACETYLASE6 (HDA6) are required for transcriptional gene silencing, but their coordinated function in gene repression requires further study. Through a genetic suppressor screen, we found that a point mutation at PKL could partially restore the developmental defects of a weak Polycomb repressive complex 1 (PRC1) mutant (ring1a-2 ring1b-3), in which RING1A expression is suppressed by a T-DNA insertion at the promoter. Compared to ring1a-2 ring1b-3, the expression of RING1A is increased, nucleosome occupancy is reduced, and the histone 3 lysine 9 acetylation (H3K9ac) level is increased at the RING1A locus in the pkl ring1a-2 ring1b-3 triple mutant. HDA6 interacts with PKL and represses RING1A expression similarly to PKL genetically and molecularly in the ring1a-2 ring1b-3 background. Furthermore, we show that PKL and HDA6 suppress the expression of a set of genes and transposable elements (TEs) by increasing nucleosome density and reducing H3K9ac. Genome-wide analysis indicated they possibly coordinately maintain DNA methylation as well. Our findings suggest that PKL and HDA6 function together to reduce H3K9ac and increase nucleosome occupancy, thereby facilitating gene/TE regulation in Arabidopsis (Arabidopsis thaliana).

Dynamic DNA methylation modifications in the cold stress response of cassava

Cassava, a crucial tropical crop, faces challenges from cold stress, necessitating an exploration of its molecular response. Here, we investigated the role of DNA methylation in moderating the response to moderate cold stress (10 °C) in cassava. Using whole-genome bisulfite sequencing, we examined DNA methylation patterns in leaf blades and petioles under control conditions, 5 h, and 48 h of cold stress. Tissue-specific responses were observed, with leaf blades exhibiting subtle changes, while petioles displayed a pronounced decrease in methylation levels under cold stress. We identified cold stress-induced differentially methylated regions (DMRs) that demonstrated both tissue and treatment specificity. Importantly, these DMRs were enriched in genes with altered expression, implying functional relevance. The cold-response transcription factor ERF105 associated with DMRs emerged as a significant and conserved regulator across tissues and treatments. Furthermore, we investigated DNA methylation dynamics in transposable elements, emphasizing the sensitivity of MITEs with bHLH binding motifs to cold stress. These findings provide insights into the epigenetic regulation of response to cold stress in cassava, contributing to an understanding of the molecular mechanisms underlying stress adaptation in this tropical plant.

Natural polymorphisms in ZMET2 encoding a DNA methyltransferase modulate the number of husk layers in maize

DNA methylation affects agronomic traits and the environmental adaptability of crops, but the natural polymorphisms in DNA methylation-related genes and their contributions to phenotypic variation in maize (Zea mays) remain elusive. Here, we show that a polymorphic 10-bp insertion/deletion variant in the 3'UTR of Zea methyltransferase2 (ZMET2) alters its transcript level and accounts for variation in the number of maize husk layers. ZMET2 encodes a chromomethylase and is required for maintaining genome-wide DNA methylation in the CHG sequence context. Disruption of ZMET2 increased the number of husk layers and resulted in thousands of differentially methylated regions, a proportion of which were also distinguishable in natural ZMET2 alleles. Population genetic analyses indicated that ZMET2 was a target of selection and might play a role in the spread of maize from tropical to temperate regions. Our results provide important insights into the natural variation of ZMET2 that confers both global and locus-specific effects on DNA methylation, which contribute to phenotypic diversity in maize.

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Multi-omic analysis of the extension of broccoli quality during storage by folic acid

INTRODUCTION

Folic acid (FA) is a critical metabolite in all living organisms and an important nutritional component of broccoli. Few studies have been conducted on the impact of an exogenous application of FA on the postharvest physiology of fruits and vegetables during storage. In this regard, the mechanism by which an exogenous application of FA extends the postharvest quality of broccoli is unclear.

OBJECTIVE

This study utilized a multicomponent analysis to investigate how an exogenous application of FA effects the postharvest quality of broccoli.

METHODS

Broccoli was soaked in 5 mg/L FA for 10 min and the effect of the treatment on the appearance and nutritional quality of broccoli was evaluated. These data were combined with transcriptomic, metabolomic, and DNA methylation data to provide insight into the potential mechanism by which FA delays senescence.

RESULTS

The FA treatment inhibited the yellowing of broccoli during storage. CHH methylation was identified as the main type of methylation that occurs in broccoli and the FA treatment was found to inhibit DNA methylation, promote the accumulation of endogenous FA and chlorophyl, and inhibit ethylene biosynthesis in stored broccoli. The FA treatment also prevented the formation of off-odors by inhibiting the degradation of glucosinolate.

CONCLUSIONS

FA treatment inhibited the loss of nutrients during the storage of broccoli, delayed its yellowing, and inhibited the generation of off-odors. Our study provides deeper insight into the mechanism by which the postharvest application of FA delays postharvest senescence in broccoli and provides the foundation for further studies of postharvest metabolism in broccoli.

Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

RNA m6A modification facilitates DNA methylation during maize kernel development

N6-methyladenosine (m6A) in mRNA and 5-methylcytosine (5mC) in DNA have critical functions for regulating gene expression and modulating plant growth and development. However, the interplay between m6A and 5mC is an elusive territory and remains unclear mechanistically in plants. We reported an occurrence of crosstalk between m6A and 5mC in maize (Zea mays) via the interaction between mRNA adenosine methylase (ZmMTA), the core component of the m6A methyltransferase complex, and decrease in DNA methylation 1 (ZmDDM1), a key chromatin-remodeling factor that regulates DNA methylation. Genes with m6A modification were coordinated with a much higher level of DNA methylation than genes without m6A modification. Dysfunction of ZmMTA caused severe arrest during maize embryogenesis and endosperm development, leading to a significant decrease in CHH methylation in the 5' region of m6A-modified genes. Instead, loss of function of ZmDDM1 had no noteworthy effects on ZmMTA-related activity. This study establishes a direct link between m6A and 5mC during maize kernel development and provides insights into the interplay between RNA modification and DNA methylation.

DNA methylome analysis reveals novel insights into active hypomethylated regulatory mechanisms of temperature-dependent flower opening in Osmanthus fragrans

Short-term ambient low temperature (ALT) stimulation is necessary for Osmanthus fragrans to facilitate continued flower opening after floral bud development reaches maturity. DNA methylation, a vital epigenetic modification, regulates various biological processes in response to temperature fluctuations. However, its role in temperature-driven flower opening remains elusive. In this study, we identified the pivotal timeframe during which O. fragrans promptly detected temperature cues. Using whole-genome bisulfite sequencing, we explored global DNA hypomethylation during this phase, with the most significant changes occurring in CHH sequence contexts. Auxin transport inhibitor (TIBA) application revealed that ALT-induced endogenous auxin accumulation promoted peduncle elongation. In our mRNA-seq analysis, we discovered that the differentially expressed genes (DEGs) with hypo-differentially methylated regions (hypo-DMRs) were mainly enriched in auxin and temperature response, RNA processing, and carbohydrate and lipid metabolism. Transcripts of three DNA demethylase genes (OfROS1a, OfDML3, OfDME) showed upregulation. Furthermore, all DNA methylase genes, except OfCMT2b, also displayed increased expression, specifically with two of them, OfCMT3a and OfCMT1, being associated with hypo-DMRs. Promoter assays showed that OfROS1a, with promoters containing low-temperature- and auxin-responsive elements, were activated by ALT and exogenous IAA at low concentrations but inhibited at high concentrations. Overexpression of OfROS1 reduced endogenous auxin levels but enhanced the expression of genes related to auxin response and spliceosome in petunia. Furthermore, OfROS1 promoted sucrose synthesis in petunia corollas. Our data characterized the rapid response of active DNA hypomethylation to ALT and suggested a possible epiregulation of temperature-dependent flower opening in O. fragrans. This study revealed the pivotal role of DNA hypomethylation in O. fragrans during the ALT-responsive phase before flower opening, involving dynamic DNA demethylation, auxin signaling modulation, and a potential feedback loop between hypomethylation and methylation.

DNA methylation in the wild: epigenetic transgenerational inheritance can mediate adaptation in clones of wild strawberry (Fragaria vesca)

Due to the accelerating climate change, it is crucial to understand how plants adapt to rapid environmental changes. Such adaptation may be mediated by epigenetic mechanisms like DNA methylation, which could heritably alter phenotypes without changing the DNA sequence, especially across clonal generations. However, we are still missing robust evidence of the adaptive potential of DNA methylation in wild clonal populations. Here, we studied genetic, epigenetic and transcriptomic variation of Fragaria vesca, a predominantly clonally reproducing herb. We examined samples from 21 natural populations across three climatically distinct geographic regions, as well as clones of the same individuals grown in a common garden. We found that epigenetic variation was partly associated with climate of origin, particularly in non-CG contexts. Importantly, a large proportion of this variation was heritable across clonal generations. Additionally, a subset of these epigenetic changes affected the expression of genes mainly involved in plant growth and responses to pathogen and abiotic stress. These findings highlight the potential influence of epigenetic changes on phenotypic traits. Our findings indicate that variation in DNA methylation, which can be environmentally inducible and heritable, may enable clonal plant populations to adjust to their environmental conditions even in the absence of genetic adaptation.

A transposable element-derived siRNAs involve DNA hypermethylation at the promoter of OsGSTZ4 for cadmium tolerance in rice

Environmental contaminants such as cadmium (Cd) pose high risks to crop production and human health. The genetic basis for regulation of Cd stress-responsive genes for plant adaptation to adverse environments remains poorly understood. In this study, we characterized a rice Zeta family glutathione-S-transferase (OsGSTZ4) gene for Cd detoxification. Heterologous expression of OsGSTZ4 in a yeast (Saccharomyces cerevisiae) conferred cellular Cd tolerance. Transgenic rice overexpressing OsGSTZ4 improved plant growth, attenuated Cd-induced toxicity, and accumulated more Cd in roots. OsGSTZ4 transcription was rapidly induced 3 h after Cd exposure and then declined to the basal level. This was followed by (days after Cd treatment) by CHH hypermethylation (by 41.2 %) at a MITE (Miniature Inverted-repeat Transposable Element) transposable element (TE) inserted in the 5'-untranscribed region (UTR) (-1,722 ∼ -1,392 bp) of OsGSTZ4. Meanwhile, three 24-nt siRNAs derived from the TE (-1,722 ∼ -1,471 bp) were detected and was also rapidly enriched under Cd stress. To validate the possibility that Cd-induced change in OsGSTZ4 expression correlates with the siRNAs-involved CHH methylation through an RdDM (RNA-directed DNA methylation) pathway, genetic analyses were performed. We found that the CHH methylation at the promoter and transcript level of OsGSTZ4 were compromised in the osdrm2 (loss of function for CHH methylation) and osrdr2i (defective in RNA-dependent RNA polymerase 2) but did not change in other types of methyltransferases such as osmet1, ossdg714 or osros1. Promoter deletion analyses confirmed that the siRNA target sequences were essential for the proper expression of OsGSTZ4. Our studies reveal an unusual feedback mechanism by which the Cd-induced rapid OsGSTZ4 expression for Cd tolerance would interplay with the late CHH hypermethylation to silence the TE through the 24-nt siRNAs- and Osdrm2-mediated RdDM pathway, and help understand the diversity of gene regulation via an epigenetic mechanism for rice adaptation to the environmental stress.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

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

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

Conceptual Framework of Epigenetic Analyses of Plant Responses to Sedentary Endoparasitic Nematodes

Epigenetic modifications including miRNA regulation, DNA methylation, and histone modifications play fundamental roles in establishing the interactions between host plants and parasitic nematodes. Over the past decade, an increasing number of studies revealed the key functions of various components of the plant epigenome in the regulation of gene expression and shaping plant responses to nematode infection. In this chapter, we provide a conceptual framework for methods used to investigate epigenetic regulation during plant-nematode interactions. We focus specifically on current and emerging methods used to study miRNA regulation and function. We also highlight various methods and analytical tools used to profile DNA methylation patterns and histone modification marks at the genome level. Our intention is simply to explain the advantages of various methods and how to overcome some limitations. With rapid development of single-cell sequencing technology and genome editing, advanced and new methodologies are expected to emerge in the near future to further improve our understanding of epigenetic regulation and function during plant-nematode interactions.

The SC10-RNase promoter displays changes in DNA methylation patterns through pistil development in self-incompatible Nicotiana alata

In Solanaceae, self-incompatibility is a genetic mechanism that prevents endogamy in plant populations. Expression of the S-determinants, S-RNase, and SLF, is tightly regulated during pistil and pollen development. However, the molecular mechanism of gene expression regulation in S-RNase-based self-incompatibility systems must be better understood. Here, we identified a 1.3 Kbp sequence upstream to the coding region of the functional SC10-RNase allele from the self-incompatible Nicotiana alata, which directs SC10-RNase expression in mature pistils. This SC10-RNase promoter includes a 300 bp region with minimal elements that sustain the SC10-RNase expression. Likewise, a fragment of a transposable element from the Gypsy family of retrotransposons is also present at the -320 bp position. Nevertheless, its presence does not affect the expression of the SC10-RNase in mature pistils. Additionally, we determined that the SC10-RNase promoter undergoes different DNA methylation states during pistil development, being the mCHH methylation context the most frequent close to the transcription start site at pistil maturity. We hypothesized that the Gypsy element at the SC10-RNase promoter might contribute to the DNA methylation remodeling on the three sequence contexts analyzed here. We propose that mCHH methylation enrichment and other regulatory elements in the S-RNase coding region regulate the specific and abundant SC10-RNase expression in mature pistils in N. alata.

miRNA-like secondary structures in maize (Zea mays) genes and transposable elements correlate with small RNAs, methylation, and expression

RNA molecules carry information in their primary sequence and also their secondary structure. Secondary structure can confer important functional information, but it is also a signal for an RNAi-like host epigenetic response mediated by small RNAs (smRNAs). In this study, we used two bioinformatic methods to predict local secondary structures across features of the maize genome, focusing on small regions that had similar folding properties to pre-miRNA loci. We found miRNA-like secondary structures to be common in genes and most, but not all, superfamilies of RNA and DNA transposable elements (TEs). The miRNA-like regions map to a higher diversity of smRNAs than regions without miRNA-like structure, explaining up to 27% of variation in smRNA mapping for some TE superfamilies. This mapping bias is more pronounced among putatively autonomous TEs relative to nonautonomous TEs. Genome-wide, miRNA-like regions are also associated with elevated methylation levels, particularly in the CHH context. Among genes, those with miRNA-like secondary structure are 1.5-fold more highly expressed, on average, than other genes. However, these genes are also more variably expressed across the 26 nested association mapping founder lines, and this variability positively correlates with the number of mapping smRNAs. We conclude that local miRNA-like structures are a nearly ubiquitous feature of expressed regions of the maize genome, that they correlate with higher smRNA mapping and methylation, and that they may represent a trade-off between functional requirements and the potentially negative consequences of smRNA production.

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

MAIN CONCLUSION

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

PSTVd infection in Nicotiana benthamiana plants has a minor yet detectable effect on CG methylation

Viroids are small circular RNAs infecting a wide range of plants. They do not code for any protein or peptide and therefore rely on their structure for their biological cycle. Observed phenotypes of viroid infected plants are thought to occur through changes at the transcriptional/translational level of the host. A mechanism involved in such changes is RNA-directed DNA methylation (RdDM). Till today, there are contradictory works about viroids interference of RdDM. In this study, we investigated the epigenetic effect of viroid infection in Nicotiana benthamiana plants. Using potato spindle tuber viroid (PSTVd) as the triggering pathogen and via bioinformatic analyses, we identified endogenous gene promoters and transposable elements targeted by 24 nt host siRNAs that differentially accumulated in PSTVd-infected and healthy plants. The methylation status of these targets was evaluated following digestion with methylation-sensitive restriction enzymes coupled with PCR amplification, and bisulfite sequencing. In addition, we used Methylation Sensitive Amplification Polymorphism (MSAP) followed by sequencing (MSAP-seq) to study genomic DNA methylation of 5-methylcytosine (5mC) in CG sites upon viroid infection. In this study we identified a limited number of target loci differentially methylated upon PSTVd infection. These results enhance our understanding of the epigenetic host changes as a result of pospiviroid infection.

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.

Methylome changes in Lolium perenne associated with long-term colonisation by the endophytic fungus Epichloë sp. LpTG-3 strain AR37

Epichloë spp. often form mutualistic interactions with cool-season grasses, such as Lolium perenne. However, the molecular mechanisms underlying this interaction remain poorly understood. In this study, we employed reduced representation bisulfite sequencing method (epiGBS) to investigate the impact of the Epichloë sp. LpTG-3 strain AR37 on the methylome of L. perenne across multiple grass generations and under drought stress conditions. Our results showed that the presence of the endophyte leads to a decrease in DNA methylation across genomic features, with differentially methylated regions primarily located in intergenic regions and CHH contexts. The presence of the endophyte was consistently associated with hypomethylation in plants across generations. This research sheds new light on the molecular mechanisms governing the mutualistic interaction between Epichloë sp. LpTG-3 strain AR37 and L. perenne. It underscores the role of methylation changes associated with endophyte infection and suggests that the observed global DNA hypomethylation in L. perenne may be influenced by factors such as the duration of the endophyte-plant association and the accumulation of genetic and epigenetic changes over time.

Rice transformation treatments leave specific epigenome changes beyond tissue culture

During transgenic plant production, tissue culture often carries epigenetic, and genetic changes that underlie somaclonal variations, leading to unpredictable phenotypes. Additionally, specific treatments for rice (Oryza sativa) transformation processes may individually or jointly contribute to somaclonal variations, but their specific impacts on rice epigenomes toward transcriptional variations remain unknown. Here, the impact of individual transformation treatments on genome-wide DNA methylation and the transcriptome were examined. In addition to activating stress-responsive genes, individual transformation components targeted different gene expression modules that were enriched in specific functional categories. The transformation treatments strongly impacted DNA methylation and expression; 75% were independent of tissue culture. Furthermore, our genome-wide analysis showed that the transformation treatments consistently resulted in global hypo-CHH methylation enriched at promoters highly associated with downregulation, particularly when the promoters were colocalized with miniature inverted-repeat transposable elements. Our results clearly highlight the specificity of impacts triggered by individual transformation treatments during rice transformation with the potential association between DNA methylation and gene expression. These changes in gene expression and DNA methylation resulting from rice transformation treatments explain a significant portion of somaclonal variations, that is, way beyond the tissue culture effect.

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

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

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

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

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

Heritable changes of epialleles in maize can be triggered in the absence of DNA methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, very little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (mediator of paramutation1) with that of their parents, wild type siblings, and backcrossed progeny in maize. Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of the CHH DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in F1 plants did require Mop1, initiation of the changes in the epigenetic state of TCM DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is not dependent on RNA-directed DNA methylation.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

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.

Dynamic DNA methylation changes reveal tissue-specific gene expression in sugarcane

DNA methylation is an important mechanism for the dynamic regulation of gene expression and silencing of transposons during plant developmental processes. Here, we analyzed genome-wide methylation patterns in sugarcane (Saccharum officinarum) leaves, roots, rinds, and piths at single-base resolution. DNA methylation patterns were similar among the different sugarcane tissues, whereas DNA methylation levels differed. We also found that DNA methylation in different genic regions or sequence contexts plays different roles in gene expression. Differences in methylation among tissues resulted in many differentially methylated regions (DMRs) between tissues, particularly CHH DMRs. Genes overlapping with DMRs tended to be differentially expressed (DEGs) between tissues, and these DMR-associated DEGs were enriched in biological pathways related to tissue function, such as photosynthesis, sucrose synthesis, stress response, transport, and metabolism. Moreover, we observed many DNA methylation valleys (DMVs), which always overlapped with transcription factors (TFs) and sucrose-related genes, such as WRKY, bZIP, WOX, SPS, and FBPase. Collectively, these findings provide significant insights into the complicated interplay between DNA methylation and gene expression and shed light on the epigenetic regulation of sucrose-related genes in sugarcane.

The evolution and function of transposons in epigenetic regulation in response to the environment

Transposable elements (TEs) make up a major proportion of plant genomes. Despite their prevalence genome-wide, TEs are often tossed aside as "junk DNA" since they rarely cause phenotypes, and epigenetic mechanisms silence TEs to prevent them from causing deleterious mutations through movement. While this bleak picture of TEs in genomes is true on average, a growing number of examples across many plant species point to TEs as drivers of phenotypic diversity and novel stress responses. Examples of TE-influenced phenotypes illustrate the many ways that novel transposition events can alter local gene expression and how this relates to potential variation in plant responses to environmental stress. Since TE families and insertions at the locus level lack evolutionary conservation, advancements in the field will require TE experts across diverse species to identify and utilize TE variation in their own systems as a means of crop improvement.

The Heterogeneity in the Landscape of Gene Dominance in Maize is Accompanied by Unique Chromatin Environments

Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.

You shall not pass! A Chromatin barrier story in plants

As in other eukaryotes, the plant genome is functionally organized in two mutually exclusive chromatin fractions, a gene-rich and transcriptionally active euchromatin, and a gene-poor, repeat-rich, and transcriptionally silent heterochromatin. In Drosophila and humans, the molecular mechanisms by which euchromatin is preserved from heterochromatin spreading have been extensively studied, leading to the identification of insulator DNA elements and associated chromatin factors (insulator proteins), which form boundaries between chromatin domains with antagonistic features. In contrast, the identity of factors assuring such a barrier function remains largely elusive in plants. Nevertheless, several genomic elements and associated protein factors have recently been shown to regulate the spreading of chromatin marks across their natural boundaries in plants. In this minireview, we focus on recent findings that describe the spreading of chromatin and propose avenues to improve the understanding of how plant chromatin architecture and transitions between different chromatin domains are defined.

Minireview: Chromatin-based regulation of iron homeostasis in plants

Plants utilize delicate mechanisms to effectively respond to changes in the availability of nutrients such as iron. The responses to iron status involve controlling gene expression at multiple levels. The regulation of iron deficiency response by a network of transcriptional regulators has been extensively studied and recent research has shed light on post-translational control of iron homeostasis. Although not as considerably investigated, an increasing number of studies suggest that histone modification and DNA methylation play critical roles during iron deficiency and contribute to fine-tuning iron homeostasis in plants. This review will focus on the current understanding of chromatin-based regulation on iron homeostasis in plants highlighting recent studies in Arabidopsis and rice. Understanding iron homeostasis in plants is vital, as it is not only relevant to fundamental biological questions, but also to agriculture, biofortification, and human health. A comprehensive overview of the effect and mechanism of chromatin-based regulation in response to iron status will ultimately provide critical insights in elucidating the complexities of iron homeostasis and contribute to improving iron nutrition in plants.

Single-base resolution methylomes of somatic embryogenesis in Theobroma cacao L. reveal epigenome modifications associated with somatic embryo abnormalities

Propagation by somatic embryogenesis in Theobroma cacao has some issues to be solved, as many morphologically abnormal somatic embryos that do not germinate into plants are frequently observed, thus hampering plant production on a commercial scale. For the first time the methylome landscape of T. cacao somatic embryogenesis was examined, using whole-genome bisulfite sequencing technique, with the aim to understand the epigenetic basis of somatic embryo abnormalities. We identified 873 differentially methylated genes (DMGs) in the CpG context between zygotic embryos, normal and abnormal somatic embryos, with important roles in development, programmed cell death, oxidative stress, and hypoxia induction, which can help to explain the morphological abnormalities of somatic embryos. We also identified the role of ethylene and its precursor 1-aminocyclopropane-1-carboxylate in several biological processes, such as hypoxia induction, cell differentiation and cell polarity, that could be associated to the development of abnormal somatic embryos. The biological processes and the hypothesis of ethylene and its precursor involvement in the somatic embryo abnormalities in cacao are discussed.

Triticale doubled haploid plant regeneration factors linked by structural equation modeling

Triticale regeneration via anther culture faces many difficulties, e.g., a low percentage of regenerated plants and the presence of albinos. Plant regeneration may be affected by abiotic stresses and by ingredients added to the induction medium. The latter influences biochemical pathways and plant regeneration efficiency. Among such ingredients, copper and silver ions acting as cofactors for enzymatic reactions are of interest. However, their role in plant tissue cultures and relationships with biochemical pathways has not been studied yet.

The study evaluated relationships between DNA methylation, changes in DNA sequence variation, and green plant regeneration efficiency influenced by copper and silver ions during triticale plant regeneration. For this purpose, a biological model based on donor plants and their regenerants, a methylation-sensitive amplified fragment length polymorphism, and structural equation modeling were employed.

The green plant regeneration efficiency varied from 0.71 to 6.06 green plants per 100 plated anthers. The values for the components of tissue culture-induced variation related to cytosine methylation in a CHH sequence context (where H is A, C, or T) were 8.65% for sequence variation, 0.76% for DNA demethylation, and 0.58% for de novo methylation. The proposed model states that copper ions affect the regeneration efficiency through cytosine methylation and may induce mutations through, e.g., oxidative processes, which may interfere with the green plant regeneration efficiency. The linear regression confirms that the plant regeneration efficiency rises with increasing copper ion concentration in the absence of Ag ions in the induction medium. The least absolute shrinkage and selection operator regression shows that de novo methylation, demethylation, and copper ions may be involved in the green plant regeneration efficiency. According to structural equation modeling, copper ions play a central role in the model determining the regeneration efficiency.

Whole-genome DNA methylation patterns of Oryza sativa (L.) and Oryza glumaepatula (Steud) genotypes associated with aluminum response

Epigenetic mechanisms in crops have emerged as a fundamental factor in plant adaptation and acclimation to biotic and abiotic stresses. Among described epigenetic mechanisms, DNA methylation has been defined as the most studied epigenetic modification involved in several developmental processes. It has been shown that contrasting methylation marks are associated with gene expression variations between cultivated and wild crop species. In this study, we analyzed single-base resolution methylome maps for Oryza sativa (a cultivated species) and Oryza glumaepatula (a wild species) genotypes grown under control conditions. Our results showed that overall, genome-wide methylation profiles are mainly conserved between both species, nevertheless, there are several differentially methylated regions with species-specific methylation patterns. In addition, we analyzed the association of identified DNA methylation marks in relation with Aluminum-tolerance levels of studied genotypes. We found several differentially methylated regions (DMRs) and DMR-associated genes (DAGs) that are linked with Al tolerance. Some of these DAGs have been previously reported as differentially expressed under Al exposure in O. sativa. Complementarily a Transposable Elements (TE) analysis revealed that specific aluminum related genes have associated-TEs potentially regulated by DNA methylation. Interestingly, the DMRs and DAGs between Al-tolerant and susceptible genotypes were different between O. sativa and O. glumaepatula, suggesting that methylation patterns related to Al responses are unique for each rice species. Our findings provide novel insights into DNA methylation patterns in wild and cultivated rice genotypes and their possible role in the regulation of plant stress responses.

Signatures of polygenic adaptation align with genome-wide methylation patterns in wild strawberry plants

Epigenetic inheritance can drive adaptive evolution independently of DNA sequence variation. However, to what extent epigenetic variation represents an autonomous evolutionary force remains largely elusive. Through gene ontology and comparative analyses of genomic and epigenomic variation of wild strawberry plants raised in distinct drought settings, we characterised genome-wide covariation between single nucleotide polymorphisms (SNPs) and differentially methylated cytosines (DMCs). Covariation between SNPs and DMCs was independent of genomic proximity, but instead associated with fitness-related processes such as stress responses, genome regulation and reproduction. We expected this functional SNP-DMC covariation to be driven by adaptive evolution canalising SNP and DMC variation, but instead observed significantly lower covariation with DMCs for adaptive rather than for neutral SNPs. Drought-induced DMCs frequently co-varied with tens of SNPs, suggesting high genomic redundancy as a broad potential basis for polygenic adaptation of gene expression. Our findings suggest that stress-responsive DMCs initially co-vary with many SNPs under increased environmental stress, and that natural selection acting upon several of these SNPs subsequently reduces standing covariation with stress-responsive DMCs. Our study supports DNA methylation profiles that represent complex quantitative traits rather than autonomous evolutionary forces. We provide a conceptual framework for polygenic regulation and adaptation shaping genome-wide methylation patterns in plants.

Genome-wide loss of CHH methylation with limited transcriptome changes in Setaria viridis DOMAINS REARRANGED METHYLTRANSFERASE (DRM) mutants

The DOMAINS REARRANGED METHYLTRANSFERASEs (DRMs) are crucial for RNA-directed DNA methylation (RdDM) in plant species. Setaria viridis is a model monocot species with a relatively compact genome that has limited transposable element (TE) content. CRISPR-based genome editing approaches were used to create loss-of-function alleles for the two putative functional DRM genes in S. viridis to probe the role of RdDM. Double mutant (drm1ab) plants exhibit some morphological abnormalities but are fully viable. Whole-genome methylation profiling provided evidence for the widespread loss of methylation in CHH sequence contexts, particularly in regions with high CHH methylation in wild-type plants. Evidence was also found for the locus-specific loss of CG and CHG methylation, even in some regions that lack CHH methylation. Transcriptome profiling identified genes with altered expression in the drm1ab mutants. However, the majority of genes with high levels of CHH methylation directly surrounding the transcription start site or in nearby promoter regions in wild-type plants do not have altered expression in the drm1ab mutant, even when this methylation is lost, suggesting limited regulation of gene expression by RdDM. Detailed analysis of the expression of TEs identified several transposons that are transcriptionally activated in drm1ab mutants. These transposons are likely to require active RdDM for the maintenance of transcriptional repression.

Comparative DNA Methylome of Phytoplasma Associated Retrograde Metamorphosis in Sesame (Sesamum indicum L.)

Phytoplasma-associated diseases such as phyllody and little leaf are critical threats to sesame cultivation worldwide. The mechanism of the dramatic conversion of flowers to leafy structures leading to yield losses and the drastic reduction in leaf size due to Phytoplasma infection remains yet to be identified. Cytosine methylation profiles of healthy and infected sesame plants studied using Whole Genome Bisulfite Sequencing (WGBS) and Quantitative analysis of DNA methylation with the real-time PCR (qAMP) technique revealed altered DNA methylation patterns upon infection. Phyllody was associated with global cytosine hypomethylation, though predominantly in the CHH (where H = A, T or C) context. Interestingly, comparable cytosine methylation levels were observed between healthy and little leaf-affected plant samples in CG, CHG and CHH contexts. Among the different genomic fractions, the highest number of differentially methylated Cytosines was found in the intergenic regions, followed by promoter, exonic and intronic regions in decreasing order. Further, most of the differentially methylated genes were hypomethylated and were mainly associated with development and defense-related processes. Loci for STOREKEEPER protein-like, a DNA-binding protein and PP2-B15, an F-Box protein, responsible for plugging sieve plates to maintain turgor pressure within the sieve tubes were found to be hypomethylated by WGBS, which was confirmed by methylation-dependent restriction digestion and qPCR. Likewise, serine/threonine-protein phosphatase-7 homolog, a positive regulator of cryptochrome signaling involved in hypocotyl and cotyledon growth and probable O-methyltransferase 3 locus were determined to be hypermethylated. Phytoplasma infection-associated global differential methylation as well as the defense and development-related loci reported here for the first time significantly elucidate the mechanism of phytoplasma-associated disease development.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Antioxidant Regulation and DNA Methylation Dynamics During Mikania micrantha Seed Germination Under Cold Stress

As a primary goal, adaptation to cold climate could expand an invasion range of exotic plants. Here, we aimed to explore the regulation strategy of M. micrantha seed development under cold stress through molecular physiology and multi-omics analysis. Significant increase of hydrogen peroxide, malondialdehyde, and electrolyte leakage observed under cold stress revealed that oxidative damage within M. micrantha seed cells was induced in the initial germination phase. Proteomic data underscored an activation of antioxidant activity to maintain redox homeostasis, with a cluster of antioxidant proteins identified. Genomic-wide transcriptome, in combination with time-series whole-genome bisulfite sequencing mining, elucidated that seven candidate genes, which were the target of DNA demethylation-dependent ROS scavenging, were possibly associated with an M. micrantha germ break. Progressive gain of CHH context DNA methylation identified in an early germination phrase suggested a role of a DNA methylation pathway, while an active DNA demethylation pathway was also initiated during late seed development, which was in line with the expression trend of methylation and demethylation-related genes verified through qRT-PCR. These data pointed out that cold-dependent DNA demethylation and an antioxidant regulatory were involved together in restoring seed germination. The expression level of total 441 genes presented an opposite trend to the methylation divergence, while the expression of total 395 genes was proved to be negatively associated with their methylation levels. These data provided new insights into molecular reprograming events during M. micrantha seed development.

Small RNAs mediate transgenerational inheritance of genome-wide trans-acting epialleles in maize

Background

Hybridization and backcrossing are commonly used in animal and plant breeding to induce heritable variation including epigenetic changes such as paramutation. However, the molecular basis for hybrid-induced epigenetic memory remains elusive.

Results

Here, we report that hybridization between the inbred parents B73 and Mo17 induces trans-acting hypermethylation and hypomethylation at thousands of loci; several hundreds (~ 3%) are transmitted through six backcrossing and three selfing generations. Notably, many transgenerational methylation patterns resemble epialleles of the nonrecurrent parent, despite > 99% of overall genomic loci are converted to the recurrent parent. These epialleles depend on 24-nt siRNAs, which are eliminated in the isogenic hybrid Mo17xB73:mop1-1 that is defective in siRNA biogenesis. This phenomenon resembles paramutation-like events and occurs in both intraspecific (Mo17xB73) and interspecific (W22xTeosinte) hybrid maize populations. Moreover, siRNA abundance and methylation levels of these epialleles can affect expression of their associated epigenes, many of which are related to stress responses.

Conclusion

Divergent siRNAs between the hybridizing parents can induce trans-acting epialleles in the hybrids, while the induced epigenetic status is maintained for transgenerational inheritance during backcross and hybrid breeding, which alters epigene expression to enhance growth and adaptation. These genetic and epigenetic principles may apply broadly from plants to animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02614-0.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Resetting of the 24-nt siRNA landscape in rice zygotes

The zygote, a totipotent stem cell, is crucial to the life cycle of sexually reproducing organisms. It is produced by the fusion of two differentiated cells-the egg and sperm, which in plants have radically different siRNA transcriptomes from each other and from multicellular embryos. Owing to technical challenges, the epigenetic changes that accompany the transition from differentiated gametes to totipotent zygote are poorly understood. Because siRNAs serve as both regulators and outputs of the epigenome, we characterized small RNA transcriptomes of zygotes from rice. Zygote small RNAs exhibit extensive maternal carryover and an apparent lack of paternal contribution, indicated by absence of sperm signature siRNAs. Zygote formation is accompanied by widespread redistribution of 24-nt siRNAs relative to gametes, such that ∼70% of the zygote siRNA loci do not overlap any egg cell siRNA loci. Newly detected siRNA loci in zygote are gene-proximal and not associated with centromeric heterochromatin, similar to canonical siRNAs, in sharp contrast to gametic siRNA loci that are gene-distal and heterochromatic. In addition, zygote but not egg siRNA loci are associated with high DNA methylation in the mature embryo. Thus, the zygote begins transitioning before the first embryonic division to an siRNA profile that is associated with future RdDM in embryogenesis. These findings indicate that, in addition to changes in gene expression, the transition to totipotency in the plant zygote is accompanied by resetting of the epigenetic reprogramming that occurred during gamete formation.