Cytosine Methylation Patterns Research Articles

Genome-wide methylome stability and parental effects in the worldwide distributed Lombardy poplar

BackgroundDespite the increasing number of epigenomic studies in plants, little is known about the forces that shape the methylome in long-lived woody perennials. The Lombardy poplar offers an ideal opportunity to investigate the impact of the individual environmental history of trees on the methylome.ResultsWe present the results of three interconnected experiments on Lombardy poplar. In the first experiment, we investigated methylome variability during a growing season and across vegetatively reproduced generations. We found that ramets collected over Europe and raised in common conditions have stable methylomes in symmetrical CG-contexts. In contrast, seasonal dynamics occurred in methylation patterns in CHH context. In the second experiment, we investigated whether methylome patterns of plants grown in a non-parental environment correlate with the parental climate. We did not observe a biological relevant pattern that significantly correlates with the parental climate. Finally, we investigated whether the parental environment has persistent carry-over effects on the vegetative offspring’s phenotype. We combined new bud set observations of three consecutive growing seasons with former published bud set data. Using a linear mixed effects analysis, we found a statistically significant but weak short-term, parental carry-over effect on the timing of bud set. However, this effect was negligible compared to the direct effects of the offspring environment.ConclusionsGenome-wide cytosine methylation patterns in symmetrical CG-context are stable in Lombardy poplar and appear to be mainly the result of random processes. In this widespread poplar clone, methylation patterns in CG-context can be used as biomarkers to infer a common ancestor and thus to investigate the recent environmental history of a specific Lombardy poplar. The Lombardy poplar shows high phenotypic plasticity in a novel environment which enabled this clonal tree to adapt and survive all over the temperate regions of the world.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation–methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular–spatial and regulatory genome diversity of the mouse brain.

Dynamics of RNA m5C modification during brain development

Post-transcriptional RNA modifications have been recognized as key regulators of neuronal differentiation and synapse development in the mammalian brain. While distinct sets of 5-methylcytosine (m5C) modified mRNAs have been detected in neuronal cells and brain tissues, no study has been performed to characterize methylated mRNA profiles in the developing brain. Here, together with regular RNA-seq, we performed transcriptome-wide bisulfite sequencing to compare RNA cytosine methylation patterns in neural stem cells (NSCs), cortical neuronal cultures, and brain tissues at three postnatal stages. Among 501 m5C sites identified, approximately 6% are consistently methylated across all five conditions. Compared to m5C sites identified in NSCs, 96% of them were hypermethylated in neurons and enriched for genes involved in positive transcriptional regulation and axon extension. In addition, brains at the early postnatal stage demonstrated substantial changes in both RNA cytosine methylation and gene expression of RNA cytosine methylation readers, writers, and erasers. Furthermore, differentially methylated transcripts were significantly enriched for genes regulating synaptic plasticity. Altogether, this study provides a brain epitranscriptomic dataset as a new resource and lays the foundation for further investigations into the role of RNA cytosine methylation during brain development.

DNA methylation aging and transcriptomic studies in horses

Cytosine methylation patterns have not yet been thoroughly studied in horses. Here, we profile n = 333 samples from 42 horse tissue types at loci that are highly conserved between mammalian species using a custom array (HorvathMammalMethylChip40). Using the blood and liver tissues from horses, we develop five epigenetic aging clocks: a multi-tissue clock, a blood clock, a liver clock and two dual-species clocks that apply to both horses and humans. In addition, using blood methylation data from three additional equid species (plains zebra, Grevy’s zebras and Somali asses), we develop another clock that applies across all equid species. Castration does not significantly impact the epigenetic aging rate of blood or liver samples from horses. Methylation and RNA data from the same tissues define the relationship between methylation and RNA expression across horse tissues. We expect that the multi-tissue atlas will become a valuable resource.

The phytochemical brazilin suppress DNMT1 expression by recruiting p53 to its promoter resulting in the epigenetic restoration of p21 in MCF7cells

BackgroundCancer is an outcome of uncontrolled cell division eventually associated with dysregulated epigenetic mechanisms, including DNA methylation. DNA methyltransferase 1 is ubiquitously expressed in the proliferating cells and is essential for the maintenance of DNA methylation. It causes the abnormal silencing of tumor suppressor genes in human cancer which is necessary for proliferation, cell cycle progression, and survival. DNMT1 is involved in tumorigenesis of several cancers, its upregulation potentially upscale the promoter level inactivation of transcription of a tumor inhibitory gene by introducing repressive methylation marks on the CpG islands. This epigenetic perturbation caused by DNMT is targeted for cancer therapeutics. PurposeTo demonstrate the proliferative inhibitory potential of brazilin in human breast cancer cell line (MCF-7) with concurrent mitigation of DNMT1 functional expression and to understand its effect on downstream targets like cell cycle inhibitor p21. Study Design/ MethodsThe impact of brazilin on the growth and proliferation of the MCF-7 cells was determined using the XTT assay. The global DNA 5-methyl cytosine methylation pattern was analyzed upon brazilin treatment. The gene and protein expression of DNMTs were determined with quantitative RTPCR and western blots respectively. The potential binding sites of transcription factors in the human DNMT1 promoter were predicted using the MatInspector tool on the Genomatix software. The chromatin immunoprecipitation (ChIP) assay was performed to demonstrate the transcription factors occupancy at the promoter. Methylation of promoter CpG islands was determined by the methylation-specific PCR (MSP) upon brazilin treatment. The molecular docking of the human DNMT1 with brazilin (ligand) was performed using the Schrödinger suite. ResultsThe heterotetracyclic compound brazilin, present in the wood of Caesalpinia sappan, inhibited the proliferation of the human breast cancer cell line (MCF-7) and reduced the DNMT1 expression with a decrease in global DNA methylation. Brazilin, by activating p38 MAPK and elevating p53 levels within the exposed cells. The elevated level of p53 enriched the occupancy at binding sites within 200 bp upstream to the transcription start site in the DNMT1 promoter, resulting in reduced DNMT1 gene expression. Furthermore, the brazilin restored the p21 levels in the exposed cells as the CpGs in the p21 promoter (-128 bp/+17 bp) were significantly demethylated as observed in the methylation-specific PCR (MSP). ConclusionHighly potential anti-proliferative molecule brazilin can modulate the DNMT1 functional expression and restore the cell cycle inhibitor p21expression. We propose that brazilin can be used in therapeutic interventions to restore the deregulated epigenetic mechanisms in cancer.

Tissue and/or developmental stage specific methylation of nrDNA in Capsicum annuum.

The nuclear ribosomal DNA (nrDNA) sequences are often used for phylogenetic analysis among organisms. Because DNA cytosine methylation and nucleolar dominancy are two common epigenetic mechanisms of nrDNA, we hypothesized that internal transcribed spacer 1 (ITS1), 5.8S rRNA and ITS2 of nrDNA sequences could be used as epigenetic biomarkers. Thus, this research was undertaken to study level and pattern of site-specific cytosine methylation of ITS1, 5.8S and ITS2 in nine tissues and/or developmental stage of pepper Capsicum annuum L. cultivar Demre Sivrisi. Tissues studied consisted of young and old roots at 30 and 90days after sowing (das), mature dry seeds and seeds at 26days of post anthesis (dpa), flowering buds at 1day before flowering, pericarps at 3, 15 and 65 dpa. Levels and patterns of DNA cytosine methylation were identified at single base resolution using bisulfite conversion sequencing. Results of this study revealed that DNA cytosine level and pattern of ITS1, 5.8S and ITS2 were different in most tissues and/or developmental stages studied. In addition, methylation levels of CG, CHG and CHH contexts were also significantly different among the regions. Based on the findings of this study, it was concluded that high level of methylation of nrDNA sequences was relatively higher as observed in transposable element and promoter. On the other hand, its tissue-specific gene expression was effective as that of gene body and promoter methylation. Overall findings revealed that methylation levels of nrDNA could be used as biomarkers for tissue identification or age estimation in plants.

Evolutionary Contribution of Duplicated Genes to Genome Evolution in the Ginseng Species Complex.

Genes duplicated by whole genome duplication (WGD) and small-scale duplication (SSD) have played important roles in adaptive evolution of all flowering plants. However, it still remains underinvestigated how the distinct models of duplication events and their contending evolutionary patterns have shaped the genome and epigenomes of extant plant species. In this study, we investigated the contribution of the WGD- and SSD-derived duplicate genes to the genome evolution of one diploid and three closely related allotetraploid Panax species based on genome, methylome, and proteome data sets. Our genome-wide comparative analyses revealed that although the ginseng species complex was recently diverged, they have evolved distinct overall patterns of nucleotide variation, cytosine methylation, and protein-level expression. In particular, genetic and epigenetic asymmetries observed in the recent WGD-derived genes are largely consistent across the ginseng species complex. In addition, our results revealed that gene duplicates generated by ancient WGD and SSD mechanisms exhibited distinct evolutionary patterns. We found the ancient WGD-derived genes (i.e., ancient collinear gene) are genetically more conserved and hypomethylated at the cytosine sites. In contrast, some of the SSD-derived genes (i.e., dispersal duplicated gene) showed hypermethylation and high variance in nucleotide variation pattern. Functional enrichment analyses of the duplicated genes indicated that adaptation-related traits (i.e., photosynthesis) created during the distant ancient WGDs are further strengthened by both the more recent WGD and SSD. Together, our findings suggest that different types of duplicated genes may have played distinct but relaying evolutionary roles in the polyploidization and speciation processes in the ginseng species complex.

Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.)

BackgroundIn fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from ‘omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments.ResultsWe retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs.ConclusionOur modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

Variations of Cytosine Methylation Patterns between Staminate and Perfect Flowers within Andromonoecious Taihangia rupestris (Rosaceae) Revealed by Methylation-Sensitive Amplification Polymorphism

DNA methylation is an epigenetic modification that plays important roles in flower development through regulation of gene expression. Taihangia rupestris is an andromonoecious plant species, which can produce staminate and perfect flowers within the same individual. In this study, we employed methylation-sensitive amplified polymorphism (MSAP) to investigate cytosine methylation patterns in staminate and perfect flowers within the andromonoecious T. rupestris. Our data showed the cytosine methylation patterns varied in staminate and perfect flowers at early and late developmental stages. The level of cytosine methylation was slightly higher in perfect flowers than in staminate flowers, and an increase of cytosine methylation levels was found alongside developmental stages in both staminate and perfect flowers. Based on MSAP profiles, a total of 116 differentially methylated fragments (DMFs) were obtained and sequenced. Of these obtained DMFs, 42 fragments showed sequence similarity to functional genes involved in a wide range of biological processes such as flower development, signal transduction, and cell proliferation, and transcription regulation. By integrative analysis of MSAP and transcriptome data, we found transcript levels correlated with cytosine methylation patterns for majority of DMF-associated genes. This study revealed that DNA methylation could play important roles in regulation of gene expression during establishment and development of staminate and perfect flower within andromonoecious T. rupestris and would be helpful to shed light on the underlying molecular mechanisms responsible for flower formation at epigenetic scale in andromonoecious sexual system.

Sugarcane mosaic virus mediated changes in cytosine methylation pattern and differentially transcribed fragments in resistance-contrasting sugarcane genotypes.

Sugarcane mosaic virus (SCMV) is the causal agent of sugarcane mosaic disease (SMD) in Brazil; it is mainly controlled by using resistant cultivars. Studies on the changes in sugarcane transcriptome provided the first insights about the molecular basis underlying the genetic resistance to SMD; nonetheless, epigenetic modifications such as cytosine methylation is also informative, considering its roles in gene expression regulation. In our previous study, differentially transcribed fragments (DTFs) were obtained using cDNA-amplified fragment length polymorphism by comparing mock- and SCMV-inoculated plants from two sugarcane cultivars with contrasting responses to SMD. In this study, the identification of unexplored DTFs was continued while the same leaf samples were used to evaluate SCMV-mediated changes in the cytosine methylation pattern by using methylation-sensitive amplification polymorphism. This analysis revealed minor changes in cytosine methylation in response to SCMV infection, but distinct changes between the cultivars with contrasting responses to SMD, with higher hypomethylation events 24 and 72 h post-inoculation in the resistant cultivar. The differentially methylated fragments (DMFs) aligned with transcripts, putative promoters, and genomic regions, with a preponderant distribution within CpG islands. The transcripts found were associated with plant immunity and other stress responses, epigenetic changes, and transposable elements. The DTFs aligned with transcripts assigned to stress responses, epigenetic changes, photosynthesis, lipid transport, and oxidoreductases, in which the transcriptional start site is located in proximity with CpG islands and tandem repeats. Real-time quantitative polymerase chain reaction results revealed significant upregulation in the resistant cultivar of aspartyl protease and VQ protein, respectively, selected from DMF and DTF alignments, suggesting their roles in genetic resistance to SMD and supporting the influence of cytosine methylation in gene expression. Thus, we identified new candidate genes for further validation and showed that the changes in cytosine methylation may regulate important mechanisms underlying the genetic resistance to SMD.

Abstract IA05: Microenvironmental signatures reveal biologic and clinical subtypes of diffuse large B-cell lymphoma

Abstract Diffuse large B-cell lymphoma (DLBCL) is a biologically and clinically heterogenous disease almost invariably treated with combinatorial chemoimmunotherapy resulting in a ~65% cure rate in first line. Further improvement of treatment outcome relies on elucidating the biology that underlies the clinical behavior of this heterogenous disease. Transcriptomic and genetic characterization of DLBCL have increased our understanding of its pathogenesis and provided potential novel therapeutic targets. Over the past few years, the role of the microenvironment in DLBCL biology has been increasingly recognized. A guiding framework is necessary to effectively and safely translate research evidence into clinical practice. Towards this goal, we characterized the DLBCL microenvironment based on the transcriptional footprint of microenvironment cells and processes in a large cohort of patients. We first developed transcriptional signatures reflecting distinct cellular subtypes of the microenvironment, biologic processes, and canonical signaling pathway activation. We then used these signatures to virtually reconstruct the lymphoma microenvironment (LME) by applying an unsupervised community detection algorithm in 4,656 DLBCL cases. As a result, we described four classes of LMEs that reflect distinct enrichment of associated signatures. 1. “Germinal center-like” (GC-like) due to the presence of signatures from cell types commonly found in germinal centers including follicular dendritic cells, lymphatic endothelial cells, total T cells, and several CD4+ T-cell subpopulations. 2. “Mesenchymal” (MS) for the abundance of signatures from stromal cells such as cancer-associated fibroblasts, fibroblastic reticular cells, vascular endothelial cells, and extracellular matrix pathways. The MS-LME class was enriched in mutations affecting antigen presentation. These lymphomas also presented with high activity of the TGFB/SMAD and HIF pathways, previously associated with favorable prognosis in DLBCL. 3. “Inflammatory” (IN) for the presence of signatures from macrophages, neutrophils, NK cells, and cytotoxic T cells. The IN-LME class was enriched for immune suppressive/pro-lymphoma cytokines, expression of the immune checkpoint molecules PD-L1 and IDO1, and higher activity of the NF-kB, JAK/STAT, and TNF signaling pathways. 4. A “depleted” (DP) form that, contrasting with the other LMEs, was characterized by an overall lower presence of microenvironment-derived signatures. We analyzed the association of LME classes with overall survival (OS) and progression-free survival (PFS) in patient cohorts of, respectively, 2,646 and 2,189 DLBCLs treated with rituximab-based chemoimmunotherapy. Analysis of OS by LME indicated significant differences in prognosis from better to poor as follows: GC-like, MS (p =0.03, vs. GC-like), IN (p = 0.006, vs. MS), and DP (p = 0.008, vs. IN) LMEs, whereas GC-like and MS have similarly favorable PFS curves (p = 0.9) followed by IN (p < 0.01) and DP (p = 0.03, vs. IN), which presented with the poorest PFS. In multivariate Cox-proportional hazard analysis controlled for C.O.O. and I.P.I., the LME subtypes remained informative with GC-like and MS LMEs associated with better outcome. When segregated by C.O.O., the DP-LME retained the poorest OS and PFS in both C.O.O. subtypes. However, in ABC-DLBCL, the LME with the best PFS and OS was MS, and in GCB-DLBCL, the best LME was GC-like, suggesting that the biologic impact of the microenvironment may be different depending on the lymphoma subtype. DHL and high-grade lymphomas (HGL) modified their outcome when LME was considered. DHL/HGL with the favorable prognosis LMEs GC-like and MS had significantly better PFS and OS than DHL/HGL with the unfavorable prognosis LMEs IN and DP. DLBCLs with a DP-LME were characterized by a minimal presence of LME signatures and tumor cells exhibited higher clonality, were enriched for mutations affecting TP53 and cell proliferation genes, as well as presented higher levels of aberrant hypermethylation of immune-related and TGFB/SMAD1 genes. Preclinical murine models and patient data demonstrated that some of these features can be pharmacologically reversed with hypomethylating agents, indicating that cytosine methylation patterning of lymphoma cells may contribute to a depleted LME. In summary, we categorized the DLBCL LME into four major transcriptionally defined subtypes with distinct biologic properties and clinical behavior, which complement genetically defined subtypes of DLBCL in guiding the development of rational therapeutic approaches for these patients. Citation Format: Leandro Cerchietti. Microenvironmental signatures reveal biologic and clinical subtypes of diffuse large B-cell lymphoma [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr IA05.

The Effect of Iodine-Containing Nano-Micelles, FS-1, on Antibiotic Resistance, Gene Expression and Epigenetic Modifications in the Genome of Multidrug Resistant MRSA Strain Staphylococcus aureus ATCC BAA-39.

Application of supplementary drugs which increase susceptibility of pathogenic bacteria to antibiotics is a promising yet unexplored approach to overcome the global problem of multidrug-resistant infections. The discovery of a new drug, an iodine-containing nano-molecular complex FS-1, which has proven to improve susceptibility to antibiotics in various pathogens, including MRSA strain Staphylococcus aureus ATCC BAA-39TM, allowed studying this phenomenon. Chromosomal DNA and total RNA samples extracted from the FS-1 treated strain (FS) and from the negative control (NC) cultures were sequenced by PacBio SMRT and Ion Torrent technologies, respectively. PacBio DNA reads were used to assemble chromosomal DNA of the NC and FS variants of S. aureus BAA-39 and to perform profiling of epigenetically modified nucleotides. Results of transcriptional profiling, variant calling and detection of epigenetic modifications in the FS variant were compared to the NC variant. Additionally, the genetic alterations caused by the treatment of S. aureus BAA-39 with FS-1 were compared to the results of a similar experiment conducted with another model organism, E. coli ATCC BAA-196. Several commonalities in responses of these phylogenetically distant microorganisms to the treatment with FS-1 were discovered, which included metabolic transition toward anaerobiosis and oxidative/osmotic stress response. S. aureus culture appeared to be more sensitive to FS-1 due to a higher penetrability of cells by iodine bound compounds, which caused carbonyl stress associated with nucleotide damaging by FS-1, abnormal epigenetic modifications and an increased rate of mutations. It was hypothesized that the disrupted pattern of adenine methylated loci within methicillin-resistance chromosome cassettes (SCCmec) may promote excision of this antibiotic resistance determinant from chromosomes while the altered pattern of cytosine methylation was behind the adaptive gene regulation in the culture FS. The selection against the antibiotic resistance in bacterial populations caused by abnormal epigenetic modifications exemplifies possible mechanisms of antibiotic resistance reversion induced by iodine-containing compounds. These finding will facilitate development of therapeutic agents against multidrug-resistant infections.

Alterations of Rice (Oryza sativa L.) DNA Methylation Patterns Associated with Gene Expression in Response to Rice Black Streaked Dwarf Virus.

Rice black-streaked dwarf virus (RBSDV) causes severe yield losses in rice (Oryza sativa L.) in China. Studies have shown that the mechanisms of DNA methylation-mediated plant defense against DNA viruses and RNA viruses are different. However, in rice its function in response to infection of RBSDV, a double-stranded RNA virus, remains unclear. In this study, high-throughput single-base resolution bisulfite sequencing (BS-Seq) was carried out to analyze the distribution pattern and characteristics of cytosine methylation in RBSDV-infected rice. Widespread differences were identified in CG and non-CG contexts between the RBSDV-infected and RBSDV-free rice. We identified a large number of differentially methylated regions (DMRs) along the genome of RBSDV-infected rice. Additionally, the transcriptome sequencing analysis obtained 1119 differentially expressed genes (DEGs). Correlation analysis of DMRs-related genes (DMGs) and DEGs filtered 102 genes with positive correlation and 71 genes with negative correlation between methylation level at promoter regions and gene expression. Key genes associated with maintaining DNA methylation in rice were analyzed by RT-qPCR and indicated that OsDMT702 might be responsible for the global increase of DNA methylation level in rice under RBSDV stress. Our results suggest important roles of rice DNA methylation in response to RBSDV and provide potential target genes for rice antiviral immunity.

Cytosine methylation patterns suggest a role of methylation in plastic and adaptive responses to temperature in European grayling (Thymallus thymallus) populations

ABSTRACT Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either ‘colder’ or ‘warmer’ spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the ‘colder’ and ‘warmer’ conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

CGGBP1-regulated cytosine methylation at CTCF-binding motifs resists stochasticity

BackgroundThe human CGGBP1 binds to GC-rich regions and interspersed repeats, maintains homeostasis of stochastic cytosine methylation and determines DNA-binding of CTCF. Interdependence between regulation of cytosine methylation and CTCF occupancy by CGGBP1 remains unknown.ResultsBy analyzing methylated DNA-sequencing data obtained from CGGBP1-depleted cells, we report that some transcription factor-binding sites, including CTCF, resist stochastic changes in cytosine methylation. By analysing CTCF-binding sites we show that cytosine methylation changes at CTCF motifs caused by CGGBP1 depletion resist stochastic changes. These CTCF-binding sites are positioned at locations where the spread of cytosine methylation in cis depends on the levels of CGGBP1.ConclusionOur findings suggest that CTCF occupancy and functions are determined by CGGBP1-regulated cytosine methylation patterns.

Inhalation exposure to cigarette smoke and inflammatory agents induces epigenetic changes in the lung

Smoking-related lung tumors are characterized by profound epigenetic changes including scrambled patterns of DNA methylation, deregulated histone acetylation, altered gene expression levels, distorted microRNA profiles, and a global loss of cytosine hydroxymethylation marks. Here, we employed an enhanced version of bisulfite sequencing (RRBS/oxRRBS) followed by next generation sequencing to separately map DNA epigenetic marks 5-methyl-dC and 5-hydroxymethyl-dC in genomic DNA isolated from lungs of A/J mice exposed whole-body to environmental cigarette smoke for 10 weeks. Exposure to cigarette smoke significantly affected the patterns of cytosine methylation and hydroxymethylation in the lungs. Differentially hydroxymethylated regions were associated with inflammatory response/disease, organismal injury, and respiratory diseases and were involved in regulation of cellular development, function, growth, and proliferation. To identify epigenetic changes in the lung associated with exposure to tobacco carcinogens and inflammation, A/J mice were intranasally treated with the tobacco carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), the inflammatory agent lipopolysaccharide (LPS), or both. NNK alone caused minimal epigenetic alterations, while exposure either to LPS or NNK/LPS in combination led to increased levels of global cytosine methylation and formylation, reduced cytosine hydroxymethylation, decreased histone acetylation, and altered expression levels of multiple genes. Our results suggest that inflammatory processes are responsible for epigenetic changes contributing to lung cancer development.

Patterns of DNA methylation changes in elite Eucalyptus clones across contrasting environments

Phenotypic plasticity refers to the ability of a single genotype to express distinct phenotypes in response to the environment, a crucial feature for sessile organisms like forest trees, especially in a scenario of global climate change. Several studies show that epigenetic regulation plays an important role in this plastic adaptation response, driving the search for associations between natural epigenetic variation with environmental cues and phenotypic traits based on patterns of cytosine methylation. Clonally propagated trees across variable sites offer a robust system to control for the confounding effect of the background genetic variation among genotypes, allowing the analysis of epigenetic modifications in response to variable environments. In this study we investigated the overall patterns of epigenetic changes by the analysis of the genome-wide DNA methylation status based on high throughput MS-DArT-seq (Methyl Sensitive DArT-seq sequencing) of reduced genome complexity representations. We compared patterns of DNA methylation of biological replicates of leaf and xylem tissue samples of four commercially planted elite Eucalyptus grandis × Eucalyptus urophylla clones and one Eucalyptus urophylla in two contrasting sites in Brazil, and the association of these methylation patterns with the environments and growth traits. DNA sequence reads were mapped against the Eucalyptus grandis reference genome, counting and annotating differentially methylated sites. A total of 90,378 MS-DArT-seq sites were identified, the majority (~70%) located in genes and 10% in transposable elements. The distribution of methylation sites showed extensive variation between the five genotypes and the environments. Sets of methylation sites exclusive to each location were identified for each clone but no consistently shared epigenetic marks for all five clones were found across environments. Multiple correspondence analysis suggests a significant contribution of the genetic background on the distribution of methylation changes. We used a gene-environment association analysis to search for association of methylation patterns with growth traits. A total of 445 methylation sites across all 11 Eucalyptus chromosomes were found significantly associated with one or more of the three measured traits (total height, estimated volume and breast high diameter). The absence of clustered differentially methylated sites is consistent with the fact that complex growth traits are governed by a large number of loci of small effect across the entire genome suggesting that this same pattern will likely hold for what regards epigenetic marks.

Silencing of Mutator Elements in Maize Involves Distinct Populations of Small RNAs and Distinct Patterns of DNA Methylation.

Transposable elements (TEs) are a ubiquitous feature of plant genomes. Because of the threat they post to genome integrity, most TEs are epigenetically silenced. However, even closely related plant species often have dramatically different populations of TEs, suggesting periodic rounds of activity and silencing. Here, we show that the process of de novo methylation of an active element in maize involves two distinct pathways, one of which is directly implicated in causing epigenetic silencing and one of which is the result of that silencing. Epigenetic changes involve changes in gene expression that can be heritably transmitted to daughter cells in the absence of changes in DNA sequence. Epigenetics has been implicated in phenomena as diverse as development, stress response, and carcinogenesis. A significant challenge facing those interested in investigating epigenetic phenomena is determining causal relationships between DNA methylation, specific classes of small RNAs, and associated changes in gene expression. Because they are the primary targets of epigenetic silencing in plants and, when active, are often targeted for de novo silencing, TEs represent a valuable source of information about these relationships. We use a naturally occurring system in which a single TE can be heritably silenced by a single derivative of that TE. By using this system it is possible to unravel causal relationships between different size classes of small RNAs, patterns of DNA methylation, and heritable silencing. Here, we show that the long terminal inverted repeats within Zea mays MuDR transposons are targeted by distinct classes of small RNAs during epigenetic silencing that are dependent on distinct silencing pathways, only one of which is associated with transcriptional silencing of the transposon. Further, these small RNAs target distinct regions of the terminal inverted repeats, resulting in different patterns of cytosine methylation with different functional consequences with respect to epigenetic silencing and the heritability of that silencing.

Epigenetic Changes in Host Ribosomal DNA Promoter Induced by an Asymptomatic Plant Virus Infection

DNA cytosine methylation is one of the main epigenetic mechanisms in higher eukaryotes and is considered to play a key role in transcriptional gene silencing. In plants, cytosine methylation can occur in all sequence contexts (CG, CHG, and CHH), and its levels are controlled by multiple pathways, including de novo methylation, maintenance methylation, and demethylation. Modulation of DNA methylation represents a potentially robust mechanism to adjust gene expression following exposure to different stresses. However, the potential involvement of epigenetics in plant-virus interactions has been scarcely explored, especially with regard to RNA viruses. Here, we studied the impact of a symptomless viral infection on the epigenetic status of the host genome. We focused our attention on the interaction between Nicotiana benthamiana and Pelargonium line pattern virus (PLPV, family Tombusviridae), and analyzed cytosine methylation in the repetitive genomic element corresponding to ribosomal DNA (rDNA). Through a combination of bisulfite sequencing and RT-qPCR, we obtained data showing that PLPV infection gives rise to a reduction in methylation at CG sites of the rDNA promoter. Such a reduction correlated with an increase and decrease, respectively, in the expression levels of some key demethylases and of MET1, the DNA methyltransferase responsible for the maintenance of CG methylation. Hypomethylation of rDNA promoter was associated with a five-fold augmentation of rRNA precursor levels. The PLPV protein p37, reported as a suppressor of post-transcriptional gene silencing, did not lead to the same effects when expressed alone and, thus, it is unlikely to act as suppressor of transcriptional gene silencing. Collectively, the results suggest that PLPV infection as a whole is able to modulate host transcriptional activity through changes in the cytosine methylation pattern arising from misregulation of methyltransferases/demethylases balance.

Codon usage bias affects α-amylase mRNA level by altering RNA stability and cytosine methylation patterns in Escherichia coli.

Codon usage bias exists in almost every organism and is reported to regulate protein translation efficiency and folding. Besides translation, the preliminary role of codon usage bias on gene transcription has also been revealed in some eukaryotes such as Neurospora crassa. In this study, we took as an example the α-amylase-coding gene (amyA) and examined the role of codon usage bias in regulating gene expression in the typical prokaryote Escherichia coli. We confirmed the higher translation efficiency on codon-optimized amyA RNAs and found that the RNA level itself was also affected by codon optimization. The decreased RNA level was caused at least in part by altered mRNA stability at the post-transcriptional level. Codon optimization also altered the number of cytosine methylation sites. Examination on dcm knockouts suggested that cytosine methylation may be a minor mechanism adopted by codon bias to regulate gene RNA levels. More studies are required to verify the global effect of codon usage and to reveal its detailed mechanism on transcription.