Role of CG and non-CG methylation in immobilization of transposons in Arabidopsis

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.

Structural basis of water-mediated cis Watson-Crick/Hoogsteen base-pair formation in non-CpG methylation

Non-CpG methylation is associated with several cellular processes, especially neuronal development and cancer, while its effect on DNA structure remains unclear. We have determined the crystal structures of DNA duplexes containing -CGCCG- regions as CCG repeat motifs that comprise a non-CpG site with or without cytosine methylation. Crystal structure analyses have revealed that the mC:G base-pair can simultaneously form two alternative conformations arising from non-CpG methylation, including a unique water-mediated cis Watson-Crick/Hoogsteen, (w)cWH, and Watson-Crick (WC) geometries, with partial occupancies of 0.1 and 0.9, respectively. NMR studies showed that an alternative conformation of methylated mC:G base-pair at non-CpG step exhibits characteristics of cWH with a syn-guanosine conformation in solution. DNA duplexes complexed with the DNA binding drug echinomycin result in increased occupancy of the (w)cWH geometry in the methylated base-pair (from 0.1 to 0.3). Our structural results demonstrated that cytosine methylation at a non-CpG step leads to an anti→syntransition of its complementary guanosine residue toward the (w)cWH geometry as a partial population of WC, in both drug-bound and naked mC:G base pairs. This particular geometry is specific to non-CpG methylated dinucleotide sites in B-form DNA. Overall, the current study provides new insights into DNA conformation during epigenetic regulation.

Mechanism of heterochromatin remodeling revealed by the DDM1 bound nucleosome structures

The SWI/SNF2 chromatin remodeling factor decreased DNA methylation 1 (DDM1) is essential for the silencing of transposable elements (TEs) in both euchromatic and heterochromatic regions. Here, we determined the cryo-EM structures of DDM1-nucleosomeH2A and DDM1-nucleosomeH2A.W complexes at near-atomic resolution in the presence of the ATP analog ADP-BeFx. The structures show that nucleosomal DNA is unwrapped more on the surface of the histone octamer containing histone H2A than that containing histone H2A.W. DDM1 embraces one DNA gyre of the nucleosome and interacts with the N-terminal tails of histone H4. Although we did not observe DDM1-H2A.W interactions in our structures, the results of the pull-down experiments suggest a direct interaction between DDM1 and the core region of histone H2A.W. Our work provides mechanistic insights into the heterochromatin remodeling process driven by DDM1 in plants.

DeepPGD: A Deep Learning Model for DNA Methylation Prediction Using Temporal Convolution, BiLSTM, and Attention Mechanism

As part of the field of DNA methylation identification, this study tackles the challenge of enhancing recognition performance by introducing a specialized deep learning framework called DeepPGD. DNA methylation, a crucial biological modification, plays a vital role in gene expression analyses, cellular differentiation, and the study of disease progression. However, accurately and efficiently identifying DNA methylation sites remains a pivotal concern in the field of bioinformatics. The issue addressed in this paper is the presence of methylation in DNA, which is a binary classification problem. To address this, our research aimed to develop a deep learning algorithm capable of more precisely identifying these sites. The DeepPGD framework combined a dual residual structure involving Temporal convolutional networks (TCNs) and bidirectional long short-term memory (BiLSTM) networks to effectively extract intricate DNA structural and sequence features. Additionally, to meet the practical requirements of DNA methylation identification, extensive experiments were conducted across a variety of biological species. The experimental results highlighted DeepPGD's exceptional performance across multiple evaluation metrics, including accuracy, Matthews' correlation coefficient (MCC), and the area under the curve (AUC). In comparison to other algorithms in the same domain, DeepPGD demonstrated superior classification and predictive capabilities across various biological species datasets. This significant advancement in algorithmic prowess not only offers substantial technical support, but also holds potential for research and practical implementation within the DNA methylation identification domain. Moreover, the DeepPGD framework shows potential for application in genomics research, biomedicine, and disease diagnostics, among other fields.

Molecular and structural basis of the chromatin remodeling activity by Arabidopsis DDM1

The histone H2A variant H2A.W occupies transposons and thus prevents access to them in Arabidopsis thaliana. H2A.W is deposited by the chromatin remodeler DDM1, which also promotes the accessibility of chromatin writers to heterochromatin by an unknown mechanism. To shed light on this question, we solve the cryo-EM structures of nucleosomes containing H2A and H2A.W, and the DDM1-H2A.W nucleosome complex. These structures show that the DNA end flexibility of the H2A nucleosome is higher than that of the H2A.W nucleosome. In the DDM1-H2A.W nucleosome complex, DDM1 binds to the N-terminal tail of H4 and the nucleosomal DNA and increases the DNA end flexibility of H2A.W nucleosomes. Based on these biochemical and structural results, we propose that DDM1 counters the low accessibility caused by nucleosomes containing H2A.W to enable the maintenance of repressive epigenetic marks on transposons and prevent their activity.

Two ARGONAUTE proteins loaded with transposon-derived small RNAs are associated with the reproductive cell lineage in Arabidopsis

In sexually propagating organisms, genetic, and epigenetic mutations are evolutionarily relevant only if they occur in the germline and are hence transmitted to the next generation. In contrast to most animals, plants are considered to lack an early segregating germline, implying that somatic cells can contribute genetic information to progeny. Here we demonstrate that 2 ARGONAUTE proteins, AGO5 and AGO9, mark cells associated with sexual reproduction in Arabidopsis (Arabidopsis thaliana) throughout development. Both AGOs are loaded with dynamically changing small RNA populations derived from highly methylated, pericentromeric, long transposons. Sequencing of single stem cell nuclei revealed that many of these transposons are co-expressed within an AGO5/9 expression domain in the shoot apical meristem (SAM). Co-occurrence of transposon expression and specific ARGONAUTE (AGO) expression in the SAM is reminiscent of germline features in animals and supports the existence of an early segregating germline in plants. Our results open the path to investigating transposon biology and epigenome dynamics at cellular resolution in the SAM stem cell niche.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

Coevolution of the CDCA7-HELLS ICF-related nucleosome remodeling complex and DNA methyltransferases

5-Methylcytosine (5mC) and DNA methyltransferases (DNMTs) are broadly conserved in eukaryotes but are also frequently lost during evolution. The mammalian SNF2 family ATPase HELLS and its plant ortholog DDM1 are critical for maintaining 5mC. Mutations in HELLS, its activator CDCA7, and the de novo DNA methyltransferase DNMT3B, cause immunodeficiency-centromeric instability-facial anomalies (ICF) syndrome, a genetic disorder associated with the loss of DNA methylation. We here examine the coevolution of CDCA7, HELLS and DNMTs. While DNMT3, the maintenance DNA methyltransferase DNMT1, HELLS, and CDCA7 are all highly conserved in vertebrates and green plants, they are frequently co-lost in other evolutionary clades. The presence-absence patterns of these genes are not random; almost all CDCA7 harboring eukaryote species also have HELLS and DNMT1 (or another maintenance methyltransferase, DNMT5). Coevolution of presence-absence patterns (CoPAP) analysis in Ecdysozoa further indicates coevolutionary linkages among CDCA7, HELLS, DNMT1 and its activator UHRF1. We hypothesize that CDCA7 becomes dispensable in species that lost HELLS or DNA methylation, and/or the loss of CDCA7 triggers the replacement of DNA methylation by other chromatin regulation mechanisms. Our study suggests that a unique specialized role of CDCA7 in HELLS-dependent DNA methylation maintenance is broadly inherited from the last eukaryotic common ancestor.

The plant response to high CO2 levels is heritable and orchestrated by DNA methylation

Plant responses to abiotic environmental challenges are known to have lasting effects on the plant beyond the initial stress exposure. Some of these lasting effects are transgenerational, affecting the next generation. The plant response to elevated carbon dioxide (CO₂ ) levels has been well studied. However, these investigations are typically limited to plants grown for a single generation in a high CO₂ environment while transgenerational studies are rare. We aimed to determine transgenerational growth responses in plants after exposure to high CO₂ by investigating the direct progeny when returned to baseline CO₂ levels. We found that both the flowering plant Arabidopsis thaliana and seedless nonvascular plant Physcomitrium patens continue to display accelerated growth rates in the progeny of plants exposed to high CO₂ . We used the model species Arabidopsis to dissect the molecular mechanism and found that DNA methylation pathways are necessary for heritability of this growth response. More specifically, the pathway of RNA-directed DNA methylation is required to initiate methylation and the proteins CMT2 and CMT3 are needed for the transgenerational propagation of this DNA methylation to the progeny plants. Together, these two DNA methylation pathways establish and then maintain a cellular memory to high CO₂ exposure.

The model diatom Phaeodactylum tricornutum provides insights into the diversity and function of microeukaryotic DNA methyltransferases

Cytosine methylation is an important epigenetic mark involved in the transcriptional control of transposable elements in mammals, plants and fungi. The Stramenopiles-Alveolate-Rhizaria (SAR) lineages are a major group of ecologically important marine microeukaryotes, including the phytoplankton groups diatoms and dinoflagellates. However, little is known about their DNA methyltransferase diversity. Here, we performed an in-silico analysis of DNA methyltransferases found in marine microeukaryotes and showed that they encode divergent DNMT3, DNMT4, DNMT5 and DNMT6 enzymes. Furthermore, we found three classes of enzymes within the DNMT5 family. Using a CRISPR/Cas9 strategy we demonstrated that the loss of the DNMT5a gene correlates with a global depletion of DNA methylation and overexpression of young transposable elements in the model diatom Phaeodactylum tricornutum. The study provides a view of the structure and function of a DNMT family in the SAR supergroup using an attractive model species.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Seminars in cell and development biology on histone variants remodelers of H2A variants associated with heterochromatin

Variants of the histone H2A occupy distinct locations in the genome. There is relatively little known about the mechanisms responsible for deposition of specific H2A variants. Notable exceptions are chromatin remodelers that control the dynamics of H2A.Z at promoters. Here we review the steps that identified the role of a specific class of chromatin remodelers, including LSH and DDM1 that deposit the variants macroH2A in mammals and H2A.W in plants, respectively. The function of these remodelers in heterochromatin is discussed together with their multiple roles in genome stability.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Parental methylation mediates how progeny respond to environments of parents and of progeny themselves

BACKGROUND AND AIMS

Environments experienced by both parents and offspring influence progeny traits, but the epigenetic mechanisms that regulate the balance of parental vs. progeny control of progeny phenotypes are not known. We tested whether DNA methylation in parents and/or progeny mediates responses to environmental cues experienced in both generations.

METHODS

Using Arabidopsis thaliana, we manipulated parental and progeny DNA methylation both chemically, via 5-azacytidine, and genetically, via mutants of methyltransferase genes, then measured progeny germination responses to simulated canopy shade in parental and progeny generations.

KEY RESULTS

We first found that germination of offspring responded to parental but not seed demethylation. We further found that parental demethylation reversed the parental effect of canopy in seeds with low (Cvi-1) to intermediate (Col) dormancy, but it obliterated the parental effect in seeds with high dormancy (Cvi-0). Demethylation did so by either suppressing germination of seeds matured under white-light (Cvi-1) or under canopy (Cvi-0), or by increasing the germination of seeds matured under canopy (Col). Disruption of parental methylation also prevented seeds from responding to their own light environment in one genotype (Cvi-0, most dormant), but it enabled seeds to respond to their own environment in another genotype (Cvi-1, least dormant). Using mutant genotypes, we found that both CG and non-CG DNA methylation were involved in parental effects on seed germination.

CONCLUSIONS

Parental methylation state influences seed germination more strongly than does the progeny's own methylation state, and it influences how seeds respond to environments of parents and progeny in a genotype-specific manner.

Heritable Epigenomic Modifications Influence Stress Resilience and Rapid Adaptations in the Brown Planthopper (Nilaparvata lugens)

DNA methylation in insects is integral to cellular differentiation, development, gene regulation, genome integrity, and phenotypic plasticity. However, its evolutionary potential and involvement in facilitating rapid adaptations in insects are enigmatic. Moreover, our understanding of these mechanisms is limited to a few insect species, of which none are pests of crops. Hence, we studied methylation patterns in the brown planthopper (BPH), a major rice pest, under pesticide and nutritional stress, across its life stages. Moreover, as the inheritance of epigenetic changes is fundamentally essential for acclimation, adaptability, and evolution, we determined the heritability and persistence of stress-induced methylation marks in BPH across generations. Our results revealed that DNA methylation pattern(s) in BPH varies/vary with environmental cues and is/are insect life-stage specific. Further, our findings provide novel insights into the heritability of stress-induced methylation marks in BPH. However, it was observed that, though heritable, these marks eventually fade in the absence of the stressors, thereby suggesting the existence of fitness cost(s) associated with the maintenance of the stressed epigenotype. Furthermore, we demonstrate how 5-azacytidine-mediated disruption of BPH methylome influences expression levels of stress-responsive genes and, thereby, highlight demethylation/methylation as a phenomenon underlying stress resilience of BPH.

DNA methylation-free Arabidopsis reveals crucial roles of DNA methylation in regulating gene expression and development

A contribution of DNA methylation to defense against invading nucleic acids and maintenance of genome integrity is uncontested; however, our understanding of the extent of involvement of this epigenetic mark in genome-wide gene regulation and plant developmental control is incomplete. Here, we knock out all five known DNA methyltransferases in Arabidopsis, generating DNA methylation-free plants. This quintuple mutant exhibits a suite of developmental defects, unequivocally demonstrating that DNA methylation is essential for multiple aspects of plant development. We show that CG methylation and non-CG methylation are required for a plethora of biological processes, including pavement cell shape, endoreduplication, cell death, flowering, trichome morphology, vasculature and meristem development, and root cell fate determination. Moreover, we find that DNA methylation has a strong dose-dependent effect on gene expression and repression of transposable elements. Taken together, our results demonstrate that DNA methylation is dispensable for Arabidopsis survival but essential for the proper regulation of multiple biological processes.

Our understanding of the extent of involvement of DNA methylation in genome-wide gene regulation and plant developmental control is incomplete. Here, the authors knock out all five known DNA methyltransferases and show the developmental and gene expression changes in the DNA methylation-free Arabidopsis plants.

N4-acetyldeoxycytosine DNA modification marks euchromatin regions in Arabidopsis thaliana

BACKGROUND

Direct analogs of chemically modified bases that carry important epigenetic information, such as 5-methylcytosine (m5C)/5-methyldeoxycytosine (5mC), 5-hydroxymethylcytosine (hm5C)/5-hydroxymethyldeoxycytosine (5hmC), and N6-methyladenosine (m6A)/N6-methyldeoxyadenosine (6mA), are detected in both RNA and DNA, respectively. The modified base N4-acetylcytosine (ac4C) is well studied in RNAs, but its presence and epigenetic roles in cellular DNA have not been explored.

RESULTS

Here, we demonstrate the existence of N4-acetyldeoxycytosine (4acC) in genomic DNA of Arabidopsis with multiple detection methods. Genome-wide profiling of 4acC modification reveals that 4acC peaks are mostly distributed in euchromatin regions and present in nearly half of the expressed protein-coding genes in Arabidopsis. 4acC is mainly located around transcription start sites and positively correlates with gene expression levels. Imbalance of 5mC does not directly affect 4acC modification. We also characterize the associations of 4acC with 5mC and histone modifications that cooperatively regulate gene expression. Moreover, 4acC is also detected in genomic DNA of rice, maize, mouse, and human by mass spectrometry.

CONCLUSIONS

Our findings reveal 4acC as a hitherto unknown DNA modification in higher eukaryotes. We identify potential interactions of this mark with other epigenetic marks in gene expression regulation.

Roles of IDM3 and SDJ1/2/3 in Establishment and/or Maintenance of DNA Methylation in Arabidopsis

Previous studies had demonstrated that in Arabidopsis, IDM3 is involved in ROS1-mediated DNA demethylation pathway, and SUVH-SDJ complex functions as a DNA methylation reader complex for enhancing gene transcription, which presumably recruits ROS1 to the promoters of target genes for DNA demethylation. Here, our analyses, however, showed that the IDM3 and SDJ1/2/3, the components of the SUVH-SDJ complex, are implicated in establishing and/or maintaining DNA methylation as well through DDR (DRD1-DMS3-RDM1) complex. idm3-3 or sdj1/2/3 mutations led to genome-wide DNA hypomethylation, and both mutants shared a large number of common hypo-DMRs (Differentially Methylated Regions) with rdm1-4 and dms3-4, suggesting that IDM3 and SDJ1/2/3 help establish and/or maintain DNA methylation, mediated by RdDM pathway, at a subset of genomic regions largely through DDR complex. IDM3 is able to strongly interact with RDM1 and DMS3, but weakly with SDJ1 and SDJ3; SDJ1 and SDJ3 is capable of interacting separately with RDM1 and DMS3. Furthermore, comparisons of DNA methylation features in idm3-3 and sdj1/2/3 indicated that idm3-3 and sdj1/2/3 mutations make differential impacts on DNA methylation levels and patterns on a genome-wide scale, indicating that they are targeted to quite distinct genomic regions for aiding in DNA methylation. Further analyses on ChIP-seq data demonstrated that RDM1, DMS3 and NRPE1 are enriched in IDM3- and SDJ1/2/3-targted regions. Altogether, our results provide clear demonstration that IDM3 and SDJ1/2/3 play a part in establishing and/or maintaining DNA methylation of a group of genomic regions, through the DDR complex.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

Comprehensive Mechanism of Gene Silencing and Its Role in Plant Growth and Development

Gene silencing is a negative feedback mechanism that regulates gene expression to define cell fate and also regulates metabolism and gene expression throughout the life of an organism. In plants, gene silencing occurs via transcriptional gene silencing (TGS) and post-transcriptional gene silencing (PTGS). TGS obscures transcription via the methylation of 5' untranslated region (5'UTR), whereas PTGS causes the methylation of a coding region to result in transcript degradation. In this review, we summarized the history and molecular mechanisms of gene silencing and underlined its specific role in plant growth and crop production.

Repression of CHROMOMETHYLASE 3 prevents epigenetic collateral damage in Arabidopsis

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

The chromatin remodeler DDM1 prevents transposon mobility through deposition of histone variant H2A.W

Mobile transposable elements (TEs) not only participate in genome evolution but also threaten genome integrity. In healthy cells, TEs that encode all of the components that are necessary for their mobility are specifically silenced, yet the precise mechanism remains unknown. Here, we characterize the mechanism used by a conserved class of chromatin remodelers that prevent TE mobility. In the Arabidopsis chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1), we identify two conserved binding domains for the histone variant H2A.W, which marks plant heterochromatin. DDM1 is necessary and sufficient for the deposition of H2A.W onto potentially mobile TEs, yet does not act on TE fragments or host protein-coding genes. DDM1-mediated H2A.W deposition changes the properties of chromatin, resulting in the silencing of TEs and, therefore, prevents their mobility. This distinct mechanism provides insights into the interplay between TEs and their host in the contexts of evolution and disease, and potentiates innovative strategies for targeted gene silencing.

Distinct Responses of Arabidopsis Telomeres and Transposable Elements to Zebularine Exposure

Involvement of epigenetic mechanisms in the regulation of telomeres and transposable elements (TEs), genomic regions with the protective and potentially detrimental function, respectively, has been frequently studied. Here, we analyzed telomere lengths in Arabidopsis thaliana plants of Columbia, Landsberg erecta and Wassilevskija ecotypes exposed repeatedly to the hypomethylation drug zebularine during germination. Shorter telomeres were detected in plants growing from seedlings germinated in the presence of zebularine with a progression in telomeric phenotype across generations, relatively high inter-individual variability, and diverse responses among ecotypes. Interestingly, the extent of telomere shortening in zebularine Columbia and Wassilevskija plants corresponded to the transcriptional activation of TEs, suggesting a correlated response of these genomic elements to the zebularine treatment. Changes in lengths of telomeres and levels of TE transcripts in leaves were not always correlated with a hypomethylation of cytosines located in these regions, indicating a cytosine methylation-independent level of their regulation. These observations, including differences among ecotypes together with distinct dynamics of the reversal of the disruption of telomere homeostasis and TEs transcriptional activation, reflect a complex involvement of epigenetic processes in the regulation of crucial genomic regions. Our results further demonstrate the ability of plant cells to cope with these changes without a critical loss of the genome stability.

RNA interference-independent reprogramming of DNA methylation in Arabidopsis

DNA methylation is important for silencing transposable elements (TEs) in diverse eukaryotes, including plants. In plant genomes, TEs are silenced by methylation of histone H3 lysine 9 (H3K9) and cytosines in both CG and non-CG contexts. The role of RNA interference (RNAi) in establishing TE-specific silent marks has been extensively studied, but the importance of RNAi-independent pathways remains largely unexplored. Here, we directly investigated transgenerational de novo DNA methylation of TEs after the loss of silent marks. Our analyses uncovered potent and precise RNAi-independent pathways for recovering non-CG methylation and H3K9 methylation in most TE genes (that is, coding regions within TEs). Characterization of a subset of TE genes without the recovery revealed the effects of H3K9 demethylation, replacement of histone H2A variants and their interaction with CG methylation, together with feedback from transcription. These chromatin components are conserved among eukaryotes and may contribute to chromatin reprogramming in a conserved manner.

SWI3B and HDA6 interact and are required for transposon silencing in Arabidopsis

Although the interplay of covalent histone acetylation/deacetylation and ATP-dependent chromatin remodelling is crucial for the regulation of chromatin structure and gene expression in eukaryotes, the underlying molecular mechanism in plants remains largely unclear. Here we show a direct interaction between Arabidopsis SWI3B, an essential subunit of the SWI/SNF chromatin-remodelling complex, and the RPD3/HDA1-type histone deacetylase HDA6 both in vitro and in vivo. Furthermore, SWI3B and HDA6 co-repress the transcription of a subset of transposons. Both SWI3B and HDA6 maintain transposon silencing by decreasing histone H3 lysine 9 acetylation, but increasing histone H3 lysine 9 di-methylation, DNA methylation and nucleosome occupancy. Our findings reveal that SWI3B and HDA6 may act in the same co-repressor complex to maintain transposon silencing in Arabidopsis.

Transgenerational acclimation to changes in ocean acidification in marine invertebrates

The rapid pace of increasing oceanic acidity poses a major threat to the fitness of the marine ecosystem, as well as the buffering capacity of the oceans. Disruption in chemical equilibrium in the ocean leads to decreased carbonate ion precipitation, resulting in calcium carbonate saturation. If these trends continue, calcifying invertebrates will experience difficultly maintaining their calcium carbonate exoskeleton and shells. Because malfunction of exoskeleton formation by calcifiers in response to ocean acidification (OA) will have non-canonical biological cascading results in the marine ecosystem, many studies have investigated the direct and indirect consequences of OA on ecosystem- and physiology-related traits of marine invertebrates. Considering that evolutionary adaptation to OA depends on the duration of OA effects, long-term exposure to OA stress over multi-generations may result in adaptive mechanisms that increase the potential fitness of marine invertebrates in response to OA. Transgenerational studies have the potential to elucidate the roles of acclimation, carryover effects, and evolutionary adaptation within and over generations in response to OA. In particular, understanding mechanisms of transgenerational responses (e.g., antioxidant responses, metabolic changes, epigenetic reprogramming) to changes in OA will enhance our understanding of marine invertebrate in response to rapid climate change.

Transposon age and non-CG methylation

Silencing of transposable elements (TEs) is established by small RNA-directed DNA methylation (RdDM). Maintenance of silencing is then based on a combination of RdDM and RNA-independent mechanisms involving DNA methyltransferase MET1 and chromodomain DNA methyltransferases (CMTs). Involvement of RdDM, according to this model should decrease with TE age but here we show a different pattern in tomato and Arabidopsis. In these species the CMTs silence long terminal repeat (LTR) transposons in the distal chromatin that are younger than those affected by RdDM. To account for these findings we propose that, after establishment of primary RdDM as in the original model, there is an RNA-independent maintenance phase involving CMTs followed by secondary RdDM. This progression of epigenetic silencing in the gene-rich distal chromatin is likely to influence the transcriptome either in cis or in trans depending on whether the mechanisms are RNA-dependent or -independent.

RNA-directed DNA methylation (RdDM) is thought to silence newly inserted transposable elements (TEs) with RNA-independent mechanisms becoming more prominent as TEs age. Here, the authors show that RdDM continues to silence the oldest intact distal TEs in tomato and Arabidopsis suggesting a second, later phase of RdDM.

CG hypomethylation leads to complex changes in DNA methylation and transpositional burst of diverse transposable elements in callus cultures of rice

CG methylation (^m CG) is essential for preserving genome stability in mammals, but this link remains obscure in plants. OsMET1-2, a major rice DNA methyltransferase, plays critical roles in maintaining ^m CG in rice. Null mutation of OsMET1-2 causes massive CG hypomethylation, rendering the mutant suitable to address the role of ^m CG in maintaining genome integrity in plants. Here, we analyzed ^m CG dynamics and genome stability in tissue cultures of OsMET1-2 homozygous (-/-) and heterozygous (+/-) mutants, and isogenic wild-type (WT). We found ^m CG levels in cultures of -/- were substantially lower than in those of WT and +/-, as expected. Unexpectedly, ^m CG levels in 1- and 3-year cultures of -/- were 77.6% and 48.7% higher, respectively, than in shoot, from which the cultures were initiated, suggesting substantial regain of ^m CG in -/- cultures, which contrasts to the general trend of ^m CG loss in all WT plant tissue cultures hitherto studied. Transpositional burst of diverse transposable elements (TEs) occurred only in -/- cultures, although no elevation of genome-wide mutation rate in the form of single nucleotide polymorphisms was detected. Altogether, our results establish an essential role of ^m CG in retaining TE immobility and hence genome stability in rice and likely in plants in general.

Dynamic Patterns of Transcript Abundance of Transposable Element Families in Maize

Transposable Elements (TEs) are mobile elements that contribute the majority of DNA sequences in the maize genome. Due to their repetitive nature, genomic studies of TEs are complicated by the difficulty of properly attributing multi-mapped short reads to specific genomic loci. Here, we utilize a method to attribute RNA-seq reads to TE families rather than particular loci in order to characterize transcript abundance for TE families in the maize genome. We applied this method to assess per-family expression of transposable elements in >800 published RNA-seq libraries representing a range of maize development, genotypes, and hybrids. While a relatively small proportion of TE families are transcribed, expression is highly dynamic with most families exhibiting tissue-specific expression. A large number of TE families were specifically detected in pollen and endosperm, consistent with reproductive dynamics that maintain silencing of TEs in the germ line. We find that B73 transcript abundance is a poor predictor of TE expression in other genotypes and that transcript levels can differ even for shared TEs. Finally, by assessing recombinant inbred line and hybrid transcriptomes, complex patterns of TE transcript abundance across genotypes emerged. Taken together, this study reveals a dynamic contribution of TEs to maize transcriptomes.

Environmental and epigenetic regulation of Rider retrotransposons in tomato

Transposable elements in crop plants are the powerful drivers of phenotypic variation that has been selected during domestication and breeding programs. In tomato, transpositions of the LTR (long terminal repeat) retrotransposon family Rider have contributed to various phenotypes of agronomical interest, such as fruit shape and colour. However, the mechanisms regulating Rider activity are largely unknown. We have developed a bioinformatics pipeline for the functional annotation of retrotransposons containing LTRs and defined all full-length Rider elements in the tomato genome. Subsequently, we showed that accumulation of Rider transcripts and transposition intermediates in the form of extrachromosomal DNA is triggered by drought stress and relies on abscisic acid signalling. We provide evidence that residual activity of Rider is controlled by epigenetic mechanisms involving siRNAs and the RNA-dependent DNA methylation pathway. Finally, we demonstrate the broad distribution of Rider-like elements in other plant species, including crops. Our work identifies Rider as an environment-responsive element and a potential source of genetic and epigenetic variation in plants.

Transposons are major constituents of plant genomes and represent a powerful source of internal genetic and epigenetic variation. For example, domestication of maize has been facilitated by a dramatic change in plant architecture, the consequence of a transposition event. Insertion of transposons near genes often confers quantitative phenotypic variation linked to changes in transcriptional patterns, as documented for blood oranges and grapes. In tomato, the most widely grown fruit crop and model for fleshy fruit biology, occurrences of several beneficial traits related to fruit shape and plant architecture are due to the activity of the transposon family Rider. While Rider represents a unique endogenous source of genetic and epigenetic variation, mechanisms regulating Rider activity remain unexplored. By achieving experimentally-controlled activation of the Rider family, we shed light on the regulation of these transposons by drought stress, signalling by phytohormones, as well as epigenetic pathways. Furthermore, we reveal the presence of Rider-like elements in other economically important crops such as rapeseed, beetroot and quinoa. This suggests that drought-inducible Rider activation could be further harnessed to generate genetic and epigenetic variation for crop breeding, and highlights the potential of transposon-directed mutagenesis for crop improvement.

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.

Old cogs, new tricks: the evolution of gene expression in a chromatin context

Sophisticated gene-regulatory mechanisms probably evolved in prokaryotes billions of years before the emergence of modern eukaryotes, which inherited the same basic enzymatic machineries. However, the epigenomic landscapes of eukaryotes are dominated by nucleosomes, which have acquired roles in genome packaging, mitotic condensation and silencing parasitic genomic elements. Although the molecular mechanisms by which nucleosomes are displaced and modified have been described, just how transcription factors, histone variants and modifications and chromatin regulators act on nucleosomes to regulate transcription is the subject of considerable ongoing study. We explore the extent to which these transcriptional regulatory components function in the context of the evolutionarily ancient role of chromatin as a barrier to processes acting on DNA and how chromatin proteins have diversified to carry out evolutionarily recent functions that accompanied the emergence of differentiation and development in multicellular eukaryotes.

UAP56 associates with DRM2 and is localized to chromatin in Arabidopsis

Repeated sequence expression and transposable element mobilization are tightly controlled by multilayer processes, which include DNA 5′‐cytosine methylation. The RNA‐directed DNA methylation (RdDM) pathway, which uses siRNAs to guide sequence‐specific directed DNA methylation, emerged specifically in plants. RdDM ensures DNA methylation maintenance on asymmetric CHH sites and specifically initiates de novo methylation in all cytosine sequence contexts through the action of DRM DNA methyltransferases, of which DRM2 is the most prominent. The RdDM pathway has been well described, but how DRM2 is recruited onto DNA targets and associates with other RdDM factors remains unknown. To address these questions, we developed biochemical approaches to allow the identification of factors that may escape genetic screens, such as proteins encoded by multigenic families. Through both conventional and affinity purification of DRM2, we identified DEAD box RNA helicases U2AF56 Associated Protein 56 (UAP56a/b), which are widespread among eukaryotes, as new DRM2 partners. We have shown that, similar to DRM2 and other RdDM actors, UAP56 has chromatin‐associated protein properties. We confirmed this association both in vitro and in vivo in reproductive tissues. In addition, our experiments also suggest that UAP56 may exhibit differential distribution in cells depending on plant organ. While originally identified for its role in splicing, our study suggests that UAP56 may also have other roles, and our findings allow us to initiate discussion about its potential role in the RdDM pathway.

CMT3 and SUVH4/KYP silence the exonic Evelknievel retroelement to allow for reconstitution of CMT1 mRNA

BACKGROUND

The Chromomethylase 1 (CMT1) has long been considered a nonessential gene because, in certain Arabidopsis ecotypes, the CMT1 gene is disrupted by the Evelknievel (EK) retroelement, inserted within exon 13, or contains frameshift mutations, resulting in a truncated, non-functional protein. In contrast to other transposable elements, no transcriptional activation of EK was observed under stress conditions (e.g., protoplasting).

RESULTS

We wanted to explore the regulatory pathway responsible for EK silencing in the Ler ecotype and its effect on CMT1 transcription. Methylome databases confirmed that EK retroelement is heavily methylated and methylation is extended toward CMT1 downstream region. Strong transcriptional activation of EK accompanied by significant reduction in non-CG methylation was found in cmt3 and kyp2, but not in ddm1 or RdDM mutants. EK activation in cmt3 and kyp2 did not interfere with upstream CMT1 expression but abolish transcription through the EK. We identified, in wild-type Ler, three spliced variants in which the entire EK is spliced out; one variant (25% of splicing incidents) facilitates proper reconstitution of an intact CMT1 mRNA. We could recover very low amount of the full-length CMT1 mRNA from WT Ler and Col, but not from cmt3 mutant.

CONCLUSIONS

Our findings highlight CMT3-SUVH4/KYP as the major pathway silencing the intragenic EK via inducing non-CG methylation. Furthermore, retroelement insertion within exons (e.g., CMT1) may not lead to a complete abolishment of the gene product when the element is kept silent. Rather the element can be spliced out to bring about reconstruction of an intact, functional mRNA and possibly retrieval of an active protein.

Ambient Temperature-Responsive Mechanisms Coordinate Regulation of Flowering Time

In plants, environmental conditions such as temperature affect survival, growth, and fitness, particularly during key stages such as seedling growth and reproduction. To survive and thrive in changing conditions, plants have evolved adaptive responses that tightly regulate developmental processes such as hypocotyl elongation and flowering time in response to environmental temperature changes. Increases in temperature, coupled with increasing fluctuations in local climate and weather, severely affect our agricultural systems; therefore, understanding the mechanisms by which plants perceive and respond to temperature is critical for agricultural sustainability. In this review, we summarize recent findings on the molecular mechanisms of ambient temperature perception as well as possible temperature sensing components in plants. Based on recent publications, we highlight several temperature response mechanisms, including the deposition and eviction of histone variants, DNA methylation, alternative splicing, protein degradation, and protein localization. We discuss roles of each proposed temperature-sensing mechanism that affects plant development, with an emphasis on flowering time. Studies of plant ambient temperature responses are advancing rapidly, and this review provides insights for future research aimed at understanding the mechanisms of temperature perception and responses in plants.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Epigenetic activation of meiotic recombination near Arabidopsis thaliana centromeres via loss of H3K9me2 and non-CG DNA methylation

Eukaryotic centromeres contain the kinetochore, which connects chromosomes to the spindle allowing segregation. During meiosis, centromeres are suppressed for inter-homolog crossover, as recombination in these regions can cause chromosome missegregation and aneuploidy. Plant centromeres are surrounded by transposon-dense pericentromeric heterochromatin that is epigenetically silenced by histone 3 lysine 9 dimethylation (H3K9me2), and DNA methylation in CG and non-CG sequence contexts. However, the role of these chromatin modifications in control of meiotic recombination in the pericentromeres is not fully understood. Here, we show that disruption of Arabidopsis thaliana H3K9me2 and non-CG DNA methylation pathways, for example, via mutation of the H3K9 methyltransferase genes KYP/SUVH4 SUVH5 SUVH6, or the CHG DNA methyltransferase gene CMT3, increases meiotic recombination in proximity to the centromeres. Using immunocytological detection of MLH1 foci and genotyping by sequencing of recombinant plants, we observe that H3K9me2 and non-CG DNA methylation pathway mutants show increased pericentromeric crossovers. Increased pericentromeric recombination in H3K9me2/non-CG mutants occurs in hybrid and inbred backgrounds and likely involves contributions from both the interfering and noninterfering crossover repair pathways. We also show that meiotic DNA double-strand breaks (DSBs) increase in H3K9me2/non-CG mutants within the pericentromeres, via purification and sequencing of SPO11-1-oligonucleotides. Therefore, H3K9me2 and non-CG DNA methylation exert a repressive effect on both meiotic DSB and crossover formation in plant pericentromeric heterochromatin. Our results may account for selection of enhancer trap Dissociation (Ds) transposons into the CMT3 gene by recombination with proximal transposon launch-pads.

Induction of epigenetic variation in Arabidopsis by over-expression of DNA METHYLTRANSFERASE1 (MET1)

Epigenetic marks such as DNA methylation and histone modification can vary among plant accessions creating epi-alleles with different levels of expression competence. Mutations in epigenetic pathway functions are powerful tools to induce epigenetic variation. As an alternative approach, we investigated the potential of over-expressing an epigenetic function, using DNA METHYLTRANSFERASE1 (MET1) for proof-of-concept. In Arabidopsis thaliana, MET1 controls maintenance of cytosine methylation at symmetrical CG positions. At some loci, which contain dense DNA methylation in CG- and non-CG context, loss of MET1 causes joint loss of all cytosines methylation marks. We find that over-expression of both catalytically active and inactive versions of MET1 stochastically generates new epi-alleles at loci encoding transposable elements, non-coding RNAs and proteins, which results for most loci in an increase in expression. Individual transformants share some common phenotypes and genes with altered gene expression. Altered expression states can be transmitted to the next generation, which does not require the continuous presence of the MET1 transgene. Long-term stability and epigenetic features differ for individual loci. Our data show that over-expression of MET1, and potentially of other genes encoding epigenetic factors, offers an alternative strategy to identify epigenetic target genes and to create novel epi-alleles.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Analysis of Global Genome Methylation Using the Cytosine-Extension Assay

DNA methylation is a reversible covalent chemical modification of DNA intended to regulate chromatin structure and gene expression in a cell- and tissue-specific manner and in response to the environment. Cytosine methylation is predominantly occurring in plants, and cytosine nucleotides in plants can be methylated at symmetrical (CpG and CpHpG) and nonsymmetrical sites. Although there exists a number of various methods for the detection of cytosine methylation, most of them are either laborious or expensive or both. Here, we describe a quick inexpensive method for the analysis of global genome methylation using a cytosine-extension assay. The assay can be used for the analysis of the total level of CpG, CpHpG, and CpHpH methylation in a given sample of plant DNA.

Analysis of somaclonal variation in transgenic and regenerated plants of Arabidopsis thaliana using methylation related metAFLP and TMD markers

We provide evidence that nucleotide sequence and methylation status changes occur in the Arabidopsis genome during in vitro tissue culture at a frequency high enough to represent an important source of variation. Somaclonal variation is a general consequence of the tissue culture process that has to be analyzed specifically when regenerated plants are obtained in any plant species. Currently, there are few studies about the variability comprising sequence changes and methylation status at the DNA level, generated by the culture of A. thaliana cells and tissues. In this work, two types of highly reproducible molecular markers, modified methylation sensitive AFLP (metAFLP) and transposon methylation display (TMD) have been used for the first time in this species to analyze the nucleotide and cytosine methylation changes induced by transformation and tissue culture protocols. We found significantly higher average methylation values (7.5%) in regenerated and transgenic plants when compared to values obtained from seed derived plants (3.2%) and that the main component of the somaclonal variation present in Arabidopsis clonal plants is genetic rather than epigenetic. However, we have found that the Arabidopsis regenerated and transgenic plants had a higher number of non-fully methylated sites flanking transposable elements than the control plants, and therefore, their mobilization can be facilitated. These data provide further evidence that changes in nucleotide sequence and methylation status occur in the Arabidopsis genome during in vitro tissue culture frequently enough to be an important source of variation in this species.

Cork Oak Young and Traumatic Periderms Show PCD Typical Chromatin Patterns but Different Chromatin-Modifying Genes Expression

Plants are subjected to adverse conditions being outer protective tissues fundamental to their survival. Tree stems are enveloped by a periderm made of cork cells, resulting from the activity of the meristem phellogen. DNA methylation and histone modifications have important roles in the regulation of plant cell differentiation. However, studies on its involvement in cork differentiation are scarce despite periderm importance. Cork oak periderm development was used as a model to study the formation and differentiation of secondary protective tissues, and their behavior after traumatic wounding (traumatic periderm). Nuclei structural changes, dynamics of DNA methylation, and posttranslational histone modifications were assessed in young and traumatic periderms, after cork harvesting. Lenticular phellogen producing atypical non-suberized cells that disaggregate and form pores was also studied, due to high impact for cork industrial uses. Immunolocalization of active and repressive marks, transcription analysis of the corresponding genes, and correlations between gene expression and cork porosity were investigated. During young periderm development, a reduction in nuclei area along with high levels of DNA methylation occurred throughout epidermis disruption. As cork cells became more differentiated, whole nuclei progressive chromatin condensation with accumulation in the nuclear periphery and increasing DNA methylation was observed. Lenticular cells nuclei were highly fragmented with faint 5-mC labeling. Phellogen nuclei were less methylated than in cork cells, and in lenticular phellogen were even lower. No significant differences were detected in H3K4me3 and H3K18ac signals between cork cells layers, although an increase in H3K4me3 signals was found from the phellogen to cork cells. Distinct gene expression patterns in young and traumatic periderms suggest that cork differentiation might be under specific silencing regulatory pathways. Significant correlations were found between QsMET1, QsMET2, and QsSUVH4 gene expression and cork porosity. This work evidences that DNA methylation and histone modifications play a role in cork differentiation and epidermis induced tension-stress. It also provides the first insights into chromatin dynamics during cork and lenticular cells differentiation pointing to a distinct type of remodeling associated with cell death.

Tracking the genome-wide outcomes of a transposable element burst over decades of amplification

Rice (Oryza sativa) has a unique combination of attributes that made it an ideal host to track the natural behavior of very active transposable elements (TEs) over generations. In this study, we have exploited its small genome and propagation by self or sibling pollination to identify and characterize two strain pairs, EG4/HEG4 and A119/A123, undergoing bursts of the nonautonomous miniature inverted repeat transposable element mPing. Comparative sequence analyses of these strains have advanced our understanding of (i) factors that contribute to sustaining a TE burst for decades, (ii) features that distinguish a natural TE burst from bursts in cell culture or mutant backgrounds, and (iii) the extent to which TEs can rapidly diversify the genome of an inbred organism.

To understand the success strategies of transposable elements (TEs) that attain high copy numbers, we analyzed two pairs of rice (Oryza sativa) strains, EG4/HEG4 and A119/A123, undergoing decades of rapid amplification (bursts) of the class 2 autonomous Ping element and the nonautonomous miniature inverted repeat transposable element (MITE) mPing. Comparative analyses of whole-genome sequences of the two strain pairs validated that each pair has been maintained for decades as inbreds since divergence from their respective last common ancestor. Strains EG4 and HEG4 differ by fewer than 160 SNPs and a total of 264 new mPing insertions. Similarly, strains A119 and A123 exhibited about half as many SNPs (277) as new mPing insertions (518). Examination of all other potentially active TEs in these genomes revealed only a single new insertion out of ∼40,000 loci surveyed. The virtual absence of any new TE insertions in these strains outside the mPing bursts demonstrates that the Ping/mPing family gradually attains high copy numbers by maintaining activity and evading host detection for dozens of generations. Evasion is possible because host recognition of mPing sequences appears to have no impact on initiation or maintenance of the burst. Ping is actively transcribed, and both Ping and mPing can transpose despite methylation of terminal sequences. This finding suggests that an important feature of MITE success is that host recognition does not lead to the silencing of the source of transposase.

Biogeographic conservation of the cytosine epigenome in the globally important marine, nitrogen-fixing cyanobacterium Trichodesmium

Cytosine methylation has been shown to regulate essential cellular processes and impact biological adaptation. Despite its evolutionary importance, only a handful of bacterial, genome-wide cytosine studies have been conducted, with none for marine bacteria. Here, we examine the genome-wide, C⁵ -Methyl-cytosine (m5C) methylome and its correlation to global transcription in the marine nitrogen-fixing cyanobacterium Trichodesmium. We characterize genome-wide methylation and highlight conserved motifs across three Trichodesmium isolates and two Trichodesmium metagenomes, thereby identifying highly conserved, novel genomic signatures of potential gene regulation in Trichodesmium. Certain gene bodies with the highest methylation levels correlate with lower expression levels. Several methylated motifs were highly conserved across spatiotemporally separated Trichodesmium isolates, thereby elucidating biogeographically conserved methylation potential. These motifs were also highly conserved in Trichodesmium metagenomic samples from natural populations suggesting them to be potential in situ markers of m5C methylation. Using these data, we highlight predicted roles of cytosine methylation in global cellular metabolism providing evidence for a 'core' m5C methylome spanning different ocean regions. These results provide important insights into the m5C methylation landscape and its biogeochemical implications in an important marine N₂ -fixer, as well as advancing evolutionary theory examining methylation influences on adaptation.

The Evolutionary Dynamics of Orthologs That Shift in Gene Body Methylation between Arabidopsis Species

DNA methylation labels a specific subset of genes in plant genomes. Recent work has shown that this gene-body methylation (gbM) is a conserved feature of orthologs, because highly methylated genes in one species tend to be highly methylated in another. In this study, we examined the exceptions to that rule by identifying genes that differ in gbM status between two plant species-Arabidopsis thaliana and Arabidopsis lyrata. Using Capsella grandiflora as an outgroup, we polarized the loss and gain of gbM for orthologs in the Arabidopsis lineage. Our survey identified a few hundred genes that differed between ingroup species, out of ∼9,000 orthologs. The estimated rate of gbM gain was ∼2 × 10-9 per gene per year for both ingroup taxa and was similar to the loss rate in A. lyrata. In contrast, A. thaliana had a ∼3-fold higher estimated rate of gbM loss per gene, consistent with a recent diminishment of genome size. As in previous studies, we found that body-methylated genes were expressed broadly across tissues and were also longer than other genic sets. Genes that differed in gbM status exhibited higher variance in expression between species than genes that were body-methylated in both species. Moreover, the gain of gbM in one lineage tended to be associated with an increase of expression in that lineage. The genes that varied in gbM status between species varied more significantly in length between species than other sets of genes; we hypothesize that length is a key feature in the transition between body-methylated and nonmethylated genes.

Epigenomic and functional analyses reveal roles of epialleles in the loss of photoperiod sensitivity during domestication of allotetraploid cottons

Background

Polyploidy is a pervasive evolutionary feature of all flowering plants and some animals, leading to genetic and epigenetic changes that affect gene expression and morphology. DNA methylation changes can produce meiotically stable epialleles, which are transmissible through selection and breeding. However, the relationship between DNA methylation and polyploid plant domestication remains elusive.

Results

We report comprehensive epigenomic and functional analyses, including ~12 million differentially methylated cytosines in domesticated allotetraploid cottons and their tetraploid and diploid relatives. Methylated genes evolve faster than unmethylated genes; DNA methylation changes between homoeologous loci are associated with homoeolog-expression bias in the allotetraploids. Significantly, methylation changes induced in the interspecific hybrids are largely maintained in the allotetraploids. Among 519 differentially methylated genes identified between wild and cultivated cottons, some contribute to domestication traits, including flowering time and seed dormancy. CONSTANS (CO) and CO-LIKE (COL) genes regulate photoperiodicity in Arabidopsis. COL2 is an epiallele in allotetraploid cottons. COL2A is hypermethylated and silenced, while COL2D is repressed in wild cottons but highly expressed due to methylation loss in all domesticated cottons tested. Inhibiting DNA methylation activates COL2 expression, and repressing COL2 in cultivated cotton delays flowering.

Conclusions

We uncover epigenomic signatures of domestication traits during cotton evolution. Demethylation of COL2 increases its expression, inducing photoperiodic flowering, which could have contributed to the suitability of cotton for cultivation worldwide. These resources should facilitate epigenetic engineering, breeding, and improvement of polyploid crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1229-8) contains supplementary material, which is available to authorized users.

Arabidopsis proteins with a transposon-related domain act in gene silencing

Transposable elements (TEs) are prevalent in most eukaryotes, and host genomes have devised silencing strategies to rein in TE activity. One of these, transcriptional silencing, is generally associated with DNA methylation and short interfering RNAs. Here we show that the Arabidopsis genes MAIL1 and MAIN define an alternative silencing pathway independent of DNA methylation and short interfering RNAs. Mutants for MAIL1 or MAIN exhibit release of silencing and appear to show impaired condensation of pericentromeric heterochromatin. Phylogenetic analysis suggests not only that MAIL1 and MAIN encode a retrotransposon-related plant mobile domain, but also that host plant mobile domains were captured by DNA transposons during plant evolution. Our results reveal a role for Arabidopsis proteins with a transposon-related domain in gene silencing.

The DNA methylation level against the background of the genome size and t-heterochromatin content in some species of the genus Secale L

Methylation of cytosine in DNA is one of the most important epigenetic modifications in eukaryotes and plays a crucial role in the regulation of gene activity and the maintenance of genomic integrity. DNA methylation and other epigenetic mechanisms affect the development, differentiation or the response of plants to biotic and abiotic stress. This study compared the level of methylation of cytosines on a global (ELISA) and genomic scale (MSAP) between the species of the genus Secale. We analyzed whether the interspecific variation of cytosine methylation was associated with the size of the genome (C-value) and the content of telomeric heterochromatin. MSAP analysis showed that S. sylvestre was the most distinct species among the studied rye taxa; however, the results clearly indicated that these differences were not statistically significant. The total methylation level of the studied loci was very similar in all taxa and ranged from 60% in S. strictum ssp. africanum to 66% in S. cereale ssp. segetale, which confirmed the lack of significant differences in the sequence methylation pattern between the pairs of rye taxa. The level of global cytosine methylation in the DNA was not significantly associated with the content of t-heterochromatin and did not overlap with the existing taxonomic rye relationships. The highest content of 5-methylcytosine was found in S. cereale ssp. segetale (83%), while very low in S. strictum ssp. strictum (53%), which was significantly different from the methylation state of all taxa, except for S. sylvestre. The other studied taxa of rye had a similar level of methylated cytosine ranging from 66.42% (S. vavilovii) to 74.41% in (S. cereale ssp. afghanicum). The results obtained in this study are evidence that the percentage of methylated cytosine cannot be inferred solely based on the genome size or t-heterochromatin. This is a significantly more complex issue.

The Combined Bisulfite Restriction Analysis (COBRA) Assay for the Analysis of Locus-Specific Changes in Methylation Patterns

DNA methylation is a heritable but reversible epigenetic mechanism of control over gene expression. The level of DNA methylation of specific genomic regions correlates with chromatin condensation, the level of gene expression, and in some cases genome stability and the frequency of homologous recombination. Here, we describe the combined bisulfite restriction analysis (COBRA) assay that allows analyzing the methylation status at a specific locus. The protocol consists of the following major steps: bisulfite conversion of non-methylated cytosines to uracils, the locus-specific PCR amplification of converted DNA, restriction digestion, the analysis of restriction patterns on the gel, and the quantification of these restriction patterns using ImageJ or a similar program.