Composite 5-methylations of cytosines modulate i-motif stability in a sequence-specific manner: Implications for DNA nanotechnology and epigenetic regulation of plant telomeric DNA

BackgroundThe i-motif is a tetrameric DNA structure based on the formation of hemiprotonated cytosine-cytosine (C⁺.C) base pairs. i-motifs are widely used in nanotechnology. In biological systems, i-motifs are involved in gene regulation and in control of genome integrity. In vivo, the i-motif forming sequences are subjects of epigenetic modifications, particularly 5-cytosine methylation. In plants, natively occurring methylation patterns lead to a complex network of C⁺.C, ^5mC⁺.C and ^5mC⁺.^5mC base-pairs in the i-motif stem. The impact of complex methylation patterns (CMPs) on i-motif formation propensity is currently unknown.MethodsWe employed CD and UV-absorption spectroscopies, native PAGE, thermal denaturation and quantum-chemical calculations to analyse the effects of native, native-like, and non-native CMPs in the i-motif stem on the i-motif stability and pK_a.ResultsCMPs have strong influence on i-motif stability and pK_a and influence these parameters in sequence-specific manner. In contrast to a general belief, i) CMPs do not invariably stabilize the i-motif, and ii) when the CMPs do stabilize the i-motif, the extent of the stabilization depends (in a complex manner) on the number and pattern of symmetric ^5mC⁺.^5mC or asymmetric ^5mC⁺.C base pairs in the i-motif stem.ConclusionsCMPs can be effectively used to fine-tune i-motif properties. Our data support the notion of epigenetic modifications as a plausible control mechanism of i-motif formation in vivo.General SignificanceOur results have implications in epigenetic regulation of telomeric DNA in plants and highlight the potential and limitations of engineered patterning of cytosine methylations on the i-motif scaffold in nanotechnological applications.     

Purification and characterization of rice DNA methyltransferase

Epigenetic modification is essential for normal development and plays important roles in gene regulation in higher plants. Multiple factors interact to regulate the establishment and maintenance of DNA methylation in plant genome. We had previously cloned and characterized DNA methyltransferase (DNA MTase) gene homologues (OsMET1) from rice. In this present study, determination of DNA MTase activity in different cellular compartments showed that DNA MTase was enriched in nuclei and the activity was remarkably increased during imbibing dry seeds. We had optimized the purification technique for DNA MTase enzyme from shoots of 10-day-old rice seedlings using the three successive chromatographic columns. The Econo-Pac Q, the Hitrap-Heparin and the Superdex-200 columns yielded a protein fraction of a specific activity of 29, 298 and 800 purification folds, compared to the original nuclear extract, respectively. The purified protein preferred hemi-methylated DNA substrate, suggesting the maintenance activity of methylation. The native rice DNA MTase was approximately 160–170 kDa and exhibited a broad pH optimum in the range of 7.6 and 8.0. The enzyme kinetics and inhibitory effects by methyl donor analogs, base analogs, cations, and cationic amines on rice DNA MTase were examined. Global cytosine methylation status of rice genome during development and in various tissue culture systems were monitored and the results suggested that the cytosine methylation level is not directly correlated with the DNA MTase activity. The purification and characterization of rice DNA MTase enzyme are expected to enhance our understanding of this enzyme function and their possible contributions in Gramineae plant development.     

i-Motif formation and spontaneous deletions in human cells

Concatemers of d(TCCC) that were first detected through their association with deletions at the RACK7 locus, are widespread throughout the human genome. Circular dichroism spectra show that d(GGGA)_n sequences form G-quadruplexes when n > 3, while i-motif structures form at d(TCCC)_n sequences at neutral pH when n ≥ 7 in vitro. In the PC3 cell line, deletions are observed only when the d(TCCC)_n variant is long enough to form significant levels of unresolved i-motif structure at neutral pH. The presence of an unresolved i-motif at a representative d(TCCC)_n element at RACK7 was suggested by experiments showing that that the region containing the d(TCCC)₉ element was susceptible to bisulfite attack in native DNA and that d(TCCC)₉ oligo formed an i-motif structure at neutral pH. This in turn suggested that that the i-motif present at this site in native DNA must be susceptible to bisulfite mediated deamination even though it is a closed structure. Bisulfite deamination of the i-motif structure in the model oligodeoxynucleotide was confirmed using mass spectrometry analysis. We conclude that while G-quadruplex formation may contribute to spontaneous mutation at these sites, deletions actually require the potential for i-motif to form and remain unresolved at neutral pH.     

Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.     

Multiplex CRISPR-Cas9 editing of DNA methyltransferases in rice uncovers a class of non-CG methylation specific for GC-rich regions

DNA methylation in the non-CG context is widespread in the plant kingdom and abundant in mammalian tissues such as the brain and pluripotent cells. Non-CG methylation in Arabidopsis thaliana is coordinately regulated by DOMAINS REARRANGED METHYLTRANSFERASE (DRM) and CHROMOMETHYLASE (CMT) proteins but has yet to be systematically studied in major crops due to difficulties in obtaining genetic materials. Here, utilizing the highly efficient multiplex CRISPR-Cas9 genome-editing system, we created single- and multiple-knockout mutants for all the nine DNA methyltransferases in rice (Oryza sativa) and profiled their whole-genome methylation status at single-nucleotide resolution. Surprisingly, the simultaneous loss of DRM2, CHROMOMETHYLASE3 (CMT2), and CMT3 functions, which completely erases all non-CG methylation in Arabidopsis, only partially reduced it in rice. The regions that remained heavily methylated in non-CG contexts in the rice Os-dcc (Osdrm2/cmt2/cmt3a) triple mutant had high GC contents. Furthermore, the residual non-CG methylation in the Os-dcc mutant was eliminated in the Os-ddccc (Osdrm2/drm3/cmt2/cmt3a/cmt3b) quintuple mutant but retained in the Os-ddcc (Osdrm2/drm3/cmt2/cmt3a) quadruple mutant, demonstrating that OsCMT3b maintains non-CG methylation in the absence of other major methyltransferases. Our results showed that OsCMT3b is subfunctionalized to accommodate a distinct cluster of non-CG-methylated sites at highly GC-rich regions in the rice genome.

Examination of knockout mutants reveals that rice methyltransferases have subfunctionalized to accommodate a distinct cluster of non-CG-methylated sites at highly GC-rich regions in the rice genome.     

An epiallele of rice AK1 affects photosynthetic capacity

Epigenetic gene variants,termed epialleles,can broaden genetic and phenotypic diversity in eukaryotes.Here,we identify a natural epiallele of OsAK1,which encodes a rice adenylate kinase.The Epiaki plants show albino in young leaf and panicle with abnormal chloroplast structures.We found that no nucleotide sequence variation but hypermethylation at promoter region caused silencing of OsAK1(Oso8go177o) in Epi-ak1 plants.OsAK1 localizes to chloroplast and many genes associated with photosynthesis processes were downregulated in Epi-ak1.Thus,the work identified a novel rice epiallele caused by DNA methylation changes,shedding light on significant roles of DNA methylation on rice development.     

Analysis of copy-number variation,insertional polymorphism,and methylation status of the tiniest class I (TRIM) and class II (MITE) transposable element families in various rice strains

Transposable elements (TEs) dominate the genetic capacity of most eukaryotes, especially plants, where they may compose up to 90% of the genome. Many studies, both in plants and animals reported that in fact non-autonomous elements that have lost their protein-coding sequences and became miniature elements were highly associated with genes, and showed a high level of transpositional activity such as mPing family in rice. In this study, we have investigated in detail the copy number, insertional polymorphism and the methylation status of the tiniest LTR retrotransposon family, termed TRIM, in nine rice strains, in comparison with mPing. While TRIM showed similar copy numbers (average of 79 insertions) in all the nine rice strains, the copy number of mPing varied dramatically (ranging from 6 to 203 insertions) in the same strains. Site-specific PCR analysis revealed that ~58% of the TRIM elements have identical insertion sites among the nine rice strains, while none of the mPing elements (100% polymorphism) have identical insertion sites in the same strains. Finally, over 65% of the TRIM insertion sites were cytosine methylated in all nine rice strains, while the level of the methylated mPing insertion sites ranged between 43 and 81.5%. The findings of this study indicate that unlike mPing, TRIM is most probably a fossil TE family in rice. In addition, the data shows that there might be a strong correlation between TE methylation and copy number.     

Imprinting in Plants and Its Underlying Mechanisms

Genomic imprinting (or imprinting) refers to an epigenetic phenomenon by which the allelic expression of a gene depends on the parent of origin. It has evolved independently in placental mammals and flowering plants. In plants, imprinting is mainly found in endosperm. Recent genome-wide surveys in Arabidopsis, rice, and maize identified hundreds of imprinted genes in endosperm. Since these genes are of diverse functions, endosperm development is regulated at different regulatory levels. The imprinted expression of only a few genes is conserved between Arabidopsis and monocots, suggesting that imprinting evolved quickly during speciation. In Arabidopsis, DEMETER (DME) mediates hypomethylation in the maternal genome at numerous loci (mainly transposons and repeats) in the central cell and results in many differentially methylated regions between parental genomes in the endosperm, and subsequent imprinted expression of some genes. In addition, histone modification mediated by Polycomb group (PcG) proteins is also involved in regulating imprinting. DME-induced hypomethylated alleles in the central cell are considered to produce small interfering RNAs (siRNAs) which are imported to the egg to reinforce DNA methylation. In parallel, the activity of DME in the vegetative cell of the male gametophyte demethylates many regions which overlap with the demethylated regions in the central cell. siRNAs from the demethylated regions are hypothesized to be also transferred into sperm to reinforce DNA methylation. Imprinting is partly the result of genome-wide epigenetic reprogramming in the central cell and vegetative cell and evolved under different selective pressures.     

Pre-conditioning the epigenetic response to high vapor pressure deficit increases the drought tolerance of Arabidopsis thaliana

Epigenetic modification of the genome via cytosine methylation is a dynamic process that responds to changes in the growing environment. This modification can also be heritable. The combination of both properties means that there is the potential for the life experiences of the parental generation to modify the methylation profiles of their offspring and so potentially to “pre-condition” them to better accommodate abiotic conditions encountered by their parents. We recently identified high vapor pressure deficit (vpd)-induced DNA methylation at 2 gene loci in the stomatal development pathway and an associated reduction in leaf stomatal frequency.¹ Here, we test whether this epigenetic modification pre-conditioned parents and their offspring to the more severe water stress of periodic drought. We found that 3 generations of high vpd-grown plants were better able to withstand periodic drought stress over 2 generations. This resistance was not directly associated with de novo methylation of the target stomata genes, but was associated with the cmt3 mutant’s inability to maintain asymmetric sequence context methylation. If our finding applies widely, it could have significant implications for evolutionary biology and breeding for stressful environments.     

Exploring the Roles of DNA Methylation in the Metal-Reducing Bacterium Shewanella oneidensis MR-1

We performed whole-genome analyses of DNA methylation in Shewanella oneidensis MR-1 to examine its possible role in regulating gene expression and other cellular processes. Single-molecule real-time (SMRT) sequencing revealed extensive methylation of adenine (N6mA) throughout the genome. These methylated bases were located in five sequence motifs, including three novel targets for type I restriction/modification enzymes. The sequence motifs targeted by putative methyltranferases were determined via SMRT sequencing of gene knockout mutants. In addition, we found that S. oneidensis MR-1 cultures grown under various culture conditions displayed different DNA methylation patterns. However, the small number of differentially methylated sites could not be directly linked to the much larger number of differentially expressed genes under these conditions, suggesting that DNA methylation is not a major regulator of gene expression in S. oneidensis MR-1. The enrichment of methylated GATC motifs in the origin of replication indicates that DNA methylation may regulate genome replication in a manner similar to that seen in Escherichia coli. Furthermore, comparative analyses suggest that many Gammaproteobacteria, including all members of the Shewanellaceae family, may also utilize DNA methylation to regulate genome replication.     

Identification of genes required for de novo DNA methylation in Arabidopsis

De novo DNA methylation in Arabidopsis thaliana is catalyzed by the methyltransferase DRM2, a homolog of the mammalian de novo methyltransferase DNMT3. DRM2 is targeted to DNA by small interfering RNAs (siRNAs) in a process known as RNA-directed DNA Methylation (RdDM). While several components of the RdDM pathway are known, a functional understanding of the underlying mechanism is far from complete. We employed both forward and reverse genetic approaches to identify factors involved in de novo methylation. We utilized the FWA transgene, which is methylated and silenced when transformed into wild-type plants, but unmethylated and expressed when transformed into de novo methylation mutants. Expression of FWA is marked by a late-flowering phenotype, which is easily scored in mutant versus wild-type plants. By reverse genetics we discovered the requirement for known RdDM effectors AGO6 and NRPE5a for efficient de novo methylation. A forward genetic approach uncovered alleles of several components of the RdDM pathway, including alleles of clsy1, ktf1 and nrpd/e2, which have not been previously shown to be required for the initial establishment of DNA methylation. Mutations were mapped and genes cloned by both traditional and whole genome sequencing approaches. The methodologies and the mutant alleles discovered will be instrumental in further studies of de novo DNA methylation.Key words: DNA methylation, Arabidopsis, de novo, genetic screen, whole-genome sequencing