DNA methylation of retrotransposons, DNA transposons and genes in sugar beet (Beta vulgaris L.)

Profiling genome-wide methylation in two maples: Fine-scale approaches to detection with nanopore technology

DNA methylation is critical to the regulation of transposable elements and gene expression and can play an important role in the adaptation of stress response mechanisms in plants. Traditional methods of methylation quantification rely on bisulfite conversion that can compromise accuracy. Recent advances in long-read sequencing technologies allow for methylation detection in real time. The associated algorithms that interpret these modifications have evolved from strictly statistical approaches to Hidden Markov Models and, recently, deep learning approaches. Much of the existing software focuses on methylation in the CG context, but methylation in other contexts is important to quantify, as it is extensively leveraged in plants. Here, we present methylation profiles for two maple species across the full range of 5mC sequence contexts using Oxford Nanopore Technologies (ONT) long-reads. Hybrid and reference-guided assemblies were generated for two new Acer accessions: Acer negundo (box elder; 65x ONT and 111X Illumina) and Acer saccharum (sugar maple; 93x ONT and 148X Illumina). The ONT reads generated for these assemblies were re-basecalled, and methylation detection was conducted in a custom pipeline with the published Acer references (PacBio assemblies) and hybrid assemblies reported herein to generate four epigenomes. Examination of the transposable element landscape revealed the dominance of LTR Copia elements and patterns of methylation associated with different classes of TEs. Methylation distributions were examined at high resolution across gene and repeat density and described within the broader angiosperm context, and more narrowly in the context of gene family dynamics and candidate nutrient stress genes.

The roles of epigenetic regulators in plant regeneration: Exploring patterns amidst complex conditions

Plants possess remarkable capability to regenerate upon tissue damage or optimal environmental stimuli. This ability not only serves as a crucial strategy for immobile plants to survive through harsh environments, but also made numerous modern plant improvements techniques possible. At the cellular level, this biological process involves dynamic changes in gene expression that redirect cell fate transitions. It is increasingly recognized that chromatin epigenetic modifications, both activating and repressive, intricately interact to regulate this process. Moreover, the outcomes of epigenetic regulation on regeneration are influenced by factors such as the differences in regenerative plant species and donor tissue types, as well as the concentration and timing of hormone treatments. In this review, we focus on several well-characterized epigenetic modifications and their regulatory roles in the expression of widely studied morphogenic regulators, aiming to enhance our understanding of the mechanisms by which epigenetic modifications govern plant regeneration.

DNA methylation regulates biosynthesis of tanshinones and phenolic acids during growth of Salvia miltiorrhiza

DNA methylation plays a crucial role in the regulation of plant growth and the biosynthesis of secondary metabolites. Danshen (Salvia miltiorrhiza) is a valuable Chinese herbal medicine commonly used to treat cardiovascular diseases; its active ingredients are tanshinones and phenolic acids, which primarily accumulate in roots. Here, we conducted a targeted metabolic analysis of S. miltiorrhiza roots at 3 distinct growth stages: 40 d old (r40), 60 d old (r60), and 90 d old (r90). The contents of tanshinones (cryptotanshinone, tanshinone I, tanshinone IIA, and rosmariquinone) and phenolic acids (rosmarinic acid and salvianolic acid B) gradually increased during plant development. Whole-genome bisulfite sequencing and transcriptome sequencing of roots at the 3 growth stages revealed an increased level of DNA methylation in the CHH context (H represents A, T, or C) context at r90 compared with r40 and r60. Increased DNA methylation levels were associated with elevated expression of various genes linked to epigenetic regulations, including CHROMOMETHYLASE2 (SmCMT2), Decrease in DNA Methylation 1 (SmDDM1), Argonaute 4 (SmAGO4), and DOMAINS REARRANGED METHYLTRANSFERASE 1 (SmDRM1). Moreover, expression levels of many genes involved in tanshinone and salvianolic acid biosynthesis, such as copalyldiphosphate synthase 5 (SmCPS5), cytochrome P450-related enzyme (SmCYP71D464), geranylgeranyl diphosphate synthase (SmGGPPS1), geranyl diphosphate synthase (SmGPPS), hydroxyphenylpyruvate reductase (SmHPPR), and hydroxyphenylpyruvate dioxygenase (SmHPPD), were altered owing to hyper-methylation, indicating that DNA methylation plays an important role in regulating tanshinone and phenolic acid accumulation. Our data shed light on the epigenetic regulation of root growth and the biosynthesis of active ingredients in S. miltiorrhiza, providing crucial clues for further improvement of active compound production via molecular breeding in S. miltiorrhiza.

Methylomes as key features for predicting recombination in some plant species

Knowing how chromosome recombination works is essential for plant breeding. It enables the design of crosses between different varieties to combine desirable traits and create new ones. This is because the meiotic crossovers between homologous chromatids are not purely random, and various strategies have been developed to describe and predict such exchange events. Recent studies have used methylation data to predict chromosomal recombination in rice using machine learning models. This approach proved successful due to the presence of a positive correlation between the CHH context cytosine methylation and recombination rates in rice chromosomes. This paper assesses the question if methylation can be used to predict recombination in four plant species: Arabidopsis, maize, sorghum, and tomato. The results indicate a positive association between CHH context methylation and recombination rates in certain plant species, with varying degrees of strength in their relationships. The CG and CHG methylation contexts show negative correlation with recombination. Methylation data was key effectively in predicting recombination in sorghum and tomato, with a mean determination coefficient of 0.65 ± 0.11 and 0.76 ± 0.05, respectively. In addition, the mean correlation values between predicted and experimental recombination rates were 0.83 ± 0.06 for sorghum and 0.90 ± 0.05 for tomato, confirming the significance of methylomes in both monocotyledonous and dicotyledonous species. The predictions for Arabidopsis and maize were not as accurate, likely due to the comparatively weaker relationships between methylation contexts and recombination, in contrast to sorghum and tomato, where stronger associations were observed. To enhance the accuracy of predictions, further evaluations using data sets closely related to each other might prove beneficial. In general, this methylome-based method holds great potential as a reliable strategy for predicting recombination rates in various plant species, offering valuable insights to breeders in their quest to develop novel and improved varieties.

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

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

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

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.

In-Plant Persistence and Systemic Transport of Nicotiana benthamiana Retrozyme RNA

Retrozymes are nonautonomous retrotransposons with hammerhead ribozymes in their long terminal repeats (LTRs). Retrozyme transcripts can be self-cleaved by the LTR ribozyme, circularized, and can undergo RNA-to-RNA replication. Here, we demonstrate that the Nicotiana benthamiana genome contains hundreds of retrozyme loci, of which nine represent full-length retrozymes. The LTR contains a promoter directing retrozyme transcription. Although retrozyme RNA is easily detected in plants, the LTR region is heavily methylated, pointing to its transcriptional silencing, which can be mediated by 24 nucleotide-long retrozyme-specific RNAs identified in N. benthamiana. A transcriptome analysis revealed that half of the retrozyme-specific RNAs in plant leaves have no exact matches to genomic retrozyme loci, containing up to 13% mismatches with the closest genomic sequences, and could arise as a result of many rounds of RNA-to-RNA replication leading to error accumulation. Using a cloned retrozyme copy, we show that retrozyme RNA is capable of replication and systemic transport in plants. The presented data suggest that retrozyme loci in the N. benthamiana genome are transcriptionally inactive, and that circular retrozyme RNA can persist in cells due to its RNA-to-RNA replication and be transported systemically, emphasizing functional and, possibly, evolutionary links of retrozymes to viroids-noncoding circular RNAs that infect plants.

Methylation in the CHH Context Allows to Predict Recombination in Rice

DNA methylation is the most studied epigenetic trait. It is considered a key factor in regulating plant development and physiology, and has been associated with the regulation of several genomic features, including transposon silencing, regulation of gene expression, and recombination rates. Nonetheless, understanding the relation between DNA methylation and recombination rates remains a challenge. This work explores the association between recombination rates and DNA methylation for two commercial rice varieties. The results show negative correlations between recombination rates and methylated cytosine counts for all contexts tested at the same time, and for CG and CHG contexts independently. In contrast, a positive correlation between recombination rates and methylated cytosine count is reported in CHH contexts. Similar behavior is observed when considering only methylated cytosines within genes, transposons, and retrotransposons. Moreover, it is shown that the centromere region strongly affects the relationship between recombination rates and methylation. Finally, machine learning regression models are applied to predict recombination using the count of methylated cytosines in the CHH context as the entrance feature. These findings shed light on the understanding of the recombination landscape of rice and represent a reference framework for future studies in rice breeding, genetics, and epigenetics.

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

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

Absence of major epigenetic and transcriptomic changes accompanying an interspecific cross between peach and almond

Hybridization has been widely used in breeding of cultivated species showing low genetic variability, such as peach (Prunus persica). The merging of two different genomes in a hybrid often triggers a so-called "genomic shock" with changes in DNA methylation and in the induction of transposable element expression and mobilization. Here, we analysed the DNA methylation and transcription levels of transposable elements and genes in leaves of Prunus persica and Prunus dulcis and in an F1 hybrid using high-throughput sequencing technologies. Contrary to the "genomic shock" expectations, we found that the overall levels of DNA methylation in the transposable elements in the hybrid are not significantly altered compared with those of the parental genomes. We also observed that the levels of transcription of the transposable elements in the hybrid are in most cases intermediate as compared with that of the parental species and we have not detected cases of higher transcription in the hybrid. We also found that the proportion of genes whose expression is altered in the hybrid compared with the parental species is low. The expression of genes potentially involved in the regulation of the activity of the transposable elements is not altered. We can conclude that the merging of the two parental genomes in this Prunus persica x Prunus dulcis hybrid does not result in a "genomic shock" with significant changes in the DNA methylation or in the transcription. The absence of major changes may facilitate using interspecific peach x almond crosses for peach improvement.

Multi-omics data integration reveals link between epigenetic modifications and gene expression in sugar beet (Beta vulgaris subsp. vulgaris) in response to cold

Background

DNA methylation is thought to influence the expression of genes, especially in response to changing environmental conditions and developmental changes. Sugar beet (Beta vulgaris ssp. vulgaris), and other biennial or perennial plants are inevitably exposed to fluctuating temperatures throughout their lifecycle and might even require such stimulus to acquire floral competence. Therefore, plants such as beets, need to fine-tune their epigenetic makeup to ensure phenotypic plasticity towards changing environmental conditions while at the same time steering essential developmental processes. Different crop species may show opposing reactions towards the same abiotic stress, or, vice versa, identical species may respond differently depending on the specific kind of stress.

Results

In this study, we investigated common effects of cold treatment on genome-wide DNA methylation and gene expression of two Beta vulgaris accessions via multi-omics data analysis. Cold exposure resulted in a pronounced reduction of DNA methylation levels, which particularly affected methylation in CHH context (and to a lesser extent CHG) and was accompanied by transcriptional downregulation of the chromomethyltransferase CMT2 and strong upregulation of several genes mediating active DNA demethylation.

Conclusion

Integration of methylomic and transcriptomic data revealed that, rather than methylation having directly influenced expression, epigenetic modifications correlated with changes in expression of known players involved in DNA (de)methylation. In particular, cold triggered upregulation of genes putatively contributing to DNA demethylation via the ROS1 pathway. Our observations suggest that these transcriptional responses precede the cold-induced global DNA-hypomethylation in non-CpG, preparing beets for additional transcriptional alterations necessary for adapting to upcoming environmental changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08312-2.

The callus formation capacity of strawberry leaf explant is modulated by DNA methylation

Shoot regeneration from leaf tissue requires de-differentiation of cells from a highly differentiated state into an active dividing state, but how this physiological transition occurs and is regulated especially at epigenetic level remains obscure. Here we have characterized the DNA methylome represented by 5-methylcytosine (5mC) in leaf and the callus tissue derived from the leaf explant of woodland strawberry Fragaria vesca. We detected an overall increase of DNA methylation and distinct 5mC enrichment patterns in the CG, CHG and CHH sequence contexts in genetic and transposable elements. Our analyses revealed an intricate relation between DNA methylation and gene expression levels in leaf or leaf-derived callus. However, when considering the genes involved in callus formation and shoot regeneration, e.g. FvePLT3/7, FveWIND3, FveWIND4, FveLOG4 and FveIAA14, their dynamic transcription levels were associated with the differentially methylated regions located in the promoters or gene bodies, indicating a regulatory role of DNA methylation in the transcriptional regulation of pluripotency acquisition in strawberry. Furthermore, application of a DNA methyltransferase inhibitor 5'-azacytidine (5'-Aza) hampered both callus formation and shoot regeneration from the leaf explant. We further showed that 5'-Aza down-regulated the genes involved in cell wall integrity, such as expansin, pectin lyase and pectin methylesterase genes, suggesting an essential role of cell wall metabolism during callus formation. This study reveals the contribution of DNA methylation in callus formation capacity and will provide a basis for developing a strategy to improve shoot regeneration for basic and applied research applications.

Dynamic changes in DNA methylation occur in TE regions and affect cell proliferation during leaf-to-callus transition in Arabidopsis

Plant somatic cells can be reprogrammed into pluripotent cell mass, called callus, through a two-step in vitro tissue culture method. Incubation on callus-inducing medium triggers active cell proliferation to form a pluripotent callus. Notably, DNA methylation is implicated during callus formation, but a detailed molecular process regulated by DNA methylation remains to be fully elucidated. Here, we compared genome-wide DNA methylation profiles between leaf and callus tissues in Arabidopsis using whole-genome bisulphite-sequencing. Global distribution of DNA methylation showed that CHG methylation was increased, whereas CHH methylation was reduced especially around transposable element (TE) regions during the leaf-to-callus transition. We further analysed differentially expressed genes around differentially methylated TEs (DMTEs) during the leaf-to-callus transition and found that genes involved in cell cycle regulation were enriched and also constituted a coexpression gene network along with pluripotency regulators. In addition, a conserved DNA sequence analysis for upstream cis-elements led us to find a putative transcription factor associated with cell fate transition. CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) was newly identified as a regulator of plant regeneration, and consistently, the cca1lhy mutant displayed altered phenotypes in callus proliferation. Overall, these results suggest that DNA methylation coordinates cell cycle regulation during callus formation, and CCA1 may act as a key upstream coordinator at least in part in the processes.

Dynamic Changes of DNA Methylation During Wild Strawberry (Fragaria nilgerrensis) Tissue Culture

Tissue culture is an important tool for asexual propagation and genetic transformation of strawberry plants. In plant tissue culture, variation of DNA methylation is a potential source of phenotypic variation in regenerated plants. However, the genome wide dynamic methylation patterns of strawberry tissue culture remain unclear. In this study, we used whole-genome bisulfite sequencing (WGBS) to study genomic DNA methylation changes of a wild strawberry Fragaria nilgerrensis at six stages: from explants of shoot tips to outplanting and acclimation. Global methylation levels showed that CG sites exhibited the highest methylation level in all stages with an average of 49.5%, followed by CHG (33.2%) and CHH (12.4%). Although CHH accounted for the lowest proportion of total cytosine methylation, it showed the most obvious methylation change and the most of these changes occurred in the transposable element regions. The overall methylation levels alternately decreased and increased during the entire tissue culture process and the distribution of DNA methylation was non-uniform among different genetic regions. Furthermore, much more differentially methylated regions (DMRs) were detected in dedifferentiation and redifferentiation stages and most of them were transposable elements, suggesting these processes involved activating or silencing of amounts of transposons. The functional enrichment of the DMR-related genes indicated that genes involved in hormone metabolic processes, plant development and the stress response changed methylation throughout the tissue culture process. Finally, the quantitative real-time PCR (qRT-PCR) was conducted to examine the association of methylation and gene expression of a set of different methylated genes. Our findings give deeper insight into the epigenetic regulation of gene expression during the plant tissue cultures process, which will be useful in the efficient control of somaclonal variations and in crop improvement.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Elevated transcription of transposable elements is accompanied by het-siRNA-driven de novo DNA methylation in grapevine embryogenic callus

Background

Somatic variation is a valuable source of trait diversity in clonally propagated crops. In grapevine, which has been clonally propagated worldwide for centuries, important phenotypes such as white berry colour are the result of genetic changes caused by transposable elements. Additionally, epiallele formation may play a role in determining geo-specific (‘terroir’) differences in grapes and thus ultimately in wine. This genomic plasticity might be co-opted for crop improvement via somatic embryogenesis, but that depends on a species-specific understanding of the epigenetic regulation of transposable element (TE) expression and silencing in these cultures. For this reason, we used whole-genome bisulphite sequencing, mRNA sequencing and small RNA sequencing to study the epigenetic status and expression of TEs in embryogenic callus, in comparison with leaf tissue.

Results

We found that compared with leaf tissue, grapevine embryogenic callus cultures accumulate relatively high genome-wide CHH methylation, particularly across heterochromatic regions. This de novo methylation is associated with an abundance of transcripts from highly replicated TE families, as well as corresponding 24 nt heterochromatic siRNAs. Methylation in the TE-specific CHG context was relatively low over TEs located within genes, and the expression of TE loci within genes was highly correlated with the expression of those genes.

Conclusions

This multi-‘omics analysis of grapevine embryogenic callus in comparison with leaf tissues reveals a high level of genome-wide transcription of TEs accompanied by RNA-dependent DNA methylation of these sequences in trans. This provides insight into the genomic conditions underlying somaclonal variation and epiallele formation in plants regenerated from embryogenic cultures, which is an important consideration when using these tissues for plant propagation and genetic improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07973-9.

Broken, silent, and in hiding: tamed endogenous pararetroviruses escape elimination from the genome of sugar beet (Beta vulgaris)

BACKGROUND AND AIMS

Endogenous pararetroviruses (EPRVs) are widespread components of plant genomes that originated from episomal DNA viruses of the Caulimoviridae family. Due to fragmentation and rearrangements, most EPRVs have lost their ability to replicate through reverse transcription and to initiate viral infection. Similar to the closely related retrotransposons, extant EPRVs were retained and often amplified in plant genomes for several million years. Here, we characterize the complete genomic EPRV fraction of the crop sugar beet (Beta vulgaris, Amaranthaceae) to understand how they shaped the beet genome and to suggest explanations for their absent virulence.

METHODS

Using next- and third-generation sequencing data and genome assembly, we reconstructed full-length in silico representatives for the three host-specific EPRVs (beetEPRVs) in the B. vulgaris genome. Focusing on the endogenous caulimovirid beetEPRV3, we investigated its chromosomal localization, abundance and distribution by fluorescent in situ and Southern hybridization.

KEY RESULTS

Full-length beetEPRVs range between 7.5 and 10.7 kb in size, are heterogeneous in structure and sequence, and occupy about 0.3 % of the beet genome. Although all three beetEPRVs were assigned to the florendoviruses, they showed variably arranged protein-coding domains, different fragmentation, and preferences for diverse sequence contexts. We observed small RNAs that specifically target the individual beetEPRVs, indicating stringent epigenetic suppression. BeetEPRV3 sequences occur along all sugar beet chromosomes, preferentially in the vicinity of each other and are associated with heterochromatic, centromeric and intercalary satellite DNAs. BeetEPRV3 members also exist in genomes of related wild species, indicating an initial beetEPRV3 integration 13.4-7.2 million years ago.

CONCLUSIONS

Our study in beet illustrates the variability of EPRV structure and sequence in a single host genome. Evidence of sequence fragmentation and epigenetic silencing implies possible plant strategies to cope with long-term persistence of EPRVs, including amplification, fixation in the heterochromatin, and containment of EPRV virulence.

Transcriptome analysis of differentially expressed genes in chrysanthemum MET1 RNA interference lines

DNA methylation is the most important epigenetic modification involved in many essential biological processes. MET1 is one of DNA methyltransferases that affect the level of methylation in the entire genome. To explore the effect of MET1 gene silencing on gene expression profile of Chrysanthemum × morifolium 'Zijingling'. The stem section and leaves at the young stage were taken for transcriptome sequencing. MET1-_RNAi leaves had 8 differentially expressed genes while 156 differentially expressed genes were observed in MET1-_RNAi stem compared with control leaves and stem. These genes encode many key proteins in plant biological processes, such as transcription factors, signal transduction mechanisms, secondary metabolite synthesis, transport and catabolism and interaction. In general, 34.58% of the differentially expressed genes in leaves and stems were affected by the reduction of the MET1 gene. The differentially expressed genes in stem and leaves of transgenic plants went through significant changes. We found adequate amount of candidate genes associated with flowering, however, the number of genes with significant differences between transgenic and control lines was not too high. Several flowering related genes were screened out for gene expression verification and all of them were obseved as consistent with transcriptome data. These candidate genes may play important role in flowering variation of chrysanthemum. This study reveals the mechanism of CmMET1 interference on the growth and development of chrysanthemum at the transcriptional level, which provides the basis for further research on the epigenetic regulation mechanism in flower induction and development.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01022-1.

Dynamic changes in transposable element and gene methylation in mulberry (Morus notabilis) in response to Botrytis cinerea

DNA methylation has been proposed to regulate plant stress resistance. However, the dynamic changes in DNA methylation in woody plants and their correlations with pathogenic responses are not fully understood. Here, we present single-base maps of the DNA methylomes of mulberry (Morus notabilis) leaves that were subjected to a mock treatment or inoculation with Botrytis cinerea. Compared with the former, the latter showed decreased mCG and mCHG levels and increased mCHH levels. DNA methylation inhibitors reduced resistance gene methylation levels and enhanced mulberry resistance, suggesting that the hypomethylation of resistance genes affects mulberry resistance to B. cinerea. Virus-induced gene silencing of MnMET1 enhanced the expression of mulberry-resistance genes, thereby increasing the plant's resistance to B. cinerea. We also found that MITEs play a dominant role in controlling DNA methylation levels. MITEs appear to be the main sources of 24-nt siRNAs that regulate gene expression through the RNA-directed DNA methylation pathway.

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.

Distinct methylome patterns contribute to ecotypic differentiation in the growth of the storage organ of a flowering plant (sacred lotus)

DNA methylation is an epigenetic modification involved in phenotypic diversity, plant development, and environmental responses. However, the mechanisms of DNA methylation underpinning the adaption of lotus (Nelumbo nucifera) ecotypes to high and low latitudes remain unsolved, especially adaptive evolution of their storage organs. Tropical and temperate lotus ecotypes have thin and enlarged rhizomes which are adapted to low and high latitudes, respectively. Here, we investigated the DNA methylomes and transcriptomes of rhizomes of the temperate and tropical lotus to address this issue. Compared with that of the tropical lotus, the DNA of the temperate lotus was significantly more hypermethylated, indicating an increase in global DNA methylation in the lotus, with rhizome enlargement. Meanwhile, genes associated with differentially methylated regions in their promoters tended to be differentially expressed between the two ecotypes. Interestingly, the genes with their expression negatively correlated with methylation levels in their promoters and genomic regions displayed significantly higher transposon coverage, while the genes showing a significant positive correlation between expression and methylation showed lesser transposon coverage. Further, we identified that DNA methylation, especially in the promoter region, was significantly correlated with the expression of many starch-biosynthetic, gibberellin-, and brassinosteroid-signalling genes associated with rhizome differentiation. Overall, our study unveiled that distinct global and local methylation patterns between the two lotus ecotypes contribute to their expression differences and adaptive phenotypic divergence of their storage organs, highlighting the role of DNA methylation in shaping the ecotypic differentiation of lotus.

DNA Modification Patterns within the Transposable Elements of the Fig (Ficus carica L.) Genome

Transposable element activity can be harmful to the host’s genome integrity, but it can also provide selective advantages. One strategy to cope with transposons is epigenetic control through DNA base modifications. We report the non-canonic DNA modification dynamics of fig (Ficus carica L.) by exploiting high-quality genome reference and related N4-methylcytosine (4mC) and N6-methyladenine (6mA) data. Overall, 1.49% of transposon nucleotides showed either 4mC or 6mA modifications: the 4mC/6mA ratio was similar in Class I and Class II transposons, with a prevalence of 4mC, which is comparable to coding genes. Different percentages of 4mC or 6mA were observed among LTR-retrotransposon lineages and sub-lineages. Furthermore, both the Copia and Gypsy retroelements showed higher modification rates in the LTR and coding regions compared with their neighbour regions. Finally, the unconventional methylation of retrotransposons is unrelated to the number of close genes, suggesting that the 4mC and 6mA frequency in LTR-retrotransposons should not be related to transcriptional repression in the adjacency of the element. In conclusion, this study highlighted unconventional DNA modification patterns in fig transposable elements. Further investigations will focus on functional implications, in regards to how modified retroelements affect the expression of neighbouring genes, and whether these epigenetic markers can spread from repeats to genes, shaping the plant phenotype.

Genome-wide DNA hypomethylation shapes nematode pattern-triggered immunity in plants

A role for DNA hypomethylation has recently been suggested in the interaction between bacteria and plants; it is unclear whether this phenomenon reflects a conserved response. Treatment of plants of monocot rice and dicot tomato with nematode-associated molecular patterns from different nematode species or bacterial pathogen-associated molecular pattern flg22 revealed global DNA hypomethylation. A similar hypomethylation response was observed during early gall induction by Meloidogyne graminicola in rice. Evidence for the causal impact of hypomethylation on immunity was revealed by a significantly reduced plant susceptibility upon treatment with DNA methylation inhibitor 5-azacytidine. Whole-genome bisulphite sequencing of young galls revealed massive hypomethylation in the CHH context, while not for CG or CHG nucleotide contexts. Further, CHH hypomethylated regions were predominantly associated with gene promoter regions, which was not correlated with activated gene expression at the same time point but, rather, was correlated with a delayed transcriptional gene activation. Finally, the relevance of CHH hypomethylation in plant defence was confirmed in rice mutants of the RNA-directed DNA methylation pathway and DECREASED DNA METHYLATION 1. We demonstrated that DNA hypomethylation is associated with reduced susceptibility in rice towards root-parasitic nematodes and is likely to be part of the basal pattern-triggered immunity response in plants.

Single-base cytosine methylation analysis in fruits of three Capsicum species

Single-base cytosine methylation analysis across fruits of Capsicum annuum, C. chinense and C. frutescens showed global average methylation ranging from 82.8-89.1%, 77.6-83.9%, and 22.4-25% at CG, CHG and CHH contexts, respectively. High gene-body methylation at CG and CHG was observed across Capsicum species. The C. annuum showed the highest proportion (>80%) of mCs at different genomic regions compared to C. chinense and C. frutescens. Cytosine methylation for transposable-elements were lower in C. frutescens compared to C. annuum and C. chinense. A total of 510,165 CG, 583112 CHG and 277,897 CHH DMRs were identified across three Capsicum species. The differentially methylated regions (DMRs) distribution analysis revealed C. frutescens as more hypo-methylated compared to C. annuum and C. chinense, and also the presence of more intergenic DMRs in Capsicum genome. At CG and CHG context, gene expression and promoter methylation showed inverse correlations. Furthermore, the observed correlation between methylation and expression of genes suggested the potential role of methylation in Capsicum fruit development/ripening.

Precise evaluation of tissue culture-induced variation during optimisation of in vitro regeneration regime in barley

KEY MESSAGE

The Taguchi method and metAFLP analysis were used to optimise barley regenerants towards maximum and minimum levels of tissue culture-induced variation. The subtle effects of symmetric and asymmetric methylation changes in regenerants were identified. Plant tissue cultures (PTCs) provide researchers with unique materials that accelerate the development of new breeding cultivars and facilitate studies on off-type regenerants. The emerging variability of regenerants derived from PTCs may have both genetic and epigenetic origins, and may be desirable or degrade the value of regenerated plants. Thus, it is crucial to determine how the PTC variation level can be controlled. The easiest way to manipulate total tissue culture-induced variation (TTCIV) is to utilise appropriate stress factors and suitable medium components. This study describes the optimisation of in vitro tissue culture-induced variation in plant regenerants derived from barley anther culture, and maximizes and minimizes regenerant variation compared with the source explants. The approach relied on methylation amplified fragment length polymorphism (metAFLP)-derived TTCIV characteristics, which were evaluated in regenerants derived under distinct tissue culture conditions and analysed via Taguchi statistics. The factors that may trigger TTCIV included CuSO4, AgNO3 and the total time spent on the induction medium. The donor plants prepared for regeneration purposes had 5.75% and 2.01% polymorphic metAFLP loci with methylation and sequence changes, respectively. The level of TTCIV (as the sum of all metAFLP characteristics analyzed) identified in optimisation and verification experiments reached 7.51 and 10.46%, respectively. In the trial designed to produce a minimum number of differences between donor and regenerant plants, CuSO4 and AgNO3 were more crucial than time, which was not a significant factor. In the trial designed to produce a maximum number of differences between donor and regenerant plants, all factors had comparable impact on variation. The Taguchi method reduced the time required for experimental trials compared with a grid method and suggested that medium modifications were required to control regenerant variation. Finally, the effects of symmetric and asymmetric methylation changes on regenerants were identified using novel aspects of the metAFLP method developed for this analysis.

Dynamics and function of DNA methylation in plants

DNA methylation is a conserved epigenetic modification that is important for gene regulation and genome stability. Aberrant patterns of DNA methylation can lead to plant developmental abnormalities. A specific DNA methylation state is an outcome of dynamic regulation by de novo methylation, maintenance of methylation and active demethylation, which are catalysed by various enzymes that are targeted by distinct regulatory pathways. In this Review, we discuss DNA methylation in plants, including methylating and demethylating enzymes and regulatory factors, and the coordination of methylation and demethylation activities by a so-called methylstat mechanism; the functions of DNA methylation in regulating transposon silencing, gene expression and chromosome interactions; the roles of DNA methylation in plant development; and the involvement of DNA methylation in plant responses to biotic and abiotic stress conditions.

Development of Beet necrotic yellow vein virus-based vectors for multiple-gene expression and guide RNA delivery in plant genome editing

Many plant viruses with monopartite or bipartite genomes have been developed as efficient expression vectors of foreign recombinant proteins. Nonetheless, due to lack of multiple insertion sites in these plant viruses, it is still a big challenge to simultaneously express multiple foreign proteins in single cells. The genome of Beet necrotic yellow vein virus (BNYVV) offers an attractive system for expression of multiple foreign proteins owning to a multipartite genome composed of five positive-stranded RNAs. Here, we have established a BNYVV full-length infectious cDNA clone under the control of the Cauliflower mosaic virus 35S promoter. We further developed a set of BNYVV-based vectors that permit efficient expression of four recombinant proteins, including some large proteins with lengths up to 880 amino acids in the model plant Nicotiana benthamiana and native host sugar beet plants. These vectors can be used to investigate the subcellular co-localization of multiple proteins in leaf, root and stem tissues of systemically infected plants. Moreover, the BNYVV-based vectors were used to deliver NbPDS guide RNAs for genome editing in transgenic plants expressing Cas9, which induced a photobleached phenotype in systemically infected leaves. Collectively, the BNYVV-based vectors will facilitate genomic research and expression of multiple proteins, in sugar beet and related crop plants.

Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

BACKGROUND

Retrotransposons are the major determinants of genome sizes and they have shaped both genes and genomes in mammalian organisms, but their overall activity, diversity, and evolution dynamics, particularly their impact on protein coding and lncRNA genes in pigs remain largely unknown.

RESULTS

In the present study, we performed de novo detection of retrotransposons in pigs by using multiple pipelines, four distinct families of pig-specific L1 s classified into 51 distinct subfamilies and representing four evolution models and three expansion waves of pig-specific SINEs represented by three distinct families were identified. ERVs were classified into 18 families and found two most "modern" subfamilies in the pig genome. The transposition activity of pig L1 was verified by experiment, the sense and antisense promoter activities of young L1 5'UTRs and ERV LTRs and expression profiles of young retrotransposons in multiple tissues and cell lines were also validated. Furthermore, retrotransposons had an extensive impact on lncRNA and protein coding genes at both the genomic and transcriptomic levels. Most protein coding and lncRNA (> 80%) genes contained retrotransposon insertions, and about half of protein coding genes (44.30%) and one-fourth (24.13%) of lncRNA genes contained the youngest retrotransposon insertions. Nearly half of protein coding genes (43.78%) could generate chimeric transcripts with retrotransposons. Significant distribution bias of retrotransposon composition, location, and orientation in lncRNA and protein coding genes, and their transcripts, were observed.

CONCLUSIONS

In the current study, we characterized the classification and evolution profile of retrotransposons in pigs, experimentally proved the transposition activity of the young pig L1 subfamily, characterized the sense and antisense expression profiles and promoter activities of young retrotransposons, and investigated their impact on lncRNA and protein coding genes by defining the mobilome landscapes at the genomic and transcriptomic levels. These findings help provide a better understanding of retrotransposon evolution in mammal and their impact on the genome and transcriptome.

Multi-omics analyses reveal epigenomics basis for cotton somatic embryogenesis through successive regeneration acclimation process

Plant regeneration via somatic embryogenesis is time-consuming and highly genotype-dependent. The plant somatic embryogenesis process provokes many epigenetics changes including DNA methylation and histone modification. Recently, an elite cotton Jin668, with an extremely high regeneration ability, was developed from its maternal inbred Y668 cultivar using a Successive Regeneration Acclimation (SRA) strategy. To reveal the underlying mechanism of SRA, we carried out a genome-wide single-base resolution methylation analysis for nonembryogenic calluses (NECs), ECs, somatic embryos (SEs) during the somatic embryogenesis procedure and the leaves of regenerated offspring plants. Jin668 (R4) regenerated plants were CHH hypomethylated compared with the R0 regenerated plants of SRA process. The increase in CHH methylation from NEC to EC was demonstrated to be associated with the RNA-dependent DNA methylation (RdDM) and the H3K9me2-dependent pathway. Intriguingly, the hypomethylated CHH differentially methylated regions (DMRs) of promoter activated some hormone-related and WUSCHEL-related homeobox genes during the somatic embryogenesis process. Inhibiting DNA methylation using zebularine treatment in NEC increased the number of embryos. Our multi-omics data provide new insights into the dynamics of DNA methylation during the plant tissue culture and regenerated offspring plants. This study also reveals that induced hypomethylation (SRA) may facilitate the higher plant regeneration ability and optimize maternal genetic cultivar.

How to Design a Whole-Genome Bisulfite Sequencing Experiment

Aside from post-translational histone modifications and small RNA populations, the epigenome of an organism is defined by the level and spectrum of DNA methylation. Methyl groups can be covalently bound to the carbon-5 of cytosines or the carbon-6 of adenine bases. DNA methylation can be found in both prokaryotes and eukaryotes. In the latter, dynamic variation is shown across species, along development, and by cell type. DNA methylation usually leads to a lower binding affinity of DNA-interacting proteins and often results in a lower expression rate of the subsequent genome region, a process also referred to as transcriptional gene silencing. We give an overview of the current state of research facilitating the planning and implementation of whole-genome bisulfite-sequencing (WGBS) experiments. We refrain from discussing alternative methods for DNA methylation analysis, such as reduced representation bisulfite sequencing (rrBS) and methylated DNA immunoprecipitation sequencing (MeDIPSeq), which have value in specific experimental contexts but are generally disadvantageous compared to WGBS.

Unravelling the epigenomic interactions between parental inbreds resulting in an altered hybrid methylome in pigeonpea

DNA methylation is an important heritable landmark conferring epigenetic changes in hybrids and has fascinated biologists and plant-breeders over the years. Although epigenetic changes have been documented in rice and maize hybrids, such investigations have not been reported in pigeonpea. Here, we report genome-wide methylation profiles of pigeonpea sterile and fertile inbred lines and their fertile F1 hybrid at single base resolution. We found that pigeonpea genome is relatively enriched in CG methylation. Identification of differentially methylated regions (DMRs) in the sterile and fertile parent revealed remarkable differences between their methylation patterns. Investigation of methylation status of parental DMRs in hybrid revealed non-additive methylation patterns resulting from trans-chromosomal methylation and trans-chromosomal demethylation events. Furthermore, we discovered several DMRs negatively associated with gene expression in the hybrid and fertile parent. Interestingly, many of those DMRs belonged to transposable elements and genes encoding pentatricopeptide repeats associated proteins, which may mediate a role in modulating the genes impacting pollen fertility. Overall, our findings provide an understanding of two parental epigenomes interacting to give rise to an altered methylome in pigeonpea hybrids, from genome-wide point of view.

DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia

Background

S. latifolia is a model organism for the study of sex chromosome evolution in plants. Its sex chromosomes include large regions in which recombination became gradually suppressed. The regions tend to expand over time resulting in the formation of evolutionary strata. Non-recombination and later accumulation of repetitive sequences is a putative cause of the size increase in the Y chromosome. Gene decay and accumulation of repetitive DNA are identified as key evolutionary events. Transposons in the X and Y chromosomes are distributed differently and there is a regulation of transposon insertion by DNA methylation of the target sequences, this points to an important role of DNA methylation during sex chromosome evolution in Silene latifolia. The aim of this study was to elucidate whether the reduced expression of the Y allele in S. latifolia is caused by genetic degeneration or if the cause is methylation triggered by transposons and repetitive sequences.

Results

Gene expression analysis in S. latifolia males has shown expression bias in both X and Y alleles. To determine whether these differences are caused by genetic degeneration or methylation spread by transposons and repetitive sequences, we selected several sex-linked genes with varying degrees of degeneration and from different evolutionary strata. Immunoprecipitation of methylated DNA (MeDIP) from promoter, exon and intron regions was used and validated through bisulfite sequencing. We found DNA methylation in males, and only in the promoter of genes of stratum I (older). The Y alleles in genes of stratum I were methylation enriched compared to X alleles. There was also abundant and high percentage methylation in the CHH context in most sequences, indicating de novo methylation through the RdDM pathway.

Conclusions

We speculate that TE accumulation and not gene decay is the cause of DNA methylation in the S. latifolia Y sex chromosome with influence on the process of heterochromatinization.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4936-y) contains supplementary material, which is available to authorized users.

Dynamic Epigenetic Changes during Plant Regeneration

Plants have the remarkable ability to drive cellular dedifferentiation and regeneration. Changes in epigenetic landscapes accompany the cell fate transition. Notably, modifications of chromatin structure occur primarily during callus formation via an in vitro tissue culture process and, thus, pluripotent callus cells have unique epigenetic signatures. Here, we highlight the latest progress in epigenetic regulation of callus formation in plants, which addresses fundamental questions related to cell fate changes and pluripotency establishment. Global and local modifications of chromatin structure underlie callus formation, and the combination and sequence of epigenetic modifications further shape intricate cell fate changes. This review illustrates how a series of chromatin marks change dynamically during callus formation and their biological relevance in plant regeneration.

Use of retrotransposon-derived genetic markers to analyse genomic variability in plants

Transposable elements (TEs) are common mobile genetic elements comprising several classes and making up the majority of eukaryotic genomes. The movement and accumulation of TEs has been a major force shaping the genes and genomes of most organisms. Most eukaryotic genomes are dominated by retrotransposons and minimal DNA transposon accumulation. The 'copy and paste' lifecycle of replicative transposition produces new genome insertions without excising the original element. Horizontal TE transfer among lineages is rare. TEs represent a reservoir of potential genomic instability and RNA-level toxicity. Many TEs appear static and nonfunctional, but some are capable of replicating and mobilising to new positions, and somatic transposition events have been observed. The overall structure of retrotransposons and the domains responsible for the phases of their replication are highly conserved in all eukaryotes. TEs are important drivers of species diversity and exhibit great variety in their structure, size and transposition mechanisms, making them important putative actors in evolution. Because TEs are abundant in plant genomes, various applications have been developed to exploit polymorphisms in TE insertion patterns, including conventional or anchored PCR, and quantitative or digital PCR with primers for the 5' or 3' junction. Alternatively, the retrotransposon junction can be mapped using high-throughput next-generation sequencing and bioinformatics. With these applications, TE insertions can be rapidly, easily and accurately identified, or new TE insertions can be found. This review provides an overview of the TE-based applications developed for plant species and assesses the contributions of TEs to the analysis of plants' genetic diversity.

siRNAs regulate DNA methylation and interfere with gene and lncRNA expression in the heterozygous polyploid switchgrass

BACKGROUND

Understanding the DNA methylome and its relationship with non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), is essential for elucidating the molecular mechanisms underlying key biological processes in plants. Few studies have examined the functional roles of the DNA methylome in grass species with highly heterozygous polyploid genomes.

RESULTS

We performed genome-wide DNA methylation profiling in the tetraploid switchgrass (Panicum virgatum L.) cultivar 'Alamo' using bisulfite sequencing. Single-base-resolution methylation patterns were observed in switchgrass leaf and root tissues, which allowed for characterization of the relationship between DNA methylation and mRNA, miRNA, and lncRNA populations. The results of this study revealed that siRNAs positively regulate DNA methylation of the mCHH sites surrounding genes, and that DNA methylation interferes with gene and lncRNA expression in switchgrass. Ninety-six genes covered by differentially methylated regions (DMRs) were annotated by GO analysis as being involved in stimulus-related processes. Functionally, 82% (79/96) of these genes were found to be hypomethylated in switchgrass root tissue. Sequencing analysis of lncRNAs identified two lncRNAs that are potential precursors of miRNAs, which are predicted to target genes that function in cellulose biosynthesis, stress regulation, and stem and root development.

CONCLUSIONS

This study characterized the DNA methylome in switchgrass and elucidated its relevance to gene and non-coding RNAs. These results provide valuable genomic resources and references that will aid further epigenetic research in this important biofuel crop.